Description Usage Arguments Details Value Function version How to cite Note Author(s) References See Also Examples

View source: R/calc_ThermalLifetime.R

The function calculates the thermal lifetime of charges for given E (in eV), s (in 1/s) and T (in deg. C.) parameters. The function can be used in two operational modes:

1 2 3 4 5 6 7 8 9 10 11 |

`E` |
numeric ( |

`s` |
numeric ( |

`T` |
numeric ( |

`output_unit` |
character ( |

`profiling` |
logical ( |

`profiling_config` |
list ( -
`n` (number of MC runs), -
`E.distribution` (distribution used for the resampling for E) and -
`s.distribution` (distribution used for the resampling for s).
Currently only the normal distribution is supported
(e.g., |

`verbose` |
logical: enables/disables verbose mode |

`plot` |
logical:
enables/disables output plot, currenlty only in combination with |

`...` |
further arguments that can be passed in combination with the plot output. Standard plot parameters are supported (plot.default) |

**Mode 1 (profiling = FALSE)**

An arbitrary set of input parameters (E, s, T) can be provided and the function calculates the thermal lifetimes using the Arrhenius equation for all possible combinations of these input parameters. An array with 3-dimensions is returned that can be used for further analyses or graphical output (see example 1)

**Mode 2 (profiling = TRUE)**

This mode tries to profile the variation of the thermal lifetime for a chosen
temperature by accounting for the provided E and s parameters and their corresponding
standard errors, e.g., `E = c(1.600, 0.001)`

The calculation based on a Monte Carlo simulation, where values are sampled from a normal
distribution (for E and s).

**Used equation (Arrhenius equation)**

*τ = 1/s exp(E/kT)*

where:
*τ* in s as the mean time an electron spends in the trap for a given *T*,
*E* trap depth in eV,
*s* the frequency factor in 1/s,
*T* the temperature in K and *k* the Boltzmann constant in eV/K (cf. Furetta, 2010).

A RLum.Results object is returned a along with a plot (for
`profiling = TRUE`

). The output object contain the following slots:

`@data`

Object | Type | Description |

`lifetimes` | array or numeric | calculated lifetimes |

`profiling_matrix` | matrix | profiling matrix used for the MC runs |

`@info`

Object | Type | Description |

`call` | `call` | the original function call |

0.1.0

Kreutzer, S., 2020. calc_ThermalLifetime(): Calculates the Thermal Lifetime using the Arrhenius equation. Function version 0.1.0. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J., 2020. Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.9.7. https://CRAN.R-project.org/package=Luminescence

The profiling is currently based on resampling from a normal distribution, this distribution assumption might be, however, not valid for given E and s paramters.

Sebastian Kreutzer, IRAMAT-CRP2A, Universite Bordeaux Montaigne (France) , RLum Developer Team

Furetta, C., 2010. Handbook of Thermoluminescence, Second Edition. ed. World Scientific.

graphics::matplot, stats::rnorm, get_RLum

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 | ```
##EXAMPLE 1
##calculation for two trap-depths with similar frequency factor for different temperatures
E <- c(1.66, 1.70)
s <- 1e+13
T <- 10:20
temp <- calc_ThermalLifetime(
E = E,
s = s,
T = T,
output_unit = "Ma"
)
contour(x = E, y = T, z = temp$lifetimes[1,,],
ylab = "Temperature [\u00B0C]",
xlab = "Trap depth [eV]",
main = "Thermal Lifetime Contour Plot"
)
mtext(side = 3, "(values quoted in Ma)")
##EXAMPLE 2
##profiling of thermal life time for E and s and their standard error
E <- c(1.600, 0.003)
s <- c(1e+13,1e+011)
T <- 20
calc_ThermalLifetime(
E = E,
s = s,
T = T,
profiling = TRUE,
output_unit = "Ma"
)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.