L'acqua è una delle sostanze più importanti sulla Terra e le sue proprietà fisiche e chimiche uniche sono essenziali per l'esistenza della vita. Fin dall'avvento della chimica computazionale, gli scienziati hanno lavorato per utilizzare modelli matematici per simulare il comportamento dell'acqua. Questi modelli non solo prevedono le proprietà fisiche dell'acqua, ma forniscono anche informazioni su come l'acqua reagisce in diversi ambienti e sul suo ruolo negli organismi viventi.
"Le proprietà dell'acqua sono strettamente correlate alla sua struttura molecolare, che può essere simulata e prevista utilizzando metodi di chimica computazionale."
I modelli idrici nella chimica computazionale vengono utilizzati principalmente per simulare aggregati di molecole d'acqua, acqua liquida e soluzioni acquose. Questi modelli si basano sulla meccanica quantistica, sulla meccanica molecolare, sui dati sperimentali o su una combinazione di questi metodi. Per imitare le proprietà specifiche delle molecole d'acqua, i ricercatori hanno sviluppato diversi tipi di modelli, che possono essere generalmente categorizzati in tre modi: (i) il numero di punti di interazione, chiamati "siti", (ii) il rigido o flessibile, e ( iii) se il modello include effetti di polarizzazione.
Nella simulazione dell'acqua, un approccio comune è quello di utilizzare un modello di solvente esplicito, cioè un modello basato su molecole specifiche. In alternativa a questi modelli espliciti, sono disponibili modelli di solvente impliciti, che utilizzano un modello continuo per trattare il comportamento dell'acqua. Un esempio di quest'area include il modello di solvente COSMO o il modello del continuo polarizzabile (PCM), o anche alcuni modelli di solventi misti.
Il modello rigido è considerato il modello più semplice per l'acqua e si basa su interazioni non legate. In questi modelli, le interazioni di legame vengono gestite implicitamente tramite vincoli globali. Le interazioni elettrostatiche sono modellate sulla base della legge di Coulomb, mentre le forze repulsive e di dispersione sono descritte utilizzando i potenziali di Lennard-Jones. Questi potenziali modelli, come TIP3P (potenziale molecolare trasferibile a tre punti) e TIP4P, sono rappresentati come:
E = ∑(kC * qi * qj / rij) + (A / rOO^12) - (B / rOO^6)
Dove kC è la costante elettrostatica, qi e qj sono le cariche parziali relative alla carica dell'elettrone e rij è la distanza tra i due atomi. In molti modelli di acqua, il termine di Lennard-Jones si applica solo alle interazioni tra atomi di ossigeno. I parametri geometrici dei vari modelli di acqua, come la distanza OH e l'angolo HOH, variano a seconda del modello.
"I modelli tridimensionali comunemente utilizzati, come TIP3P, funzionano bene nel calcolo delle prestazioni termiche specifiche."
Ad esempio, il modello SPC/E aggiunge una correzione di polarizzazione alla funzione di energia potenziale, rendendo la densità e la diffusione dell'acqua risultanti costanti migliori rispetto al modello SPC. Il modello TIP3P è ampiamente utilizzato nel campo di forza CHARMM e al modello originale sono state apportate lievi modifiche per renderlo più adatto alla simulazione di molecole biologiche.
Modelli flessibili vs. rigidiIl modello flessibile SPC dell'acqua è un modello tridimensionale riparametrizzato dell'acqua. A differenza del modello rigido SPC, il modello flessibile può descrivere correttamente la densità e la costante dielettrica dell'acqua nelle simulazioni di dinamica molecolare. Questo modello è implementato in diversi programmi di calcolo, come MDynaMix e Abalone.
Il modello a quattro siti migliora la distribuzione della carica delle molecole d'acqua aggiungendo un atomo fittizio vicino all'atomo di ossigeno del modello a tre siti. Il primo modello di questo tipo può essere fatto risalire al modello di Bernal-Fowler del 1933. Sebbene il modello fosse storicamente importante, non riproduceva molto bene le proprietà chiave dell'acqua.
Il modello TIP4P è ampiamente utilizzato nei software di chimica computazionale e svolge un ruolo fondamentale nella simulazione dei sistemi biomolecolari, mentre nuovi modelli dell'acqua come il modello OPC possono descrivere con maggiore accuratezza le proprietà elettriche dell'acqua.
Sebbene il modello a cinque bit abbia un elevato costo computazionale, negli ultimi anni ha fatto gradualmente progressi con l'introduzione del modello TIP5P. Il modello a cinque bit riproduce meglio la geometria del dimero d'acqua ed è in grado di catturare con precisione i dati sperimentali. Il modello a sei bit incorpora nei dati tutte le caratteristiche dei modelli precedenti ed è progettato specificamente per studiare i sistemi di acqua e ghiaccio.
"Nella chimica computazionale, la simulazione dell'acqua non è solo una sfida tecnica, ma anche una chiave per comprendere il funzionamento della vita."
Il costo computazionale del modello idrico aumenta con il numero di siti. Nelle simulazioni di dinamica molecolare, all'aumentare del numero di siti aumenta anche il numero di distanze interatomiche da calcolare. Tuttavia, lo sviluppo di questi modelli non è solo una narrazione matematica, ma un microcosmo di come l'acqua si comporta realmente in natura. Con il progresso della tecnologia, saremo in grado di trovare modelli che sveleranno altri misteri sull'acqua nel prossimo futuro?