Gradient dynamics model for chemically driven running drops

1. Il Problema

Il lavoro affronta la dinamica delle gocce liquide parzialmente bagnanti su un substrato solido, che sono in grado di auto-propulsarsi a causa di reazioni chimiche. Nello specifico, il focus è sulle gocce che trasportano una sostanza chimica che, attraverso reazioni reversibili, modifica l'umettabilità (wettability) del substrato sottostante.

• Contesto: Esistono modelli precedenti per gocce "in movimento" (running drops), ma spesso si basano su accoppiamenti chemomeccanici non termodinamicamente consistenti (ad esempio, assumendo un accoppiamento "unilaterale" dove la concentrazione influisce sulla pressione di disgiunzione, ma non viceversa sulla cinetica di reazione o diffusione).
• Limitazioni dei modelli esistenti: I modelli per il caso di copertura irreversibile del substrato (dove la sostanza chimica rimane permanentemente) sono stati resi termodinamicamente consistenti in lavori precedenti. Tuttavia, il caso reversibile (dove la sostanza chimica può desorbire nel fluido circostante) presenta sfide maggiori. I modelli precedenti per il caso reversibile non garantivano una consistenza termodinamica completa e non potevano descrivere correttamente la guida fuori dall'equilibrio necessaria per un moto sostenuto indefinito.
• Obiettivo: Sviluppare un modello termodinamicamente coerente per gocce in movimento guidate chimicamente in regime di adsorbimento reversibile, capace di descrivere la transizione da uno stato di riposo a uno stato di auto-propulsione sostenuta.

2. Metodologia

Gli autori adottano un approccio basato sulla dinamica dei gradienti (gradient dynamics), una struttura matematica che garantisce la consistenza termodinamica per sistemi dissipativi.

• Sistema Fisico: Si considerano gocce sessili confinate in uno spazio verticale stretto (gap) tra due piastre parallele. Questo confinamento è cruciale per derivare una descrizione coerente per un sistema chiuso, permettendo di tracciare le concentrazioni chimiche sia all'interno della goccia, sul substrato che nel fluido circostante.
• Variabili di Campo: Il modello utilizza quattro variabili di campo scalari:
  1. $h(\vec{x}, t)$: altezza locale del film liquido.
  2. $a(\vec{x}, t)$: densità delle particelle nella goccia (specie A).
  3. $b(\vec{x}, t)$: densità delle particelle adsorbite sul substrato (specie B).
  4. $c(\vec{x}, t)$: densità delle particelle disciolte nel fluido circostante (specie C).
• Reazioni Chimiche: Si modellano due reazioni reversibili monomolecolari:
  • $A \rightleftharpoons B$: Adsorbimento dalla goccia al substrato.
  • $B \rightleftharpoons C$: Desorbimento dal substrato nel fluido circostante.
  • Le velocità di reazione dipendono dall'altezza del film ($h$) per garantire che l'adsorbimento avvenga solo sotto la goccia e il desorbimento solo nelle zone libere.
• Energia Libera e Dinamica: Viene definito un funzionale di energia libera $F$ che include termini di bagnabilità (energia di bagnatura mesoscopica), tensione superficiale, entropia di miscelazione e gradienti di concentrazione. L'evoluzione temporale è descritta da equazioni di tipo gradiente che includono:
  • Trasporto spaziale (flusso di massa guidato da gradienti di potenziale chimico).
  • Scambi locali (reazioni chimiche).
  • Vincoli di conservazione della massa e del volume.
• Guida Fuori dall'Equilibrio (Chemostats): Per ottenere un moto sostenuto, il sistema chiuso viene trasformato in un sistema aperto trattando la goccia e il fluido circostante come chemostats (serbatoi chimici). Si impongono potenziali chimici esterni costanti e distinti ($\mu_a$ per la goccia, $\mu_c$ per il fluido).
  • Se $\mu_a \neq \mu_c$, il sistema è spinto fuori dall'equilibrio termodinamico.
  • La differenza di potenziale crea un flusso netto di particelle: adsorbimento sotto la goccia e desorbimento lontano da essa, mantenendo un gradiente di umettabilità locale.

