Estoy tratando de fabricar lo que se llamaTrama de Arrheniusen diagramas de páginas. Básicamente es un gráfico semilogarítmico (eje y) donde se traza una función dependiente de la temperatura frente a 1/kB*T (siendo kB elConstante de Boltzmann) Esto funciona hasta ahora, pero, por ejemplo, 1600°C corresponden aproximadamente a 6,2 1/eV (recíproco).electronvoltio)
Científicamente eso es correcto, pero nadie puede darle un significado físico a ese número (no puedo configurar mi horno a 6,2 1/eV o en una charla, donde las diapositivas solo son visibles durante unos segundos, podría confundir al oyente). Así que quería poner un segundo eje x en el gráfico que mostrara la temperatura en grados Celsius. Pensé en la capacidad de pgfplots para transformar las coordenadas pero no lo resuelvo.
A continuación hay un MWE de lo que probé. Debe haber algún error en mi proceso de pensamiento sobre la transformación de coordenadas, pero no puedo resolverlo. Los puntos finales de la trama coinciden, pero entre ellos el mapeo está ligeramente desviado. ¿Alguien puede sugerir cómo corregir esto?
\documentclass[tikz=true]{standalone}
\usepackage[T1]{fontenc}
\usepackage[utf8]{inputenc}
\usepackage[english]{babel}
\usepackage{mathptmx}
\usepackage{siunitx}
\usepackage{pgfplots}
\pgfplotsset{compat=newest}
\begin{document}
\begin{tikzpicture}
\begin{semilogyaxis}[
axis x line* = top,
axis y line* = left,
xlabel={$\frac{1}{k_B \cdot T} / \si[per-mode=reciprocal]{\per\electronvolt}$},
ylabel={$k$},
]
\addplot+[x = kb, y = k] table {
kb T k
7.61876 1250 3.98E-17
7.37661 1300 1.68E-17
7.23857 1330 1.00E-17
7.14938 1350 7.09E-18
7.10560 1360 5.97E-18
7.06236 1370 5.03E-18
7.01964 1380 4.23E-18
6.97743 1390 3.56E-18
6.93573 1400 3.00E-18
6.89452 1410 2.52E-18
6.85380 1420 2.12E-18
6.83362 1425 1.95E-18
6.77379 1440 1.50E-18
6.69562 1460 1.07E-18
6.61924 1480 7.55E-19
6.54458 1500 5.35E-19
6.45359 1525 3.48E-19
6.36509 1550 2.26E-19
6.19519 1600 9.54E-20
};
\end{semilogyaxis}
% % extra x axis
\begin{semilogyaxis}[
x dir = reverse,
axis y line* = right,
axis x line* = bottom,
xlabel={$T / \si{\celsius}$},
x coord trafo/.code={\pgfmathparse{(1/(#1*8.617333262145E-05)-273.15)}},
]
\addplot+[x = T, y = k] table {
kb T k
7.61876 1250 3.98E-17
7.37661 1300 1.68E-17
7.23857 1330 1.00E-17
7.14938 1350 7.09E-18
7.10560 1360 5.97E-18
7.06236 1370 5.03E-18
7.01964 1380 4.23E-18
6.97743 1390 3.56E-18
6.93573 1400 3.00E-18
6.89452 1410 2.52E-18
6.85380 1420 2.12E-18
6.83362 1425 1.95E-18
6.77379 1440 1.50E-18
6.69562 1460 1.07E-18
6.61924 1480 7.55E-19
6.54458 1500 5.35E-19
6.45359 1525 3.48E-19
6.36509 1550 2.26E-19
6.19519 1600 9.54E-20
};
\end{semilogyaxis}
%
\end{tikzpicture}
\end{document}
La tabla son datos ficticios, que se calculan: se da T, kB se calcula a partir de 1/kB*T y k es una constante de velocidad imaginaria calculada a partir de una ecuación de Arrhenius.
Respuesta1
No es 100% lo que estás buscando, pero eso es al menos algo quesaberes correcto. Supongo que usted mismo podrá hacer las modificaciones necesarias desde aquí.
% used PGFPlots v1.16
\documentclass[border=5pt]{standalone}
\usepackage{xfp} % <-- needed for accuracy
\usepackage{siunitx}
\usepackage{pgfplotstable}
\pgfplotsset{compat=1.16}
\pgfplotstableread{
kb T k
7.61876 1250 3.98E-17
7.37661 1300 1.68E-17
7.23857 1330 1.00E-17
7.14938 1350 7.09E-18
7.10560 1360 5.97E-18
7.06236 1370 5.03E-18
7.01964 1380 4.23E-18
6.97743 1390 3.56E-18
6.93573 1400 3.00E-18
6.89452 1410 2.52E-18
6.85380 1420 2.12E-18
6.83362 1425 1.95E-18
6.77379 1440 1.50E-18
6.69562 1460 1.07E-18
6.61924 1480 7.55E-19
6.54458 1500 5.35E-19
6.45359 1525 3.48E-19
6.36509 1550 2.26E-19
6.19519 1600 9.54E-20
}{\data}
\begin{document}
\begin{tikzpicture}
\begin{semilogyaxis}[
axis lines*=left,
xlabel={$\frac{\SI{1000}{\kelvin}}{T}$},
ylabel={$k$},
]
% please note that table options have to be given *after* the
% table keyword. Otherwise they will not be applied
\addplot table [
x expr={1000/\thisrow{T}},
y=k,
] {\data};
\end{semilogyaxis}
% extra x axis
\begin{semilogyaxis}[
axis lines*=right,
xlabel={$T / \si{\kelvin}$},
% avoid scaling (to avoid unexpected result)
scaled ticks=false,
% state values which should be shown on axis
xtick={1250,1300,1350,1400,1500,1600},
% perform coordinate transformation
% to ensure right results (accuracy), use the `\fpeval` command
% from the `xfp` package
% (compare results when not using it. This will be way off for high T)
x coord trafo/.code={\pgfmathparse{\fpeval{1/(#1)}}},
x coord inv trafo/.code={\pgfmathparse{\fpeval{1/(#1)}}},
]
\addplot+ [red] table [
x=T,
y=k,
] {\data};
\end{semilogyaxis}
\end{tikzpicture}
\end{document}