El trabajo va bastante por libre: las investigaciones vecinas tratan del viento electrónico en el grafeno, de la limpieza ultravioleta del grafeno, de vórtices en luz líquida, y también de gotas sobre un cabello y de hidrógeno en el acero. Casi no hay vecinos directos en el tema de la disociación del agua sobre defectos del grafeno; el más cercano está a una distancia de 0.575. Esto significa que el trabajo llena un hueco: hasta ahora la disociación del agua en una vacante de grafeno solo se había estudiado para una molécula aislada en fase gaseosa, y la influencia del disolvente seguía siendo desconocida.
Los autores entrenaron un potencial equivariante MACE con datos de DFT híbrido con polarización de espín (revPBE0-D3). El conjunto de entrenamiento, de unas 9000 estructuras, incluía vacantes en fase gaseosa, configuraciones parcialmente solvatadas y configuraciones totalmente solvatadas. Para acelerar los sucesos raros usaron metadinámica con múltiples caminantes (multiple-walker well-tempered metadynamics) en el paquete de programas LAMMPS+PLUMED. Como variables colectivas emplearon los números de coordinación carbono–hidrógeno y carbono–oxígeno, que siguen el enlace de los fragmentos de agua con la vacante. La fase gaseosa se modeló con una molécula de agua en la celda y el sistema solvatado, con unas 60 moléculas de agua. El tiempo total de metadinámica fue de 6 ns para la fase gaseosa y de 40 ns para la disolución. Además se realizaron 2 ns de dinámica molecular sin sesgo para analizar los tiempos de retención del agua en los defectos.
En fase gaseosa el agua se disocia en la vacante a través de un único estado de transición concertado con una barrera de 28.4 kcal/mol, y forma el intermedio Hydroxy-SV. En disolución se abren dos canales competidores. El primero, el básico, lleva a la formación de SV-H y de un ion hidróxido en disolución con una barrera de 21.9 kcal/mol. El segundo, el ácido, da SV-OH y un ion hidronio con una barrera de 23.7 kcal/mol. Ambas barreras son más bajas que la de fase gaseosa. El intermedio ácido es más estable que el básico en unos 20 kcal/mol, lo que indica una preferencia termodinámica por la ruta ácida. El producto final SV-Ketone queda 74 kcal/mol por debajo de los reactivos en disolución, frente a 78 kcal/mol en fase gaseosa. El análisis de los tiempos de retención del agua mostró que Hydroxy-SV retiene el agua 55 ps y SV-Ketone 27 ps, mientras que la vacante desnuda y SV-Pyran son comparables al grafeno puro (unidades de picosegundos).
El trabajo se distingue de los vecinos en que por primera vez compara cuantitativamente los mecanismos de disociación del agua sobre un mismo defecto en fase gaseosa y en disolvente explícito a nivel de energía libre, usando un potencial entrenado por aprendizaje automático con la precisión del DFT híbrido. Los estudios vecinos o bien tratan otras propiedades del grafeno o bien no abordan la reactividad de los defectos en medio acuoso. Aquí, en cambio, se muestra que la solvatación cambia cualitativamente la ruta de reacción: en lugar de un único canal concertado surgen dos rutas iónicas estabilizadas por el disolvente. Esto subraya la necesidad de tener en cuenta explícitamente el disolvente al modelar reacciones en las interfaces de materiales 2D.
El tema puede desarrollarse hacia el estudio de efectos cooperativos entre vacantes vecinas: si reaccionan de forma independiente o si muestran una reactividad acoplada que pueda aprovecharse para química de superficie selectiva. También es prometedor extender el enfoque a otros materiales 2D defectuosos, por ejemplo el nitruro de boro hexagonal monocapa con fuertes efectos de carga superficial. La metodología de potenciales entrenados por aprendizaje automático con DFT híbrido y metadinámica puede aplicarse a una amplia clase de interfaces sólido–líquido reactivas.
Los resultados influirán en la comprensión de la corrosión de los electrodos de carbono en las pilas de combustible, en el diseño de membranas y dispositivos nanofluídicos basados en grafeno, y también en las estrategias de funcionalización de materiales carbonosos mediante química controlada de defectos.
El siguiente paso será estudiar la interacción de varias vacantes y su influencia conjunta en la disociación del agua. También es importante comprobar si los mecanismos encontrados se mantienen para otros defectos, como los bordes del grafeno o las divacantes.
El trabajo está directamente relacionado con el problema no resuelto del control de la mojabilidad y la reactividad de los materiales carbonosos en medio acuoso. Hasta ahora no estaba claro cómo influyen los defectos en la disociación del agua en condiciones de disolvente explícito, y este estudio cubre ese hueco para la vacante única de grafeno.
Una vacante única de grafeno es un birradical con dos electrones desapareados: uno en un orbital σ localizado y otro en un orbital π, lo que la hace parecida a un carbeno molecular. Esa naturaleza «de dos caras» permite que la vacante actúe a la vez como aceptora y como donante de densidad electrónica, y abre dos rutas de reacción distintas con el agua.
En la novela «La era del diamante», de Neal Stephenson, la nanotecnología utiliza estructuras de carbono para controlar procesos químicos a escala atómica. Aunque allí se habla de sistemas más complejos, la idea de los defectos puntuales como centros activos de transformaciones químicas conecta con el papel de las vacantes de grafeno en la disociación del agua.