Los subneptunos son la clase de exoplanetas más común en la Galaxia, pero su estructura interna sigue siendo tema de debate. Tradicionalmente se pensaba que bajo la atmósfera de hidrógeno había un límite nítido con el magma de silicato. Sin embargo, cálculos recientes mostraron que a las presiones y temperaturas características del interior de los subneptunos, el silicato y el hidrógeno pueden mezclarse por completo. Esto plantea la necesidad de describir las propiedades termodinámicas de dicho fluido supercrítico, sin lo cual no es posible construir modelos correctos de la estructura interna de los planetas. El trabajo se enmarca en las investigaciones sobre equilibrios de fases y ecuaciones de estado para sistemas de silicato-hidrógeno, que se han desarrollado activamente en los últimos años.
Los autores realizaron cálculos de dinámica molecular de primeros principios en el marco de la teoría del funcional de la densidad, utilizando el paquete de software VASP y el funcional de intercambio-correlación PBEsol. Se modeló una celda de 243 átomos (27 Mg, 27 Si, 81 O, 108 H), lo que corresponde a la composición MgSiO3H4. Primero se fundió MgSiO3 cristalino y luego se añadió hidrógeno. Los cálculos se llevaron a cabo a temperaturas de 4000, 6000 y 8000 K y presiones de hasta ~600 GPa. Los datos obtenidos de presión y energía interna se ajustaron con la energía libre de Helmholtz, lo que permitió derivar la densidad, la entropía, la capacidad calorífica, la expansión térmica, el módulo de compresión volumétrica y el parámetro de Grüneisen a partir de una única superficie termodinámica.
El fluido supercrítico MgSiO3H4 resultó ser menos denso que el fundido de silicato seco en todas las condiciones estudiadas. La desviación de la mezcla ideal de volúmenes de MgSiO3 y H2 alcanza varios por ciento y aumenta con la presión. El análisis estructural mostró que, al aumentar la presión, el número de coordinación Si–H aumenta, mientras que la cantidad de enlaces H–H estables cae bruscamente: a 4000 K, la fracción de átomos de hidrógeno que participan en contactos H–H de larga duración disminuye por debajo de 0,1 a presiones de alrededor de 340 GPa. Esto indica una transición de hidrógeno molecular a una forma más atómica, integrada en la red de silicato. La aplicación de la ecuación de estado a modelos de subneptunos de 6 masas terrestres con el mismo contenido total de hidrógeno del 4% mostró que, dependiendo de la presión en el límite atmósfera-interior (1 GPa frente a 5 GPa), el radio del planeta cambia en un 5,7%.
El trabajo proporciona por primera vez una ecuación de estado autoconsistente para el fluido supercrítico de silicato-hidrógeno, lo que cubre un vacío importante en el modelado de subneptunos. A diferencia de trabajos vecinos que principalmente evaluaban equilibrios de fases o mezclas a nivel de modelos termodinámicos, aquí se obtienen datos directos de primeros principios sobre la densidad y las derivadas termodinámicas. Esto permite tener en cuenta cuantitativamente el hidrógeno disuelto en el interior y su influencia en el radio y la evolución térmica de los planetas.
La continuación natural es ampliar la composición con hierro, que puede influir fuertemente en la densidad, las condiciones redox y las transiciones de fase. También se necesitan verificaciones experimentales a altas presiones y temperaturas, que por ahora son extremadamente escasas. Los esfuerzos conjuntos de teóricos y experimentalistas permitirán construir ecuaciones de estado para el rango completo de composiciones, presiones y temperaturas relevantes para los subneptunos.
Los resultados afectarán directamente a los modelos de estructura interna y evolución de los subneptunos, así como a la interpretación de los datos de masas y radios de exoplanetas.
Incluir hierro en los cálculos y realizar mediciones experimentales de densidad y compresibilidad de fluidos supercríticos de silicato-hidrógeno.
El trabajo aborda un problema fundamental de la planetología: cómo determinar la composición y la estructura interna de los exoplanetas a partir de datos observacionales limitados. Muestra que ignorar la no idealidad de la mezcla de silicato e hidrógeno puede llevar a errores sistemáticos en las estimaciones del contenido de hidrógeno y del estado térmico del interior.
Esta sección está dirigida al autor del trabajo y a sus colegas. Los trabajos vecinos los encontró la búsqueda por significado en nuestro archivo, no por coincidencia de palabras. No evaluamos: mostramos el campo alrededor.
Qué hace el trabajo
Los autores realizaron cálculos de dinámica molecular de primeros principios para la mezcla supercrítica MgSiO3H4 a presiones de hasta 600 GPa y temperaturas de 4000–8000 K, obtuvieron la ecuación de estado y mostraron que la densidad de la mezcla es menor que la del silicato seco, y que el comportamiento se desvía de la mezcla ideal. Al aplicar la ecuación de estado a modelos de subneptunos de 6 masas terrestres, encontraron que la distribución de hidrógeno entre la atmósfera y el interior afecta el radio del planeta (diferencia del 5,7% para presiones límite de 1 y 5 GPa).
Dónde está su fuerza
La fortaleza del trabajo radica en obtener una ecuación de estado autoconsistente para el fluido supercrítico de silicato-hidrógeno desde primeros principios, lo que llena un vacío en el modelado de subneptunos y permite considerar la distribución real de hidrógeno en el interior.
Hacia dónde puede avanzar
Vale la pena probar
Por qué puede importar
Tener en cuenta la mezcla supercrítica de silicato e hidrógeno puede cambiar las estimaciones del radio de los subneptunos en varios por ciento, lo cual es importante para interpretar la relación observada entre radio y masa y para identificar la estructura interna de los exoplanetas.
Trabajos cercanos en nuestro archivo
El trabajo más cercano del archivo está a 0.65 por significado: el nivel habitual de nuestro corpus.
A una presión de unos 340 GPa y una temperatura de 4000 K, solo uno de cada diez átomos de hidrógeno mantiene un enlace con otro hidrógeno durante más de 100 femtosegundos; los demás ya están integrados en la red de silicato.