Le risposte precedenti di mpprogram6771 e MSalters lo hanno risolto. Vorrei aggiungere che, poiché $ \ ce {CO2} $ è una molecola molto piccola, puoi, con un po 'di sforzo, impostare un piccolo valore numerico sperimentare per rispondere alla tua domanda e persino ottenere cariche parziali approssimative in ciascun atomo e momento di dipolo dell'intera molecola, utilizzando solo software gratuito / open source.
Innanzitutto, devi installare un software di modellazione molecolare in la tua macchina. Quello che mi piace di più è Avogadro. Ha una straordinaria usabilità e molte funzionalità per progettare e visualizzare i tuoi composti. Anche Ghemical era buono, ma sembra che non sia più mantenuto da anni. Non riuscivo più a farlo funzionare correttamente.
Nella mia macchina utilizzo Ubuntu MATE 18.04 (una variante GNU / Linux) come sistema operativo. Lì sono in grado di installare Avogadro con un semplice comando nel terminale:
sudo apt-get install avogadro
Con Avogadro puoi assemblare il $ \ ce {CO2} $ , che unisce l'atomo di carbonio ed entrambi gli atomi di ossigeno con doppi legami. Oltre all'editor molecolare, avrai bisogno di un altro software, in grado di prendere i dati sulla molecola che hai assemblato e fare una serie di calcoli di meccanica quantistica su di essa, per darti una risposta approssimativa alle tue domande.
Esiste una grande varietà di software di meccanica quantistica, come mostra questa pagina su Wikipedia. Sfortunatamente, IMHO il panorama degli strumenti gratuiti / open source in questo campo è frammentato e la maggior parte è molto indietro rispetto ad Avogadro in termini di usabilità, bloccato nel livello medio di facilità d'uso degli anni '80 (a volte a livello di compilazione da soli ), e le alternative proprietarie hanno licenze restrittive e / o sono allettanti e costose, fuori dalla portata di persone senza affiliazione istituzionale. L'Accademia tratta male i suoi creatori di strumenti volontari, come alcune persone straordinarie di matematica possono dirti in prima persona. Prima o poi dovremo rimediare. Abbiamo bisogno di un William Stein in chimica computazionale. Spero solo che riceva un trattamento migliore dopo aver intensificato il compito.
Tuttavia, tra i numerosi pacchetti supportati dal generatore di input di Avogadro, la mia raccomandazione è Psi4, per un principiante. È facile da installare come Avogadro, se utilizzi Ubuntu o qualsiasi distribuzione basata su Debian.
sudo apt-get install psi4
Hanno un sito ben documentato, con una sezione dedicata all'istruzione con progetti semplici e bacheche di messaggi amichevoli. La versione disponibile nel repository di Ubuntu è funzionale, ma abbastanza obsoleta, 1.1.5, a partire da marzo 2020. Se si vuole seriamente impararla, il mio consiglio è di scaricarla direttamente dal loro sito. L'ultima versione stabile di marzo 2020 è la 1.3.2. Ma per il bene di questa risposta, il valore predefinito del repository è sufficiente.
Dopo aver assemblato la tua molecola e aver eseguito alcune ottimizzazioni preliminari della geometria all'interno di Avogadro, puoi generare un file di testo di input preliminare con il suo plugin Psi4 nel menu Extra → PSI4 . La mia versione preliminare è iniziata in questo modo:
set basis aug-cc-pVDZmolecule {0 1 C -3.47367 0.73246 0.22361 O -2.43476 1.12414 -0.22175 O -4.51237 0.34053 0.66926} ottimizza ('B3LYP-D ')
Il plugin Avogadro per Psi4 è molto semplice, quindi dovremo mettere a punto il modello a mano. Una serie di buoni modelli che puoi modificare in base alle tue esigenze è un'ottima cosa da avere quando impari a usare un nuovo pacchetto. Dovremmo averne di più. Ma prima di tutto, vediamo cosa abbiamo sul nostro proto-input. Ha tre sezioni. La prima sezione specifica un set di base, aug-cc-pVDZ (i chimici computazionali amano banchettare con la zuppa di alfabeto). In breve, un set di base è un insieme di funzioni matematiche facili da calcolare, utilizzato per emulare gli orbitali atomici e molecolari reali e difficili da calcolare, in questo modo:
La seconda sezione ha le coordinate x, y, z di ogni atomo della molecola, e anche la sua carica complessiva (in questo caso 0) e la molteplicità (in questo caso 1, poiché tutti gli elettroni sono accoppiati). La terza sezione dice che tipo di informazioni vogliamo calcolare dalle nostre informazioni iniziali, in questo caso, la geometria ottimale della molecola (ottimizzare), e il meccanismo algoritmico scelto per elaborarla, in questo caso, B3LYP-D (un'altra porzione of alphabet soup), una variante della teoria del funzionale della densità (DFT).
Ho cambiato il modello generato da Avogadro come segue:
memory 4 Gbset base aug-cc-pVTZmolecola {0 1 C -3.47367 0.73246 0.22361 O -2.43476 1.12414 -0.22175 O -4.51237 0.34053 0.66926} ottimizza ('B3LYP-D') E, wfn = energia ('B3LYP-D', return_wfn = True ) oeprop (wfn, "MULLIKEN_CHARGES", "DIPOLE", title = "CO2 B3LYP-D")
Ho facoltativamente aumentato il limite sulla memoria di sistema a 4 GB, dal valore predefinito del sistema, poiché la mia macchina ha una buona quantità di memoria. Poiché la molecola è piccola e l'impatto sul runtime sarà probabilmente accettabile, ho anche cambiato il precedente set di base, aug-cc-pVDZ, con uno più dettagliato, aug-cc-pVTZ. Aggiunta anche una sezione che chiede a Psi4 di restituire un oggetto funzione d'onda (wfn) per il sistema, oltre alla sua energia (E). Infine, seguendo la guida sul manuale di Psi4 qui, ho aggiunto una sezione che richiede le nostre informazioni di interesse, le cariche parziali stimate su ciascun atomo, fornite da analisi Mulliken, e il momento di dipolo stimato sulla molecola $ \ ce {CO2} $ .
Ora possiamo salvare il file di testo con i nostri dati di input ed eseguire Psi4 nel terminale:
psi4 carbon_dioxide.in
Dopo un po 'di tempo, Psi4 finirà la corsa e restituirà i suoi risultati in un file di output denominato carbon_dioxide.out che contiene un'enorme quantità di informazioni. Ma la sezione di maggiore interesse per la tua domanda è proprio alla fine:
Proprietà calcolate utilizzando la matrice di densità CO2 B3LYP-D Momento del dipolo nucleare: (au) X: -0.0000 Y: 0.0000 Z : 0.0000 Momento dipolo elettronico: (au) X: 0.0000 Y: 0.0000 Z: -0.0000 Momento dipolo: (au) X: 0.0000 Y: 0.0000 Z: -0.0000 Totale: 0.0000 Momento dipolo: (Debye) X: 0.0000 Y: 0.0000 Z: -0.0001 Totale: 0.0001 Costi Mulliken: (au) Simbolo centrale Alpha Beta Spin Totale 1 C 2.80993 2.80993 0.00000 0.38015 2 O 4.09503 4.09503 0.00000 -0.19006 3 O 4.09504 4.09504 0.00000 -0.19008 Totale alfa = 11.00000, Totale beta = 11.00000, Totale carica = 0,00000 *** Psi4 terminato correttamente. Compra una birra a uno sviluppatore!
I risultati indicano esattamente la situazione che hai previsto intuitivamente, con entrambi gli atomi di ossigeno che allontanano la densità elettronica dall'atomo di carbonio centrale e l'atomo di carbonio che diventa leggermente positivo e gli atomi di ossigeno leggermente negativi. In effetti, siamo stati in grado di utilizzare il computer come una sorta di armatura potenziata per la mente.
All'inizio, la tua intuizione poteva solo fornire una guida vaga nella direzione del trasferimento della densità elettronica, dall'ossigeno al carbonio. Ora possiamo confermarlo e aumentare la nostra intuizione con stime numeriche, una perdita media di 0,38 elettroni nell'atomo di carbonio e un guadagno medio di 0,19 elettroni in ogni atomo di ossigeno. Meraviglioso.
Nonostante la separazione di carica, i risultati del nostro piccolo esperimento numerico indicano anche un momento di dipolo prossimo allo zero, come vediamo. Non ci dice esplicitamente perché. Ma la nostra intuizione geometrica suggerisce una via d'uscita. Poiché ci sono due atomi di ossigeno, l'effetto della separazione della carica su entrambi può annullarsi. L'output di Psi4 conferma che, poiché la carica parziale su ciascun atomo di ossigeno è la stessa entro quattro cifre decimali, ed entrambi prendono posizioni opposte in una geometria lineare.
C'è una molecola simile, ma senza la possibilità di cancellazione della separazione di carica, $ \ ce {CO} $ , monossido di carbonio, con un solo ossigeno. Per fare un confronto, ho creato il file di input equivalente.
memory 4 Gbset base aug-cc-pVTZmolecule {0 1 C -3.99710 1.44942 0.00000 O -2.86898 1.44942 0.00000} ottimizzare (' B3LYP-D ') E, wfn = energy (' B3LYP-D ', return_wfn = True) oeprop (wfn, "MULLIKEN_CHARGES", "DIPOLE", title = "CO B3LYP-D")
E lo ha eseguito.
psi4 carbon_monoxide.in
Ancora una volta i risultati indicano una certa misura della separazione di carica.
Proprietà calcolate utilizzando la matrice di densità CO B3LYP-D Momento di dipolo nucleare: (au) X: 0,0000 Y: 0,0000 Z: 0,0023
Momento di dipolo elettronico: (au) X: 0.0000 Y: 0.0000 Z: 0.0348 Momento di dipolo: (au) X: 0.0000 Y: 0.0000 Z: 0.0371 Totale: 0.0371 Momento di dipolo: (Debye) X: 0.0000 Y: 0.0000 Z: 0.0944 Totale : 0,0944 Costi Mulliken: (au) Simbolo centrale Alpha Beta Spin Totale 1 C 2.95397 2.95397 0.00000 0.09206 2 O 4.04603 4.04603 0.00000 -0.09206 Alpha totale = 7.00000, Beta totale = 7.00000, Carica totale = 0.00000 *** Psi4 uscito correttamente. Compra una birra a uno sviluppatore!
Ma questa volta il dipolo era diverso da zero, con un valore stimato intorno a 0,094 debye. L ' articolo di Wikipedia sul monossido di carbonio ci fornisce un valore misurato di 0,122 debye. Quindi abbiamo una stima inferiore di circa il 23% rispetto al valore reale. La differenza può sorgere o come una limitazione intrinseca del nostro modello (la scienza contro l'ingegneria), o perché ho armeggiato da qualche parte nell'input che ho dato a Psi4 o nelle mie ipotesi per trattare il problema (sempre molto probabile).
Sarebbe interessante controllare la letteratura sull'argomento, se si vuole andare più a fondo. Ad ogni modo, il contrasto nei risultati tra $ \ ce {CO2} $ e $ \ ce {CO} $ span > indica chiaramente la cancellazione reciproca per spiegare la mancanza di un dipolo in $ \ ce {CO2} $ .