Deuterium atom exposure of self-damaged tungsten#
This validation case is NRA measurements of D in self-damaged W performed by Markelj et al. [25].
The experimental procedure included three main phases. A high-purity polycrystalline W sample of 0.8 mm thickness was pre-damaged with 20 MeV W ions. The pre-damaged sample was then exposed to low-energy (~0.3 eV) D atomic flux of \(5.8\times 10^{18}\,\textrm{m}^{-2}\textrm{s}^{-1}\) at 600 K. The exposure continued until the D fluence of \(1\times 10^{24}\,\mathrm{m}^{-2}\) was reached. Finally, an isothermal desorption of D for 43 h. at 600 K was conducted.
The FESTIM model is based on the approach of Hodille et al. [26]. Only isothermal D exposure and desorption phases are simulated omitting intermediate cooling/re-heating steps. For the surface processes, adsorption of low-energy atoms, desorption of molecules (Langmuir-Hinshelwood recombination), and recombination of an adsorbed atom with an incident atom (Eley-Rideal recombination) are considered.
The D diffusivity in W is defined by scaling the corresponding value for H (Fernandez et al. [24]) by a factor of \(1/\sqrt{2}\). Five types of trapping sites are included to reproduce the experimental data: two intrinsic traps and three extrinsic traps with sigmoidal distribution (\(f\)) within the damaged layer:
\[
f(x)=\dfrac{1}{1+\exp\left(\dfrac{x-x_0}{\Delta x}\right)},
\]
where \(x_0=2.2\,\mu\textrm{m}\), \(\Delta x=0.154\,\mu\textrm{m}\).
The FESTIM results are compared to the experimental data and the results of MHIMS simulation, both taken from [26].
FESTIM model#
Show code cell source
import festim as F
import fenics as f
import numpy as np
import sympy as sp
import h_transport_materials as htm
################### PARAMETERS ###################
# Exposure conditions
Gamma_atom = 5.8e18
T_exposure = 600
t_exposure = 1e24 / Gamma_atom
t_des = 52 * 3600
final_time = t_exposure + t_des
# Sample
L = 0.8e-3 # half thickness, m
# W properties
rho_W = 6.3e28 # W atomic concentration, m^-3
n_IS = 6 * rho_W # concentration of interstitial sites, m^-3
n_surf = 6.9 * rho_W ** (2 / 3) # concentration of adsorption sites, m^-2
nu0 = 1e13 # attempt frequency, s^-1
SP = 0.19
D_H = htm.diffusivities.filter(material=htm.Tungsten, author="fernandez")[0]
D0 = D_H.pre_exp.magnitude / np.sqrt(2) # diffusivity pre-factor, m^2 s^-1
E_diff = D_H.act_energy.magnitude # diffusion activation energy, eV
lambda_IS = 110e-12 # distance between 2 IS sites, m
sigma_exc = 1.7e-21 # Cross-section for the direct abstraction, m^2
lambda_des = 1 / np.sqrt(n_surf)
# Transitions
E_bs = E_diff # energy barrier from bulk to surface, eV
E_sb = 1.545
E_des = 0.87
################### FUNCTIONS ###############
def k_sb(T, surf_conc, t):
return nu0 * f.exp(-E_sb / F.k_B / T)
def k_bs(T, surf_conc, t):
return nu0 * f.exp(-E_bs / F.k_B / T)
def J_vs_left(T, surf_conc, t):
G_atom = Gamma_atom * f.conditional(t <= t_exposure, 1, 0)
phi_atom = SP * G_atom * (1 - surf_conc / n_surf)
phi_exc = G_atom * sigma_exc * surf_conc
phi_des = 2 * nu0 * (lambda_des * surf_conc) ** 2 * f.exp(-2 * E_des / F.k_B / T)
return phi_atom - phi_exc - phi_des
def J_vs_right(T, surf_conc, t):
phi_des = 2 * nu0 * (lambda_des * surf_conc) ** 2 * f.exp(-2 * E_des / F.k_B / T)
return -phi_des
################### MODEL ###################
W_model = F.Simulation(log_level=40)
# Mesh
vertices = np.concatenate(
[
np.linspace(0, 5e-8, num=100),
np.linspace(5e-8, 5e-6, num=400),
np.linspace(5e-6, L, num=500),
]
)
W_model.mesh = F.MeshFromVertices(vertices)
# Materials
tungsten = F.Material(id=1, D_0=D0, E_D=E_diff)
W_model.materials = tungsten
distr = 1 / (1 + sp.exp((F.x - 2.2e-6) / 1.54e-7))
traps = F.Traps(
[
F.Trap(
k_0=D0 / (n_IS * lambda_IS**2),
E_k=E_diff,
p_0=nu0,
E_p=0.85,
density=1e-4 * rho_W,
materials=tungsten,
),
F.Trap(
k_0=D0 / (n_IS * lambda_IS**2),
E_k=E_diff,
p_0=nu0,
E_p=1.00,
density=1e-4 * rho_W,
materials=tungsten,
),
F.Trap(
k_0=D0 / (n_IS * lambda_IS**2),
E_k=E_diff,
p_0=nu0,
E_p=1.65,
density=0.19e-2 * rho_W * distr,
materials=tungsten,
),
F.Trap(
k_0=D0 / (n_IS * lambda_IS**2),
E_k=E_diff,
p_0=nu0,
E_p=1.85,
density=0.16e-2 * rho_W * distr,
materials=tungsten,
),
F.Trap(
k_0=D0 / (n_IS * lambda_IS**2),
E_k=E_diff,
p_0=nu0,
E_p=2.06,
density=0.02e-2 * rho_W * distr,
materials=tungsten,
),
]
)
W_model.traps = traps
W_model.T = T_exposure
BC_left = F.SurfaceKinetics(
k_sb=k_sb,
k_bs=k_bs,
lambda_IS=lambda_IS,
n_surf=n_surf,
n_IS=n_IS,
J_vs=J_vs_left,
surfaces=1,
initial_condition=0,
t=F.t,
)
BC_right = F.SurfaceKinetics(
k_sb=k_sb,
k_bs=k_bs,
lambda_IS=lambda_IS,
n_surf=n_surf,
n_IS=n_IS,
J_vs=J_vs_right,
surfaces=2,
initial_condition=0,
t=F.t,
)
W_model.boundary_conditions = [BC_left, BC_right]
# Exports
export_fluences = [
5.22e22,
1.25e23,
4.8e23,
6.3e23,
1e24,
]
export_times = [fluence / Gamma_atom for fluence in export_fluences]
export_times += [export_times[-1] + 20 * 3600, export_times[-1] + 52 * 3600]
derived_quantities = F.DerivedQuantities(
[
F.TotalVolume(field="retention", volume=1),
F.AdsorbedHydrogen(surface=1),
F.AdsorbedHydrogen(surface=2),
],
show_units=True,
)
TXT = F.TXTExport(field="retention", filename="./FESTIM_sim.txt", times=export_times)
W_model.exports = [derived_quantities] + [TXT]
W_model.dt = F.Stepsize(
initial_value=1e-5,
stepsize_change_ratio=1.25,
max_stepsize=500,
dt_min=1e-6,
milestones=export_times,
)
W_model.settings = F.Settings(
absolute_tolerance=1e7,
relative_tolerance=1e-13,
final_time=final_time,
traps_element_type="DG",
)
W_model.initialise()
W_model.run()
Show code cell output
Calling FFC just-in-time (JIT) compiler, this may take some time.
Defining initial values
Defining variational problem
Defining source terms
Defining boundary conditions
Time stepping...
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
0.0 % 1.0e-05 s Elapsed time so far: 2.7 s
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Calling FFC just-in-time (JIT) compiler, this may take some time.
0.0 % 2.3e-05 s Elapsed time so far: 4.5 s
0.0 % 3.8e-05 s Elapsed time so far: 4.5 s
0.0 % 5.8e-05 s Elapsed time so far: 4.6 s
0.0 % 8.2e-05 s Elapsed time so far: 4.6 s
0.0 % 1.1e-04 s Elapsed time so far: 4.7 s
0.0 % 1.5e-04 s Elapsed time so far: 4.8 s
0.0 % 2.0e-04 s Elapsed time so far: 4.9 s
0.0 % 2.6e-04 s Elapsed time so far: 5.0 s
0.0 % 3.3e-04 s Elapsed time so far: 5.0 s
0.0 % 4.3e-04 s Elapsed time so far: 5.1 s
0.0 % 5.4e-04 s Elapsed time so far: 5.2 s
0.0 % 6.9e-04 s Elapsed time so far: 5.3 s
0.0 % 8.7e-04 s Elapsed time so far: 5.4 s
0.0 % 1.1e-03 s Elapsed time so far: 5.5 s
0.0 % 1.4e-03 s Elapsed time so far: 5.5 s
0.0 % 1.7e-03 s Elapsed time so far: 5.6 s
0.0 % 2.2e-03 s Elapsed time so far: 5.7 s
0.0 % 2.7e-03 s Elapsed time so far: 5.8 s
0.0 % 3.4e-03 s Elapsed time so far: 5.9 s
0.0 % 4.3e-03 s Elapsed time so far: 6.0 s
0.0 % 5.4e-03 s Elapsed time so far: 6.0 s
0.0 % 6.7e-03 s Elapsed time so far: 6.1 s
0.0 % 8.4e-03 s Elapsed time so far: 6.2 s
0.0 % 1.1e-02 s Elapsed time so far: 6.3 s
0.0 % 1.3e-02 s Elapsed time so far: 6.4 s
0.0 % 1.7e-02 s Elapsed time so far: 6.5 s
0.0 % 2.1e-02 s Elapsed time so far: 6.5 s
0.0 % 2.6e-02 s Elapsed time so far: 6.6 s
0.0 % 3.2e-02 s Elapsed time so far: 6.7 s
0.0 % 4.0e-02 s Elapsed time so far: 6.8 s
0.0 % 5.0e-02 s Elapsed time so far: 6.9 s
0.0 % 6.3e-02 s Elapsed time so far: 7.0 s
0.0 % 7.9e-02 s Elapsed time so far: 7.0 s
0.0 % 9.9e-02 s Elapsed time so far: 7.1 s
0.0 % 1.2e-01 s Elapsed time so far: 7.2 s
0.0 % 1.5e-01 s Elapsed time so far: 7.3 s
0.0 % 1.9e-01 s Elapsed time so far: 7.4 s
0.0 % 2.4e-01 s Elapsed time so far: 7.4 s
0.0 % 3.0e-01 s Elapsed time so far: 7.5 s
0.0 % 3.8e-01 s Elapsed time so far: 7.6 s
0.0 % 4.7e-01 s Elapsed time so far: 7.7 s
0.0 % 5.9e-01 s Elapsed time so far: 7.8 s
0.0 % 7.3e-01 s Elapsed time so far: 7.9 s
0.0 % 9.2e-01 s Elapsed time so far: 8.0 s
0.0 % 1.1e+00 s Elapsed time so far: 8.1 s
0.0 % 1.4e+00 s Elapsed time so far: 8.2 s
0.0 % 1.8e+00 s Elapsed time so far: 8.3 s
0.0 % 2.2e+00 s Elapsed time so far: 8.4 s
0.0 % 2.8e+00 s Elapsed time so far: 8.5 s
0.0 % 3.5e+00 s Elapsed time so far: 8.6 s
0.0 % 4.4e+00 s Elapsed time so far: 8.8 s
0.0 % 5.5e+00 s Elapsed time so far: 8.9 s
0.0 % 6.8e+00 s Elapsed time so far: 9.0 s
0.0 % 8.6e+00 s Elapsed time so far: 9.1 s
0.0 % 1.1e+01 s Elapsed time so far: 9.2 s
0.0 % 1.3e+01 s Elapsed time so far: 9.3 s
0.0 % 1.7e+01 s Elapsed time so far: 9.4 s
0.0 % 2.1e+01 s Elapsed time so far: 9.6 s
0.0 % 2.6e+01 s Elapsed time so far: 9.7 s
0.0 % 3.3e+01 s Elapsed time so far: 9.8 s
0.0 % 4.1e+01 s Elapsed time so far: 10.0 s
0.0 % 5.1e+01 s Elapsed time so far: 10.1 s
0.0 % 6.4e+01 s Elapsed time so far: 10.3 s
0.0 % 8.0e+01 s Elapsed time so far: 10.5 s
0.0 % 1.0e+02 s Elapsed time so far: 10.6 s
0.0 % 1.2e+02 s Elapsed time so far: 10.7 s
0.0 % 1.6e+02 s Elapsed time so far: 10.9 s
0.1 % 1.9e+02 s Elapsed time so far: 11.1 s
0.1 % 2.3e+02 s Elapsed time so far: 11.2 s
0.1 % 2.6e+02 s Elapsed time so far: 11.3 s
0.1 % 3.1e+02 s Elapsed time so far: 11.7 s
0.1 % 3.5e+02 s Elapsed time so far: 12.0 s
0.1 % 4.0e+02 s Elapsed time so far: 12.1 s
0.1 % 4.6e+02 s Elapsed time so far: 12.3 s
0.1 % 5.4e+02 s Elapsed time so far: 12.5 s
0.2 % 6.0e+02 s Elapsed time so far: 12.7 s
0.2 % 6.7e+02 s Elapsed time so far: 12.9 s
0.2 % 7.7e+02 s Elapsed time so far: 13.0 s
0.2 % 8.4e+02 s Elapsed time so far: 13.2 s
0.3 % 9.4e+02 s Elapsed time so far: 13.5 s
0.3 % 1.1e+03 s Elapsed time so far: 13.8 s
0.3 % 1.2e+03 s Elapsed time so far: 14.1 s
0.3 % 1.3e+03 s Elapsed time so far: 14.4 s
0.4 % 1.4e+03 s Elapsed time so far: 14.5 s
0.4 % 1.5e+03 s Elapsed time so far: 14.6 s
0.5 % 1.6e+03 s Elapsed time so far: 14.8 s
0.5 % 1.8e+03 s Elapsed time so far: 15.0 s
0.5 % 1.9e+03 s Elapsed time so far: 15.1 s
0.6 % 2.1e+03 s Elapsed time so far: 15.3 s
0.6 % 2.2e+03 s Elapsed time so far: 15.4 s
0.7 % 2.4e+03 s Elapsed time so far: 15.6 s
0.7 % 2.6e+03 s Elapsed time so far: 15.7 s
0.8 % 2.8e+03 s Elapsed time so far: 15.9 s
0.8 % 3.0e+03 s Elapsed time so far: 16.0 s
0.9 % 3.2e+03 s Elapsed time so far: 16.2 s
0.9 % 3.4e+03 s Elapsed time so far: 16.3 s
1.0 % 3.7e+03 s Elapsed time so far: 16.5 s
1.1 % 3.9e+03 s Elapsed time so far: 16.6 s
1.2 % 4.2e+03 s Elapsed time so far: 16.8 s
1.2 % 4.4e+03 s Elapsed time so far: 16.9 s
1.3 % 4.7e+03 s Elapsed time so far: 17.1 s
1.4 % 5.1e+03 s Elapsed time so far: 17.2 s
1.5 % 5.4e+03 s Elapsed time so far: 17.4 s
1.6 % 5.8e+03 s Elapsed time so far: 17.5 s
1.7 % 6.0e+03 s Elapsed time so far: 17.7 s
1.8 % 6.4e+03 s Elapsed time so far: 17.8 s
1.9 % 6.9e+03 s Elapsed time so far: 18.0 s
2.0 % 7.2e+03 s Elapsed time so far: 18.1 s
2.1 % 7.7e+03 s Elapsed time so far: 18.3 s
2.2 % 8.0e+03 s Elapsed time so far: 18.4 s
2.4 % 8.5e+03 s Elapsed time so far: 18.6 s
2.5 % 8.8e+03 s Elapsed time so far: 18.7 s
2.5 % 9.0e+03 s Elapsed time so far: 18.8 s
Calling FFC just-in-time (JIT) compiler, this may take some time.
2.6 % 9.2e+03 s Elapsed time so far: 19.4 s
2.6 % 9.4e+03 s Elapsed time so far: 19.6 s
2.7 % 9.7e+03 s Elapsed time so far: 19.7 s
2.8 % 1.0e+04 s Elapsed time so far: 19.8 s
2.9 % 1.1e+04 s Elapsed time so far: 20.0 s
3.1 % 1.1e+04 s Elapsed time so far: 20.1 s
3.2 % 1.2e+04 s Elapsed time so far: 20.3 s
3.4 % 1.2e+04 s Elapsed time so far: 20.4 s
3.5 % 1.3e+04 s Elapsed time so far: 20.5 s
3.6 % 1.3e+04 s Elapsed time so far: 20.7 s
3.8 % 1.4e+04 s Elapsed time so far: 20.8 s
3.9 % 1.4e+04 s Elapsed time so far: 20.9 s
4.0 % 1.5e+04 s Elapsed time so far: 21.1 s
4.2 % 1.5e+04 s Elapsed time so far: 21.2 s
4.3 % 1.6e+04 s Elapsed time so far: 21.4 s
4.5 % 1.6e+04 s Elapsed time so far: 21.5 s
4.6 % 1.7e+04 s Elapsed time so far: 21.6 s
4.7 % 1.7e+04 s Elapsed time so far: 21.8 s
4.9 % 1.8e+04 s Elapsed time so far: 21.9 s
5.0 % 1.8e+04 s Elapsed time so far: 22.0 s
5.2 % 1.9e+04 s Elapsed time so far: 22.2 s
5.3 % 1.9e+04 s Elapsed time so far: 22.3 s
5.4 % 2.0e+04 s Elapsed time so far: 22.5 s
5.6 % 2.0e+04 s Elapsed time so far: 22.6 s
5.7 % 2.1e+04 s Elapsed time so far: 22.7 s
5.9 % 2.1e+04 s Elapsed time so far: 22.9 s
6.0 % 2.2e+04 s Elapsed time so far: 23.0 s
6.1 % 2.2e+04 s Elapsed time so far: 23.2 s
6.3 % 2.3e+04 s Elapsed time so far: 23.3 s
6.4 % 2.3e+04 s Elapsed time so far: 23.4 s
6.5 % 2.4e+04 s Elapsed time so far: 23.6 s
6.7 % 2.4e+04 s Elapsed time so far: 23.7 s
6.8 % 2.5e+04 s Elapsed time so far: 23.8 s
7.0 % 2.5e+04 s Elapsed time so far: 24.0 s
7.1 % 2.6e+04 s Elapsed time so far: 24.1 s
7.2 % 2.6e+04 s Elapsed time so far: 24.2 s
7.4 % 2.7e+04 s Elapsed time so far: 24.4 s
7.5 % 2.7e+04 s Elapsed time so far: 24.5 s
7.7 % 2.8e+04 s Elapsed time so far: 24.6 s
7.8 % 2.8e+04 s Elapsed time so far: 24.8 s
7.9 % 2.9e+04 s Elapsed time so far: 24.9 s
8.1 % 2.9e+04 s Elapsed time so far: 25.1 s
8.2 % 3.0e+04 s Elapsed time so far: 25.2 s
8.4 % 3.0e+04 s Elapsed time so far: 25.3 s
8.5 % 3.1e+04 s Elapsed time so far: 25.5 s
8.6 % 3.1e+04 s Elapsed time so far: 25.6 s
8.8 % 3.2e+04 s Elapsed time so far: 25.8 s
8.9 % 3.2e+04 s Elapsed time so far: 25.9 s
9.1 % 3.3e+04 s Elapsed time so far: 26.0 s
9.2 % 3.3e+04 s Elapsed time so far: 26.2 s
9.3 % 3.4e+04 s Elapsed time so far: 26.3 s
9.5 % 3.4e+04 s Elapsed time so far: 26.4 s
9.6 % 3.5e+04 s Elapsed time so far: 26.6 s
9.8 % 3.5e+04 s Elapsed time so far: 26.7 s
9.9 % 3.6e+04 s Elapsed time so far: 26.8 s
10.0 % 3.6e+04 s Elapsed time so far: 27.0 s
10.2 % 3.7e+04 s Elapsed time so far: 27.1 s
10.3 % 3.7e+04 s Elapsed time so far: 27.2 s
10.4 % 3.8e+04 s Elapsed time so far: 27.4 s
10.6 % 3.8e+04 s Elapsed time so far: 27.5 s
10.7 % 3.9e+04 s Elapsed time so far: 27.7 s
10.9 % 3.9e+04 s Elapsed time so far: 27.8 s
11.0 % 4.0e+04 s Elapsed time so far: 28.0 s
11.1 % 4.0e+04 s Elapsed time so far: 28.1 s
11.3 % 4.1e+04 s Elapsed time so far: 28.2 s
11.4 % 4.1e+04 s Elapsed time so far: 28.4 s
11.6 % 4.2e+04 s Elapsed time so far: 28.5 s
11.7 % 4.2e+04 s Elapsed time so far: 28.6 s
11.8 % 4.3e+04 s Elapsed time so far: 28.8 s
12.0 % 4.3e+04 s Elapsed time so far: 28.9 s
12.1 % 4.4e+04 s Elapsed time so far: 29.0 s
12.2 % 4.4e+04 s Elapsed time so far: 29.2 s
12.4 % 4.5e+04 s Elapsed time so far: 29.3 s
12.5 % 4.5e+04 s Elapsed time so far: 29.5 s
12.7 % 4.6e+04 s Elapsed time so far: 29.6 s
12.8 % 4.6e+04 s Elapsed time so far: 29.7 s
12.9 % 4.7e+04 s Elapsed time so far: 29.9 s
13.1 % 4.7e+04 s Elapsed time so far: 30.0 s
13.2 % 4.8e+04 s Elapsed time so far: 30.1 s
13.4 % 4.8e+04 s Elapsed time so far: 30.3 s
13.5 % 4.9e+04 s Elapsed time so far: 30.4 s
13.6 % 4.9e+04 s Elapsed time so far: 30.6 s
13.8 % 5.0e+04 s Elapsed time so far: 30.7 s
13.9 % 5.0e+04 s Elapsed time so far: 30.8 s
14.1 % 5.1e+04 s Elapsed time so far: 31.0 s
14.2 % 5.1e+04 s Elapsed time so far: 31.1 s
14.3 % 5.2e+04 s Elapsed time so far: 31.2 s
14.5 % 5.2e+04 s Elapsed time so far: 31.4 s
14.6 % 5.3e+04 s Elapsed time so far: 31.5 s
14.8 % 5.3e+04 s Elapsed time so far: 31.6 s
14.9 % 5.4e+04 s Elapsed time so far: 31.7 s
15.0 % 5.4e+04 s Elapsed time so far: 31.8 s
15.2 % 5.5e+04 s Elapsed time so far: 32.0 s
15.3 % 5.5e+04 s Elapsed time so far: 32.1 s
15.4 % 5.6e+04 s Elapsed time so far: 32.2 s
15.6 % 5.6e+04 s Elapsed time so far: 32.3 s
15.7 % 5.7e+04 s Elapsed time so far: 32.4 s
15.9 % 5.7e+04 s Elapsed time so far: 32.5 s
16.0 % 5.8e+04 s Elapsed time so far: 32.6 s
16.1 % 5.8e+04 s Elapsed time so far: 32.7 s
16.3 % 5.9e+04 s Elapsed time so far: 32.8 s
16.4 % 5.9e+04 s Elapsed time so far: 32.9 s
16.6 % 6.0e+04 s Elapsed time so far: 33.0 s
16.7 % 6.0e+04 s Elapsed time so far: 33.1 s
16.8 % 6.1e+04 s Elapsed time so far: 33.2 s
17.0 % 6.1e+04 s Elapsed time so far: 33.4 s
17.1 % 6.2e+04 s Elapsed time so far: 33.5 s
17.3 % 6.2e+04 s Elapsed time so far: 33.6 s
17.4 % 6.3e+04 s Elapsed time so far: 33.7 s
17.5 % 6.3e+04 s Elapsed time so far: 33.8 s
17.7 % 6.4e+04 s Elapsed time so far: 33.9 s
17.8 % 6.4e+04 s Elapsed time so far: 34.0 s
17.9 % 6.5e+04 s Elapsed time so far: 34.1 s
18.1 % 6.5e+04 s Elapsed time so far: 34.2 s
18.2 % 6.6e+04 s Elapsed time so far: 34.3 s
18.4 % 6.6e+04 s Elapsed time so far: 34.5 s
18.5 % 6.7e+04 s Elapsed time so far: 34.6 s
18.6 % 6.7e+04 s Elapsed time so far: 34.7 s
18.8 % 6.8e+04 s Elapsed time so far: 34.8 s
18.9 % 6.8e+04 s Elapsed time so far: 34.9 s
19.1 % 6.9e+04 s Elapsed time so far: 35.0 s
19.2 % 6.9e+04 s Elapsed time so far: 35.1 s
19.3 % 7.0e+04 s Elapsed time so far: 35.2 s
19.5 % 7.0e+04 s Elapsed time so far: 35.3 s
19.6 % 7.1e+04 s Elapsed time so far: 35.4 s
19.8 % 7.1e+04 s Elapsed time so far: 35.5 s
19.9 % 7.2e+04 s Elapsed time so far: 35.7 s
20.0 % 7.2e+04 s Elapsed time so far: 35.8 s
20.2 % 7.3e+04 s Elapsed time so far: 35.9 s
20.3 % 7.3e+04 s Elapsed time so far: 36.0 s
20.4 % 7.4e+04 s Elapsed time so far: 36.1 s
20.6 % 7.4e+04 s Elapsed time so far: 36.2 s
20.7 % 7.5e+04 s Elapsed time so far: 36.3 s
20.9 % 7.5e+04 s Elapsed time so far: 36.4 s
21.0 % 7.6e+04 s Elapsed time so far: 36.5 s
21.1 % 7.6e+04 s Elapsed time so far: 36.6 s
21.3 % 7.7e+04 s Elapsed time so far: 36.7 s
21.4 % 7.7e+04 s Elapsed time so far: 36.9 s
21.6 % 7.8e+04 s Elapsed time so far: 37.0 s
21.7 % 7.8e+04 s Elapsed time so far: 37.1 s
21.8 % 7.9e+04 s Elapsed time so far: 37.2 s
22.0 % 7.9e+04 s Elapsed time so far: 37.3 s
22.1 % 8.0e+04 s Elapsed time so far: 37.4 s
22.3 % 8.0e+04 s Elapsed time so far: 37.5 s
22.4 % 8.1e+04 s Elapsed time so far: 37.6 s
22.5 % 8.1e+04 s Elapsed time so far: 37.7 s
22.7 % 8.2e+04 s Elapsed time so far: 37.8 s
22.8 % 8.2e+04 s Elapsed time so far: 38.0 s
23.0 % 8.3e+04 s Elapsed time so far: 38.1 s
23.0 % 8.3e+04 s Elapsed time so far: 38.2 s
23.1 % 8.3e+04 s Elapsed time so far: 38.3 s
23.2 % 8.3e+04 s Elapsed time so far: 38.4 s
23.3 % 8.4e+04 s Elapsed time so far: 38.5 s
23.4 % 8.4e+04 s Elapsed time so far: 38.6 s
23.6 % 8.5e+04 s Elapsed time so far: 38.7 s
23.7 % 8.5e+04 s Elapsed time so far: 38.8 s
23.8 % 8.6e+04 s Elapsed time so far: 38.9 s
24.0 % 8.6e+04 s Elapsed time so far: 39.1 s
24.1 % 8.7e+04 s Elapsed time so far: 39.2 s
24.3 % 8.7e+04 s Elapsed time so far: 39.3 s
24.4 % 8.8e+04 s Elapsed time so far: 39.4 s
24.5 % 8.8e+04 s Elapsed time so far: 39.5 s
24.7 % 8.9e+04 s Elapsed time so far: 39.6 s
24.8 % 8.9e+04 s Elapsed time so far: 39.7 s
25.0 % 9.0e+04 s Elapsed time so far: 39.8 s
25.1 % 9.0e+04 s Elapsed time so far: 39.9 s
25.2 % 9.1e+04 s Elapsed time so far: 40.0 s
25.4 % 9.1e+04 s Elapsed time so far: 40.1 s
25.5 % 9.2e+04 s Elapsed time so far: 40.2 s
25.6 % 9.2e+04 s Elapsed time so far: 40.4 s
25.8 % 9.3e+04 s Elapsed time so far: 40.5 s
25.9 % 9.3e+04 s Elapsed time so far: 40.6 s
26.1 % 9.4e+04 s Elapsed time so far: 40.7 s
26.2 % 9.4e+04 s Elapsed time so far: 40.8 s
26.4 % 9.5e+04 s Elapsed time so far: 40.9 s
26.5 % 9.5e+04 s Elapsed time so far: 41.0 s
26.6 % 9.6e+04 s Elapsed time so far: 41.1 s
26.8 % 9.6e+04 s Elapsed time so far: 41.2 s
26.9 % 9.7e+04 s Elapsed time so far: 41.3 s
27.0 % 9.7e+04 s Elapsed time so far: 41.4 s
27.2 % 9.8e+04 s Elapsed time so far: 41.5 s
27.3 % 9.8e+04 s Elapsed time so far: 41.7 s
27.5 % 9.9e+04 s Elapsed time so far: 41.8 s
27.6 % 9.9e+04 s Elapsed time so far: 41.9 s
27.7 % 1.0e+05 s Elapsed time so far: 42.0 s
27.9 % 1.0e+05 s Elapsed time so far: 42.1 s
28.0 % 1.0e+05 s Elapsed time so far: 42.2 s
28.1 % 1.0e+05 s Elapsed time so far: 42.3 s
28.3 % 1.0e+05 s Elapsed time so far: 42.4 s
28.4 % 1.0e+05 s Elapsed time so far: 42.5 s
28.6 % 1.0e+05 s Elapsed time so far: 42.6 s
28.7 % 1.0e+05 s Elapsed time so far: 42.7 s
28.9 % 1.0e+05 s Elapsed time so far: 42.9 s
29.0 % 1.0e+05 s Elapsed time so far: 43.0 s
29.1 % 1.0e+05 s Elapsed time so far: 43.1 s
29.3 % 1.1e+05 s Elapsed time so far: 43.2 s
29.4 % 1.1e+05 s Elapsed time so far: 43.3 s
29.5 % 1.1e+05 s Elapsed time so far: 43.4 s
29.7 % 1.1e+05 s Elapsed time so far: 43.5 s
29.8 % 1.1e+05 s Elapsed time so far: 43.6 s
30.0 % 1.1e+05 s Elapsed time so far: 43.7 s
30.1 % 1.1e+05 s Elapsed time so far: 43.8 s
30.2 % 1.1e+05 s Elapsed time so far: 43.9 s
30.3 % 1.1e+05 s Elapsed time so far: 44.1 s
30.5 % 1.1e+05 s Elapsed time so far: 44.2 s
30.6 % 1.1e+05 s Elapsed time so far: 44.3 s
30.8 % 1.1e+05 s Elapsed time so far: 44.4 s
30.9 % 1.1e+05 s Elapsed time so far: 44.5 s
31.0 % 1.1e+05 s Elapsed time so far: 44.6 s
31.2 % 1.1e+05 s Elapsed time so far: 44.7 s
31.3 % 1.1e+05 s Elapsed time so far: 44.8 s
31.4 % 1.1e+05 s Elapsed time so far: 44.9 s
31.6 % 1.1e+05 s Elapsed time so far: 45.0 s
31.7 % 1.1e+05 s Elapsed time so far: 45.1 s
31.9 % 1.1e+05 s Elapsed time so far: 45.2 s
32.0 % 1.2e+05 s Elapsed time so far: 45.4 s
32.1 % 1.2e+05 s Elapsed time so far: 45.5 s
32.3 % 1.2e+05 s Elapsed time so far: 45.6 s
32.4 % 1.2e+05 s Elapsed time so far: 45.7 s
32.6 % 1.2e+05 s Elapsed time so far: 45.8 s
32.7 % 1.2e+05 s Elapsed time so far: 45.9 s
32.8 % 1.2e+05 s Elapsed time so far: 46.0 s
33.0 % 1.2e+05 s Elapsed time so far: 46.1 s
33.1 % 1.2e+05 s Elapsed time so far: 46.2 s
33.3 % 1.2e+05 s Elapsed time so far: 46.3 s
33.4 % 1.2e+05 s Elapsed time so far: 46.4 s
33.5 % 1.2e+05 s Elapsed time so far: 46.6 s
33.7 % 1.2e+05 s Elapsed time so far: 46.7 s
33.8 % 1.2e+05 s Elapsed time so far: 46.8 s
34.0 % 1.2e+05 s Elapsed time so far: 46.9 s
34.1 % 1.2e+05 s Elapsed time so far: 47.0 s
34.2 % 1.2e+05 s Elapsed time so far: 47.1 s
34.4 % 1.2e+05 s Elapsed time so far: 47.2 s
34.5 % 1.2e+05 s Elapsed time so far: 47.3 s
34.6 % 1.2e+05 s Elapsed time so far: 47.4 s
34.8 % 1.3e+05 s Elapsed time so far: 47.5 s
34.9 % 1.3e+05 s Elapsed time so far: 47.6 s
35.1 % 1.3e+05 s Elapsed time so far: 47.7 s
35.2 % 1.3e+05 s Elapsed time so far: 47.9 s
35.3 % 1.3e+05 s Elapsed time so far: 48.0 s
35.5 % 1.3e+05 s Elapsed time so far: 48.1 s
35.6 % 1.3e+05 s Elapsed time so far: 48.2 s
35.8 % 1.3e+05 s Elapsed time so far: 48.3 s
35.9 % 1.3e+05 s Elapsed time so far: 48.4 s
36.0 % 1.3e+05 s Elapsed time so far: 48.5 s
36.2 % 1.3e+05 s Elapsed time so far: 48.6 s
36.3 % 1.3e+05 s Elapsed time so far: 48.7 s
36.5 % 1.3e+05 s Elapsed time so far: 48.8 s
36.6 % 1.3e+05 s Elapsed time so far: 48.9 s
36.7 % 1.3e+05 s Elapsed time so far: 49.0 s
36.9 % 1.3e+05 s Elapsed time so far: 49.1 s
37.0 % 1.3e+05 s Elapsed time so far: 49.2 s
37.1 % 1.3e+05 s Elapsed time so far: 49.4 s
37.3 % 1.3e+05 s Elapsed time so far: 49.5 s
37.4 % 1.3e+05 s Elapsed time so far: 49.6 s
37.6 % 1.4e+05 s Elapsed time so far: 49.7 s
37.7 % 1.4e+05 s Elapsed time so far: 49.8 s
37.8 % 1.4e+05 s Elapsed time so far: 49.9 s
38.0 % 1.4e+05 s Elapsed time so far: 50.0 s
38.1 % 1.4e+05 s Elapsed time so far: 50.1 s
38.3 % 1.4e+05 s Elapsed time so far: 50.2 s
38.4 % 1.4e+05 s Elapsed time so far: 50.3 s
38.5 % 1.4e+05 s Elapsed time so far: 50.4 s
38.7 % 1.4e+05 s Elapsed time so far: 50.5 s
38.8 % 1.4e+05 s Elapsed time so far: 50.7 s
39.0 % 1.4e+05 s Elapsed time so far: 50.8 s
39.1 % 1.4e+05 s Elapsed time so far: 50.9 s
39.2 % 1.4e+05 s Elapsed time so far: 51.0 s
39.4 % 1.4e+05 s Elapsed time so far: 51.1 s
39.5 % 1.4e+05 s Elapsed time so far: 51.2 s
39.6 % 1.4e+05 s Elapsed time so far: 51.3 s
39.8 % 1.4e+05 s Elapsed time so far: 51.4 s
39.9 % 1.4e+05 s Elapsed time so far: 51.5 s
40.1 % 1.4e+05 s Elapsed time so far: 51.6 s
40.2 % 1.4e+05 s Elapsed time so far: 51.7 s
40.4 % 1.5e+05 s Elapsed time so far: 51.9 s
40.5 % 1.5e+05 s Elapsed time so far: 52.0 s
40.6 % 1.5e+05 s Elapsed time so far: 52.1 s
40.8 % 1.5e+05 s Elapsed time so far: 52.2 s
40.9 % 1.5e+05 s Elapsed time so far: 52.3 s
41.0 % 1.5e+05 s Elapsed time so far: 52.4 s
41.2 % 1.5e+05 s Elapsed time so far: 52.5 s
41.3 % 1.5e+05 s Elapsed time so far: 52.6 s
41.5 % 1.5e+05 s Elapsed time so far: 52.7 s
41.6 % 1.5e+05 s Elapsed time so far: 52.8 s
41.7 % 1.5e+05 s Elapsed time so far: 52.9 s
41.9 % 1.5e+05 s Elapsed time so far: 53.0 s
42.0 % 1.5e+05 s Elapsed time so far: 53.2 s
42.1 % 1.5e+05 s Elapsed time so far: 53.3 s
42.3 % 1.5e+05 s Elapsed time so far: 53.4 s
42.4 % 1.5e+05 s Elapsed time so far: 53.5 s
42.6 % 1.5e+05 s Elapsed time so far: 53.6 s
42.7 % 1.5e+05 s Elapsed time so far: 53.7 s
42.9 % 1.5e+05 s Elapsed time so far: 53.8 s
43.0 % 1.5e+05 s Elapsed time so far: 53.9 s
43.1 % 1.6e+05 s Elapsed time so far: 54.0 s
43.3 % 1.6e+05 s Elapsed time so far: 54.1 s
43.4 % 1.6e+05 s Elapsed time so far: 54.2 s
43.5 % 1.6e+05 s Elapsed time so far: 54.4 s
43.7 % 1.6e+05 s Elapsed time so far: 54.5 s
43.8 % 1.6e+05 s Elapsed time so far: 54.6 s
44.0 % 1.6e+05 s Elapsed time so far: 54.7 s
44.1 % 1.6e+05 s Elapsed time so far: 54.8 s
44.2 % 1.6e+05 s Elapsed time so far: 54.9 s
44.4 % 1.6e+05 s Elapsed time so far: 55.0 s
44.5 % 1.6e+05 s Elapsed time so far: 55.1 s
44.7 % 1.6e+05 s Elapsed time so far: 55.2 s
44.8 % 1.6e+05 s Elapsed time so far: 55.4 s
44.9 % 1.6e+05 s Elapsed time so far: 55.5 s
45.1 % 1.6e+05 s Elapsed time so far: 55.6 s
45.2 % 1.6e+05 s Elapsed time so far: 55.7 s
45.4 % 1.6e+05 s Elapsed time so far: 55.8 s
45.5 % 1.6e+05 s Elapsed time so far: 55.9 s
45.6 % 1.6e+05 s Elapsed time so far: 56.0 s
45.8 % 1.6e+05 s Elapsed time so far: 56.1 s
45.9 % 1.7e+05 s Elapsed time so far: 56.2 s
46.0 % 1.7e+05 s Elapsed time so far: 56.3 s
46.2 % 1.7e+05 s Elapsed time so far: 56.5 s
46.3 % 1.7e+05 s Elapsed time so far: 56.6 s
46.5 % 1.7e+05 s Elapsed time so far: 56.7 s
46.6 % 1.7e+05 s Elapsed time so far: 56.8 s
46.7 % 1.7e+05 s Elapsed time so far: 56.9 s
46.9 % 1.7e+05 s Elapsed time so far: 57.0 s
47.0 % 1.7e+05 s Elapsed time so far: 57.1 s
47.2 % 1.7e+05 s Elapsed time so far: 57.2 s
47.3 % 1.7e+05 s Elapsed time so far: 57.3 s
47.4 % 1.7e+05 s Elapsed time so far: 57.5 s
47.6 % 1.7e+05 s Elapsed time so far: 57.6 s
47.7 % 1.7e+05 s Elapsed time so far: 57.7 s
47.9 % 1.7e+05 s Elapsed time so far: 57.8 s
47.9 % 1.7e+05 s Elapsed time so far: 57.9 s
48.1 % 1.7e+05 s Elapsed time so far: 58.1 s
48.1 % 1.7e+05 s Elapsed time so far: 58.2 s
48.3 % 1.7e+05 s Elapsed time so far: 58.4 s
48.4 % 1.7e+05 s Elapsed time so far: 58.5 s
48.5 % 1.7e+05 s Elapsed time so far: 58.6 s
48.7 % 1.8e+05 s Elapsed time so far: 58.7 s
48.8 % 1.8e+05 s Elapsed time so far: 58.8 s
49.0 % 1.8e+05 s Elapsed time so far: 59.0 s
49.1 % 1.8e+05 s Elapsed time so far: 59.1 s
49.2 % 1.8e+05 s Elapsed time so far: 59.2 s
49.4 % 1.8e+05 s Elapsed time so far: 59.3 s
49.5 % 1.8e+05 s Elapsed time so far: 59.4 s
49.6 % 1.8e+05 s Elapsed time so far: 59.5 s
49.8 % 1.8e+05 s Elapsed time so far: 59.6 s
49.9 % 1.8e+05 s Elapsed time so far: 59.7 s
50.1 % 1.8e+05 s Elapsed time so far: 59.9 s
50.2 % 1.8e+05 s Elapsed time so far: 60.0 s
50.3 % 1.8e+05 s Elapsed time so far: 60.1 s
50.5 % 1.8e+05 s Elapsed time so far: 60.2 s
50.6 % 1.8e+05 s Elapsed time so far: 60.3 s
50.8 % 1.8e+05 s Elapsed time so far: 60.4 s
50.9 % 1.8e+05 s Elapsed time so far: 60.5 s
51.0 % 1.8e+05 s Elapsed time so far: 60.7 s
51.2 % 1.8e+05 s Elapsed time so far: 60.8 s
51.3 % 1.8e+05 s Elapsed time so far: 60.9 s
51.5 % 1.9e+05 s Elapsed time so far: 61.1 s
51.6 % 1.9e+05 s Elapsed time so far: 61.2 s
51.7 % 1.9e+05 s Elapsed time so far: 61.3 s
51.9 % 1.9e+05 s Elapsed time so far: 61.4 s
52.0 % 1.9e+05 s Elapsed time so far: 61.5 s
52.1 % 1.9e+05 s Elapsed time so far: 61.6 s
52.3 % 1.9e+05 s Elapsed time so far: 61.7 s
52.4 % 1.9e+05 s Elapsed time so far: 61.8 s
52.6 % 1.9e+05 s Elapsed time so far: 61.9 s
52.7 % 1.9e+05 s Elapsed time so far: 62.1 s
52.9 % 1.9e+05 s Elapsed time so far: 62.2 s
53.0 % 1.9e+05 s Elapsed time so far: 62.3 s
53.1 % 1.9e+05 s Elapsed time so far: 62.4 s
53.3 % 1.9e+05 s Elapsed time so far: 62.5 s
53.4 % 1.9e+05 s Elapsed time so far: 62.6 s
53.5 % 1.9e+05 s Elapsed time so far: 62.7 s
53.7 % 1.9e+05 s Elapsed time so far: 62.8 s
53.8 % 1.9e+05 s Elapsed time so far: 62.9 s
54.0 % 1.9e+05 s Elapsed time so far: 63.0 s
54.1 % 1.9e+05 s Elapsed time so far: 63.1 s
54.2 % 2.0e+05 s Elapsed time so far: 63.2 s
54.4 % 2.0e+05 s Elapsed time so far: 63.4 s
54.5 % 2.0e+05 s Elapsed time so far: 63.5 s
54.6 % 2.0e+05 s Elapsed time so far: 63.6 s
54.8 % 2.0e+05 s Elapsed time so far: 63.7 s
54.9 % 2.0e+05 s Elapsed time so far: 63.8 s
55.1 % 2.0e+05 s Elapsed time so far: 63.9 s
55.2 % 2.0e+05 s Elapsed time so far: 64.0 s
55.4 % 2.0e+05 s Elapsed time so far: 64.1 s
55.5 % 2.0e+05 s Elapsed time so far: 64.2 s
55.6 % 2.0e+05 s Elapsed time so far: 64.3 s
55.8 % 2.0e+05 s Elapsed time so far: 64.4 s
55.9 % 2.0e+05 s Elapsed time so far: 64.5 s
56.0 % 2.0e+05 s Elapsed time so far: 64.7 s
56.2 % 2.0e+05 s Elapsed time so far: 64.8 s
56.3 % 2.0e+05 s Elapsed time so far: 64.9 s
56.5 % 2.0e+05 s Elapsed time so far: 65.0 s
56.6 % 2.0e+05 s Elapsed time so far: 65.1 s
56.7 % 2.0e+05 s Elapsed time so far: 65.2 s
56.9 % 2.0e+05 s Elapsed time so far: 65.3 s
57.0 % 2.1e+05 s Elapsed time so far: 65.4 s
57.2 % 2.1e+05 s Elapsed time so far: 65.5 s
57.3 % 2.1e+05 s Elapsed time so far: 65.7 s
57.4 % 2.1e+05 s Elapsed time so far: 65.8 s
57.6 % 2.1e+05 s Elapsed time so far: 65.9 s
57.7 % 2.1e+05 s Elapsed time so far: 66.0 s
57.9 % 2.1e+05 s Elapsed time so far: 66.1 s
58.0 % 2.1e+05 s Elapsed time so far: 66.2 s
58.1 % 2.1e+05 s Elapsed time so far: 66.3 s
58.3 % 2.1e+05 s Elapsed time so far: 66.4 s
58.4 % 2.1e+05 s Elapsed time so far: 66.5 s
58.5 % 2.1e+05 s Elapsed time so far: 66.6 s
58.7 % 2.1e+05 s Elapsed time so far: 66.7 s
58.8 % 2.1e+05 s Elapsed time so far: 66.8 s
59.0 % 2.1e+05 s Elapsed time so far: 66.9 s
59.1 % 2.1e+05 s Elapsed time so far: 67.0 s
59.2 % 2.1e+05 s Elapsed time so far: 67.2 s
59.4 % 2.1e+05 s Elapsed time so far: 67.3 s
59.5 % 2.1e+05 s Elapsed time so far: 67.4 s
59.7 % 2.1e+05 s Elapsed time so far: 67.5 s
59.8 % 2.2e+05 s Elapsed time so far: 67.6 s
59.9 % 2.2e+05 s Elapsed time so far: 67.7 s
60.1 % 2.2e+05 s Elapsed time so far: 67.8 s
60.2 % 2.2e+05 s Elapsed time so far: 67.9 s
60.4 % 2.2e+05 s Elapsed time so far: 68.0 s
60.5 % 2.2e+05 s Elapsed time so far: 68.1 s
60.6 % 2.2e+05 s Elapsed time so far: 68.2 s
60.8 % 2.2e+05 s Elapsed time so far: 68.3 s
60.9 % 2.2e+05 s Elapsed time so far: 68.4 s
61.0 % 2.2e+05 s Elapsed time so far: 68.4 s
61.2 % 2.2e+05 s Elapsed time so far: 68.5 s
61.3 % 2.2e+05 s Elapsed time so far: 68.6 s
61.5 % 2.2e+05 s Elapsed time so far: 68.7 s
61.6 % 2.2e+05 s Elapsed time so far: 68.8 s
61.8 % 2.2e+05 s Elapsed time so far: 68.8 s
61.9 % 2.2e+05 s Elapsed time so far: 68.9 s
62.0 % 2.2e+05 s Elapsed time so far: 69.0 s
62.2 % 2.2e+05 s Elapsed time so far: 69.1 s
62.3 % 2.2e+05 s Elapsed time so far: 69.2 s
62.4 % 2.2e+05 s Elapsed time so far: 69.2 s
62.6 % 2.3e+05 s Elapsed time so far: 69.3 s
62.7 % 2.3e+05 s Elapsed time so far: 69.4 s
62.9 % 2.3e+05 s Elapsed time so far: 69.5 s
63.0 % 2.3e+05 s Elapsed time so far: 69.6 s
63.1 % 2.3e+05 s Elapsed time so far: 69.7 s
63.3 % 2.3e+05 s Elapsed time so far: 69.8 s
63.4 % 2.3e+05 s Elapsed time so far: 69.8 s
63.5 % 2.3e+05 s Elapsed time so far: 69.9 s
63.7 % 2.3e+05 s Elapsed time so far: 70.0 s
63.8 % 2.3e+05 s Elapsed time so far: 70.1 s
64.0 % 2.3e+05 s Elapsed time so far: 70.2 s
64.1 % 2.3e+05 s Elapsed time so far: 70.3 s
64.2 % 2.3e+05 s Elapsed time so far: 70.3 s
64.4 % 2.3e+05 s Elapsed time so far: 70.4 s
64.5 % 2.3e+05 s Elapsed time so far: 70.5 s
64.7 % 2.3e+05 s Elapsed time so far: 70.6 s
64.8 % 2.3e+05 s Elapsed time so far: 70.7 s
64.9 % 2.3e+05 s Elapsed time so far: 70.7 s
65.1 % 2.3e+05 s Elapsed time so far: 70.8 s
65.2 % 2.3e+05 s Elapsed time so far: 70.9 s
65.4 % 2.4e+05 s Elapsed time so far: 71.0 s
65.5 % 2.4e+05 s Elapsed time so far: 71.1 s
65.6 % 2.4e+05 s Elapsed time so far: 71.2 s
65.8 % 2.4e+05 s Elapsed time so far: 71.2 s
65.9 % 2.4e+05 s Elapsed time so far: 71.3 s
66.1 % 2.4e+05 s Elapsed time so far: 71.4 s
66.2 % 2.4e+05 s Elapsed time so far: 71.5 s
66.3 % 2.4e+05 s Elapsed time so far: 71.6 s
66.5 % 2.4e+05 s Elapsed time so far: 71.7 s
66.6 % 2.4e+05 s Elapsed time so far: 71.7 s
66.8 % 2.4e+05 s Elapsed time so far: 71.8 s
66.9 % 2.4e+05 s Elapsed time so far: 71.9 s
67.0 % 2.4e+05 s Elapsed time so far: 72.0 s
67.2 % 2.4e+05 s Elapsed time so far: 72.1 s
67.3 % 2.4e+05 s Elapsed time so far: 72.2 s
67.5 % 2.4e+05 s Elapsed time so far: 72.2 s
67.6 % 2.4e+05 s Elapsed time so far: 72.3 s
67.7 % 2.4e+05 s Elapsed time so far: 72.4 s
67.9 % 2.4e+05 s Elapsed time so far: 72.5 s
68.0 % 2.4e+05 s Elapsed time so far: 72.6 s
68.1 % 2.4e+05 s Elapsed time so far: 72.7 s
68.2 % 2.5e+05 s Elapsed time so far: 72.7 s
68.4 % 2.5e+05 s Elapsed time so far: 72.8 s
68.5 % 2.5e+05 s Elapsed time so far: 72.9 s
68.7 % 2.5e+05 s Elapsed time so far: 73.0 s
68.8 % 2.5e+05 s Elapsed time so far: 73.1 s
68.9 % 2.5e+05 s Elapsed time so far: 73.1 s
69.1 % 2.5e+05 s Elapsed time so far: 73.2 s
69.2 % 2.5e+05 s Elapsed time so far: 73.3 s
69.3 % 2.5e+05 s Elapsed time so far: 73.4 s
69.5 % 2.5e+05 s Elapsed time so far: 73.5 s
69.6 % 2.5e+05 s Elapsed time so far: 73.6 s
69.8 % 2.5e+05 s Elapsed time so far: 73.6 s
69.9 % 2.5e+05 s Elapsed time so far: 73.7 s
70.0 % 2.5e+05 s Elapsed time so far: 73.8 s
70.2 % 2.5e+05 s Elapsed time so far: 73.9 s
70.3 % 2.5e+05 s Elapsed time so far: 74.0 s
70.5 % 2.5e+05 s Elapsed time so far: 74.0 s
70.6 % 2.5e+05 s Elapsed time so far: 74.1 s
70.7 % 2.5e+05 s Elapsed time so far: 74.2 s
70.9 % 2.5e+05 s Elapsed time so far: 74.3 s
71.0 % 2.6e+05 s Elapsed time so far: 74.4 s
71.2 % 2.6e+05 s Elapsed time so far: 74.5 s
71.3 % 2.6e+05 s Elapsed time so far: 74.5 s
71.4 % 2.6e+05 s Elapsed time so far: 74.6 s
71.6 % 2.6e+05 s Elapsed time so far: 74.7 s
71.7 % 2.6e+05 s Elapsed time so far: 74.8 s
71.8 % 2.6e+05 s Elapsed time so far: 74.9 s
72.0 % 2.6e+05 s Elapsed time so far: 74.9 s
72.1 % 2.6e+05 s Elapsed time so far: 75.0 s
72.3 % 2.6e+05 s Elapsed time so far: 75.1 s
72.4 % 2.6e+05 s Elapsed time so far: 75.2 s
72.5 % 2.6e+05 s Elapsed time so far: 75.3 s
72.7 % 2.6e+05 s Elapsed time so far: 75.4 s
72.8 % 2.6e+05 s Elapsed time so far: 75.4 s
73.0 % 2.6e+05 s Elapsed time so far: 75.5 s
73.1 % 2.6e+05 s Elapsed time so far: 75.6 s
73.2 % 2.6e+05 s Elapsed time so far: 75.7 s
73.4 % 2.6e+05 s Elapsed time so far: 75.8 s
73.5 % 2.6e+05 s Elapsed time so far: 75.9 s
73.7 % 2.6e+05 s Elapsed time so far: 75.9 s
73.8 % 2.7e+05 s Elapsed time so far: 76.0 s
73.9 % 2.7e+05 s Elapsed time so far: 76.1 s
74.1 % 2.7e+05 s Elapsed time so far: 76.2 s
74.2 % 2.7e+05 s Elapsed time so far: 76.3 s
74.3 % 2.7e+05 s Elapsed time so far: 76.3 s
74.5 % 2.7e+05 s Elapsed time so far: 76.4 s
74.6 % 2.7e+05 s Elapsed time so far: 76.5 s
74.8 % 2.7e+05 s Elapsed time so far: 76.6 s
74.9 % 2.7e+05 s Elapsed time so far: 76.7 s
75.0 % 2.7e+05 s Elapsed time so far: 76.7 s
75.2 % 2.7e+05 s Elapsed time so far: 76.8 s
75.3 % 2.7e+05 s Elapsed time so far: 76.9 s
75.5 % 2.7e+05 s Elapsed time so far: 77.0 s
75.6 % 2.7e+05 s Elapsed time so far: 77.1 s
75.7 % 2.7e+05 s Elapsed time so far: 77.2 s
75.9 % 2.7e+05 s Elapsed time so far: 77.2 s
76.0 % 2.7e+05 s Elapsed time so far: 77.3 s
76.2 % 2.7e+05 s Elapsed time so far: 77.4 s
76.3 % 2.7e+05 s Elapsed time so far: 77.5 s
76.4 % 2.7e+05 s Elapsed time so far: 77.6 s
76.6 % 2.8e+05 s Elapsed time so far: 77.7 s
76.7 % 2.8e+05 s Elapsed time so far: 77.8 s
76.8 % 2.8e+05 s Elapsed time so far: 77.8 s
77.0 % 2.8e+05 s Elapsed time so far: 77.9 s
77.1 % 2.8e+05 s Elapsed time so far: 78.0 s
77.3 % 2.8e+05 s Elapsed time so far: 78.1 s
77.4 % 2.8e+05 s Elapsed time so far: 78.2 s
77.5 % 2.8e+05 s Elapsed time so far: 78.2 s
77.7 % 2.8e+05 s Elapsed time so far: 78.3 s
77.8 % 2.8e+05 s Elapsed time so far: 78.4 s
78.0 % 2.8e+05 s Elapsed time so far: 78.5 s
78.1 % 2.8e+05 s Elapsed time so far: 78.6 s
78.2 % 2.8e+05 s Elapsed time so far: 78.7 s
78.4 % 2.8e+05 s Elapsed time so far: 78.7 s
78.5 % 2.8e+05 s Elapsed time so far: 78.8 s
78.7 % 2.8e+05 s Elapsed time so far: 78.9 s
78.8 % 2.8e+05 s Elapsed time so far: 79.0 s
78.9 % 2.8e+05 s Elapsed time so far: 79.1 s
79.1 % 2.8e+05 s Elapsed time so far: 79.1 s
79.2 % 2.8e+05 s Elapsed time so far: 79.2 s
79.4 % 2.9e+05 s Elapsed time so far: 79.3 s
79.5 % 2.9e+05 s Elapsed time so far: 79.4 s
79.6 % 2.9e+05 s Elapsed time so far: 79.5 s
79.8 % 2.9e+05 s Elapsed time so far: 79.6 s
79.9 % 2.9e+05 s Elapsed time so far: 79.7 s
80.0 % 2.9e+05 s Elapsed time so far: 79.7 s
80.2 % 2.9e+05 s Elapsed time so far: 79.8 s
80.3 % 2.9e+05 s Elapsed time so far: 79.9 s
80.5 % 2.9e+05 s Elapsed time so far: 80.0 s
80.6 % 2.9e+05 s Elapsed time so far: 80.1 s
80.8 % 2.9e+05 s Elapsed time so far: 80.1 s
80.9 % 2.9e+05 s Elapsed time so far: 80.2 s
81.0 % 2.9e+05 s Elapsed time so far: 80.3 s
81.2 % 2.9e+05 s Elapsed time so far: 80.4 s
81.3 % 2.9e+05 s Elapsed time so far: 80.5 s
81.4 % 2.9e+05 s Elapsed time so far: 80.6 s
81.6 % 2.9e+05 s Elapsed time so far: 80.6 s
81.7 % 2.9e+05 s Elapsed time so far: 80.7 s
81.9 % 2.9e+05 s Elapsed time so far: 80.8 s
82.0 % 2.9e+05 s Elapsed time so far: 80.9 s
82.1 % 3.0e+05 s Elapsed time so far: 81.0 s
82.3 % 3.0e+05 s Elapsed time so far: 81.0 s
82.4 % 3.0e+05 s Elapsed time so far: 81.1 s
82.5 % 3.0e+05 s Elapsed time so far: 81.2 s
82.7 % 3.0e+05 s Elapsed time so far: 81.3 s
82.8 % 3.0e+05 s Elapsed time so far: 81.4 s
83.0 % 3.0e+05 s Elapsed time so far: 81.5 s
83.1 % 3.0e+05 s Elapsed time so far: 81.5 s
83.2 % 3.0e+05 s Elapsed time so far: 81.6 s
83.4 % 3.0e+05 s Elapsed time so far: 81.7 s
83.5 % 3.0e+05 s Elapsed time so far: 81.8 s
83.7 % 3.0e+05 s Elapsed time so far: 81.9 s
83.8 % 3.0e+05 s Elapsed time so far: 82.0 s
83.9 % 3.0e+05 s Elapsed time so far: 82.0 s
84.1 % 3.0e+05 s Elapsed time so far: 82.1 s
84.2 % 3.0e+05 s Elapsed time so far: 82.2 s
84.4 % 3.0e+05 s Elapsed time so far: 82.3 s
84.5 % 3.0e+05 s Elapsed time so far: 82.4 s
84.6 % 3.0e+05 s Elapsed time so far: 82.5 s
84.8 % 3.0e+05 s Elapsed time so far: 82.5 s
84.9 % 3.1e+05 s Elapsed time so far: 82.6 s
85.1 % 3.1e+05 s Elapsed time so far: 82.7 s
85.2 % 3.1e+05 s Elapsed time so far: 82.8 s
85.3 % 3.1e+05 s Elapsed time so far: 82.9 s
85.5 % 3.1e+05 s Elapsed time so far: 82.9 s
85.6 % 3.1e+05 s Elapsed time so far: 83.0 s
85.8 % 3.1e+05 s Elapsed time so far: 83.1 s
85.9 % 3.1e+05 s Elapsed time so far: 83.2 s
86.0 % 3.1e+05 s Elapsed time so far: 83.3 s
86.2 % 3.1e+05 s Elapsed time so far: 83.4 s
86.3 % 3.1e+05 s Elapsed time so far: 83.4 s
86.5 % 3.1e+05 s Elapsed time so far: 83.5 s
86.6 % 3.1e+05 s Elapsed time so far: 83.6 s
86.7 % 3.1e+05 s Elapsed time so far: 83.7 s
86.9 % 3.1e+05 s Elapsed time so far: 83.8 s
87.0 % 3.1e+05 s Elapsed time so far: 83.8 s
87.1 % 3.1e+05 s Elapsed time so far: 83.9 s
87.3 % 3.1e+05 s Elapsed time so far: 84.0 s
87.4 % 3.1e+05 s Elapsed time so far: 84.1 s
87.6 % 3.1e+05 s Elapsed time so far: 84.2 s
87.7 % 3.2e+05 s Elapsed time so far: 84.2 s
87.8 % 3.2e+05 s Elapsed time so far: 84.3 s
88.0 % 3.2e+05 s Elapsed time so far: 84.4 s
88.1 % 3.2e+05 s Elapsed time so far: 84.5 s
88.2 % 3.2e+05 s Elapsed time so far: 84.6 s
88.4 % 3.2e+05 s Elapsed time so far: 84.7 s
88.5 % 3.2e+05 s Elapsed time so far: 84.7 s
88.7 % 3.2e+05 s Elapsed time so far: 84.8 s
88.8 % 3.2e+05 s Elapsed time so far: 84.9 s
89.0 % 3.2e+05 s Elapsed time so far: 85.0 s
89.1 % 3.2e+05 s Elapsed time so far: 85.1 s
89.2 % 3.2e+05 s Elapsed time so far: 85.1 s
89.4 % 3.2e+05 s Elapsed time so far: 85.2 s
89.5 % 3.2e+05 s Elapsed time so far: 85.3 s
89.6 % 3.2e+05 s Elapsed time so far: 85.4 s
89.8 % 3.2e+05 s Elapsed time so far: 85.5 s
89.9 % 3.2e+05 s Elapsed time so far: 85.6 s
90.1 % 3.2e+05 s Elapsed time so far: 85.6 s
90.2 % 3.2e+05 s Elapsed time so far: 85.7 s
90.3 % 3.2e+05 s Elapsed time so far: 85.8 s
90.5 % 3.3e+05 s Elapsed time so far: 85.9 s
90.6 % 3.3e+05 s Elapsed time so far: 86.0 s
90.8 % 3.3e+05 s Elapsed time so far: 86.1 s
90.9 % 3.3e+05 s Elapsed time so far: 86.1 s
91.0 % 3.3e+05 s Elapsed time so far: 86.2 s
91.2 % 3.3e+05 s Elapsed time so far: 86.3 s
91.3 % 3.3e+05 s Elapsed time so far: 86.4 s
91.5 % 3.3e+05 s Elapsed time so far: 86.5 s
91.6 % 3.3e+05 s Elapsed time so far: 86.5 s
91.7 % 3.3e+05 s Elapsed time so far: 86.6 s
91.9 % 3.3e+05 s Elapsed time so far: 86.7 s
92.0 % 3.3e+05 s Elapsed time so far: 86.8 s
92.2 % 3.3e+05 s Elapsed time so far: 86.9 s
92.3 % 3.3e+05 s Elapsed time so far: 87.0 s
92.4 % 3.3e+05 s Elapsed time so far: 87.0 s
92.6 % 3.3e+05 s Elapsed time so far: 87.1 s
92.7 % 3.3e+05 s Elapsed time so far: 87.2 s
92.8 % 3.3e+05 s Elapsed time so far: 87.3 s
93.0 % 3.3e+05 s Elapsed time so far: 87.4 s
93.1 % 3.3e+05 s Elapsed time so far: 87.5 s
93.3 % 3.4e+05 s Elapsed time so far: 87.5 s
93.4 % 3.4e+05 s Elapsed time so far: 87.6 s
93.5 % 3.4e+05 s Elapsed time so far: 87.7 s
93.7 % 3.4e+05 s Elapsed time so far: 87.8 s
93.8 % 3.4e+05 s Elapsed time so far: 87.9 s
94.0 % 3.4e+05 s Elapsed time so far: 88.0 s
94.1 % 3.4e+05 s Elapsed time so far: 88.0 s
94.2 % 3.4e+05 s Elapsed time so far: 88.1 s
94.4 % 3.4e+05 s Elapsed time so far: 88.2 s
94.5 % 3.4e+05 s Elapsed time so far: 88.3 s
94.7 % 3.4e+05 s Elapsed time so far: 88.4 s
94.8 % 3.4e+05 s Elapsed time so far: 88.4 s
94.9 % 3.4e+05 s Elapsed time so far: 88.5 s
95.1 % 3.4e+05 s Elapsed time so far: 88.6 s
95.2 % 3.4e+05 s Elapsed time so far: 88.7 s
95.3 % 3.4e+05 s Elapsed time so far: 88.8 s
95.5 % 3.4e+05 s Elapsed time so far: 88.9 s
95.6 % 3.4e+05 s Elapsed time so far: 88.9 s
95.8 % 3.4e+05 s Elapsed time so far: 89.0 s
95.9 % 3.4e+05 s Elapsed time so far: 89.1 s
96.0 % 3.5e+05 s Elapsed time so far: 89.2 s
96.2 % 3.5e+05 s Elapsed time so far: 89.3 s
96.3 % 3.5e+05 s Elapsed time so far: 89.3 s
96.5 % 3.5e+05 s Elapsed time so far: 89.4 s
96.6 % 3.5e+05 s Elapsed time so far: 89.5 s
96.7 % 3.5e+05 s Elapsed time so far: 89.6 s
96.9 % 3.5e+05 s Elapsed time so far: 89.7 s
97.0 % 3.5e+05 s Elapsed time so far: 89.8 s
97.2 % 3.5e+05 s Elapsed time so far: 89.9 s
97.3 % 3.5e+05 s Elapsed time so far: 89.9 s
97.4 % 3.5e+05 s Elapsed time so far: 90.0 s
97.6 % 3.5e+05 s Elapsed time so far: 90.1 s
97.7 % 3.5e+05 s Elapsed time so far: 90.2 s
97.8 % 3.5e+05 s Elapsed time so far: 90.3 s
98.0 % 3.5e+05 s Elapsed time so far: 90.3 s
98.1 % 3.5e+05 s Elapsed time so far: 90.4 s
98.3 % 3.5e+05 s Elapsed time so far: 90.5 s
98.4 % 3.5e+05 s Elapsed time so far: 90.6 s
98.5 % 3.5e+05 s Elapsed time so far: 90.7 s
98.7 % 3.5e+05 s Elapsed time so far: 90.7 s
98.8 % 3.6e+05 s Elapsed time so far: 90.8 s
99.0 % 3.6e+05 s Elapsed time so far: 90.9 s
99.1 % 3.6e+05 s Elapsed time so far: 91.0 s
99.2 % 3.6e+05 s Elapsed time so far: 91.1 s
99.4 % 3.6e+05 s Elapsed time so far: 91.2 s
99.5 % 3.6e+05 s Elapsed time so far: 91.2 s
99.7 % 3.6e+05 s Elapsed time so far: 91.3 s
99.8 % 3.6e+05 s Elapsed time so far: 91.4 s
99.9 % 3.6e+05 s Elapsed time so far: 91.5 s
100.0 % 3.6e+05 s Elapsed time so far: 91.6 s
Comparison with experimental data and MHIMS: D content#
The results produced by FESTIM are in good agreement with the experimental data and correlate perfectly with MHIMS.
Show code cell source
import plotly.graph_objects as go
import plotly.express as px
from scipy.interpolate import interp1d
def RMSPE(x_sim, x_exp):
error = np.sqrt(np.mean((x_sim - x_exp) ** 2)) / np.mean(x_exp)
return error
retention = (
np.array(derived_quantities[0].data)
+ np.array(derived_quantities[1].data)
+ np.array(derived_quantities[2].data)
)
t = np.array(derived_quantities.t)
exp_ret = np.loadtxt("./reference_data/exp_ret.csv", delimiter=",", skiprows=1)
interp_ret = interp1d(t / 3600, retention, fill_value="extrapolate")
error = RMSPE(
interp_ret(
exp_ret[:, 0],
),
exp_ret[:, 1],
)
print(f"RMSPE between FESTIM and experimental data is {error*100:.2f}%")
MHIMS_ret = np.loadtxt("./reference_data/MHIMS_ret.csv", delimiter=",", skiprows=1)
fig = go.Figure()
fig.add_trace(
go.Scatter(
x=t / 3600,
y=retention / 1e20,
mode="lines",
line=dict(width=4, color=px.colors.qualitative.Plotly[1]),
name="FESTIM",
)
)
fig.add_trace(
go.Scatter(
x=MHIMS_ret[::100, 0],
y=MHIMS_ret[::100, 1] / 1e20,
mode="markers",
marker_symbol="square",
marker=dict(size=10, color=px.colors.qualitative.Plotly[2], opacity=0.6),
name="MHIMS",
)
)
fig.add_trace(
go.Scatter(
x=exp_ret[:, 0],
y=exp_ret[:, 1] / 1e20,
mode="markers",
marker=dict(size=10, color=px.colors.qualitative.Plotly[0], opacity=0.6),
name="Exp.",
)
)
fig.update_yaxes(
title_text="D inventory, 10<sup>20</sup> m<sup>-2</sup>",
range=[0, 5],
tick0=0,
dtick=1,
)
fig.update_xaxes(title_text="Time, h", range=[0, 100], tick0=0, dtick=10)
fig.update_layout(template="simple_white", height=600)
# The writing-reading block below is needed to avoid the issue with compatibility
# of Plotly plots and dollarmath syntax extension in Jupyter Book
# For mode details, see https://github.com/jupyter-book/jupyter-book/issues/1528
fig.write_html("./markelj_comparison_ret.html")
from IPython.display import HTML, display
display(HTML("./markelj_comparison_ret.html"))
RMSPE between FESTIM and experimental data is 7.30%
Comparison with experimental data: D depth distribution#
FESTIM reproduces well the experimental NRA measurements.
Show code cell source
fig = go.Figure()
FESTIM_profiles = np.genfromtxt("./FESTIM_sim.txt", names=True, delimiter=",")
NRA_exp = np.loadtxt("./reference_data/exp_NRA.csv", delimiter=",", skiprows=1)
NRA_sto = np.loadtxt("./reference_data/sto_NRA.csv", delimiter=",", skiprows=1)
def create_slider(fig):
steps = []
for i in range(0, len(fig.data), 2):
step = dict(
method="update",
args=[
{"visible": [False] * len(fig.data)},
{"title": "Time: " + f"{export_times[int(i/2)]/3600:.2f}"},
], # layout attribute
label=f"{export_times[int(i/2)]/3600:.2f} h",
)
step["args"][0]["visible"][i] = True # Toggle i'th trace to "visible"
step["args"][0]["visible"][i + 1] = True # Toggle i+1'th trace to "visible"
steps.append(step)
sliders = [
dict(active=0, currentvalue={"prefix": "Time: "}, pad={"t": 50}, steps=steps)
]
return sliders
for i, t in enumerate(export_times):
x = FESTIM_profiles["x"]
y = FESTIM_profiles[f"t{t:.2e}s".replace(".", "").replace("+", "")]
# order y by x
x, y = zip(*sorted(zip(x, y)))
color = px.colors.qualitative.Plotly[i]
fig.add_trace(
go.Scatter(
x=np.array(x) / 1e-6,
y=np.array(y) / rho_W * 100,
mode="lines",
line=dict(width=3.5, color=px.colors.qualitative.Plotly[i]),
name="FESTIM",
visible=False,
)
)
if i <= 4:
fig.add_trace(
go.Scatter(
x=NRA_exp[:, 0],
y=NRA_exp[:, i + 1],
mode="lines",
line=dict(
width=3.5, color=px.colors.qualitative.Plotly[i], dash="dash"
),
name="NRA",
visible=False,
)
)
else:
fig.add_trace(
go.Scatter(
x=NRA_sto[:, 0],
y=NRA_sto[:, i - 4],
mode="lines",
line=dict(
width=3.5, color=px.colors.qualitative.Plotly[i], dash="dash"
),
name="NRA",
visible=False,
)
)
fig.data[0].visible = True
fig.data[1].visible = True
fig.update_yaxes(title_text="D concentration, at.%", range=[0, 0.5], tick0=0, dtick=0.1)
fig.update_xaxes(title_text="Depth, µm", range=[0, 5], tick0=0, dtick=1)
fig.update_layout(template="simple_white", sliders=create_slider(fig), height=600)
# The writing-reading block below is needed to avoid the issue with compatibility
# of Plotly plots and dollarmath syntax extension in Jupyter Book
# For mode details, see https://github.com/jupyter-book/jupyter-book/issues/1528
fig.write_html("./markelj_profiles_exp.html")
display(HTML("./markelj_profiles_exp.html"))
Comparison with MHIMS: D depth distribution#
FESTIM agrees with MHIMS. Slight differences are due to the use of the precise value of Fernandez’s diffusivity pre-factor. FESTIM uses the value from HTM divided by \(\sqrt{2}\): \(1.93 \times 10^{-7} / \sqrt{2}\), whereas the presented diffusivity pre-factor in [26] is \(1.9 \times 10^{-7} / \sqrt{2}\).
Show code cell source
fig = go.Figure()
FESTIM_profiles = np.genfromtxt("./FESTIM_sim.txt", names=True, delimiter=",")
MHIMS_profiles = np.loadtxt(
"./reference_data/MHIMS_profiles.csv", delimiter=",", skiprows=1
)
for i, t in enumerate(export_times):
x = FESTIM_profiles["x"]
y = FESTIM_profiles[f"t{t:.2e}s".replace(".", "").replace("+", "")]
# order y by x
x, y = zip(*sorted(zip(x, y)))
color = px.colors.qualitative.Plotly[i]
fig.add_trace(
go.Scatter(
x=np.array(x) / 1e-6,
y=np.array(y) / rho_W * 100,
mode="lines",
line=dict(width=3.5, color=px.colors.qualitative.Plotly[i]),
name="FESTIM",
visible=False,
)
)
fig.add_trace(
go.Scatter(
x=MHIMS_profiles[:, 0],
y=MHIMS_profiles[:, i + 1],
mode="markers",
marker=dict(size=10, color=px.colors.qualitative.Plotly[i], opacity=0.4),
name="MHIMS",
visible=False,
)
)
fig.data[0].visible = True
fig.data[1].visible = True
fig.update_yaxes(title_text="D concentration, at.%", range=[0, 0.5], tick0=0, dtick=0.1)
fig.update_xaxes(title_text="Depth, µm", range=[0, 5], tick0=0, dtick=1)
fig.update_layout(template="simple_white", sliders=create_slider(fig), height=600)
# The writing-reading block below is needed to avoid the issue with compatibility
# of Plotly plots and dollarmath syntax extension in Jupyter Book
# For mode details, see https://github.com/jupyter-book/jupyter-book/issues/1528
fig.write_html("./markelj_profiles_MHIMS.html")
display(HTML("./markelj_profiles_MHIMS.html"))