Arquitectura matemática

Este capítulo expone la arquitectura formal de la Schubart Master Formula: qué se integra, sobre qué dominio y por qué la expresión resultante es dimensionalmente consistente. Tras haber motivado el marco físicamente en los capítulos anteriores, el objetivo aquí es hacer plenamente explícito su esqueleto matemático, de modo que todo trabajo numérico o analítico posterior descanse sobre una definición inequívoca.

La representación física expandida

La Schubart Master Formula canónica es la forma compacta de volumen $P(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV$. Para hacer plenamente explícita su arquitectura matemática, este capítulo trabaja con su representación física expandida,

$$P(t) = \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}\;dE\,d\theta .$$

Aquí $P(t)$ es la potencia eléctrica utilizable instantánea [W]; $\Phi(E,\theta,t)$ es la densidad de flujo de radiación; $\sigma(E)$ es la sección eficaz de dispersión; $\Delta p$ es la transferencia de momento dirigido; $v_{\text{ph}}$ es la velocidad de los fonones; y $\eta_{\text{conv}}$ es la eficiencia de conversión. La integral doble se extiende sobre la variable de energía $E$ y la variable angular $\theta$. Esta es la representación física expandida; la misma expresión, junto con el glosario estándar de símbolos, se reproduce en la ecuación y en la referencia de notación.

Leer el integrando como un producto de etapas

El integrando se factoriza en una cadena de etapas físicas, y esta factorización es la idea organizadora de todo el marco. Leyendo de izquierda a derecha: $\Phi$ suministra el recurso incidente; $\sigma$ fija con qué intensidad se acopla ese recurso al material; $\Delta p$ convierte el acoplamiento en momento dirigido; $v_{\text{ph}}$ transporta ese momento; y $\eta_{\text{conv}}$ lo rectifica en potencia utilizable. Dado que aparecen como un producto, todo el integrando se anula si cualquier factor se anula — la expresión matemática del argumento de la «cadena» desarrollado en conversión de energía en no equilibrio. La forma de producto significa también que los factores son, a primer orden, ajustables de forma independiente: mejorar el transporte no cambia el acoplamiento, y viceversa, que es lo que hace al marco útil para aislar cuellos de botella.

La factorización puede hacerse precisa. Solo dos factores portan las variables de integración — el flujo $\Phi(E,\theta,t)$ depende de $E$ y $\theta$, y la sección eficaz $\sigma(E)$ depende de $E$ — mientras que $\Delta p$, $v_{\text{ph}}$ y $\eta_{\text{conv}}$ se tratan, en la forma básica, como constantes efectivas representativas del convertidor. Bajo esa suposición los factores constantes salen de la integral,

$$P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \iint \Phi(E,\theta,t)\,\sigma(E)\;dE\,d\theta ,$$

de modo que todo el contenido espectral y angular se concentra en la integral de solapamiento $\iint \Phi\,\sigma\,dE\,d\theta$. Esta separación es lo que autoriza la metodología de «aislar el cuello de botella»: puede estudiarse el solapamiento entorno-material independientemente del prefactor de transporte-rectificación. Cuando alguno de los tres prefactores depende de hecho de la energía o del ángulo — por ejemplo $\Delta p(E,\theta)$ a partir del promedio microscópico — debe permanecer dentro de la integral, y la separación anterior solo vale de forma aproximada. El marco enuncia la forma general con todos los factores dentro de la integral precisamente para que esta generalización no requiera cambio alguno de estructura.

Diferenciar la forma factorizada da la sensibilidad de la salida a cada etapa. Dado que la dependencia es multiplicativa, las sensibilidades relativas simplemente se suman:

$$\frac{dP}{P} = \frac{d(\Delta p)}{\Delta p} + \frac{dv_{\text{ph}}}{v_{\text{ph}}} + \frac{d\eta_{\text{conv}}}{\eta_{\text{conv}}} + \frac{d\!\left(\iint \Phi\sigma\right)}{\iint \Phi\sigma} .$$

Una mejora del diez por ciento en cualquier factor individual produce una mejora del diez por ciento en la salida, y el peor factor domina la ganancia alcanzable. Esta estructura de derivada logarítmica es la columna vertebral cuantitativa del argumento del cuello de botella.

