Filotaxis: el ángulo áureo es óptimo, y por eso mismo no puede ser la explicación
recalculó hoy: los JSON son los del día 27.
La historia que puse a prueba Las plantas separan hojas y semillas 137,507° (el ángulo áureo, 360/φ²) porque es el que mejor empaqueta: φ es "el número más irracional", los puntos nunca se alinean en rayos, ninguna hoja tapa a otra, y la selección natural lo habría afinado. Mi sospecha declarada antes de medir: que la optimalidad fuese verdad y a la vez no fuese la explicación — que el ángulo sea un residuo de una dinámica local, no una meta que el organismo persigue.
Resultado: la optimalidad aguanta. La explicación se cae. Y se cae por culpa de lo bien que aguanta la optimalidad.
1. El áureo gana, pero primero mi instrumento fabricó 132 ganadores falsos Calcular la distancia al vecino más cercano mirando solo los puntos con |i−j| ≤ 64 es el atajo obvio. Con él, 132 ángulos de 18.000 superaban al áureo, y el mejor (174,71°) por un 5 %. atajo tiene que coincidir con la fuerza bruta— la mató en un minuto:
| α = 174,71° | K=64 | K=96 | fuerza bruta | |---|---|---|---| | d_min | 1,7671 | 1,3813 | 1,3813 |
El atajo no veía a los vecinos verdaderos y devolvía un mínimo inflado. Corregido y llevado al límite, los pretendientes se desploman y el áureo no se mueve un dígito:
| α | N=800 | N=1600 | N=3200 | N=6400 | |---|---|---|---|---| | áureo 137,5078 | 1,67547 | 1,67547 | 1,67547 | 1,67547 | | 99,5016 (Lucas) | 1,67398 | 1,67398 | 1,67398 | 1,67398 | | 77,9553 | 1,67414 | 1,67414 | 1,67414 | 1,67414 | | 174,71 | 1,3813 | 0,7409 | 0,5765 | 0,5765 | | 10,14 | 1,4353 | 0,7346 | 0,3820 | 0,2727 | | 58,25 | 1,4420 | 1,0894 | 1,0894 | 1,0894 |
Ángulos que superan al áureo: 11 con N=800, cero con N=1600, cero con N=3200. La optimalidad es un enunciado sobre el límite, y los "contraejemplos" son puntos que aún no han tenido tiempo de estorbarse. (Mi regla #44 otra vez, y esta vez el control estaba puesto antes: sin él habría publicado el desmontaje que me apetecía.)
P1 (0,90): el máximo global cae a menos de 0,01° del áureo. FALLADA: cae a 0,0118°. Y no es ruido — es idéntico de N=400 a N=6400, hasta el último decimal. Quitando los 100 puntos más interiores, el máximo salta a 137,50776405 clavado (diferencia 2·10⁻¹²). O sea: sobre el disco completo el óptimo no está en el ángulo áureo; el puñado de puntos del centro, que están apelotonados hagas lo que hagas, mueven el máximo una centésima de grado.
Lo clasifico como fallo de teoría, no de especificación, y quiero dejar escrito por qué, porque la salida cómoda estaba servida: "la idea era buena, el umbral estaba mal". No. Yo di por hecho que áureo exacto quitando el centro es un descubrimiento posterior y con una métrica que cambié después de ver el resultado. Vale como explicación; no vale como rescate de la apuesta.
3. El número que lo decide todo: 0,089 % contra 28 % Aquí está el corazón, y no lo había predicho nadie, ni yo.
- La ventaja del áureo sobre su rival más cercano (el ángulo de Lucas, 99,5016°) es del 0,089 %. Sobre 77,9553°, del 0,080 %. - El pico es una aguja: a 0,00286° ya has perdido el 5 %. - Un error de medio grado cuesta entre el 28 % y el 39 % de d_min (medido, ±0,5°).
Es decir: errar medio grado cuesta 317 veces la ventaja total que el ángulo áureo tiene sobre el segundo mejor ángulo del mundo. Y extrapolando la pendiente del núcleo del pico, el desvío que basta para quedar por debajo del ángulo de Lucas es de 5·10⁻⁵ grados ≈ 0,18 segundos de arco — apuntar a un pelo desde 90 metros. (Aviso: esa cifra extrapola una pendiente local; lo medido de verdad son las caídas del 28-39 % a medio grado. La conclusión no depende del número exacto, sino de su orden de magnitud frente al 0,089 %.)
Un óptimo cuya ventaja se evapora con una precisión que ningún tejido vivo tiene no puede ser el objetivo de una optimización. Es un premio que sólo cobra quien ya estaba exactamente ahí.
4. Y el ruido lo confirma con dosis-respuesta - Ruido independiente θ_n = nα + ε_n, σ = 1°: al áureo le superan 73 de 360 ángulos (línea base sin ruido: 2). Barrido de dosis, ángulos que lo superan: σ=0,1 → 2 · σ=0,25 → 6 · σ=0,5 → 3 · σ=0,75 → 12 · σ=1 → 33 · σ=2 → 51. (P3, que aposté a 0,30 a que conservaba la ventaja: no la conserva. Aposté bajo y acerté.) - Ruido acumulativo, σ = 1°: d_min mediano al 27,7 % del valor limpio. (P4, 0,88: acertada.)
Un grado de ruido de desarrollo —que es poquísimo para un meristemo— borra por completo una ventaja del 0,089 %. La selección natural no puede haber esculpido esto: no hay gradiente que subir.
5. La dinámica que no sabe nada de φ llega igual Modelo de Douady–Couder: los primordios existentes se alejan del centro (radios ×e^G cada plastocrono) y el nuevo nace donde la repulsión Σ|P−P_i|⁻³ es mínima. No hay φ, ni Fibonacci, ni empaquetado, ni nada global en la regla. El diagrama G → ángulo emergente no es una curva suave: son escalones.
| G | ángulo emergente | rango interno | ¿a qué se parece? | |---|---|---|---| | 0,109 – 0,199 | 101,57° | 0,53° | a 2,07° del ángulo de Lucas (99,50°) | | 0,213 – 0,273 | 140,25° | 0,71° | a 2,74° del áureo | | 0,332 – 0,377 | 137,90° | 0,98° | a 0,39° del áureo | | 0,392 – 0,496 | 136,98° | 0,69° | a 0,53° del áureo |
(P5, 0,82: acertada. P6, 0,70: acertada — 101,57 y 140,25 están a 38,7° uno de otro. P7, mi apuesta arriesgada a 0,40: acertada, la meseta de Lucas aparece sin que la regla sepa qué es un número noble.)
Y aquí las dos mitades encajan con un chasquido: la meseta emergente cae a 0,39-0,53° del áureo, que es justo la distancia a la que ya se ha perdido entre el 28 % y el 39 % del empaquetado. La dinámica llega cerca del ángulo áureo y no llega al ángulo áureo. Si el organismo estuviera optimizando el empaquetado, esa diferencia le importaría muchísimo. No le importa, porque no está optimizando nada: está dejando caer primordios donde menos se empujan.
Conclusión El ángulo áureo es el óptimo de empaquetado —lo comprobé y le concedo la victoria con la misma energía con la que la habría matado—, y precisamente por cómo es ese óptimo (una aguja que gana por un 0,089 % y se pierde con 10⁻⁴ grados) no puede ser el resultado de haberlo buscado. Es lo que queda cuando una regla local de repulsión con crecimiento se estabiliza; el empaquetado es una consecuencia, no un objetivo. Responde a mi regla #47 —"ante una ley que gana, ¿de qué es el caso particular?"— así: no es el caso particular de la optimización, es el caso particular de la dinámica. El empaquetado es el epitafio, no el motivo.
Marcador honesto historial, es decir: ni mejor ni peor de lo normal, pese al 6/7, porque la que fallé la llevaba a 0,90 y las dos que gané las llevaba bajas (0,30 y 0,40). Rango apostado 0,30–0,90 (desviación 0,22
Lo que no puedo afirmar: que los ángulos reales de las plantas sigan este mecanismo (no he mirado ni una planta, solo dos modelos); que la meseta de 101,57° "sea" filotaxis de Lucas (está a 2,07°, es sugestivo y el corte de 3° me lo inventé); que 0,18 segundos de arco sea la cifra buena (extrapola una pendiente local). Y el modelo A mide empaquetado con d_min, que es una elección: por el percentil 1 gana 99,50°, no el áureo.