3. Contributi Chiave

1. Coerenza Termodinamica: Il modello è il primo a fornire una descrizione termodinamicamente consistente per il caso di copertura reversibile del substrato, rispettando il bilancio dell'entropia e le relazioni di Onsager.
2. Meccanismo di Auto-Propulsione Sostenuta: Dimostrano che, a differenza del caso irreversibile (dove il moto è transitorio e la goccia si ferma), la combinazione di due chemostats con potenziali diversi permette un moto perpetuo della goccia su un dominio periodico.
3. Struttura Matematica Unificata: Derivano un modello ridotto che descrive l'evoluzione del profilo della goccia e dell'adsorbato, includendo un termine di reazione aggiuntivo rispetto ai modelli precedenti che cattura lo scambio con il fluido circostante.
4. Analisi delle Biforcazioni: Identificano il meccanismo matematico attraverso cui nasce il moto: una biforcazione di tipo drift-pitchfork (drift-pitchfork bifurcation).

4. Risultati Numerici e Analisi

Gli autori hanno simulato le equazioni derivate utilizzando il metodo degli elementi finiti (libreria oomph-lib) su un dominio periodico 1D.

• Dinamica Passiva vs. Attiva:
  • Caso Passivo ($\mu_a = \mu_c$): Il sistema evolve verso l'equilibrio termodinamico. La goccia si ferma, il gradiente di concentrazione iniziale si smorza e si forma un profilo di adsorbato simmetrico.
  • Caso Attivo ($\mu_a > \mu_c$): Il sistema rimane fuori dall'equilibrio. Il gradiente di concentrazione sotto la goccia viene amplificato: le particelle si accumulano sotto la goccia (rendendo il substrato meno bagnabile) e vengono rimosse dietro di essa. Questo crea un gradiente di umettabilità asimmetrico che spinge la goccia in avanti.
• Soglia di Movimento: Il moto inizia solo quando la differenza tra i potenziali chimici supera una soglia critica. Questo corrisponde a una differenza finita di umettabilità tra la zona sotto la goccia e quella esterna. A differenza dei gradienti di umettabilità imposti esternamente (dove anche gradienti infinitesimi causano moto), qui il gradiente è dinamico e richiede una soglia finita per rompere la simmetria.
• Biforcazioni:
  • L'emergere del moto è caratterizzato da una biforcazione di drift-pitchfork supercritica. La velocità della goccia ($v$) cresce come la radice quadrata della distanza dalla soglia di instabilità ($v \propto |\mu - \mu_{crit}|^{1/2}$).
  • L'analisi della dipendenza dalle velocità di reazione ($r_1, r_2$) mostra che il moto esiste solo in un intervallo intermedio: se le reazioni sono troppo lente o troppo veloci, la goccia rimane ferma.
  • È stato osservato che variando la diffusività dell'adsorbato ($D_b$), la biforcazione può diventare subcritica, portando a bistabilità (coesistenza di stati di riposo e di moto).

5. Significato e Implicazioni

• Avanzamento Teorico: Il lavoro colma un divario fondamentale nella modellazione dell'attività chimica nelle gocce, fornendo un quadro termodinamico rigoroso che mancava per i sistemi reversibili. Questo è essenziale per distinguere tra rilassamento verso l'equilibrio e stati stazionari fuori equilibrio sostenuti.
• Rilevanza Sperimentale: Il modello offre predizioni quantitative (come le soglie di potenziale chimico e le velocità) che possono essere confrontate con esperimenti su gocce catalitiche o sistemi biologici (es. condensati biomolecolari).
• Prospettive Future: Il modello è estensibile per includere fluidi con viscosità comparabili (accoppiamento idrodinamico completo), substrati liquidi, effetti Marangoni solutali e reazioni chimiche più complesse (es. reazioni di ordine superiore o cooperatività), aprendo la strada allo studio di pattern di moto complessi e caotici.

In sintesi, il paper stabilisce un nuovo standard per la modellazione di gocce attive chimicamente, dimostrando come la rottura controllata dell'equilibrio termodinamico tramite chemostats possa generare moto auto-sostenuto attraverso gradienti di umettabilità dinamici.