El dominio de integración

El dominio de integración merece cuidado. La integral de energía se extiende sobre el rango de energías presentes en el espectro ambiental y al que el convertidor puede acoplarse; formalmente puede tomarse desde cero hasta infinito, en el entendimiento de que el integrando es despreciable fuera del soporte de $\Phi(E,\theta,t)\,\sigma(E)$. La integral angular se extiende sobre el ángulo sólido subtendido por el flujo incidente en la superficie del convertidor. En muchos tratamientos $\theta$ es un ángulo polar y una integración de ángulo sólido completa introduciría la medida habitual $\sin\theta$; el marco mantiene $\theta$ explícito precisamente para que pueda insertarse, sin cambiar la estructura global, cualquier medida angular apropiada a una geometría dada. El punto formal importante es que los límites los fija el soporte físico del integrando, no una elección arbitraria.

Dependencia temporal

La salida porta un argumento temporal explícito porque el flujo lo hace. El tiempo entra a través de $\Phi(E,\theta,t)$; los factores de acoplamiento y rectificación se tratan como cuasi-estáticos, lo que significa que responden de manera efectivamente instantánea a los cambios de flujo en las escalas de tiempo de interés. En consecuencia, $P(t)$ hereda su dependencia temporal del entorno, y el marco describe naturalmente un convertidor cuya salida sigue un campo ambiental fluctuante. Si alguno de los parámetros del material fuera él mismo dependiente del tiempo — por ejemplo, si $\eta_{\text{conv}}$ derivara con la temperatura — también podría dársele un argumento temporal sin alterar la forma de la integral.

La suposición cuasi-estática puede enunciarse como una separación de escalas de tiempo. Sea $\tau_{\text{env}}$ el tiempo característico en el que varía el flujo ambiental y $\tau_{\text{resp}}$ el tiempo de respuesta interno del convertidor (fijado por la dispersión y la recolección). La forma instantánea de la fórmula maestra es válida cuando

$$\tau_{\text{resp}} \ll \tau_{\text{env}} ,$$

de modo que en cada instante el convertidor ha alcanzado efectivamente su respuesta estacionaria al flujo actual. Cuando esto se cumple, puede promediarse la salida sobre una ventana de observación $T$ para obtener una potencia media,

$$\langle P \rangle = \frac{1}{T}\int_0^T P(t)\;dt ,$$

que es la magnitud que la mayoría de los experimentos reportan en realidad. Si la separación fallara — un flujo que varía más rápido de lo que el convertidor puede responder — el marco requeriría una convolución con un núcleo de respuesta en lugar del producto instantáneo, una extensión que deja intacta la estructura espectral–angular pero reemplaza los factores cuasi-estáticos por funciones de memoria.

Orden de integración y separabilidad

Dado que el integrando es un producto de una parte espectral–angular y un conjunto de factores cuasi-constantes, el teorema de Fubini garantiza que la integral doble puede evaluarse en cualquiera de los dos órdenes, $\int dE \int d\theta$ o $\int d\theta \int dE$, siempre que el integrando sea absolutamente integrable — lo cual es, pues $\Phi$ y $\sigma$ son no negativos y el integrando tiene soporte finito. Esta libertad es prácticamente útil: puede integrarse primero sobre el ángulo para obtener un flujo espectral efectivo $\Phi_{\text{eff}}(E)=\int \Phi(E,\theta,t)\,d\theta$ y luego realizar una única integral de solapamiento espectral,

$$P(t) = \Delta p\, v_{\text{ph}}\, \eta_{\text{conv}} \int_0^\infty \Phi_{\text{eff}}(E,t)\,\sigma(E)\;dE ,$$

o, a la inversa, integrar primero sobre la energía para obtener un patrón angular de potencia. Las dos rutas deben coincidir, y comprobar que lo hacen es una prueba de autoconsistencia simple pero eficaz para cualquier implementación numérica. Cuando un factor acopla genuinamente las dos variables — por ejemplo una sección eficaz dependiente del ángulo $\sigma(E,\theta)$ — se pierde la separabilidad y debe conservarse la integral doble completa, pero el intercambio del orden de integración sigue siendo válido. Esta es otra razón por la que el marco se escribe con ambas integrales explícitas en lugar de colapsadas de antemano en una forma unidimensional efectiva.

