# Cuaderno de preguntas — proyecto cell
> Registro de las preguntas conceptuales que van surgiendo durante el proyecto y sus respuestas, condensadas para revisitar rápido. Los documentos grandes (`ADN_del_atomo_al_cromosoma.md`, `ADN_simulador_fisica_atomica.md`) son el manual; esto es la libreta de lo hablado. Cada entrada lleva su fecha; las nuevas se añaden al final.
---
## Índice
- [P1 — ¿Qué es la esfera de van der Waals? ¿Es como un campo magnético?](#p1)
- [P2 — ¿Por qué ver el 1BNA (F1) antes de validar la física (F2)?](#p2)
- [P3 — ¿Qué son PME, SHAKE y los «diedros bsc1» que le faltan a nuestro motor?](#p3)
- [P4 — Las cinco fuerzas del sistema, puestas en la misma regla (kJ/mol)](#p4)
- [P5 — ¿Qué es la temperatura para átomos y moléculas, y por qué «evapora» la gota del sandbox?](#p5)
- [P6 — ¿Los átomos «orbitan» como los planetas, por el balance cinética/potencial?](#p6)
- [P7 — Choques, disipación… ¿y si todos los átomos se movieran alineados, sería temperatura?](#p7)
- [P8 — ¿De dónde sale el puente de hidrógeno? ¿Por qué no hay «puentes de hierro»?](#p8)
- [P9 — ¿El universo es EM + reglas cuánticas + un «primer movimiento»? ¿La gravedad es emergente como la temperatura?](#p9)
- [P10 — ¿El CERN acelera átomos para superar la barrera de Pauli y ver qué pasa?](#p10)
- [P11 — ¿Una gota de agua en el espacio, sin gravedad ni atmósfera, mantendría su forma de bola?](#p11)
- [P12 — ¿Los átomos «decaen» emitiendo luz hasta desaparecer? ¿Qué es entonces el fondo de microondas?](#p12)
- [P13 — ¿Por qué las bases son rígidas en F3? ¿Podría su forma mantenerse sola, por las fuerzas?](#p13)
- [P14 — ¿Qué significa que las distancias atómicas «se miden en cristales»?](#p14)
- [P15 — Estados de la materia, romper enlaces con calor… y qué es kT (que no son teslas)](#p15)
- [P16 — ¿La simulación F5 se ve en vivo o se graba? ¿Los fotogramas se conservan o se pierden al cerrar?](#p16)
- [P17 — ¿Cuánto ocuparía grabar los 100 ns con vibraciones de femtosegundos? (la cámara de alta velocidad)](#p17)
- [P18 — En F3, ¿a qué velocidad vemos la realidad y con qué muestreo? (y la GPU que era la equivocada)](#p18)
- [P19 — Las cinco leyes que hay en el código (y cómo se mapean a las cinco fuerzas de la química); ¿el muro de Pauli se deforma según la molécula?](#p19)
- [P20 — ¿Podemos simular sal disolviéndose, ósmosis con paredes virtuales y «aplicar electricidad»? ¿Dónde hace falta la cuántica para no hacer un mock?](#p20)
- [P21 — La presión: su origen atómico (paredes y gravedad), por qué es estadística como la T, y dónde vive en nuestro simulador](#p21)
- [P22 — La gota «parece viva»: ¿una ameba usa ese movimiento a su antojo?](#p22)
- [P23 — ¿Qué significa que la nube «no es esférica ni se repolariza»? ¿Qué hacen AMOEBA y Drude?](#p23)
- [P24 — ¿Los enlaces están «hardcoded»? ¿No los rompemos nunca? ¿Y entonces qué hace el muelle angular?](#p24)
- [P25 — ¿Qué es MD? (y qué quería decir «en MD la presión es uniforme»)](#p25)
- [P26 — ¿Puedo reproducir los experimentos de F7 en casa? ¿Es peligroso lo que el modelo no sabe simular?](#p26)
- [P27 — ¿La bomba sodio-potasio contrae el músculo? ¿Se puede simular en directo, átomo a átomo?](#p27)
- [P28 — ¿Por qué importan los nanotubos? ¿Qué es el ATP y qué papel juega (en la bomba, en el músculo y en la ameba)?](#p28)
- [P29 — ¿El ADN está en el núcleo así, en esta agua con sal y plegado, tal como lo modelamos en F5?](#p29)
- [P30 — F7a–F7e entregadas: lo que se midió en `soluciones.html` (presión, sal, ósmosis, campo, Debye) y lo que enseñó construirlo](#p30)
- [P31 — F9a entregada: la física de F5 corriendo mientras miras, y por qué el mismo campo eléctrico da respuestas opuestas en una gota y en una caja periódica](#p31)
- [P32 — F9b–F9e entregadas: un canal, una bomba de juguete, la bomba de verdad y un músculo. Qué distingue a un canal de una bomba, y las tres trampas que casi arruinan el experimento](#p32)
- [P33 — F10: romper un enlace de verdad. Qué *es* un enlace covalente, y por qué la cuántica cabe en el navegador cuando la molécula tiene un solo electrón](#p33)
- [P34 — Hacia F11: qué se puede *ver emerger* del átomo (niveles, orbitales, espectro, capas, radiactividad) y qué hay que declarar (el espín y el signo de Pauli). Las librerías cuánticas que sí se instalan aquí, con sus controles](#p34)
- [P35 — F11a entregada: por qué el electrón solo puede tener ciertas energías (la onda que se dispara), qué es un paquete de onda, por qué se ensancha, y cómo atraviesa una barrera que no puede saltar](#p35)
---
## P1 (31/08/2026) — ¿Qué es la esfera de van der Waals? ¿Es como un campo magnético?
**No es magnetismo.** Un átomo no tiene superficie: es una nube de electrones que se desvanece sin borde. Pero al juntar dos átomos aparece a cierta distancia un **muro de repulsión brutal**. El *radio de van der Waals* es "dónde está el muro": el radio de contacto efectivo, **medido** en cristales (C 1,70 Å · N 1,55 · O 1,52 · P 1,80 · H 1,20 — los del visor).
El muro tiene dos orígenes, y ninguno es magnético:
1. Una parte menor de **repulsión eléctrica** entre nubes (ambas negativas).
2. La parte dominante: el **[principio de exclusión de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)** (cuántica pura). Si dos nubes se solapan, sus electrones tendrían que compartir estados — prohibido — y algunos deben saltar a energías mucho mayores. Ese coste crece de forma explosiva y *se comporta* como fuerza repulsiva. En el [campo de fuerzas](Del_atomo_al_cromosoma.html#g-campo-de-fuerzas)[(1)](Del_atomo_al_cromosoma.html#campo)[(2)](#p19) es el término (σ/r)¹² de [Lennard-Jones](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](#p19)[(2)](sandbox.html).
Es también la respuesta a por qué la materia es sólida al tacto: nunca "tocas" nada — sientes el muro de Pauli.
**El magnetismo casi no pinta nada aquí**: los electrones van emparejados con espines opuestos y se cancelan (solo los radicales libres son magnéticos, y dañan el ADN por química, no por magnetismo). El magnetismo nuclear es tan débil que solo sirve para observar (RMN).
**No confundir las dos "auras" de un átomo** (serán capas distintas del simulador):
- **Esfera vdW** (corto alcance): el volumen duro de contacto — lo que se ve en modo esferas.
- **Campo eléctrico** de sus cargas parciales (largo alcance, apantallado por el medio): lo que pintaremos en F4 con [isosuperficies](Del_atomo_al_cromosoma.html#g-isosuperficie)[(1)](campos.html) y líneas de campo.
Remate: el muro (r⁻¹²), el pegamento (dispersión r⁻⁶, la que sostiene el [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice) de las bases) y el campo eléctrico salen todos de la **misma nube de electrones**, comportándose de tres maneras según la distancia.
---
## P2 (31/08/2026) — ¿Por qué ver el 1BNA (F1) antes de validar la física (F2)?
Porque **F1 no contiene física**: es *abrir una foto* — leer coordenadas medidas por rayos X y pintar esferas. F1 construye la **ventana** (parseo, cámara, render); F2+ construye el mundo.
El principio de ingeniería: **nunca depurar dos incógnitas a la vez.** Con F1 primero, cualquier fallo visual es del visor (no hay física que pueda fallar); al llegar a F2, cualquier rareza es física (la ventana ya está verificada). Cada fase introduce una sola incógnita nueva.
Además, cada fase lleva su **control experimental** (la validación que preguntabas): F2 el agua (densidad, ~3 puentes de H por molécula), F3 la geometría y energía del par A·T, F5 el RMSD contra el cristal [1BNA](visor.html)[(1)](Del_atomo_al_cromosoma.html#helice)[(2)](trayectoria.html), F7 las curvas de fusión ([Tm](Del_atomo_al_cromosoma.html#g-tm)[(1)](Del_atomo_al_cromosoma.html#medio)) reales. La cadena completa está en la sección 10 de `ADN_simulador_fisica_atomica.md`.
---
## P3 (31/08/2026) — ¿Qué son PME, SHAKE y los «diedros bsc1» que le faltan a nuestro motor?
Las tres piezas por las que en F5 usaremos OpenMM en vez del sandbox:
- **[PME](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared) (Particle Mesh Ewald)** — la fuerza de Coulomb decae como 1/r: tan despacio que no se puede "cortar a 1 nm" sin romper la física (cargas que aparecen/desaparecen en el borde → artefactos), y sumar todos los pares es O(N²). PME parte el campo en dos: lo cercano se suma directo, y lo lejano se resuelve **sobre una rejilla de voxels con transformadas de Fourier** (como calcular la marea con una malla, no molécula a molécula). Nuestro sandbox esquiva el problema porque la gota es pequeña y sin bordes periódicos.
- **SHAKE** — la vibración O–H/N–H (período ~10 fs) obliga a pasos de ~1 fs o el integrador explota. Pero esas vibraciones están **congeladas cuánticamente** a 37 °C (su cuanto de energía ≫ kT: en la realidad ni oscilan). SHAKE las elimina fijando las longitudes X–H como bielas rígidas → paso de 2 fs, mitad de coste, y más fiel, no menos. Nuestro sandbox usa muelles flexibles con dt = 0,4 fs: más simple de entender, más caro de calcular.
- **Diedros bsc1/OL15** — los parámetros de torsión del esqueleto del ADN (α…ζ, χ) son química destilada en números. Con los de los años 90, las simulaciones de más de ~10 ns **se degradaban solas** (giros espurios α/γ deformaban la hélice B). parmbsc0 (2007), parmbsc1 (2015) y OL15 los recalibraron contra cuántica + RMN + cristalografía. Usar OpenMM es heredar décadas de esa calibración.
- **La receta de F5 completa**: [1BNA](visor.html)[(1)](Del_atomo_al_cromosoma.html#helice)[(2)](trayectoria.html) + hidrógenos → caja de ~7.000 aguas + 22 Na⁺ + NaCl 150 mM → minimizar → calentar a 310 K → 100+ ns de producción (≈5×10⁷ pasos; una noche de GPU) guardando foto cada 10 ps → 10.000 fotogramas DCD/XTC → reproducidos en nuestro visor con línea de tiempo → medir RMSD contra el cristal y ver si aparece sola la espina de hidratación.
---
## P4 (31/08/2026) — Las cinco fuerzas del sistema, puestas en la misma regla (kJ/mol)
La regla de medir: la **agitación térmica a 37 °C vale kT ≈ 2,6 kJ/mol**. Lo que ronda ese valor se rompe y reforma sin parar; lo que está muy por encima es permanente salvo enzima mediante.
| Fuerza | Energía típica | Qué es | Papel |
|---|---|---|---|
| **Covalente** (C–C, C–N, P–O, O–H) | 300–600 | compartición de electrones | esqueletos y enlaces: irrompibles a T corporal |
| **Puente de hidrógeno** | 4–30 (neto en agua, mucho menos) | dipolo–dipolo direccional D–H···A | especificidad: A·T (2), G·C (3), red del agua |
| **Apilamiento π + van der Waals** | 2–8 por escalón | dispersión de London (dipolos inducidos) | **el principal estabilizador de la hélice** |
| **Electrostática iónica** | variable, apantallada (Debye ~0,8 nm) | Coulomb entre cargas netas | fosfatos(−) se repelen; histonas(+) los atrapan |
| **Efecto hidrofóbico** | ~3 por CH₂ | entropía del agua | esconde las bases dentro de la hélice |
En el **sandbox F2** estas fuerzas son literalmente cinco términos de código en `fuerzas()`: LJ (Pauli + dispersión), Coulomb (cargas parciales), muelle de enlace, muelle angular y pared. El [efecto hidrofóbico](Del_atomo_al_cromosoma.html#fuerzas)[(1)](#p19) no es un término: **emerge** de los otros (y se vio: la gota se cohesiona sola).
La gran moraleja: **la vida funciona por separación de escalas** — información en enlaces covalentes (estables), estructura en interacciones débiles (dinámicas, reversibles, sensibles a T).
---
## P5 (01/09/2026) — ¿Qué es la temperatura para átomos y moléculas, y por qué «evapora» la gota del sandbox?
### La temperatura no es una cosa: es una estadística
Un átomo **no tiene temperatura**; tiene velocidad. La temperatura es una propiedad del *colectivo*: mide la **anchura de la distribución de velocidades[(1)](Del_atomo_al_cromosoma.html#motor)[(1)](Del_atomo_al_cromosoma.html#motor)[(1)](Del_atomo_al_cromosoma.html#motor)[(1)](Del_atomo_al_cromosoma.html#motor)** (distribución de Maxwell-Boltzmann). La definición exacta es equipartición: la energía cinética media por cada dirección de movimiento vale ½kT. En el sandbox esto es literal — la función `temperatura()` suma la energía cinética de los 300 átomos y despeja T de esa fórmula. Por eso T fluctúa unos grados alrededor del objetivo: es una media sobre solo 100 moléculas.
Números para calibrar la intuición: a 300 K una molécula de agua se mueve a ~**650 m/s de media** (velocidad de bala) — pero recorre menos de 1 Å entre choque y choque, por eso no lo parece: es el *temblor browniano* que ves en el sandbox. Subir de 300 a 500 K solo multiplica las velocidades típicas por √(500/300) ≈ 1,3.
### Dónde vive T en el código
En dos sitios, ninguno es "un calefactor":
1. **Velocidades iniciales**: se sortean de la distribución de Maxwell a la T objetivo.
2. **El [termostato de Langevin](Del_atomo_al_cromosoma.html#g-langevin)[(1)](Del_atomo_al_cromosoma.html#motor)**: a cada átomo, en cada paso, fricción (saca energía) + empujones aleatorios (la meten), calibrados por el teorema de fluctuación-disipación para que el balance se estacione exactamente en T objetivo. Pulsar ↑ no "calienta" nada directamente: **recalibra el tamaño de los empujones**, y la energía fluye hasta el nuevo equilibrio. (El termostato representa el baño térmico infinito que no simulamos — es parte de la "pared".)
### El punto conceptual clave: la temperatura NO cambia las fuerzas
Ni un término de `fuerzas()` contiene T. A 50 K y a 500 K, dos moléculas a la misma distancia sienten **exactamente la misma atracción**. Lo que cambia es el **resultado de la competición**:
> **energía cinética (agitación, ~kT) contra energía potencial (los pozos que atrapan, fijos)**
- Un puente de hidrógeno del sandbox es un pozo efectivo del orden de ~8 kJ/mol.
- kT vale: 0,4 (a 50 K) · 2,1 (250 K) · **2,5 (300 K)** · 3,1 (373 K) · 4,2 (500 K) kJ/mol.
- La probabilidad de que una molécula reúna energía para saltar fuera de un pozo de profundidad ΔE va como el **factor de Boltzmann, e^(−ΔE/kT)** — exponencial. Para ΔE = 8 kJ/mol: e^(−8/2,5) ≈ 4% de los intentos a 300 K, e^(−8/4,2) ≈ 15% a 500 K. Casi 4× más escapes… y hay miles de "intentos" (choques) por nanosegundo. Un cambio moderado de T produce un cambio drástico de comportamiento **porque está dentro de una exponencial**.
### La evaporación, vista de cerca
La distribución de Maxwell tiene **colas**: siempre hay alguna molécula anormalmente rápida. Si está en la superficie de la gota y su velocidad apunta hacia fuera, escapa de la atracción colectiva: eso **es** evaporar — la cola de la distribución goteando hacia fuera. Por eso evaporar enfría (se van las rápidas y la media baja… pero nuestro termostato repone). En el sandbox la pared las devuelve, así que a T alta se forma un equilibrio gota↔vapor dentro de la burbuja — exactamente lo que vimos a 500 K.
Y el hielo es el otro extremo: a T baja casi ninguna molécula puede saltar de su pozo → la red de puentes se congela y solo queda vibración. Frío = gana el potencial (orden); calor = gana la cinética (libertad). Los **cambios de fase** no están programados: emergen de esta competición.
Hay una segunda cara, más fina: T también decide **cuánto pesa el desorden** (la entropía) en el balance (F = E − TS). El estado evaporado tiene muchísimas más configuraciones posibles; a T alta ese peso gana aunque cueste energía. El calor no rompe la red "empujando": hace que el desorden salga a cuenta.
### Por qué la gota nunca se "descompone" en H y O
El mismo factor de Boltzmann con la tabla de P4: romper un enlace covalente O–H (~450 kJ/mol) a 500 K tiene probabilidad e^(−450/4,2) ≈ 10⁻⁴⁷ por intento — no ocurre ni en la edad del universo. **La separación de escalas de P4 en acción**: la agitación térmica deshace puentes (≈kT) pero jamás toca lo covalente (≈200·kT). Por eso hervir agua da vapor de H₂O y no hidrógeno y oxígeno.
### Datos medidos en nuestro propio sandbox (31/08/2026)
| T objetivo | T medida | Puentes de H (100 H₂O) | Estado |
|---|---|---|---|
| 300 K | 294–296 K | 145–149 (~2,9/molécula) | gota líquida cohesionada |
| 500 K | 500 K | 100 | expansión, evaporación parcial |
| 250 K | 246 K | 162 | recondensación, red densa |
Experimentos sugeridos: baja a 50 K (↓) y verás la red inmóvil, solo vibrando (vidrio/hielo amorfo); a 300 K cronometra un puente cian — vive ~1 ps antes de recablearse. Aviso honesto: con 100 moléculas, pared y modelo flexible, las T de transición no caerán exactamente en 273/373 K — el *comportamiento* es el real, los números finos no tienen por qué.
---
## P6 (01/09/2026) — ¿Los átomos «orbitan» como los planetas, por el balance cinética/potencial?
El instinto es correcto: el equilibrio cinética-contra-pozo que hace orbitar planetas es el mismo balance que estabiliza la materia. Tanto, que **Bohr (1913) modeló el átomo exactamente así** (electrón-planeta, Coulomb en vez de gravedad)… y falla por la cuántica: el electrón no es una bolita, y una carga orbitando radiaría y caería al núcleo en ~10⁻¹¹ s. Lo que sostiene de verdad al átomo: comprimir la nube electrónica eleva su energía cinética *cuántica* (Heisenberg: más confinado = más agitado); el tamaño del átomo es el mínimo del tira-y-afloja **presión de Heisenberg vs atracción de Coulomb**. El balance existe — la órbita no.
¿Qué «orbita» entonces en el mundo molecular?
- **Diferencia clave de pozos**: la gravedad es atracción 1/r **sin suelo** → la única estabilidad posible es orbitar. El pozo interatómico (LJ) **tiene suelo** ([muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)) → la estabilidad es **oscilar alrededor de la distancia de equilibrio**: la *vibración*. Un enlace O–H es un sistema planetario con suelo: no orbita, bota.
- La **rotación de la molécula entera** (los H barriendo alrededor del O en un agua que gira) es lo más parecido a un mini-sistema orbital.
- Existen «resonancias orbitantes» en colisiones: dos átomos que se capturan transitoriamente y dan unas vueltas antes de soltarse.
---
## P7 (01/09/2026) — Choques, disipación… ¿y si todos los átomos se movieran alineados, sería temperatura?
**Lo confirmado del modelo propuesto** (que era casi todo correcto):
- Los átomos no se paran solos: en el fondo **no hay fricción**; fricción y viscosidad son *emergentes* (energía ordenada fugándose a movimiento desordenado).
- «Entregar velocidad a los átomos externos» = **conducción del calor**: intercambio choque a choque hasta igualar la distribución. (Canal extra que el sandbox no tiene: **radiación** infrarroja — las cargas aceleradas emiten; así enfría lo que está en el vacío y así ve una cámara térmica.)
- Las moléculas guardan energía en **traslación + rotación + vibración**; equipartición: ½kT por modo. (Capacidad calorífica = contar modos.)
**El matiz corregido**: a nivel atómico los choques son **elásticos** — la energía total se conserva siempre, convirtiéndose cinética↔potencial contra el [muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones). «Inelástico» no existe en el fondo: un choque inelástico macroscópico es energía ordenada desordenándose. La inelasticidad *es* el calentamiento visto desde arriba.
**La respuesta a la pregunta estrella**: NO — y la pista estaba en la propia pregunta («sumar y ver el *resto*»). La temperatura se mide **tras restar el movimiento colectivo** (en el marco del centro de masas): la dispersión *residual* es lo térmico. Cien moléculas alineadas a 650 m/s tienen **T = 0**: eso es viento, no calor (la ISS va a 7,7 km/s y no está caliente por ello). Movimiento **ordenado** = energía mecánica; **desordenado** = temperatura.
El tráfico entre ambas es de sentido único:
- Orden → desorden: facilísimo (el haz choca y se aleatoriza). Por eso arde la reentrada, queman los frenos, calienta martillear.
- Desorden → orden: espontáneamente jamás — no lo prohíbe Newton, lo prohíbe la probabilidad (10²³ velocidades no se alinean solas). Ese «jamás estadístico» es la **segunda ley de la termodinámica**; la entropía no es una fuerza, es contar casos.
**Conexión con nuestra simulación**: el artefacto famoso *flying ice cube* — termostatos antiguos fugaban energía al movimiento colectivo y el sistema acababa volando/girando como un todo, congelado por dentro, mientras el medidor (que no restaba la media) juraba 300 K. Nuestra `temperatura()` del sandbox tiene ese defecto teórico (no resta el centro de masas); nos salva que la fricción de [Langevin](Del_atomo_al_cromosoma.html#g-langevin)[(1)](Del_atomo_al_cromosoma.html#motor) frena la deriva. Pendiente F3: restar el centro de masas y mostrar en el HUD «velocidad colectiva vs T interna».
---
## P8 (01/09/2026) — ¿De dónde sale el puente de hidrógeno? ¿Por qué no hay «puentes de hierro»?
**No es una fuerza fundamental — y el sandbox F2 lo demostró**: en `fuerzas()` no hay término «puente de H» (solo Coulomb + LJ + muelles) y aun así emergieron 145. Es una **circunstancia** de la única fuerza de la química (electromagnetismo) bajo dos reglas cuánticas (Pauli + reparto de electrones en orbitales/electronegatividad). Nada nuevo que creerse.
**Por qué el H y solo el H**:
1. Un único electrón y **ninguna capa interna debajo**: enlazado a un átomo muy electronegativo (O 3,44 · N 3,04 · F 3,98 vs H 2,20) queda por la cara opuesta como **protón casi desnudo**. Cualquier otro elemento desplumado conserva capas internas que lo mantienen grande y apantallado.
2. **Diminuto** → el par solitario del aceptor se arrima a ~1,8–2 Å; con Coulomb ~q/r², poca distancia = mucha fuerza (10–40 kJ/mol ≈ kT×4–15: rompible pero pegajoso, la escala de la vida).
3. **Direccional**: el aceptor ofrece un par solitario (orbital que apunta) → geometría casi lineal D–H···A → especificidad (A·T, G·C).
**Por qué no el Fe**: es electropositivo (1,83) — en Fe–H la polaridad se invierte (el H queda −, hidruro); un metal nunca queda como protón desnudo. Su versión del puente es la **química de coordinación** (Fe²⁺ agarrando O₂ en la hemoglobina; Mg²⁺ coordinando fosfatos en las polimerasas).
**Sí existen primos en la tabla periódica** (prueba de que es circunstancia, no fuerza): enlace de **halógeno** (agujero σ positivo en Cl/Br/I; IUPAC 2013, diseño de fármacos), enlaces de **calcógeno/pnictógeno/tetrel** (misma jugada con S/Se, P/As, C/Si), **puentes salinos** (cargas completas). El de H es el campeón porque solo él carece de capas internas.
**La pizca cuántica de verdad**: descompuesto con SAPT, ~2/3 es electrostática clásica (por eso emerge con cargas puntuales), más polarización y una pizca de **transferencia de carga** (el par solitario se derrama en el σ* del O–H) que afina la direccionalidad y hace los puentes **cooperativos** (cadenas que se refuerzan — el hielo). Los modelos de cargas fijas pierden parte de esa cooperatividad: límite honesto del sandbox.
**Moraleja**: en química hay UNA fuerza y DOS reglas; covalente, vdW, puente de H, halógeno, coordinación… son regímenes de la misma fuerza.
---
## P9 (01/09/2026) — ¿El universo es EM + reglas cuánticas + un «primer movimiento»? ¿La gravedad es emergente como la temperatura?
**El inventario completo son CUATRO interacciones** (EM, nuclear fuerte, nuclear débil, gravedad). Que a la química le baste la EM es separación de escalas otra vez: la **fuerte** pega los núcleos (sin ella, solo hidrógeno — la tabla periódica es regalo suyo) a ~10⁸ kJ/mol, un millón de veces la escala química → el núcleo es un «dado congelado», como lo covalente para lo térmico. La **débil** solo transmuta (radiactividad, fusión solar). Y «orbitales/electronegatividad» no es regla aparte: se deriva de la mecánica cuántica (Schrödinger + Coulomb + Pauli). Inventario limpio de la química: **EM + mecánica cuántica (con Pauli) + estadística**, núcleos de regalo.
**La gravedad, en la física establecida, NO es un resto del EM**: es interacción fundamental aparte — y ni siquiera «fuerza»: curvatura del espacio-tiempo (relatividad general). A escala atómica es ~10³⁹ veces más débil que la eléctrica (el sandbox la ignora con razón); domina el cosmos porque el EM **se autoanula** (+/− se cancelan, la materia es neutra) y la gravedad solo suma. Reparto: EM gobierna lo pequeño, gravedad lo grande; la frontera es la neutralidad.
**PERO** la intuición «emergente y estadística como la temperatura» es una hipótesis seria de frontera: la **gravedad entrópica** (Jacobson 1995, Verlinde 2010) — minoritaria, no demostrada, con problemas abiertos. No hay aún teoría cuántica de la gravedad que lo zanje.
**El big bang, contado al revés**: no metió «primer movimiento» en átomos preexistentes — no había átomos, y el universo nació máximamente caliente (el movimiento venía de serie). La expansión lo **enfrió**, y la historia cósmica es la competición de Boltzmann de P5 en cascada — cada ligadura congela sus estructuras cuando kT cae por debajo de su ΔE:
- ~10⁻⁵ s: quarks → protones/neutrones (fuerte)
- ~3 min: núcleos de H y He
- ~380.000 años: **átomos** (kT < 13,6 eV; la luz liberada es el fondo cósmico de microondas)
- ~10⁸ años: la gravedad amontona gas → estrellas, que forjan C/N/O/P (y Fe en supernovas)
- planetas a ~300 K: covalente estable + puentes rompibles-pegajosos → química, agua, ADN, vida
La biología vive en el último escalón frío de la cascada; la gota del sandbox a 300 K es un fósil del enfriamiento del universo.
**¿De dónde sale el orden si el desorden gana (P7)?** De la gravedad + expansión: agruparse y radiar calor al vacío sube la entropía total aunque nazcan estructuras locales. Las estrellas quedan como fuentes de energía ordenada; la vida procesa ese chorro (orden solar entra, infrarrojo sale). La segunda ley no se viola: se surfea.
**Inventario final del universo (honesto)**: 4 interacciones + mecánica cuántica (Pauli incluido) + estadística + condiciones iniciales (energía, expansión, fluctuaciones primordiales). La química del proyecto es el sector EM de ese todo.
---
## P10 (01/09/2026) — ¿El CERN acelera átomos para superar la barrera de Pauli y ver qué pasa?
La intuición («velocidades brutales para superar una barrera y ver qué pasa») es correcta — pero la barrera no es la de Pauli. **La escalera de energías por barrera**:
| Energía por partícula | Barrera que se supera | Quién lo hace |
|---|---|---|
| ~0,03 eV (kT a 300 K) | puentes de H | tu sandbox |
| 1–10 eV | **[muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)** entre nubes y enlaces covalentes | una llama, cualquier química |
| 10 eV – keV | arrancar electrones (ionizar) | plasmas, el Sol |
| ~MeV | **barrera de Coulomb entre núcleos** (los + se repelen) | fusión: estrellas, tokamaks |
| TeV (10¹² eV) | ninguna: **crear materia desde energía** (E=mc²) | LHC del CERN |
O sea: el muro de Pauli lo «supera» una vela — es escala química. Los aceleradores usan partículas cargadas (protones, iones, electrones — no átomos neutros: se aceleran con campos EM) para (1) saltar la barrera de Coulomb nuclear y transmutar núcleos, (2) a lo TeV, **convertir energía cinética en partículas nuevas** (así se encontró el Higgs), y (3) usar longitudes de onda diminutas como microscopio (más momento = más resolución). En colisiones Pb-Pb, el LHC llega a **fundir protones en plasma de quarks y gluones**: recrea el primer escalón de la cascada de P9. Un acelerador es la cascada de enfriamiento del universo **corriendo hacia atrás, localmente**.
---
## P11 (01/09/2026) — ¿Una gota de agua en el espacio, sin gravedad ni atmósfera, mantendría su forma de bola?
**Sí — y el sandbox F2 casi lo es.** La cohesión (puentes de H) mantiene el líquido unido y la **tensión superficial** (emergente, la vimos) minimiza la superficie → esfera. Es lo que hacen las bolas de agua flotantes de los astronautas en la ISS. Tres matices importantes:
1. **La deriva del centro de masas del sandbox es un artefacto del termostato**: los empujones de [Langevin](Del_atomo_al_cromosoma.html#g-langevin)[(1)](Del_atomo_al_cromosoma.html#motor) no conservan el momento (son un baño externo). En el espacio real, aislada, el momento **se conserva exactamente**: los choques internos no pueden mover el centro de masas (fuerzas internas se cancelan por pares) → la gota viaja en línea recta a velocidad constante (Newton I), sin paseo aleatorio.
2. **En vacío real la gota se evapora sin retorno**: no hay pared que devuelva a las de la cola de Maxwell (P5). Y al fugarse las rápidas, el resto se enfría → la gota **hierve y se congela a la vez**: acaba siendo una bolita de hielo menguante. Nuestra pared mock existe precisamente para evitar esto.
3. **A lo grande, cambia el escultor de esferas**: en gotas manda la tensión superficial; en masas enormes manda la **autogravedad** — que también hace esferas (por eso los planetas son redondos). Dos fuerzas distintas, misma forma; el relevo ocurre a escalas de cientos de km (lunas heladas).
---
## P12 (01/09/2026) — ¿Los átomos «decaen» emitiendo luz hasta desaparecer? ¿Qué es entonces el fondo de microondas?
Eran **dos cosas distintas** que convenía separar:
**1. La radiación térmica no gasta átomos: gasta calor.** Un átomo/molécula agitado emite infrarrojo porque sus cargas se aceleran (choques, vibraciones) — al emitir pierde energía *cinética*, no sustancia: **se enfría, no se consume**. Y la cuántica pone el suelo: al llegar al estado fundamental ya no puede radiar más (la solución al problema de Bohr de P6 — el estado fundamental es estable). La materia no desaparece radiando: la radiación es un canal de enfriamiento (el que le falta a nuestro sandbox, que usa [Langevin](Del_atomo_al_cromosoma.html#g-langevin)[(1)](Del_atomo_al_cromosoma.html#motor) en su lugar).
**2. El fondo cósmico de microondas (CMB) no es átomos muriendo: es la primera luz liberada, estirada.** Antes de los 380.000 años el universo era un plasma opaco: la luz chocaba sin parar con los electrones libres (niebla luminosa). Al enfriarse lo bastante para que los electrones se quedaran en los átomos (recombinación, P9), la niebla se levantó de golpe: el universo se volvió **transparente**, y los fotones que volaban en ese instante quedaron libres — y siguen volando desde entonces, sin que nada los haya absorbido. Eran luz visible/naranja de un brillo a ~3.000 K.
**¿Por qué microondas hoy?** La **expansión estira las longitudes de onda** (corrimiento al rojo cosmológico): factor ~1.100 desde entonces → de ~1 µm a ~1 mm; el «gas de fotones» pasó de 3.000 K a 2,7 K. Los fotones tampoco desaparecieron: se estiraron con el espacio que atraviesan.
**Nada desaparece en ninguna de las dos historias**: la materia se volvió gas transparente (y luego estrellas, y tú); la luz siguió su viaje, cada vez más roja. El CMB es la foto más antigua posible del universo — el fogonazo del momento en que se levantó la niebla, llegándonos desde todas las direcciones a la vez.
---
## P13 (01/09/2026) — ¿Por qué las bases son rígidas en F3? ¿Podría su forma mantenerse sola, por las fuerzas?
**La escalera de tres niveles para «tener forma»:**
1. **Rígido (F3)**: la forma es un dato de entrada. Justificación física: los anillos aromáticos SON placas tiesas — sus vibraciones internas tienen amplitudes de ~0,05 Å (covalente: 300–600 kJ/mol vs kT 2,5) y las X–H están congeladas cuánticamente (P3). Lo que se pierde honestamente: deformaciones inducidas pequeñas pero reales (el *propeller twist* de ~15° de los pares en el ADN, la respiración del anillo, el giro del grupo amino).
2. **Flexible con muelles** — «que las fuerzas hagan el trabajo» *es exactamente esto y es el estándar de la MD*: enlaces + ángulos + **diedros** + **impropios** (los términos de planaridad) como muelles parametrizados, y la forma se mantiene sola… si los muelles están bien calibrados. Nuestra agua TIP3P ya es de este nivel. El coste computacional es casi nulo (los términos enlazados son baratos; el coste sigue en los pares no enlazados); **el coste real es la calidad de los parámetros**: sin torsiones/impropios bien afinados la molécula adopta formas incorrectas — la historia de décadas de bsc1 (P3). En F5, OpenMM nos dará las bases flexibles con los muelles refinados oficiales (y para un sandbox F3b se podrían importar parámetros GAFF).
3. **Forma verdaderamente emergente, sin muelles**: imposible en clásica — el enlace covalente es compartición de electrones (cuántica: intercambio, deslocalización) y **no se puede representar con fuerzas clásicas por pares**; con solo LJ+Coulomb la molécula se desarma o colapsa. Para que la forma emerja de verdad hay que resolver electrones: DFT/dinámica ab initio (~10⁶× más cara: minutos por femtosegundo para 30 átomos) o el camino intermedio moderno, **potenciales de aprendizaje automático** (ANI/MACE): redes entrenadas con millones de cálculos cuánticos que mantienen las moléculas unidas emergentemente — hasta pueden romper enlaces — a ~100–1.000× el coste clásico. OpenMM 8 los integra: es viable para una base en GPU.
**Sobre «¿los enlaces internos no se ven afectados por el resto de fuerzas?»**: en el nivel 2 sí responden — un tirón externo estira ligeramente los muelles — pero el efecto es diminuto porque la rigidez covalente aplasta a las fuerzas intermoleculares (200:1). Por eso el nivel 1 es tan buena aproximación… hasta que te importan los detalles finos.
**Experimento de la competencia (mismo día)**: `bases.html?agua=40` — en vacío A·T se capturaron a t=163 ps y jamás se soltaron (E=−55 kJ/mol); con 40 aguas, a t=154 ps seguían sin un solo puente A·T (llegaron a 7,3 Å a t=114 y ¡rebotaron!), cada base abrigada por su capa de solvatación (E base·agua = −364 kJ/mol, ~62 puentes con agua en la gota). **En agua, emparejarse casi no renta**: el agua paga puentes equivalentes. Es la razón profunda de que la biología necesite el andamio del ADN (esqueleto + [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice)) para hacer fiable el apareamiento — y de que la doble hélice pueda abrirse y cerrarse.
---
## P14 (01/09/2026) — ¿Qué significa que las distancias atómicas «se miden en cristales»?
**Un cristal** es un sólido cuyos átomos/moléculas se repiten en una red 3D perfecta (naranjas apiladas). Se forma cuando el enfriamiento da tiempo a que cada molécula encuentre su encaje de mínima energía — el potencial ganando a kT *con orden*. (Enfriado brusco → vidrio: sólido desordenado, un líquido congelado; como nuestra gota de F2 a 50 K.)
**Cristalografía de rayos X**: los rayos X tienen longitud de onda ~1 Å — la escala interatómica — y al atravesar la red periódica difractan: el cristal actúa de rejilla de difracción 3D y la interferencia constructiva (Bragg) concentra la luz dispersada en puntos discretos. Del patrón de puntos e intensidades se reconstruye el mapa de densidad electrónica → posiciones atómicas con precisión de centésimas de Å. **El cristal es el amplificador**: una sola molécula dispersa demasiado poco, pero 10¹⁵ copias perfectamente alineadas suman su señal en fase. Así se midió el [1BNA](visor.html)[(1)](Del_atomo_al_cromosoma.html#helice)[(2)](trayectoria.html) (Dickerson 1981, «rayos X, 1,9 Å» del visor F1) y casi todo el PDB. (Lo moderno: crio-EM promedia millones de imágenes de moléculas sueltas — sin cristal.)
**Los radios de van der Waals «medidos en cristales»**: se revisan miles de estructuras cristalinas y se hace el histograma de las distancias de contacto entre átomos *no enlazados* de moléculas vecinas. Dos carbonos jamás se acercan a menos de ~3,4 Å → r(C) = 1,7 Å. Es la posición del [muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones), medida estadísticamente sobre la naturaleza.
**La conexión con la temperatura (preciosa)**: los átomos del cristal siguen vibrando (kT) y eso emborrona la foto. Los cristalógrafos miden ese emborronamiento átomo a átomo: es el **factor B** — y está en nuestro propio fichero: el último número de cada línea ATOM de `data/1BNA.pdb` (p. ej. el 64.35 del O5' inicial). La temperatura, literalmente, desenfoca la foto; por eso hoy se cristaliza en frío (crio, 100 K).
---
## P15 (01/09/2026) — Estados de la materia, romper enlaces con calor… y qué es kT
**Corrección primero**: la temperatura NO es una distancia media — es la **anchura de la [distribución de velocidades](#p5)[(1)](Del_atomo_al_cromosoma.html#motor)** (½m⟨v²⟩ por eje = ½kT; P5/P7). La densidad (distancias) es otra variable independiente: hay gas caliente diluido y líquido frío denso.
**¿Se va un átomo y queda la misma T?** T es **intensiva**: no depende de cuántos átomos midas — la mitad de átomos con la misma distribución = misma T. PERO la evaporación no se lleva átomos al azar: **se lleva los más rápidos** (la cola de Maxwell) → la distribución restante queda más estrecha → **la T baja**. Es el *enfriamiento evaporativo*: por eso sudar enfría y por eso soplas la sopa — y es la técnica con que los laboratorios llegan a nanokelvin (evaporar lo más caliente de una trampa de átomos). En F2 no se ve porque el termostato repone lo evaporado energéticamente.
**La escalera de estados = la competición kT contra pozos, por niveles**:
| Estado | Competición | Agua |
|---|---|---|
| **Sólido** | kT ≪ pozos intermoleculares: presos en sus sitios, solo vibran; con orden de largo alcance = **cristal** | hielo, <273 K |
| **Líquido** | kT ≈ pozos: sueltan y recapturan vecinos sin parar; cohesión sin sitios fijos → fluye | 273–373 K |
| **Gas** | kT > cohesión: libres; interacción solo en choques breves | vapor |
| **(Disociación)** | los choques rompen **enlaces covalentes** de la molécula | 2.500–3.500 K: H₂O → H + OH |
| **Plasma** | kT ~ energías de ionización: los choques arrancan **electrones** | >10⁴ K; el Sol; el universo pre-CMB (P12) |
(Las fronteras dependen también de la presión/densidad; y sí: el estado no es propiedad de la molécula sino del colectivo y su competición.)
**¿El calor puede «disolver» enlaces covalentes? SÍ — y el razonamiento propuesto es la ecuación de Arrhenius re-derivada**: probabilidad por choque ~ e^(−Ea/kT) (Boltzmann), multiplicada por la frecuencia de intentos (~10¹³ vibraciones/s) y por un factor de orientación (hace falta chocar en la dirección adecuada, con los enlaces en el estado vibracional propicio). Eso es exactamente k = A·e^(−Ea/kT), la fórmula maestra de la cinética química. Números para el O–H (~460 kJ/mol): a 300 K, e^(−460/2,5) ≈ 10⁻⁸⁰ → jamás; a 3.000 K, kT=25 → e^(−18) ≈ 10⁻⁸ × 10¹³ intentos/s → se rompe en microsegundos. La termólisis del agua es real a partir de ~2.200 K. Por eso a 300 K lo covalente es «solo lectura» (la información genética, estable) — y por eso la química viva no sube la T sino que **baja la Ea**: catalizadores y enzimas.
**Y kT no son kiloteslas** 😄 — el tesla (T, campo magnético) no pinta nada aquí. En kT: **k = constante de Boltzmann** (k_B = 0,0083145 kJ/mol/K — «cuánta energía por grado») y **T = temperatura absoluta en kelvin** (0 K = agitación cero; 273 K = 0 °C; 300 K ≈ 27 °C; 310 K = tu cuerpo). El producto **kT es una ENERGÍA**: el tamaño típico del empujón térmico — 2,5 kJ/mol a 300 K. Es el conversor de divisas entre temperatura y energía, y por eso toda comparación del proyecto es «X contra kT». (Los físicos lo escriben por molécula en julios, 4,1×10⁻²¹ J; los químicos por mol, ×Avogadro. Misma cosa.)
---
## P16 (02/09/2026) — ¿La simulación F5 se ve en vivo o se graba? ¿Los fotogramas se conservan o se pierden al cerrar?
**Confirmación primero, porque el planteamiento era exacto**: sí — en F5 solo metemos las leyes (Coulomb, Pauli/LJ, muelles calibrados, [PME](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared), termostato) y las condiciones del medio (agua, sal, 310 K, 1 bar); que la doble hélice **siga plegada** no está impuesto en ningún sitio: es la *consecuencia* — puentes de H, [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice), apantallamiento — y el **RMSD contra el cristal** es el juez que lo certifica (≈1,6 Å en la prueba: se sostiene).
**Los dos modos de ver simulación que tiene el proyecto**:
1. **En directo** (sandboxes F2–F4): la física se calcula mientras miras, a ~0,5–1 ps por segundo de reloj — unas 10¹¹ veces más lento que la realidad. Viable porque son sistemas pequeños (cientos de átomos).
2. **Grabar y reproducir** (F5): con 18.416 átomos, la GPU produce ~50 ns/día — «verla en directo» sería mirar 2 días de pantalla para presenciar 100 ns. Así que la GPU trabaja sola, guarda un fotograma cada 50 ps, y después la película se reproduce **a la velocidad que se quiera**, adelante y atrás, con pausa y con zoom de nave. No es renunciar al simulador: es el microscopio con máquina del tiempo.
**Los fotogramas son FICHEROS y no caducan**: `trayectorias/f5_visor.bin` (para nuestro reproductor) y `f5_full.dcd` (formato estándar, VMD/ChimeraX). Generados una vez, se reproducen para siempre; cerrar el navegador o apagar el PC no pierde nada. Matiz: `trayectorias/` está fuera del git (pesa mucho y es regenerable) — si esas 2 jornadas de GPU te son valiosas, copia la carpeta a un disco de respaldo.
**Lo que NO guarda la película**: entre fotograma y fotograma (50 ps = 25.000 pasos de 2 fs) la física se calcula entera pero no se almacena — las vibraciones de femtosegundos viven y mueren entre fotos. El reproductor interpola entre fotogramas para que el movimiento se vea continuo; si algún fenómeno rápido interesara, se baja el `--stride_ps` a cambio de ficheros más gordos.
**E interrumpir/reanudar (mejora de este día)**: el checkpoint (posiciones, velocidades, hasta el estado de los generadores aleatorios) se guarda tras cada fotograma, y la *definición* del sistema queda congelada en disco (`f5_system.xml` + `f5_sistema.pdb` + `f5_ref.npz`) — sin eso, reconstruir colocaría las aguas al azar y el checkpoint no encajaría. Cortar cuando haga falta (Ctrl+C, apagar) arriesga ~1,5 min de GPU como mucho; `python f5_simular.py --resume --ns 100` continúa exactamente donde quedó, tantas veces como toque.
---
## P17 (02/09/2026) — ¿Cuánto ocuparía grabar los 100 ns con vibraciones de femtosegundos? (la cámara de alta velocidad)
**La estimación** (812 átomos en el binario del visor; 9.752 bytes por fotograma):
| Foto cada… | Fotogramas en 100 ns | Tamaño | Comentario |
|---|---|---|---|
| 50 ps (defecto) | 2.000 | **19,5 MB** | película larga: respiración, iones, extremos |
| 10 ps | 10.000 | 97,5 MB | movimiento fino |
| 1 ps | 100.000 | 975 MB | límite práctico: el reproductor carga todo en RAM |
| 100 fs | 10⁶ | 9,7 GB | ya no cabe en el navegador |
| 4 fs (vibraciones) | 2,5×10⁷ | **244 GB** | absurdo |
| 2 fs (cada paso) | 5×10⁷ | 487 GB | ni el disco ni nadie |
(El DCD completo, aparte: ~223 MB a su ritmo de 100 ps.)
**La solución correcta: cámara de alta velocidad.** Las vibraciones son estadísticamente idénticas en todo momento — no hay que grabarlas 100 ns; basta un **tramo corto a resolución de femtosegundos**: 4 ps a 4 fs/foto = 1.000 fotogramas = **9,8 MB**. Implementado como modo aparte que parte del checkpoint sin tocar la película larga:
```
python f5_simular.py --camara vibra --stride_ps 0.004 --ps 4
```
y se ve con `trayectoria.html?traj=f5_cam_vibra.bin`. Dos películas del mismo universo: la larga (50 ps/foto, los procesos lentos) y la lupa temporal (4 fs/foto, el temblor).
**Matiz honesto**: en F5, SHAKE congela a propósito las vibraciones X–H (P3: están congeladas cuánticamente en la realidad). Lo que vibra en la cámara rápida son los enlaces pesados (C–C, C–N: períodos de 20–40 fs, bien muestreados a 4 fs), los ángulos, los azúcares fruncéndose y los fosfatos basculando. Y el checkpoint no se guarda por fotograma en este modo (sería carísimo a este ritmo): cada 15 segundos de reloj.
---
## P18 (02/09/2026) — En F3, ¿a qué velocidad vemos la realidad y con qué muestreo? (y la GPU que era la equivocada)
**Velocidad de visionado** (sandbox F3/F4; paso de física Δt = 0,4 fs, con ~30–60 fotogramas de pantalla por segundo):
- **En vacío** (80 pasos/fotograma): ~1–1,9 ps de realidad por segundo de pantalla → vemos el mundo **~10¹² veces (un billón) más lento** que la realidad.
- **Con 40 aguas** (20 pasos/fotograma): ~0,4–0,5 ps/s → ~2×10¹² más lento.
- Las teclas **,/.** cambian los pasos por fotograma (hasta 200 ≈ 5 ps/s): más realidad por segundo a cambio de ver el temblor menos fluido. El HUD siempre muestra la cifra medida real («X ps/s»).
**El muestreo tiene tres capas distintas** (¡no confundirlas!):
1. **La física no se salta nada**: se integra cada 0,4 fs, todos los pasos, siempre — a diferencia de la película F5, aquí no hay huecos entre lo calculado.
2. **La pantalla enseña 1 de cada N estados**: con 80 pasos/fotograma, una foto cada 32 fs (vacío) o cada 8 fs (agua). Entre fotos la física ocurre pero no se pinta.
3. **Los medidores van a su ritmo**: detector de puentes cada 5 fotogramas (~160 fs vacío), HUD cada 10.
Comparado con F5: el sandbox es «en directo, sin huecos, mientras miras»; F5 es «grabado con huecos de 50 ps, reproducible a cualquier velocidad» — y su cámara `--camara` graba tramos cada pocos fs (P17).
**Y la GPU (cazada por el usuario en el Administrador de tareas)**: OpenMM elegía por defecto la **Intel integrada** en OpenCL. Sondeo con el sistema F5 real (18.548 átomos): Intel Graphics 52,8 ns/día · **GTX 1650 125,6 ns/día** — nuestras ejecuciones iban a 49,7: eran la Intel. Arreglado en `f5_simular.py` (enumera dispositivos, lee `DeviceName`, prefiere la discreta): **×2,5 — los 100 ns pasan de ~2 días a ~19 horas.** Moraleja doble: (1) en un portátil/equipo con dos GPUs, nunca asumas cuál trabaja — míralo en el Administrador de tareas; (2) los números «sospechosamente bajos» son pistas, no ruido.
---
## P19 (02/09/2026) — Las cinco leyes que hay en el código (y cómo se mapean a las cinco fuerzas de la química); ¿el muro de Pauli se deforma según la molécula?
Hay **dos listas de cinco** que conviene no mezclar: las cinco *fuerzas de la química* (P4: covalente, puente de H, [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice)/vdW, iónica, hidrofóbica) y los cinco *términos de `fuerzas()`* del sandbox, que son los que el ordenador calcula de verdad. Todo lo de la primera lista sale de la segunda.
**Los cinco términos del código** (`sandbox.html`, `bases.html`; F5 hereda los mismos vía AMBER más diedros):
| # | Término | Fórmula | Qué es físicamente | Ejemplo para la intuición |
|---|---|---|---|---|
| 1a | **[Lennard-Jones](Del_atomo_al_cromosoma.html#g-lennard-jones)(1)[(2)](sandbox.html), muro** | +4ε(σ/r)¹² | **Pauli**: dos nubes electrónicas no pueden solaparse | por qué la mesa aguanta la mano: nunca tocas nada, sientes el muro |
| 1b | **Lennard-Jones(1)[(2)](sandbox.html)(1)[(2)](sandbox.html)(1)[(2)](sandbox.html), pegamento** | −4ε(σ/r)⁶ | **dispersión de London** (vdW): dipolos instantáneos que se inducen mutuamente | las patas del geco pegándose al cristal; el apilamiento de las bases |
| 2 | **Coulomb** | k·qᵢqⱼ/r | atracción/repulsión entre **cargas parciales** fijas (O −0,834e, H +0,417e) | el chorro de agua del grifo se curva hacia un peine frotado |
| 3 | **Muelle de enlace** | k_b(r−r₀)² | el enlace covalente que no se rompe: oscila alrededor de r₀ | la vibración O–H a ~10 fs; el infrarrojo que absorbe el agua |
| 4 | **Muelle angular** | k_θ(θ−θ₀)² | los 104,5° del agua (y en F5, diedros para la planaridad y las torsiones) | la geometría tetraédrica/trigonal de tu orgánica, convertida en muelle |
| 5 | **Pared (mock)** | ½K(r−R)² solo fuera | la frontera del mundo simulado; en F5 la sustituye la [caja periódica](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared) (PBC) | la burbuja marrón del sandbox |
Y **dos ingredientes que no son fuerzas pero sin los que no hay nada**: (a) **Newton**, F=ma integrado con [velocity-Verlet](Del_atomo_al_cromosoma.html#g-velocity-verlet)[(1)](Del_atomo_al_cromosoma.html#motor) cada 0,4 fs (2 fs en F5 con SHAKE); (b) **Boltzmann**, en forma de [termostato de Langevin](Del_atomo_al_cromosoma.html#g-langevin)[(1)](Del_atomo_al_cromosoma.html#motor): fricción + empujones aleatorios calibrados a kT (ese es el «Bol…» de la memoria). El primero mueve; el segundo fija a qué temperatura se juega.
**Mapa química → código** (lo importante: tres de las cinco fuerzas químicas *no están programadas*):
| Fuerza química (P4) | De dónde sale en el código |
|---|---|
| Covalente | términos 3 y 4 (muelles); la topología es dato de entrada |
| Puente de H | **emerge** de 2 + 1a (Coulomb da la atracción, Pauli fija la distancia) |
| Apilamiento / vdW | término 1b (decenas de pares C···C a 3,4 Å) |
| Iónica (fosfatos, Na⁺, Cl⁻) | término 2 con cargas netas ±1; el medio la apantalla (F6, Debye) |
| Hidrofóbica | **emerge**: no es un término, es la entropía del agua que prefiere no ordenarse alrededor de lo apolar |
**Sí: el muro de Pauli[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones) «se deforma» según con quién esté enlazado el átomo.** La nube electrónica de un H unido a O o N se la lleva el vecino electronegativo; el H queda casi desnudo y su muro retrocede. En `bases.html` está literalmente en la tabla de tipos: **σ(H sobre N/O) = 1,07 Å frente a σ(H sobre C) = 2,60 Å**; O carbonilo 2,96 frente a O del agua 3,15. Por eso un puente de H llega a 1,8–2 Å (P8). Lo que el [campo de fuerzas](Del_atomo_al_cromosoma.html#g-campo-de-fuerzas)[(1)](Del_atomo_al_cromosoma.html#campo)(2) clásico *no* captura: que la nube no es esférica (anillos aromáticos más finos de canto que de plano, pares solitarios que apuntan) ni se repolariza al acercarse otra molécula (para eso, [AMOEBA](#p23)[(1)](Del_atomo_al_cromosoma.html#existe)/Drude). El muro es donde está la nube, y la nube la esculpen los enlaces.
---
## P20 (02/09/2026) — ¿Podemos simular sal disolviéndose, ósmosis con paredes virtuales y «aplicar electricidad»? ¿Dónde hace falta la cuántica para no hacer un mock?
**Respuesta corta: las tres cosas se hacen con exactamente las leyes de P19 y son práctica estándar de investigación.** La diferencia entre *mock* y *simulación* no es el nivel de detalle: es que en un mock el resultado está escrito y en la simulación solo están escritas las leyes, el resultado emerge y **se contrasta con un número medido**.
**1. Tirar sal al agua.** Na⁺ y Cl⁻ son dos líneas de parámetros cada uno (masa, carga ±1, σ, ε; los estándar son los de Joung–Cheatham 2008 para TIP3P). Con solo eso se ve: el cristal pierde primero los iones de las esquinas y aristas (menos vecinos que retenerlos), cada Na⁺ se rodea de ~6 aguas con el O hacia dentro, cada Cl⁻ de ~6–7 con los H hacia dentro, y el cristal se redondea y se deshace. [Lanaro & Patey 2015](https://pubs.acs.org/doi/abs/10.1021/jp512358s) lo hicieron con ~2.400 iones y vieron tres etapas (esquinas → disolución a ritmo fijo → final). Detalle precioso: energéticamente la disolución de NaCl es casi un empate (red ≈ 787 kJ/mol contra hidratación ≈ 769; ΔH_disolución ≈ +4 kJ/mol): **disuelve por entropía**, y el simulador lo captura porque la entropía es emergente (contar configuraciones). Límite honesto: los modelos clásicos subestiman la solubilidad (~la mitad de los 6 mol/kg reales en los mejores), y en el sandbox JS (~300 átomos, O(N²)) solo veríamos las primeras esquinas caer; la disolución entera es tarea de OpenMM y del reproductor `trayectoria.html`.
**2. Ósmosis con paredes virtuales.** La idea de «paredes que simulan la membrana» es literalmente el método de [Luo & Roux 2010](https://research.westernu.edu/luo-lab/publications/simulation-osmotic-pressure-concentrated-aqueous-salt-solutions): un potencial semiarmónico que **solo actúa sobre los iones** (el agua lo atraviesa libremente) divide la caja en compartimento salado y agua pura; la fuerza media que los iones ejercen contra la pared, dividida por su área, **es la presión osmótica[(1)](soluciones.html?exp=osmosis)[(1)](soluciones.html?exp=osmosis)[(1)](soluciones.html?exp=osmosis)[(1)](soluciones.html?exp=osmosis)**. Validación: van 't Hoff π = i·c·R·T da ~7 bar para NaCl 150 mM a 310 K (la del plasma sanguíneo: por eso un glóbulo rojo en agua pura revienta), ~27 bar para el agua de mar (por eso la ósmosis inversa necesita bombas de >30 bar).
**3. «Aplicar electricidad».** Son tres fenómenos distintos:
| Fenómeno | Qué es | ¿Clásico basta? |
|---|---|---|
| **Electroforesis / electrodiálisis** | un campo E empuja los iones (Na⁺ hacia un lado, Cl⁻ al otro); con paredes selectivas → desalación | **sí**: en OpenMM es una `CustomExternalForce` con energía −q·E·x ([cookbook](https://openmm.github.io/openmm-cookbook/dev/notebooks/cookbook/Applying%20a%20Constant%20External%20Force.html)) |
| **Electroósmosis** | los iones arrastran agua junto a paredes cargadas | **sí**, misma receta |
| **Electrólisis** | en los electrodos: H₂, Cl₂, O₂ — transferencia de electrones, enlaces que se rompen | **no**: hace falta [QM/MM](Del_atomo_al_cromosoma.html#g-qm)[(1)](#p24), potenciales reactivos o ML. Frontera del proyecto |
Matiz de escala: para que un ion cruce una caja de nanómetros en nanosegundos hay que aplicar ~0,1 V/nm, un millón de veces el campo de un gel de electroforesis (~10 V/cm). Es legítimo porque la respuesta es lineal y lo que se mide (movilidad, conductividad) coincide con lo experimental. Y una frontera cuántica útil: el modelo clásico jamás rompería el agua por muy fuerte que sea el campo; la dinámica ab initio dice que se disocia a partir de ~3,5 V/nm.
**Lo que sí exige cuántica (y por qué)**: reacciones redox en electrodos, el salto de protones (H₃O⁺ por mecanismo de Grotthuss, el pH), romper/formar enlaces covalentes, repolarización de nubes. Todo lo demás de esta pregunta — disolución, hidratación, ósmosis, migración iónica, apantallamiento — es electromagnetismo + Pauli + estadística, y **lo clásico es el estándar validado**. Candidatos de fase: **F7b** sal disolviéndose (OpenMM, reproducida en `trayectoria.html`), **F7c** paredes de Luo–Roux midiendo π contra van 't Hoff, **F7d** campo E con iones migrando.
---
## P21 (02/09/2026) — La presión: su origen atómico (paredes y gravedad), por qué es estadística como la T, y dónde vive en nuestro simulador
**Las dos intuiciones propuestas son correctas y son la misma cosa.** Presión = **flujo de momento a través de una superficie**: cuántos átomos golpean cada segundo por unidad de área, multiplicado por cuánto momento entrega cada golpe. Da igual si la superficie es una pared real, una pared imaginaria dentro del fluido, o la «tapa» que forma el peso de la columna de arriba.
- **Pared**: cada rebote transfiere 2mv⊥ a la pared; sumados, F/A = P. En el sandbox, la fuerza media de la burbuja sobre los átomos dividida por 4πR² **es** la presión de la gota (no está en el HUD: candidato a añadir).
- **Gravedad / atmósfera**: cada capa de aire debe soportar el peso de todas las de arriba, y lo hace *con choques*: la de abajo está más densa. Fórmula barométrica: n(h) ∝ e^(−mgh/kT), altura característica kT/(mg) ≈ 8 km. En una caja de nanómetros mgh/kT ≈ 10⁻¹¹: **la gravedad no pinta nada en MD** (P9 otra vez), y la presión es uniforme en toda la caja.
- **Gas ideal**: P = n·k·T. Dos mandos independientes: **densidad** (cuántos) y **temperatura** (cómo de rápidos). Temperatura = anchura de la [distribución de velocidades](#p5)[(1)](Del_atomo_al_cromosoma.html#motor) (P5); presión = esa anchura *por* la densidad, medida contra una superficie. Ambas son estadísticas del colectivo: un átomo no tiene ni T ni P.
- **Líquidos**: la fórmula real (virial) suma dos partes: la cinética nkT y la de las **interacciones** (Σ r·F/3V): atracciones tiran hacia dentro (negativa), el [muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones) empuja hacia fuera. En el agua las dos son enormes (~miles de bar) y casi se cancelan: el 1 bar de la atmósfera es el *resto*. Por eso el agua es casi incompresible (módulo ≈ 2,2 GPa: comprimirla un 1 % cuesta ~220 bar) y por eso la estructura del ADN apenas nota la presión: hacen falta cientos o miles de bar para mover [Tm](Del_atomo_al_cromosoma.html#g-tm)[(1)](Del_atomo_al_cromosoma.html#medio) unos grados.
**Dónde está la presión en el proyecto (más de lo que parecía)**:
1. **F5 la controla desde el principio**: `f5_simular.py` añade un `MonteCarloBarostat(1 bar, 310 K)`: cada pocos pasos propone cambiar el volumen de la caja y acepta con Boltzmann. El resultado es la columna *Density* de `trayectorias/f5_log.csv`, estable en ~1,02 g/mL: el barostato manteniendo el 1 bar. Cambiar `1*bar` por `1000*bar` es todo lo que hace falta para simular el fondo del mar.
2. **Sandbox**: la burbuja es de radio fijo, luego lo que se fija es el volumen y la presión es la que resulta. Extensiones directas: mostrar P (pared y virial) en el HUD y teclas para encoger/agrandar R_PARED y ver cómo se dispara P al comprimir un líquido, o cómo cae a cero al expandir y aparecer el vapor (P5).
**Remate**: temperatura y presión son las dos estadísticas de Maxwell-Boltzmann leídas en dos sitios: la T en el histograma de velocidades, la P en la frontera. La ecuación de los gases PV = nRT de bachillerato es exactamente esa contabilidad.
---
## P22 (02/09/2026) — La gota «parece viva»: ¿una ameba usa ese movimiento a su antojo?
**Casi: no lo manda, lo rectifica.** La gota del sandbox es agitación térmica *en equilibrio*: intensa, pero sin dirección, sin memoria y sin objetivo; a escala de nanómetros el temblor es 10¹² veces más rápido que cualquier gesto de una ameba. Una ameba no puede ordenar a sus moléculas que se muevan hacia la presa. Lo que hace es lo que hacen todas las máquinas de la vida (doc F0, sección 8): funcionar como **trinquetes brownianos[(1)](#p28)[(1)](#p28)[(1)](#p28)[(1)](#p28)**. En el pseudópodo, la membrana fluctúa térmicamente; cuando un temblor abre un hueco entre la punta de un filamento de actina y la membrana, se cuela un monómero nuevo y la membrana **ya no puede volver atrás**: el azar empuja en todas direcciones, la célula solo deja pasar los empujones que le convienen y paga con [ATP](#p28)[(1)](musculo.html) el cerrojo (modelo de Peskin, Odell y Oster, 1993). Lo mismo con la miosina que contrae la cola. Y el temblor es, además, el transporte gratuito de la célula: la difusión reparte ATP, iones y proteínas sin que nadie las lleve.
**La diferencia entre la gota y la ameba no es el movimiento (es el mismo) sino tres cosas que la gota no tiene**: información (ADN: qué máquinas construir), un flujo de energía que la mantiene fuera del equilibrio (ATP, el chorro solar de P9) y esos cerrojos que convierten azar en dirección. Sin ellas, «vivo» es solo Boltzmann; con ellas, Boltzmann es el motor.
---
## P23 (02/09/2026) — ¿Qué significa que la nube «no es esférica ni se repolariza»? ¿Qué hacen AMOEBA y Drude?
**Lo que es un átomo en nuestro modelo clásico**: una bola dura (la σ de [Lennard-Jones](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](#p19)[(2)](sandbox.html), el [muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones)) con una pegatina de carga fija en el centro (la carga parcial). Dos simplificaciones respecto al átomo real, que es una nube de electrones:
**1. La nube tiene forma, la bola no.** Los enlaces esculpen la nube: el O del agua tiene dos **pares solitarios** que apuntan en direcciones concretas (por eso acepta puentes de H por ahí y no por cualquier lado); un anillo aromático tiene su nube π encima y debajo del plano; un Cl unido a un C tiene un casquete *positivo* en la punta del enlace (el «agujero σ» de P8). Una carga puntual en el núcleo no puede representar nada de eso. Parches clásicos: el modelo de agua TIP4P añade un cuarto punto de carga fuera del O para imitar los pares solitarios; los campos modernos ponen «sitios virtuales» extra. AMOEBA[(1)](Del_atomo_al_cromosoma.html#existe)[(1)](Del_atomo_al_cromosoma.html#existe)[(1)](Del_atomo_al_cromosoma.html#existe)[(1)](Del_atomo_al_cromosoma.html#existe) lo hace en serio: cada átomo lleva carga + **dipolo + cuadrupolo** (multipolos), que describen la forma de la nube, no solo su cantidad.
**2. La nube se deforma, la pegatina no.** Acerca un Na⁺ a un agua y los electrones del agua se corren hacia el catión: su dipolo crece. Es la **polarización**, y es reversible: cuando el ion se va, la nube vuelve. Números que lo miden: el dipolo del agua aislada es 1,85 D; en el líquido, rodeada de otras aguas que la polarizan, ~2,5–3 D. Nuestro TIP3P lleva un dipolo fijo de 2,35 D: está **pre-polarizado en promedio** para el agua líquida, y funciona mientras el entorno se parezca al agua líquida. Falla cuando no: junto a un Mg²⁺ (que polariza brutalmente y en clásico se une «mal» a los fosfatos), en una interfase, dentro de una proteína, en la interacción catión–π. Analogía: la nube real es un globo que se aplasta cuando lo aprietas; la nuestra es una bola de billar con una flecha pintada.
**Cómo lo arreglan los campos polarizables** (los dos los soporta OpenMM):
- **Drude**: a cada átomo pesado se le cuelga una partícula cargada con un muelle; el campo de los vecinos la desplaza y ese desplazamiento *es* la nube deformándose. Coste ×3–5.
- **AMOEBA**: multipolos (forma) + dipolos inducibles que en cada paso se recalculan hasta autoconsistencia (cada dipolo inducido cambia el campo que induce a los demás; hay que iterar). Coste ~×10.
**Qué veríamos distinto con ellos**: [capas de hidratación](#p20)[(1)](soluciones.html?exp=sal) de iones más fieles, el dipolo del agua creciendo junto al Na⁺ *en vivo*, pares iónicos y solubilidad más correctos, Mg²⁺ en los fosfatos bien. Para el ADN en NaCl a 310 K, el fijo aguanta (el RMSD de F5 lo demuestra); para el «campo que se deforma por los vecinos» que pediste en F0b, haría falta esto.
---
## P24 (02/09/2026) — ¿Los enlaces están «hardcoded»? ¿No los rompemos nunca? ¿Y entonces qué hace el muelle angular?
**Sí y no, y la clave es distinguir dos cosas que el lenguaje mezcla**:
| | Qué es | En nuestro modelo |
|---|---|---|
| **Conectividad** (topología) | quién está enlazado con quién, con qué orden de enlace | **dato de entrada, fijo para siempre**: la lista de enlaces del SDF/PDB |
| **Conformación** (forma) | dónde está cada átomo: longitudes, ángulos, torsiones, plegado | **emerge** de los muelles + las fuerzas no enlazadas, en cada paso |
Analogía: una cadena de bicicleta. Qué eslabón engancha con cuál está fijado en fábrica (conectividad). La forma que toma la cadena la deciden las fuerzas (conformación). El muelle angular, el de enlace y (en F5) los diedros son los eslabones: fijan longitudes, ángulos y torsiones preferidos, con cierta elasticidad. Con ellos y con Coulomb + LJ, la molécula **se pliega sola**: la doble hélice de F5 se mantiene en forma B (RMSD 1,6–2,5 Å) porque los diedros bsc1/OL15 + [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice) + fosfatos + agua la sostienen, no porque esté fijada. (La excepción fue F3: las bases como cuerpos rígidos, forma *hardcoded*, P13. En F5 ya son flexibles.) Un ejemplo célebre: en 2011 se plegaron proteínas enteras desde el hilo estirado con exactamente esta física, conectividad fija y forma emergente.
**Por qué el muelle no se rompe nunca**: E = k(r−r₀)² crece sin límite. Cuanto más estiras, más fuerte tira: energía de enlace infinita. Un enlace real sigue la **curva de Morse**: sube, se aplana y a partir de la energía de disociación (~450 kJ/mol para O–H) suelta. A 300 K solo se visita el fondo del pozo (kT es 1/200 de eso, P5/P15), y ahí las dos curvas coinciden: el muelle es exacto donde importa. Es un modelo «solo lectura» de la química: las moléculas conservan su identidad.
**Por qué romper enlaces no es «poner un Morse»**: se podría dejar que el muelle suelte, pero al romperse un enlace los electrones se redistribuyen: cambian las cargas parciales, cambian los tipos de átomo, pueden formarse enlaces nuevos con otro vecino. Sin electrones, el modelo clásico no tiene con qué decidir nada de eso. Consecuencias concretas en nuestro laboratorio: el agua **no puede autoionizarse** (no hay H₃O⁺ ni OH⁻: no existe el pH), un protón no puede saltar de un agua a otra (mecanismo de Grotthuss), y en un electrodo no puede pasar un electrón del metal al agua (electrólisis). Las salidas, de menor a mayor: **ReaxFF** ([campo de fuerzas](Del_atomo_al_cromosoma.html#g-campo-de-fuerzas)[(1)](Del_atomo_al_cromosoma.html#campo)[(2)](#p19) con órdenes de enlace variables que forma y rompe enlaces con reglas), **[QM/MM](Del_atomo_al_cromosoma.html#g-qm)(1)** (una región pequeña cuántica, el resto clásico), **potenciales ML** (ANI, MACE: redes entrenadas con millones de cálculos cuánticos, forman y rompen enlaces a ~100–1.000× el coste clásico; OpenMM 8 los integra).
---
## P25 (02/09/2026) — ¿Qué es MD? (y qué quería decir «en MD la presión es uniforme»)
**MD = *Molecular Dynamics*, dinámica molecular[(1)](Del_atomo_al_cromosoma.html#g-md)[(2)](Del_atomo_al_cromosoma.html#traduccion)[(1)](Del_atomo_al_cromosoma.html#g-md)[(2)](Del_atomo_al_cromosoma.html#traduccion)[(1)](Del_atomo_al_cromosoma.html#g-md)[(2)](Del_atomo_al_cromosoma.html#traduccion)[(1)](Del_atomo_al_cromosoma.html#g-md)[(2)](Del_atomo_al_cromosoma.html#traduccion)**: la técnica de simular materia átomo a átomo aplicando la segunda ley de Newton a cada átomo con las fuerzas de un [campo de fuerzas](Del_atomo_al_cromosoma.html#g-campo-de-fuerzas)[(1)](Del_atomo_al_cromosoma.html#campo)[(2)](#p19) (P19), a pasos de femtosegundos, y grabando la trayectoria. Es **todo lo que hacemos desde F2**: `sandbox.html` es un motor de MD escrito por nosotros; OpenMM (F5) es un motor de MD profesional; [oxDNA](Del_atomo_al_cromosoma.html#g-oxdna) es MD con bolas en vez de átomos (grano grueso). El bucle es siempre el mismo (doc F0b, sección 5): fuerzas → mover → termostato → foto. Lo que NO es MD: la mecánica cuántica (resuelve electrones, no mueve núcleos con Newton), el Monte Carlo (el barostato de F5 lo usa: prueba cambios al azar y acepta con Boltzmann, sin tiempo), y la animación dirigida (Drew Berry: coreografía, no física).
**«En MD la presión es uniforme»** quería decir: en una caja de simulación de nanómetros, la gravedad es despreciable (mgh/kT ≈ 10⁻¹¹, P21), así que **no hay gradiente hidrostático**: no hay «abajo» más comprimido que «arriba», como sí lo hay en la atmósfera o en el mar. La presión vale lo mismo en cualquier punto de la caja, salvo fluctuaciones térmicas. Matiz: en una interfase (superficie de la gota, junto a una pared) la presión deja de ser igual en todas las direcciones, y esa diferencia es exactamente la tensión superficial que vimos en F2.
---
## P26 (02/09/2026) — ¿Puedo reproducir los experimentos de F7 en casa? ¿Es peligroso lo que el modelo no sabe simular?
**Sí a lo primero, y la segunda pregunta mezcla dos cosas.** Que la simulación necesite «las salidas de mayor ambición» (cuántica, potenciales reactivos) para un fenómeno es una **limitación del modelo**; el peligro de un experimento es una **propiedad de la realidad**. No van juntos: la electrólisis del agua salada no la podemos simular *y además* es la única de la lista con riesgo (cloro en los electrodos); el pH y el hielo tampoco los podemos simular y son inofensivos; la disolución y la ósmosis las simulamos perfectamente y son igual de inofensivas. **Lo que decide es la tabla de seguridad, no el modelo.**
**Lo que sí se puede hacer en la cocina, y qué número se compara con el sandbox** (detalle en `plan_F7_soluciones.md`, secciones 3, 5 y 8):
- **Sal y azúcar disolviéndose**: termómetro digital: la salmuera baja 1–3 °C (ΔH ≈ +4 kJ/mol: endotérmica pero espontánea, gana la entropía). El NaCl disuelve casi lo mismo en frío que en caliente porque su ΔH ≈ 0 (van 't Hoff); el azúcar, no. Cronómetro: fría vs caliente (Arrhenius).
- **Ósmosis**: rodajas de patata o pepino pesadas antes y después de 1–2 h en agua pura y en salmueras de 1, 2 y 3 cucharadas por litro; ganan masa en agua pura, la pierden en sal, y hay una concentración isotónica donde no cambian. Huevo sin cáscara (48 h en vinagre) en agua vs sirope. Curva masa vs concentración contra la π(c) del sandbox.
- **Conductividad**: multímetro en ohmios entre dos clavos (sin riesgo) o LED + pila de 9 V (brilla más con más sal). **Segundos, ventilado, recipiente abierto**: con NaCl los electrodos hacen un poco de cloro e hipoclorito. Bicarbonato o sal de Epsom evitan el cloro. Nunca red eléctrica.
- **Presión**: jeringa tapada: el aire cede, el agua no.
- **Vapor con sal**: la salmuera saturada hierve a ~106–108 °C (cuidado con las salpicaduras).
- **Colorante en agua fría y caliente**: difusión vs temperatura.
**Y lo que el sandbox no podrá enseñar aunque la cocina sí**: el hielo con sal (TIP3P funde a ~146 K y nuclear hielo tarda microsegundos), las burbujas de los electrodos, el pH. Límite honesto, declarado en el plan (sección 7).
---
## P27 (02/09/2026) — ¿La bomba sodio-potasio contrae el músculo? ¿Se puede simular en directo, átomo a átomo?
**No contrae: recarga.** La cadena real de «excitar → contraer → relajar» es: canales de Na⁺ que se abren (potencial de acción, ~1 ms) → el impulso abre el retículo y sale **Ca²⁺** → el Ca²⁺ destapa la actina → las cabezas de **miosina** tiran gastando [ATP](#p28)[(1)](musculo.html) (contracción) → **SERCA** bombea el Ca²⁺ de vuelta (relajación). La bomba Na⁺/K⁺ es el paso de mantenimiento: saca 3 Na⁺ y mete 2 K⁺ por cada ATP, ~100 ciclos por segundo, siempre, para reponer lo que cada impulso gasta; sin ella el músculo deja de ser excitable en minutos (ouabaína, digoxina). SERCA y la Na⁺/K⁺ son de la misma familia (ATPasas de tipo P): mismo aspartato fosforilable, mismo ciclo de [acceso alternante](#p32)[(1)](directo.html)[(2)](bomba.html) E1/E2.
**¿En directo y átomo a átomo? La bomba entera, no; el principio, sí.** La escalera de escalas en la GTX 1650 (F5 mide ~100 ns/día con 18.000 átomos):
| Sistema | Átomos | Lo que dura el fenómeno | Tiempo de reloj |
|---|---|---|---|
| electrolito bajo campo | 1.000–5.000 | ps–ns | segundos–minutos, en directo |
| poro de juguete con voltaje | 5.000–10.000 | un ion cada 1–10 ns | minutos por ion |
| canal real en membrana | 25.000–100.000 | un ion cada 10–100 ns | horas–días por ion (película) |
| **bomba Na⁺/K⁺ en membrana** | ~250.000 | ciclo de **10 ms**; transiciones µs–ms; ATP = química | 10⁷ ns a ~10 ns/día ≈ **miles de años**; y la fosforilación no es clásica (P24) |
| sarcómero | > 10¹² | ms | fuera de toda escala atómica |
La investigación simula la bomba durante cientos de nanosegundos alrededor de sus estructuras cristalinas (E1 con 3 Na⁺, PDB 3WGU; E2 con 2 K⁺, PDB 2ZXE) para ver cómo entran los iones, y la transición entre estados hay que **empujarla** con dinámica dirigida porque su duración excede lo que la MD sin sesgo alcanza (Castillo, Roux et al., 2015). **Nadie ve un ciclo completo emerger de la física atómica.**
**Lo que sí hará F9** (plan en `plan_F9_directo.md`): la tubería «OpenMM en directo» estrenada con un electrolito excitable (campo on/off), un poro con voltaje (corriente, curva I–V, selectividad: un canal nunca sube el gradiente), una **bomba de juguete** con acceso alternante coreografiado sobre iones y agua reales (dos compuertas que nunca están abiertas a la vez + un sitio que cambia de preferencia; el gradiente y el potencial de membrana emergen; el ciclo se declara como coreografía, igual que el ATP lo coreografía en la realidad), la visita a la bomba real en el visor, y el músculo excitable a su escala (modelo eléctrico + Ca²⁺ + puentes cruzados, en JS, no atómico).
---
## P28 (02/09/2026) — ¿Por qué importan los nanotubos? ¿Qué es el ATP y qué papel juega (en la bomba, en el músculo y en la ameba)?
### Los nanotubos: el canal más simple que existe
Un **nanotubo de carbono** es una lámina de grafeno (carbonos en panal) enrollada en cilindro: diámetro de 0,7 a 2 nm, longitud de micras, rígido, liso a escala atómica e hidrófobo por dentro (Iijima, 1991). Para nosotros importa porque es **el canal ideal de juguete**: geometría exacta, sin partes móviles, parámetros de un solo tipo de átomo. Y porque la física que ocurre dentro es la misma que en los canales biológicos:
- **Agua en fila india**: en 2001 (Hummer, Rasaiah y Noworyta, *Nature*) la MD mostró que un nanotubo de 0,8 nm se llena solo de agua y la conduce como una cadena de una molécula de ancho, a pulsos: exactamente como la gramicidina y las aquaporinas. La red de puentes de H de F2, forzada a una dimensión.
- **Selectividad por deshidratación**: un ion que entra en un tubo estrecho tiene que soltar parte de sus ~6 aguas (F7b); eso cuesta decenas de kJ/mol y crea una barrera que depende del tamaño y la carga. Es el principio del filtro de selectividad del canal de K⁺ (sus carbonilos imitan las aguas que el K⁺ pierde) y del «borde cargado» del poro de F9b.
- **Aplicaciones reales**: nanotubos cortos insertados en bicapas lipídicas funcionan como canales artificiales (más permeables al agua que la aquaporina); membranas de nanotubos alineados para desalar (agua rápida, iones rechazados: la ósmosis inversa de F7c hecha material); nanoporos para secuenciar ADN (los secuenciadores de Oxford Nanopore leen la hebra por la corriente iónica que bloquea al pasar por un poro proteico: nuestro F9b es la caricatura de ese detector).
Y la biología tiene sus propios «nanotubos»: los **microtúbulos** (tubulina, 25 nm), los raíles por donde kinesina y dineína cargan vesículas y el huso que reparte los cromosomas en la mitosis (gastando ATP[(1)](musculo.html)[(1)](musculo.html)[(1)](musculo.html)[(1)](musculo.html), claro). El ejemplo `nanotube.xml` de NanoVer es un nanotubo de carbono por el que se empuja un metano en realidad virtual: la misma idea que F9b.
### El ATP: la moneda, no la pila
**Adenosina trifosfato**: adenina + ribosa + tres fosfatos en fila (α, β, γ). Fíjate: es **el mismo nucleótido del ADN** (el dATP del doc F0, sección 3.2, con ribosa en vez de desoxirribosa). Con cuatro cargas negativas a pH 7 casi siempre va agarrado a un Mg²⁺ (MgATP²⁻, el Mg²⁺ de P23).
**De dónde sale su energía**: romper el fosfato γ (ATP → ADP + Pᵢ) libera ~30 kJ/mol en condiciones estándar y **~50–60 kJ/mol en la célula**, unos **20 kT**. Tres razones: los fosfatos cargados se repelen y separarlos alivia (Coulomb, P19); los productos se estabilizan mejor por resonancia e hidratación; y, sobre todo, **la célula mantiene la proporción ATP/ADP lejísimos del equilibrio** (unos diez órdenes de magnitud). El ATP no es una pila con energía «dentro»: es un desequilibrio químico mantenido a propósito, como la burbuja mantenida a 300 K por el termostato. Quien lo recarga es la **ATP sintasa**, un motor rotatorio en la mitocondria movido por un gradiente de protones (que sale de quemar glucosa con oxígeno): un ser humano recicla del orden de su propio peso en ATP cada día.
**Por qué 20 kT es el número mágico**: el factor de Boltzmann de P5 dice que un paso que cuesta 20 kT se deshace solo con probabilidad e^(−20) ≈ 10⁻⁹. Gastar un ATP en un paso lo vuelve **irreversible en la práctica**: es el **cerrojo** del trinquete de P22. Toda máquina de la vida hace lo mismo: deja que el azar térmico proponga (gratis, a 650 m/s) y usa un ATP para que no pueda volver atrás.
**Cómo lo usa cada máquina que hemos nombrado**:
- **Bomba Na⁺/K⁺ y SERCA** (ATPasas de tipo P): el ATP no «empuja»; **pega su fosfato γ a un aspartato** de la bomba (E1 → E1P). Ese añadido covalente cambia la forma de la proteína (abre una puerta y cierra la otra: [acceso alternante](#p32)[(1)](directo.html)[(2)](bomba.html)), y quitar el fosfato después la devuelve (E2P → E2). Es una llave que se pone y se quita, y por eso el ciclo es química y no lo puede hacer nuestra MD clásica (P24).
- **Miosina**: al unirse un ATP la cabeza **suelta** la actina; al hidrolizarlo se «arma» como un resorte; al agarrar la actina y soltar el fosfato da el golpe de fuerza (~10 nm); suelta ADP; espera otro ATP. Sin ATP las cabezas se quedan agarradas: el *rigor mortis* es literalmente eso.
- **La ameba** (corrigiendo la frase «el ATP previene la deformación»): el ATP **no impide** que la membrana se deforme, la membrana se deforma sola, térmicamente, sin parar. Lo que hace el ATP es que una deformación **no pueda deshacerse**: cada monómero de actina lleva un ATP; la forma con ATP es la que se añade al extremo del filamento cuando un temblor abre un hueco, y una vez añadida la membrana no vuelve; más tarde la actina hidroliza su ATP, se marca como «vieja» y se desmonta por el otro extremo para reciclarse (*treadmilling*). El coste energético se paga sobre todo en mantener muchos monómeros cargados con ATP listos para pegarse (unos 5 kT de empuje por monómero, ~1 pN por filamento). La **dirección** la pone otra química: la ameba activa los nucleadores de actina (Arp2/3, forminas) solo en el lado donde huele comida, y la miosina aprieta la cola. Azar térmico para moverse, ATP para no retroceder, señalización para elegir dónde.
**Moraleja**: temperatura (P5), presión (P21) y ATP son la misma historia contada tres veces: un desequilibrio mantenido a propósito, medido en kT. La vida es el arte de gastar 20 kT en el momento justo.
---
## P29 (03/09/2026) — ¿El ADN está en el núcleo así, en esta agua con sal y plegado, tal como lo modelamos en F5?
**Sí en lo que importa a esta escala; no en casi todo lo demás.** Lo simulado es, con precisión, **ADN en un tubo de ensayo con suero fisiológico** — el sistema de referencia de media bioquímica (PCR, secuenciación, cristalografía) y la razón de que existan números con los que comparar.
**Lo que es real y por eso el modelo vale:**
- **La doble hélice B** local: escalera, 32 puentes, [apilamiento](Del_atomo_al_cromosoma.html#g-stacking)[(1)](Del_atomo_al_cromosoma.html#helice), surcos — a escala de nm el ADN nuclear *es* esto.
- **Agua líquida a 310 K, pH ~7,3, iones apantallando**: hidratación, espina del surco menor, atmósfera iónica en los fosfatos — física idéntica a la nuclear. La fuerza iónica ~150 mM y su λ_D ≈ 8 Å (F6) son las correctas en magnitud.
**Lo que es aproximación de laboratorio, no de núcleo:**
1. **Los iones son los del tubo, no los de la célula.** NaCl 150 mM es el «suero fisiológico» estándar; dentro de la célula manda el **K⁺ (~140 mM)**, el Na⁺ es minoritario (~10 mM) y el Cl⁻ escaso (los aniones son fosfatos, proteínas, glutamato). Faltan el **Mg²⁺** (0,5–1 mM libre: divalente, puentea dos fosfatos, el ion que más quiere al ADN) y las **poliaminas** (espermina⁴⁺, espermidina³⁺: cationes largos que se tumban en los surcos y compactan).
2. **Nuestro ADN va desnudo; el nuclear va vestido.** El 75–90 % está enrollado en **[nucleosoma](Del_atomo_al_cromosoma.html#g-histona)[(1)](Del_atomo_al_cromosoma.html#cromatina)s** (147 pb en torno a 8 histonas, radio 4 nm, argininas en el surco menor cada 10 pb). Nuestros 12 pb equivalen a un trozo de *linker* suelto — o al oligonucleótido sintético que el [1BNA](visor.html)[(1)](Del_atomo_al_cromosoma.html#helice)[(2)](trayectoria.html) fue en realidad.
3. **Caja diluida vs multitud.** El núcleo tiene 100–300 g/L de macromoléculas: más gel que disolución, con condensados. Nuestra [caja periódica](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared) es una disolución diluida de ADN; en cromatina densa hay ~100 mg/mL.
4. **Plegado: solo el primer nivel.** La hélice sí; nucleosoma → collar → bucles → TADs → cromosoma, no. El ladrillo está perfecto; el edificio, ausente.
5. **Los extremos que se deshilachan son artefacto del corte.** Un fragmento tiene dos bordes sin vecino que apile; en un cromosoma de 50–250 millones de pb no hay extremos libres (salvo telómeros). El **núcleo de la hélice** (8 pares centrales, RMSD ~1–1,5 Å) es lo representativo; los picos de RMSD a 3 Å son el peaje del fragmento.
**Escalera para acercarlo al núcleo**: K⁺/Mg²⁺ en vez de NaCl (un parámetro) → poliaminas (fácil) → nucleosoma 1AOI con histonas (~300.000 átomos con agua; factible en la GTX a días por decena de ns) → hacinamiento (ya grano grueso).
---
## P30 (02/09/2026) — F7a–F7e entregadas: lo que se midió en `soluciones.html` y lo que enseñó construirlo
Un solo visor nuevo, `soluciones.html?exp=presion|sal|osmosis|campo|debye`, con el motor de F2 generalizado a iones (parámetros de Na⁺, Cl⁻, K⁺, Mg²⁺, Ca²⁺ copiados del `amber14/tip3p.xml` de OpenMM: los mismos de F5), un soluto neutro y una carga fija; `campo7.js` es la capa de campos de F4/F6 con el término −E·x. Rendimiento medido: 100 pasos de 462 átomos en 0,18 s → **0,2–0,4 ps de realidad por segundo de pantalla** (sube con `.`). Los números de la primera tarde:
**F7a, presión** (100 aguas; presión = fuerza media de la pared / área, media móvil de 5 ps):
| Burbuja | R (Å) | Req gota (Å) | «densidad» (g/cm³) | P (bar) |
|---|---|---|---|---|
| holgada (de partida) | 11,1 | 9,3 | 0,88 | 7 |
| pegada a la gota (tecla 3) | 10,1 | 9,2 | 0,93 | 61 |
| −6 % | 9,5 | 9,1 | 0,95 | 155 |
| −12 % | 9,0 | 8,9 | 1,03 | 393 |
| −18 % | 8,4 | 8,6 | 1,14 | 481 |
| −24 % | 7,9 | 8,5 | 1,16 | 907 |
La lección se ve: de ~0 a ~900 bar apretando, con la gota cediendo poco. Honestidad: la pared es blanda (armónica, K = 10 kJ/mol/Ų: los átomos la penetran) y Req es una estimación de esfera uniforme, así que el módulo efectivo sale ~0,4 GPa frente a los 2,2 GPa del agua real: el comportamiento es el correcto, el número fino no. Con la burbuja holgada la gota vive a P ≈ 0: un líquido a su propia presión de vapor.
**F7b, la sal** (cristal 2×2×3 de NaCl, 12 iones, 150 aguas): a 300 K el cristal aguanta decenas de ps (a los 8 ps, 0/12 «disueltos»; las esquinas ya llevan 2–3 aguas). A **400 K**, 3/12 a 2,4 ps, 5/12 a 7 ps, **6/12 a 9,6 ps** y 8/12 a 10 ps en la captura: las esquinas y aristas primero, como en Lanaro y Patey. Hidratación de los iones libres ≈ 5–6 aguas (Na⁺, radio 3,2 Å) y 6–8 (Cl⁻, 3,9 Å). El «calor del baño» no resuelve los +4 kJ/mol por par: en una gota de 150 aguas el termostato mueve cientos de kJ/mol por ps solo con la relajación del agua; queda como indicador, no como medida.
**F7c, ósmosis** (160 aguas, 4 pares NaCl en el lado x<0, pared iónica en x=0): en los primeros 10 ps **π medida = 0** porque ningún ion había tocado aún la pared (difunden a ~0,1 Ų/ps): la media de Luo–Roux necesita cientos de ps. La demostración rápida es el **pistón**: libre y sin presión externa avanzó de 0 a **+1,7 Å hacia el agua pura en ~25 ps** (el lado salado crece: ósmosis); con **300 bar** externos retrocedió a 0,3 Å, el lado salado se concentró (2,05 → 2,35 M) y π medida subió a **~90 bar** frente a ~60 de van 't Hoff (φ ≈ 1,5, ruidoso: media de 10 ps con 8 partículas). Comportamiento correcto; número aún grueso.
**F7d, campo eléctrico** (150 aguas, 4 pares, E = 1 V/nm): la sorpresa útil. Aplicando E a **todas** las cargas, la deriva era ~0: la gota se polariza como un dieléctrico en vacío y el campo interior cae a 3E/(ε+2) ≈ E/27 (electrostática de bachillerato, comprobada sin querer). Por eso el campo actúa por defecto **solo sobre los iones** (E es entonces el campo dentro del medio, como en una [caja periódica](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared)) y la tecla 3 enseña el apantallamiento. Con eso, en modo circuito: **Na⁺ ≈ +1,0 Å/ps, Cl⁻ ≈ −0,7 Å/ps** ⇒ movilidades ~1·10⁻⁷ y ~7·10⁻⁸ m²/V·s (reales a dilución infinita: 5,2·10⁻⁸ y 7,9·10⁻⁸): el orden de magnitud correcto en una gota de 4 pares y sin termostato sobre los iones. Tres cargas elementales completaron el circuito en 8 ps.
**F7e, Debye** (carga fija −2 e, 3 pares): un Na⁺ se pega en contacto a ~2,2 Å (par iónico: con −2 e y un sitio pequeño, condensación, el «Mg²⁺» de P23 al revés); λD teórica 4,4 Å para la fuerza iónica de la gota. El perfil radial acumulado necesita cientos de ps (10–30 min de pantalla) para dibujar la exponencial: se deja como experimento de sobremesa.
**Lo que enseñó construirlo**: (1) la pared iónica y el pistón son fáciles; lo caro es *promediar*; (2) una gota en vacío no es un electrolito entre electrodos: el dieléctrico manda; (3) el «electrodo» tiene que estar en la superficie del líquido, porque un ion hidratado no abandona el agua; (4) 500 átomos es el techo cómodo del JS: para más, la vía de rigor o WebGPU.
---
## P31 (03/09/2026) — F9a entregada: la física de F5 corriendo mientras miras, y por qué el mismo campo da respuestas opuestas en una gota y en una caja periódica
Hasta ahora el laboratorio tenía dos vías: el motor JS, que corre en el navegador y se entiende entero pero simplifica la física, y OpenMM, que hace física de investigación pero en diferido (F5: dos días de GPU y luego la película). **F9a es la tercera**: OpenMM corriendo en la GPU *mientras miras*. Un proceso Python mantiene la simulación viva y sirve el último fotograma por HTTP; `directo.html` lo pide 20 veces por segundo y le devuelve las teclas. Es lo que hace NanoVer con sus clientes de realidad virtual, aquí con una tubería de tres rutas (`/info`, `/frame`, `/orden`) y sin instalar nada nuevo.
**El experimento**: NaCl 1 M en agua, caja cúbica de 3 nm periódica con [PME](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared), AMBER + TIP3P, 2 fs con SHAKE — **la misma física de F5**, 2.597 átomos. Un campo eléctrico conmutable desde el teclado. Pulsar `→` **excita** (los iones derivan y nace una corriente); pulsar `0` **relaja** (el temblor térmico borra la deriva en picosegundos). La gráfica lo enseña como una onda cuadrada naranja (el campo) con la conductividad verde siguiéndola.
**Velocidad medida** en la GTX 1650 (con la Intel integrada descartada a mano, P18 otra vez): 6,5 ps/s a solas y 4–4,5 ps/s con el navegador pidiendo fotogramas. Frente a los 0,2–0,4 ps/s del motor JS con 500 átomos (P30), es **entre 15 y 25 veces más rápido con cinco veces más átomos**. Leer el estado de la GPU cuesta 0,2 ms: la tubería no es el cuello de botella.
**Lo medido** (concentración real 0,98 M):
| Campo | Promediado | Deriva Na⁺ | Deriva Cl⁻ | Conductividad | Agua alineada |
|---|---|---|---|---|---|
| 1,0 V/nm | 112 ps | +0,041 nm/ps | −0,054 nm/ps | 9,0 S/m | 0,71 |
| 0,5 V/nm | 300 ps | +0,022 | −0,029 | 9,6 S/m | 0,58 |
| **0,1 V/nm** (lineal) | **9,9 ns** | **+0,0045** | **−0,0062** | **10,1 S/m** | 0,22 |
| 0 | — | 0,002 (ruido) | 0,002 (ruido) | — | 0,00 |
Contra el laboratorio: la conductividad del NaCl 1 M a 25 °C es **8,5 S/m**; medimos 10,1, un 19 % por encima. Las movilidades salen 4,5·10⁻⁸ (Na⁺) y 6,2·10⁻⁸ (Cl⁻) frente a 5,2 y 7,9·10⁻⁸ **a dilución infinita** — y deben salir más bajas, porque a 1 M los iones se estorban entre sí. La cuenta ideal (Nernst–Einstein con las movilidades de dilución infinita) daría 12,6 S/m: el experimento baja a 8,5 y nosotros a 10,1, o sea que **capturamos parte del estorbo entre iones, pero menos del que hay**. A eso se suma que el agua TIP3P es demasiado fluida (su viscosidad es un tercio de la real), lo que empuja las movilidades hacia arriba. Un 20 % es exactamente el margen que la literatura da a los modelos sin polarización.
**Y la lección que no esperaba**: en P30 (F7d) descubrimos que aplicar el campo a *todas* las cargas mataba la deriva, porque la gota se polariza como un dieléctrico y el campo interior cae 27 veces; allí había que empujar solo a los iones. **Aquí, en la caja periódica, hay que hacer justo lo contrario**: el campo se aplica a todas las cargas, y es correcto. La razón es la frontera. Una gota en vacío tiene superficie, y en esa superficie se acumula carga de polarización que crea un campo despolarizante hacia dentro. Una caja periódica no tiene superficie: se repite hasta el infinito, así que no hay dónde acumular esa carga y el campo aplicado *es* el campo dentro del medio. Misma pregunta física, respuestas opuestas, y quien decide es la frontera — la misma «pared mock» de la Fase 0b, que resulta no ser un detalle técnico sino parte de la física del experimento.
**Dos detalles de oficio**, por si vuelven a aparecer: (1º) la corriente hay que medirla **solo sobre los iones**. El campo sí empuja a las cargas parciales del agua (eso es física), pero el agua no *transporta* carga: sus dipolos se orientan y se quedan quietos. Su contribución media es cero y su ruido es unas cien veces la señal de los iones, porque cada hidrógeno va a 3 nm/ps. Con el agua dentro, la conductividad saltaba entre 15 y 38 S/m; sin ella, se clava en 9–10. (2) A 0,5 V/nm el agua ya está alineada al 58 %: por encima de 0,1 V/nm esto **deja de ser respuesta lineal**, y por eso la medida buena es la de 0,1 con nanosegundos de promedio, aunque la vistosa sea la de 1 V/nm.
---
## P32 (05/09/2026) — F9b–F9e entregadas: un canal, una bomba de juguete, la bomba de verdad y un músculo
Cuatro piezas que, juntas, contestan la pregunta con la que empezó todo esto ([P27](#p27)): ¿se puede ver una bomba sodio-potasio funcionando? La respuesta honesta resultó ser «la de verdad no, pero su principio sí, y con física de verdad debajo».
### F9b · un canal: `directo.html` con `--exp poro`
Una membrana de partículas **fijas** (masa cero en OpenMM: no se mueven nunca, así que son puro [muro de Pauli](Del_atomo_al_cromosoma.html#g-lennard-jones)[(1)](#p19)) con un agujero de 1,4 nm de diámetro, agua y sal a los dos lados, y un voltaje. Medido a 1,8 V: unos **35 cruces en 1,2 ns**, con 2 iones y 35 aguas dentro del poro en todo momento. Los cationes van en un sentido y los aniones en el otro; los dos llevan corriente en la misma dirección.
Y lo que **no** pasa, que es la lección: por mucho que esperes, los dos lados se igualan y ahí se quedan. **Un canal nunca sube un gradiente.** Para subirlo hace falta gastar energía, y eso es una bomba.
### F9c · una bomba de juguete: `--exp bomba`
Dos membranas separan dos compartimentos de verdad, y el poro tiene **dos compuertas que nunca se abren a la vez** y un **sitio de unión** que cambia de preferencia entre Na⁺ y K⁺ *y se desplaza hacia el lado que está abierto*. Eso es el **acceso alternante[(1)](directo.html)[(2)](bomba.html)[(1)](directo.html)[(2)](bomba.html)[(1)](directo.html)[(2)](bomba.html)[(1)](directo.html)[(2)](bomba.html)**, y basta: medido, unos **2 Na⁺ hacia fuera y 2 K⁺ hacia dentro cada 3 ciclos**, a unos 45 segundos de reloj por ciclo. El gradiente **sube**.
Honestidad, dicha también en el propio visor: el ciclo está **coreografiado**. En la célula lo dirige el [ATP](#p28)[(1)](musculo.html) fosforilando un aspartato, química que un [campo de fuerzas](Del_atomo_al_cromosoma.html#g-campo-de-fuerzas)[(1)](Del_atomo_al_cromosoma.html#campo)[(2)](#p19) clásico no puede hacer ([P24](#p24), P28). Los iones, el agua, el gradiente y el potencial de membrana sí emergen solos. Y la de verdad mueve 3 Na⁺ por cada 2 K⁺ — por eso es electrogénica y genera corriente; la nuestra va uno por uno.
### F9d · la bomba de verdad: `bomba.html`
No es una simulación: son **coordenadas medidas** por rayos X, la misma clase de dato que el [1BNA de F1](visor.html). Dos fotos del mismo motor: **3WGU** (E1, con 3 Na⁺ ocluidos mirando al citoplasma) y **2ZXE** (E2, con 2 K⁺ mirando al exterior). El visor no las «sabe» de antemano: busca en el propio dato qué oxígenos y nitrógenos agarran a cada ion a menos de 3,2 Å, y sale el sitio de unión solo. Lo que encuentra:
| | Coordinación | Residuos que forman la jaula |
|---|---|---|
| Na⁺ (3WGU) | 5–6 vecinos a **2,29–2,49 Å** | Tyr771, Thr774, Ser775, Asp808, Gln923, Asn776, Asp804, Glu779 |
| K⁺ (2ZXE) | 4–6 vecinos a **2,69–2,76 Å** | Asp811, Ser782, Thr779, Asn783, Glu786 |
Son los residuos que la literatura da para esos sitios, y el dato bonito es la diferencia de distancias: **el K⁺ es más grande y se coordina más lejos**. La selectividad de la bomba está ahí, en unas décimas de ångström. Un ion nunca viaja solo: viaja en una jaula de oxígenos que le sustituyen el agua que traía ([P20](#p20) sobre la capa de hidratación).
### F9e · el músculo, a su escala: `musculo.html`
Aquí no hay átomos, y se dice: un sarcómero tiene más de 10¹² y su contracción dura milisegundos. Este es el nivel al que un músculo **sí** se puede simular y comparar con uno real: Hodgkin-Huxley para los canales, calcio con su bomba SERCA, y puentes cruzados de miosina. La cadena completa, con los números saliendo donde deben: reposo **−85 mV**, calcio de 0,05 a **9 µM** y vuelta, una sacudida suelta al **51 %** de fuerza, y a 40 Hz los impulsos **se funden en un tétanos**.
Y el experimento que cierra el círculo con F9c: pulsa `O` (ouabaína, que es lo que hace la digital del corazón) y la bomba se para. El sodio interno sube de 12 a 40 mM, el potencial de Nernst del sodio cae de 67 a 34 mV, y el músculo deja de responder. **La bomba no contrae nada: recarga la batería que cada impulso gasta** (P27).
### Las tres trampas que casi arruinan el experimento
Valen para cualquiera que monte algo así, y las tres se descubrieron mirando números raros:
1. **Una sola membrana en una [caja periódica](Del_atomo_al_cromosoma.html#g-pbc)[(1)](Del_atomo_al_cromosoma.html#pared) no separa nada.** Los dos «lados» se tocan al dar la vuelta por el borde (PBC). Para el canal eso está *bien* — el camino de vuelta cierra el circuito y el poro queda en serie, así que hay corriente estacionaria. Para la bomba es fatal: si los iones pueden volver por detrás, nunca se acumula gradiente. Por eso la bomba lleva una segunda membrana, sin agujero, pegada a la frontera.
2. **Un poro estrecho e hidrófobo se seca solo.** Salía con cero aguas dentro, y sin agua no pasa ningún ion: entrar significaría desnudarse entera la capa de hidratación. No es un fallo del programa, es un fenómeno real con nombre (*gating* hidrofóbico). Se arregló como lo arregla un canal de verdad: tapizando el borde de grupos polares. Pero con cuidado — a ±0,25 e el poro se moja igual de bien y **cada sitio negativo se convierte en un pozo de más de 100 kJ/mol** que atrapa al catión dentro y la bomba dejaba de soltarlo. Mojar sí, atrapar no: ±0,12 e.
3. **Contar cruces por cambio de signo no vale.** Un ion parado *dentro* del poro tiembla alrededor del plano de la membrana y sumaba un cruce en cada fotograma: salieron **342 cruces en 574 ps**, físicamente imposible. Hay que contar con histéresis: el ion tiene que salir del todo por un lado y aparecer del todo por el otro.
Hay una cuarta, de oficio: el visor se quedaba congelado sin avisar porque una excepción dentro del bucle de dibujo impedía pedir el siguiente fotograma. Ahora el bucle captura sus propios fallos y sigue.
---
## P33 (05/09/2026) — F10: romper un enlace de verdad, y qué *es* un enlace covalente
Desde F2 hasta F9, en todo el laboratorio, un enlace covalente ha sido **un muelle que no se rompe nunca** ([P24](#p24)): su energía crece sin límite, así que cuanto más estiras, más tira. F10 tira de uno hasta partirlo, con tres modelos dibujados sobre la misma curva de energía, en `enlace.html`.
**1. El muelle.** Tiras y no pasa nada: el enlace apenas cede y la fuerza sube y sube. Medido: a 14 ps de tirón la distancia había pasado de 1,32 a 1,62 Å y la fuerza ya iba por 3,9 nN, subiendo. No se romperá jamás, por definición.
**2. Morse.** E = Dₑ(1−e^(−a(r−r₀)))². Sí cede, y a la fuerza correcta. Lo que **no** sabe es por qué: al romperse deja a cada trozo con la carga parcial que llevaba pintada, y no tiene forma de saber que los fragmentos deberían salir neutros.
**3. El cuántico, que es de verdad.** Aquí está lo bueno. Para la molécula más simple que existe —el **H₂⁺**: dos protones y **un solo electrón**— la [ecuación de Schrödinger](Del_atomo_al_cromosoma.html#electrones)[(1)](#p23) se resuelve. Se construye el orbital como suma de los dos 1s (aproximación LCAO), y las tres integrales que hacen falta tienen forma cerrada:
- solapamiento: S = e^(−R)(1 + R + R²/3)
- atracción del otro núcleo: J = −1/R + (1 + 1/R)e^(−2R)
- y la de intercambio, que es la que crea el enlace: K = −e^(−R)(1 + R)
La energía sale E = (J + K)/(1 + S) + 1/R, en unidades atómicas, **sin un solo parámetro ajustado**. Son quince líneas de JavaScript, y comprobado contra el libro: da r₀ = **1,319 Å** y Dₑ = **170 kJ/mol**, frente a los 2,5 a₀ = 1,323 Å y 1,76 eV = 170 kJ/mol que es el resultado LCAO clásico.
Y entonces se ve **qué es un enlace covalente**, que es lo que llevo explicando con palabras desde P8: la nube del visor muestra la densidad electrónica **acumulada entre los dos núcleos**. Eso es todo. El electrón está a la vez cerca de los dos protones, y esa densidad de en medio los pega. Pulsa `1` y pasas al orbital **antienlazante**: la densidad en el punto medio cae a **cero exacto** —hay un nodo— y el «enlace» se va a 724 kJ/mol por encima: no hay pozo, no hay molécula. Enlace y antienlace, con la misma fórmula y un signo cambiado.
**El detalle que más me gusta**: el muelle y el Morse no se sacan de tablas, se **ajustan a la curva cuántica** (su mínimo, su profundidad y su curvatura). Por eso los tres coinciden en el fondo del pozo por construcción, y se ve exactamente dónde se separan. Y el fondo del pozo es donde vive la agitación térmica ([P5](#p5)): por eso el muelle de F2 y de F5 es *correcto donde importa* y solo miente cuando estiras de verdad.
**Las fuerzas de ruptura**, midiendo con la pinza contra la teórica a·Dₑ/2:
| Enlace | Medida | Teórica |
|---|---|---|
| H₂⁺ | 1,81 nN | 1,86 nN |
| H₂ | 6,86 | 7,37 |
| O–H (el del agua) | 8,71 | 9,35 |
| C–C | 6,06 | 6,15 |
Siempre un poco por debajo, y está bien que así sea: el temblor térmico ayuda a saltar la barrera antes de llegar al máximo. Es el efecto de Bell y Evans, y es la misma razón por la que el microscopio de fuerza atómica midió 2,0 nN en 1999 tirando 10²⁰ veces más despacio que nosotros.
### Sobre la cuántica «de verdad», y hacia dónde sigue
Fui a montar el escalón con un potencial cuántico externo y resultó que **en esta máquina no hay ninguno instalable con pip**: `tblite`, `xtb` y `pyscf` no publican ruedas para Windows en ninguna versión de Python. Solo `torchani` (que arrastra PyTorch, 2,5 GB). Así que el escalón cuántico se hizo calculándolo nosotros — y encaja mucho mejor con el proyecto, que va justo de eso.
Y esto abre la puerta a lo siguiente, porque el mapa es más amable de lo que parece:
| Qué | Se resuelve | Aquí |
|---|---|---|
| un electrón (H, He⁺, H₂⁺) | **exacto** | ya está hecho; los orbitales del hidrógeno tienen fórmula cerrada |
| un electrón en cualquier potencial | **exacto numéricamente** | integrando la ecuación radial (método de Numerov) y buscando autovalores por tiro: ~150 líneas |
| varios electrones (He, Li, C) | aproximado, campo autoconsistente | asequible sobre la misma malla; ahí aparece el apantallamiento y **con él la tabla periódica** |
| precisión de investigación | HF/DFT/post-HF | PySCF, Psi4, ORCA: hace falta WSL o conda |
O sea que lo más instructivo se puede hacer entero en el navegador, sin instalar nada, igual que F2 y F10. Y el último escalón es el que explicaría **por qué el sodio cede un electrón y el cloro lo coge**, que es exactamente donde empezó [F7b](soluciones.html?exp=sal).
---
## P34 (05/09/2026) — Hacia F11: qué se puede ver emerger del átomo y qué hay que declarar
La pregunta fue: quiero entender el átomo —la radiación, los electrones y su comportamiento, el núcleo, Pauli y su porqué, los espines, un protón solo— con lo estudiado en el instituto mostrado en simulaciones con librerías de cuántica. El plan completo está en `plan_F11_atomo.md`; aquí va lo conceptual, que es lo que cambia la manera de mirar.
**La regla del laboratorio, aplicada al átomo.** En F2 declaramos cinco leyes y miramos emerger los [puentes de hidrógeno](Del_atomo_al_cromosoma.html#fuerzas)[(1)](sandbox.html)[(2)](#p8). En el átomo se declara aún menos: la [ecuación de Schrödinger](Del_atomo_al_cromosoma.html#electrones)[(1)](#p23) (una contabilidad de energía para una onda) y el Coulomb del núcleo. Y emerge casi todo lo que el instituto enseña como lista: que solo ciertas energías tienen solución (**la cuantización sale de una condición de contorno**: si eliges una energía cualquiera la onda se dispara a infinito, y solo en los niveles se aplana), los orbitales con sus nodos, el radio de Bohr, los −13,6 eV, el espectro con las rayas de Balmer en su color, las vidas medias de los estados excitados y las reglas de selección. Hoy lo comprobé: un solver radial de cien líneas encuentra −13,604 / −3,401 / −1,512 / −0,850 eV en 0,1 segundos, sin que nadie le diga dónde buscar.
**Lo que no emerge y hay que declarar**, dicho con la misma honestidad que el ciclo de la bomba en [P32](#p32): **el espín**, que no sale de Schrödinger sino de Dirac (relatividad), y **el signo menos de los fermiones**, que es el teorema espín-estadística de Pauli de 1940 y exige teoría cuántica de campos. Son las dos únicas cosas. Todo lo que viene después de ellas —la exclusión, las capas, la tabla periódica y **el [muro de Pauli](#p19)[(1)](Del_atomo_al_cromosoma.html#g-lennard-jones) que el laboratorio lleva pintado a mano desde F2**— se calcula.
**Por qué Pauli, en dos líneas.** Dos electrones son indistinguibles de verdad, así que intercambiarlos no puede cambiar nada observable, y eso solo deja dos opciones: la función conjunta es simétrica o antisimétrica. Si es antisimétrica, al poner los dos en el mismo sitio se anula: **ψ(x,x) = 0**. La exclusión no es un postulado, es una consecuencia. Y el muro de F2 es lo mismo visto desde fuera: dos nubes que intentan solaparse obligan a la onda a arrugarse, y arrugarse cuesta energía cinética. En F11e se calculará con dos átomos de helio y se dibujará encima del r⁻¹² de siempre.
**La radiación son dos cosas.** La luz de los átomos: una superposición de dos estados hace oscilar la nube a la frecuencia (E₂−E₁)/h, y una carga que oscila radia —es el electromagnetismo del instituto en una antena de un ángstrom—. Y la radiactividad del núcleo: la α es **efecto túnel** (Gamow, 1928), y por eso un factor 2 en energía separa los 0,3 µs del polonio-212 de los 4.470 millones de años del uranio-238; la β es la fuerza débil convirtiendo un neutrón en protón; la γ es la luz del núcleo, en MeV.
**El núcleo**: si el átomo fuera un campo de fútbol, el protón sería una lenteja en el centro. Lo pega la fuerza fuerte, que solo actúa a menos de 2 fm; la energía de enlace es masa que falta, y su curva tiene máximo en el hierro: a un lado se gana fusionando (el Sol), al otro partiendo.
**Un protón solo** tiene tamaño (0,84 fm), espín, un imán con tres quarks dentro, y en agua **nunca está solo**: salta de molécula en molécula por los puentes de hidrógeno de F2, que es justo lo que F2 no puede simular ([P24](#p24)). Cuando captura un electrón nace el hidrógeno y salen 13,6 eV de luz; y sus dos imanes, el del protón y el del electrón, dan la raya de 21 cm con la que se cartografía la Galaxia, calculable con la densidad del electrón justo encima del protón.
**Las librerías**, por fin buenas noticias tras las de F10: **QuTiP** (espines: hoy dio 1006/994 en Stern–Gerlach y la vuelta de Rabi exacta), **qmsolve** (Schrödinger en malla con evolución temporal), **SymPy** (el hidrógeno exacto en fórmula cerrada) y **pyqint** (Hartree–Fock de verdad: H₂ a −1,1168 hartree frente a −1,1167 de libro, y el diente de sierra de la tabla periódica asomando ya con la base mínima). Las cuatro se instalan con `pip`, sin WSL. Lo único que no conviene es la malla 3D bruta: 27 minutos para el hidrógeno y encima desafinado (el 1s sale a −9,82 eV, porque una malla de 70³ puntos no resuelve el pico del Coulomb en el núcleo), mientras que la ecuación radial lo clava en 0,1 s. Se resuelve en radial y el 3D solo se pinta.
Y dos escalas nuevas para el laboratorio: el electrón vive en attosegundos y el núcleo entre microsegundos y miles de millones de años. Con F11 cubriremos 33 órdenes de magnitud en tiempo.
---
## P35 (05/09/2026) — F11a: por qué el electrón solo puede tener ciertas energías
`onda.html`, la primera fase dentro del átomo. Un electrón en una dimensión, cinco pozos (la caja, el muelle, el pozo doble, una caricatura del Coulomb y una barrera) y tres modos.
**El tiro: la cuantización es una condición de contorno.** La [ecuación de Schrödinger](Del_atomo_al_cromosoma.html#electrones)[(1)](#p23) estacionaria es una contabilidad de energía para una onda que no cambia de forma: lo que la onda tiene de arrugada (cinética) más el pozo (potencial) igual a E. Eliges una E cualquiera con las flechas e integras de izquierda a derecha. La onda sale del pozo y **se dispara** hacia +∞ o −∞. Una onda que se dispara no puede ser una probabilidad, porque no suma 1. Solo para ciertas E exactas la cola muere suavemente, y ahí el visor dice «¡encaja!». Abajo se dibuja lo que vale la onda en el borde según la energía: **cada cruce por cero es un nivel**. Nadie le dijo dónde estaban. Al lado, la bolita clásica de la misma energía, para la que *cualquier* E vale: esa es toda la diferencia entre una partícula y una onda.
Y los números que salieron del tiro son los del libro: la caja de 1 nm da 0,37603 · 1,50412 · 3,38427 · 6,01648 eV, que es n²·h²/(8mL²) con **cinco decimales**; el muelle con ħω = 1 eV da 0,5 · 1,5 · 2,5 …, equiespaciados y con el fondo a ½ħω, no a cero (energía del punto cero: ni a temperatura cero se para). El nivel n tiene siempre n − 1 nodos. Y el Coulomb aplanado da −18,2 · −7,5 · −4,1 · −2,5 eV: los saltos se acortan hacia arriba, como en el hidrógeno que viene en F11b.
**El pozo doble** es el H₂⁺ de F10 en una dimensión: dos pozos iguales y cada nivel se **desdobla** en dos, el simétrico y el antisimétrico, separados 0,43 eV por debajo de la barrera. Ese desdoblamiento *es* el túnel: un paquete metido en el pozo izquierdo se cuela al derecho y vuelve cada h/ΔE = 9,6 fs. Por encima de la barrera los niveles dejan de ir en pareja: ya no la notan.
**El tiempo: un paquete de onda.** La ecuación con tiempo, resuelta por Crank–Nicolson (la norma se queda en 1,000000 y la energía en 3,4901 eV con cuatro decimales durante los 11 fs enteros). Un paquete estrecho **se ensancha**: Δx·Δp arranca en 0,500 ħ, el mínimo exacto que permite la incertidumbre, y ya vale 0,9 ħ a los 0,17 fs y 1,3 ħ a los 0,34; la anchura medida (0,145 nm a 0,17 fs) sigue la fórmula del paquete libre (0,146) hasta que toca la pared. En la caja el paquete rebota, se deshace en un revoltijo… y **renace** entero cada h/E₁ = 11 fs: medido, a los 5,5 fs reaparece en +0,12 nm (había salido de −0,10) con casi su anchura inicial, y a los 11 fs vuelve a −0,13 nm. Solo pasa porque los niveles van exactamente como n², y es un fenómeno real que se ha visto en átomos de Rydberg. El paquete gaussiano del muelle, en cambio, oscila sin ensancharse: el estado coherente, lo más parecido a una bolita que la cuántica permite. Y en el pozo doble, el estado (ψ₁+ψ₂)/√2 metido a la izquierda pasa al pozo derecho y vuelve: 1 % → 47 % → **92 %** → 47 % → 1 % en los 9,6 fs que dicta h/ΔE.
**El túnel.** Un electrón de 2 eV contra una barrera de 3 eV y 0,30 nm. Clásicamente rebota siempre, y la bolita naranja lo hace. Cuánticamente, a los 4 fs el **18,9 %** de la probabilidad está al otro lado. La fórmula exacta de la barrera rectangular, promediada sobre las energías que lleva el paquete, dice 18,6 %. (Para la energía media sola dice 15,3 %: el paquete lleva componentes más rápidas que pasan más, y la T depende exponencialmente de la energía, así que promediar importa.) Con `+` la barrera se engorda 0,05 nm y la T se hunde: cae como e^(−2κw). Cuando el pozo sea el núcleo y la partícula una α, esto será la radiactividad, y ese exponencial explicará por qué un factor 2 en energía separa microsegundos de miles de millones de años ([P34](#p34)).
**El carril de rigor**: `f11_onda.py` calcula los mismos niveles con qmsolve y con SciPy, y el visor los enseña como anillos verdes sobre los puntos del tiro. Coinciden a la milésima de eV; qmsolve se desvía un 0,2 % en la caja porque pone muros de 10⁵ eV donde nosotros ponemos infinito, y eso también se dice.
---
*Cuaderno vivo: cada nueva pregunta conceptual del chat se añade aquí con su fecha.*