Consistencia dimensional

Un marco de este tipo solo es creíble si es dimensionalmente sólido, de modo que conviene comprobarlo. La densidad de flujo $\Phi$ tiene unidades de energía por área, por tiempo, por ángulo sólido, por energía — esto es, $\mathrm{W\,m^{-2}\,sr^{-1}\,J^{-1}}$ cuando la densidad de energía se expresa por julio. La sección eficaz $\sigma$ aporta $\mathrm{m^{2}}$, la transferencia de momento $\Delta p$ aporta $\mathrm{kg\,m\,s^{-1}}$, la velocidad de los fonones $v_{\text{ph}}$ aporta $\mathrm{m\,s^{-1}}$, y $\eta_{\text{conv}}$ es adimensional. La medida de integración $dE\,d\theta$ aporta energía por ángulo sólido, $\mathrm{J\,sr}$. En cualquier realización concreta del marco, el producto de estas unidades, junto con la constante de proporcionalidad que relacione los factores geométricos con la potencia, debe reducirse a vatios. El marco exige que esta reducción se cumpla; una aplicación específica debe suministrar las constantes que hagan que la contabilidad cierre exactamente. Esta disciplina — insistir en que las unidades cuadren antes de confiar en cualquier número — es lo que impide que el formalismo se deslice hacia un terreno no falsable.

Positividad y acotación

Dos propiedades estructurales nos tranquilizan en cuanto a que la fórmula se comporta bien. Primero, la positividad: el flujo $\Phi$ es no negativo por definición, la sección eficaz $\sigma$ es no negativa, $v_{\text{ph}}$ es positiva, y $\eta_{\text{conv}}$ está en $[0,1]$. El único factor que puede llevar signo es la transferencia de momento dirigido $\Delta p$, que es positiva cuando el empuje neto va hacia el eje de recolección y sería negativa para un empuje en contra de él. En un convertidor correctamente orientado $\Delta p \ge 0$, de modo que el integrando es no negativo y $P(t) \ge 0$: el marco no puede reportar potencia entregada negativa, como exige el sentido físico. Segundo, la acotación: dado que $\eta_{\text{conv}} \le 1$, la potencia entregada nunca puede exceder la potencia total rectificable que cruza la frontera,

$$P(t) \;\le\; \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\;dE\,d\theta ,$$

que a su vez está acotada por el flujo incidente. No hay configuración de los parámetros que permita a la salida exceder el recurso suministrado por el entorno. La positividad y la acotación juntas son las garantías matemáticas que corresponden, respectivamente, a «ninguna corriente inversa espontánea» y «ninguna salida más allá del recurso disponible» — las dos cosas que un modelo creíble de conversión de energía nunca debe violar.

De la estructura a la evaluación

Con la arquitectura fijada, el trabajo restante es de interpretación y evaluación. El significado y las unidades de cada factor se catalogan en definiciones de variables; la superficie de control a través de la cual se traza el balance se especifica en fronteras del sistema; la construcción completa paso a paso se da en derivación completa; y las relaciones rectoras que se siguen de la estructura se recopilan en relaciones rectoras.

Lecturas adicionales

• R. Courant y D. Hilbert, Methods of Mathematical Physics, Vol. 1 (Interscience, 1953). Referencia estándar para los operadores integrales y los argumentos de separabilidad usados aquí.
• G. B. Arfken, H. J. Weber y F. E. Harris, Mathematical Methods for Physicists, 7.ª ed. (Academic Press, 2013). Integración multidimensional, medidas y razonamiento dimensional.
• H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2.ª ed. (Wiley, 1985). Sobre la consistencia dimensional y las unidades de las expresiones de balance de energía.
• J. M. Ziman, Electrons and Phonons (Oxford University Press, 1960). Fundamento físico del prefactor de transporte $v_{\text{ph}}$.