Base locale e derivata covariante (prima parte)

Come è noto un sistema di riferimento cartesiano è formato da n rette ortogonali che si intersecano in un punto O detto origine, ognuna delle rette è orientata e riporta una unità di misura: in questo modo è possibile identificare qualsiasi punto dello spazio euclideo Rⁿ con una n-upla di numeri reali (x¹, x²,..., xⁿ) in modo univoco (vedi Wikipedia).

Generalizzando è possibile costruire geometricamente, a partire da un sistema di riferimento cartesiano, un altro riferimento qualsiasi detto curvilineo, che avrà lo stesso numero di coordinate ma nel quale le linee coordinate sono generalmente delle curve (vedi Wikipedia).

Ad esempio, come avevamo già visto nel post "Trasformazioni di basi, vettori e... co-vettori!", consideriamo un sistema di coordinate cartesiano bidimensionale (x¹,x²) e un nuovo sistema di coordinate curvilinee (x¹,x²) che sono note in funzione delle prime*:

x¹=x¹(x¹,x²)   ,   x²=x²(x¹,x²)

e dove vale anche la trasformazione inversa:

x¹=x¹(x¹,x²)   ,   x²=x²(x¹,x²)

ed inoltre tali funzioni sono per definizione differenziabili (funzioni lisce).
Nota: le coordinate devono essere indipendenti e quindi ∂xⁱ/∂x^j=δ_ij cioè ∂xⁱ/∂x^j=0 se i≠j e ∂xⁱ/∂x^j=1 solo se i=j (δ_ij è la delta di Kronecker). 

Consideriamo ad esempio il riferimento cartesiano rappresentato in figura dove è indicata la retta x¹ e il relativo vettore di base unitario e₁ con origine nel punto O da cui parte una linea curva coordinata x¹ con vettore di base unitario e₁ ad essa tangente come illustrato in figura:

Ora si osservi che tra il tratto infinitesimo dx¹ della retta x¹ e il tratto infinitesimo dx¹ tangente alla linea curva x¹ (che quindi approssima la linea in quel punto) esiste la seguente relazione trigonometrica:

dx¹=dx¹cosα

e quindi la proiezione del vettore di base e₁ sulla retta x¹ è pari a

|e₁|cosα=|e₁|(dx¹/dx¹)

dove |e₁|=1 è il modulo unitario di e₁ mentre α è l'angolo tra dx¹ e dx¹.
In modo equivalente la proiezione del vettore di base e₁ sulla retta x² è:

|e₁|sinα=|e₁|(dx²/dx¹)

valendo di nuovo la relazione trigonometrica:

dx²=dx¹sinα.

Perciò la nuova base locale e₁ espressa in funzione delle basi cartesiane e₁ ed e₂ (che ricordiamo sono per semplicità tutti vettori di modulo 1) è:

e₁=(|e₁|cosα)e₁+(|e₁|sinα)e₂

e quindi utilizzando le relazioni ricavate sopra

e₁=(∂x¹/∂x¹)e₁+(∂x²/∂x¹)e₂

avendo posto |e₁|=|e₂|=1 e ∂x¹/∂x¹=cosα , ∂x²/∂x¹=sinα.
Nota: abbiamo indicato le derivate parziali dato che le coordinate sono funzioni di più variabili.

Si osservi però che dal rapporto trigonometrico di due infinitesimi come dx¹/dx¹ e dx²/dx¹ siamo passati alle derivate parziali ∂x¹/∂x¹ e ∂x²/∂x¹ supponendo (correttamente) che la variazione di x¹ rispetto a x¹ e quella di x² sempre rispetto a x¹ siano rispettivamente pari a cosα e sinα.
Nota: ciò è vero poiché x¹ approssima (nell'origine) la coordinata di un riferimento ruotato di un angolo α rispetto a quello cartesiano (x¹, x²)**.

In modo analogo per la base e₂ si ha (anche se non è mostrato in figura):

e₂=(∂x¹/∂x²)e₁+(∂x²/∂x²)e₂

dove ricordiamo che le basi sono state tutte normalizzate:

|e₁|=|e₂|=1   e   |e₁|=|e₂|=1.

Nota: per costruzione geometrica le basi (e₁,e₂) sono tangenti alle linee coordinate (x¹,x²) e ciò vale in generale per più coordinate.

Se viceversa volessimo derivare le basi cartesiane e₁ e e₂ a partire da quelle curvilinee e₁ e e₂ un ragionamento analogo ci porterebbe ad ottenere:

e₁=(∂x¹/∂x¹)e₁+(∂x²/∂x¹)e₂

e₂=(∂x¹/∂x²)e₁+(∂x²/∂x²)e₂.

che rappresentano le relazioni inverse di quelle prima ottenute.

Quindi riscrivendo quanto abbiamo ottenuto sopra in forma più generale (applicando la notazione di Einstein sugli indici ripetuti) si ha:

e_i=(∂x^j/∂xⁱ)e_j   ,   e_i=(∂x^j/∂xⁱ)e_j

dove gli indici i,j=1,..., n indicano il numero di coordinate e le relative basi.
Nota: si osservi che gli elementi (∂x^j/∂xⁱ) e (∂x^j/∂xⁱ) definiscono rispettivamente la matrice jacobiana e la sua inversa.

Ma facciamo subito un esempio introducendo le coordinate curvilinee (r,θ) (ponendo cioè x¹=r e x²=θ) dove r≥0 è la distanza dall'origine (polo) mentre 0≤θ≤2π è l'angolo tra r e l'asse X (coordinate polari).

Per calcolare le nuove basi e₁=e_r e e₂=e_θ è utile definire le coordinate cartesiane (x,y) in funzione delle nuove coordinate (r,θ):

x(r,θ)=rcosθ  ,   y(r,θ)=rsinθ.

Possiamo quindi ottenere le nuove basi e_r e e_θ utilizzando le equazioni alle derivate parziali ottenute sopra (dove x¹=r e x²=θ):

e_r=(∂x/∂r)e_x+(∂y/∂r)e_y

e_θ=(∂x/∂θ)e_x+(∂y/∂θ)e_y.

Perciò le basi del nuovo riferimento curvilineo di coordinate (r,θ) espresse in funzione delle basi cartesiane note e_x ed e_y sono:

e_r=cosθe_x+sinθe_y

e_θ=-rsinθe_x+rcosθe_y

ed inoltre essendo perpendicolari il loro prodotto scalare è nullo:

<e_r,e_θ>=-rcosθsinθ+rsinθcosθ=0.

Nota: la base e_θ non è unitaria poiché risulta |e_θ|=r tuttavia possiamo porre come base unitaria θ=-sinθe_x+cosθe_y e quindi e_θ=rθ.

È importante osservare che le basi e_r e e_θ dipendono dalle coordinate (r,θ) come in effetti capita generalmente per le basi curvilinee: viceversa le basi cartesiane e_x ed e_y sono sempre le stesse in ogni punto dello spazio.

Come vedremo nel prossimo post le relazioni ricavate sopra saranno utili per definire la trasformazione di un vettore affinché resti invariato quando passa da un riferimento ad un altro e definiremo la sua derivata covariante.

(*) Ricordiamo che gli apici indicano entità che si trasformano in modo controvariante (come le componenti di un vettore) mentre i pedici indicano entità che si trasformano in modo covariante (come ad esempio le relative basi) come descritto in "Trasformazioni di basi, vettori e... co-vettori!".
(**) In due dimensioni la trasformazione delle coordinate per una rotazione degli assi è: x¹=x¹cosα-x²sinα ,  x²=x¹sinα+x²cosα da cui segue subito ∂x¹/∂x¹=cosα , ∂x²/∂x¹=sinα come già derivato sopra.

[Una ottima esposizione di questi concetti si trova nella Playlist Video di Dermot Green - Queen's University Belfast]

La Sfera di Bloch

È noto che lo stato quantistico di un qubit (quantum bit) è così definito:

|Ψ>=α|0>+β|1>

dove α e β sono numeri complessi che devono soddisfare la condizione:

|α|²+|β|²=1

per la normalizzazione a 1 della probabilità complessiva P=|α|²+|β|².
Nota: i coefficienti |α|² e |β|² rappresentano rispettivamente la probabilità che si verifichi lo stato |0> oppure |1>.

Ricordiamo che valendo la formula di Eulero e^iø=cosø+isinø, in generale per un numero complesso x+iy si può scrivere:

x+iy=r(cosø+isinø)=re^iø

dove r=(x²+y²)^1/2 è il modulo e ø=arctan(y/x) è l'angolo tra r e la sua proiezione sull'asse X (pari a rcosø). Inoltre ricordiamo che vale la relazione:

|x+iy|²=|re^iø|²=r²(e^iøe^-iø)=r²

poiché |e^iø|²=(e^iøe^-iø)=1 essendo e^-iø il complesso coniugato di e^iø.
Nota: è noto che il modulo al quadrato di un numero complesso è dato dal prodotto di quel numero per il suo coniugato.

Quindi possiamo porre per i coefficienti complessi α e β:

α=r₁e^iø₁   e   β=r₂e^iø2

dove ø₁ e ø₂ sono due angoli qualsiasi compresi tra 0 e 2π.
Se inoltre poniamo r₁=rcosθ e r₂=rsinθ con r=1 si ha come richiesto:

|α|²+|β|²=|r₁e^iø₁|²+|r₂e^iø₂|²=cosθ²+sinθ²=1

poiché |e^iø₁|²=|e^iø₂|²=1 come visto sopra.

Quindi per lo stato di un qubit vale la seguente relazione:

|Ψ>=α|0>+β|1>=cosθ(e^iø₁)|0>+sinθ(e^iø₂)|1>

definita in funzione degli angoli θ, ø₁ e ø₂ (di cui daremo di seguito una rappresentazione geometrica).

Si osservi però che moltiplicando lo stato |Ψ> per e^-iø₁ si ottiene un nuovo stato |Ψ'>=e^-iø₁|Ψ>:

|Ψ'>=e^-iø₁|Ψ>=e^-iø₁(α|0>+β|1>)=cosθ|0>+sinθ(e^iø)|1>

dove ø=ø₂-ø₁. Tuttavia la probabilità degli stati |0> e |1> resta invariata infatti:

|e^-iø₁α|²=(e^-iø₁α)(e^-iø₁α)*=(e^-iø₁e^iø₁)(αα*)=|α|²

e lo stesso vale per |e^-iø₁β|²=|β|².
Nota: qui il simbolo (*) indica il valore coniugato del numero complesso.

Perciò possiamo omettere il termine e^-iø₁ dallo stato |Ψ'> poiché non ha effetti osservabili sperimentalmente e possiamo riscrivere:

|Ψ'> => |Ψ>=cosθ|0>+sinθ(e^iø)|1>.

Infine per ottenere una rappresentazione geometrica riscriviamo |Ψ> in funzione delle variabili x, iy e z ponendo:

cosθ=z

sinθ(e^iø)=sinθcosø+isinθsinø=x+iy

o in modo equivalente

x=sinθcosø
y=sinθsinø
z=cosθ.

Ciò significa che se x, y e z vengono interpretate come coordinate, esse rappresentano le coordinate polari di una sfera di raggio unitario* |Ψ|=1:

Perciò lo stato |Ψ> in funzione delle coordinate x, iy e z diventa:

|Ψ>=cosθ|0>+sinθ(e^iø)|1>=z|0>+(x+iy)|1>

quindi le coordinate di |Ψ> al quadrato, cioè z² e |(x+iy)|², rappresentano rispettivamente la probabilità che si verifichi lo stato |0> oppure |1>.
Nota: si ricordi che in generale risulta |(x+iy)|²=(x+iy)(x-iy)=x²+y².

Si noti che quando θ=0 risulta |Ψ>=|0> mentre se θ=π/2 si ha** |Ψ>=|1> quindi per una descrizione completa della sfera possiamo porre:

|Ψ>=cos(θ/2)|0>+sin(θ/2)(e^iø)|1>

in modo che possiamo far variare 0≤θ≤π e 0≤ø<2π.

Questa rappresentazione geometrica*** dello spazio degli stati puri di un sistema quantistico a 2 stati è detta sfera di Bloch (vedi Wikipedia).

(*) Risulta |Ψ|²=<Ψ|Ψ>=(<0|α*+<1|β*)|(α|0>+β|1>)=α*α+β*β=1 essendo <0|0>=1, <0|1>=0, <1|0>=0 e <1|1>=1 per l'ortogonalità di |0> e |1>.
(**) Se poniamo θ=π/2 allora |Ψ>=e^iø|1> che è equivalente a |Ψ>=|1> essendo |e^iø|² =1 (cioè la probabilità dello stato |1> non cambia).
(***)  Correttamente risulta per le probabilità dei due stati:
|cos(θ/2)|²+|sin(θ/2)(e^iø)|²=cos²(θ/2)+sin²(θ/2)=1 essendo |e^iø|² =1.  

È possibile un modello esteso dell'elettrone?

È noto che nel Modello Standard le particelle elementari vengono considerate a tutti gli effetti come puntiformi, tuttavia ciò porta a valori infiniti per alcune quantità che le caratterizzano, legate alle mutue interazioni tra particelle e campi.
Nota: con l'aggettivo puntiforme si intende una particella elementare priva di qualsiasi struttura interna.

Vediamo infatti come si definisce l'energia di una particella carica di massa a riposo m nella teoria classica (relativistica), considerando l'energia del campo elettrostatico da essa generato oltre alla sua energia di massa (vedi Wikipedia):

mc²=m₀c²+(1/2)∫ ε₀E²dV     (1.1)

dove m₀ è la massa nuda della particella priva di campo elettrico* mentre E=F/e=e/4πε₀R² è il campo elettrico a distanza R dalla particella, F è la forza di Coulomb ed ε₀ è la permittività elettrica del vuoto.
Nota: in tutti i post useremo, anche se non dichiarate, le unità di misura del sistema metrico internazionale SI.

Se, ad esempio, integriamo tra il raggio r ipotetico di un elettrone e il volume V=(4/3)πR³ che lo circonda all'infinito si ha:

(1/2)∫ ε₀E²4πR²dR=(1/2)∫ (e²/4πε₀R²)dR     (1.2)

posto dV=4πR²dR ed essendo E=e/4πε₀R².
Quindi calcolando l'integrale tra il raggio r ed infinito si ottiene:

-(e²/8πε₀R)|_r^∞=e²/8πε₀r     (1.3)

per cui l'energia complessiva dell'elettrone con il suo campo è (secondo la eq.1.1):

mc²=m₀c²+e²/8πε₀r     (1.4)

ciò significa che se l'elettrone è puntiforme (cioè r -> 0) allora la sua energia mc² tende come anticipato ad infinito!

Nelle teorie di campo quantistiche l'approccio è più complesso di quanto esposto qui, tuttavia si devono comunque usare procedure matematiche di rinormalizzazione per eliminare le divergenze che insorgono nei calcoli, come ad esempio quello che determina l'anomalia magnetica dell'elettrone.

In un blog dedicato mostreremo come sia possibile definire, secondo una proposta dell'autore, un modello esteso dell'elettrone dove tutte le quantità che lo definiscono sono finite.

 

Dovremo però tenere conto, nello sviluppo del modello e.m. esteso,  che le attuali misure sperimentali indicano una dimensione della carica elettrica non superiore a circa 10^-19 metri, che è la migliore risoluzione degli odierni acceleratori di particelle (tale risoluzione è dell'ordine della lunghezza d'onda λ delle particelle-sonda)**. 

Nota: come vedremo in un Blog dedicato la carica elettrica del modello esteso è assunta come puntiforme mentre è la massa ad essere distribuita in modo esteso sulla superficie del modello.

(*) La massa nuda m₀ è un parametro libero della teoria che non è possibile misurare direttamente: non possiamo in effetti separare una particella carica dal campo elettrico che essa stessa genera.

(**) Dato che le particelle circolano negli acceleratori a velocità v prossime a c si ha E=hc²/λv≈hc/λ da cui λ≈hc/E e per E=14 TeV (energia massima oggi raggiungibile) si ha λ≈10^-19 m (si ricordi che 1 eV≈1,6x10^-19 Joule).
[Ricordiamo che per una particella di massa m si ha E=mc² e p=mv da cui p=Ev/c² ed essendo p=h/λ segue E=hc²/λv]

ATTENZIONE
Per il seguito di questo post vedi il Blog: Electron Extended Model dove verrà proposto dall'autore un modello dell'elettrone non puntiforme.

INDICE DEI POST

I Teoremi di Gödel, l'IA e... un'ipotesi di lavoro! (seconda parte)

Nel post precedente (a cui rimandiamo) abbiamo introdotto i Teoremi di incompletezza di Gödel con l'obiettivo di mostrare poi (cioè in questo post) come essi ci pongano davanti ad una disgiunzione: o la mente umana è equivalente ad una macchina di Turing per quanto complessa oppure siamo in presenza di un fenomeno completamente nuovo, mai studiato prima.
Nota: vedi anche l'articolo "La disgiunzione di Gödel" di F. Beccuti.

Introduciamo quindi quella che è stata definita Macchina di Turing: in pratica questo termine indica uno qualsiasi degli attuali computer poiché essi sono realizzazioni fisiche di questa macchina ideale e universale in grado di eseguire qualsiasi algoritmo si possa formalizzare.

Come è noto si è dimostrata la perfetta equivalenza tra ogni sistema formale assiomatizzabile S e la macchina ideale di Turing: cioè è possibile programmare un computer che produca tutti e soli i teoremi di un dato sistema S e, viceversa, qualsiasi programmazione di un computer che produce formule, può essere rappresentata da un sistema formale S che derivi gli stessi risultati.

Quindi la scommessa dell'intelligenza artificiale è proprio quella di supporre che l'insieme delle capacità cognitive del nostro cervello, in particolare il processo del pensiero razionale, possa essere completamente riprodotto ed espresso da un programma evoluto per computer.

L'obiezione più nota a questo programma di ricerca è quella del filosofo Lucas nel celebre articolo "Menti, Macchine e Gödel" (del 1961):
"Data qualsiasi macchina che sia coerente e capace di fare semplice aritmetica, c'è una formula che essa è incapace di produrre come vera - cioè la formula è indimostrabile nel sistema - tuttavia noi la possiamo vedere come vera. Perciò nessuna macchina può essere un modello completo o adeguato della mente, le menti sono essenzialmente differenti dalle macchine"*.

Questa tesi segue proprio dall'argomento di incompletezza di Gödel, in particolare dal primo teorema (vedi il precedente post), ed è confermata dal Teorema di indefinibilità di Tarski (del 1936) che afferma che non è possibile definire la nozione di verità all'interno di un sistema formale.
Nota: si può definire la nozione di verità solo facendo una meta-analisi al di fuori del sistema, ad esempio usando la logica del secondo ordine.

Quindi sembrerebbe stabilita la tesi di Lucas secondo cui le nostre capacità cognitive, in particolare quelle che determinano il pensiero razionale, sono di certo superiori a quelle di una qualsiasi macchina o computer.

Tuttavia dobbiamo ricordare che il teorema di Gödel fa in effetti una affermazione che è del tutto condizionale:

"Se S è coerente allora G non è dimostrabile".

Ma la nostra mente è veramente in grado di riconoscere se un qualunque sistema formale è coerente dato che questa proprietà non può essere provata all'interno di un qualsiasi sistema?
Nota: se S non è coerente si può dimostrare G (ma anche non-G) quindi la mente potrebbe essere un sistema incoerente e dimostrare che G è vera.

Inoltre ciò dovrebbe valere per qualsiasi sitema formale (come ad esempio sistemi più complessi che includono gli assiomi dell'infinito), perciò non è detto che la mente umana riesca sempre a riconoscere che un sistema è coerente.
Nota: la mente umana potrebbe essere un sistema coerente che non può dimostrare G (e quindi è incompleta) ma che non sa di essere coerente.

Lo stesso Gödel, che non era proprio un meccanicista, affermò nella Gibbs Lecture (del 1951), che potrebbe essere che "la mente umana (nel regno della matematica pura) [...] sia dunque equivalente ad una macchina finita che è incapace di comprendere interamente il suo funzionamento".

In definitiva, chi si occupa di intelligenza artificiale o di processi cognitivi e apprendimento, è costretto a fare una ben definita scelta o ipotesi di lavoro:
a) la mente umana non è riducibile ad una macchina di Turing che computa, quindi dobbiamo studiare le sue capacità cognitive in modo del tutto nuovo, poiché non possiamo trattarla come se fosse un oggetto computazionale**;
oppure
b) il nostro cervello funziona come un computer per quanto evoluto, tuttavia se la nostra mente è coerente, siamo costretti ad accettare che ci siano dei problemi matematici irresolubili, come ad esempio dimostrare la propria coerenza interna***.
Nota: per approfondire l'interessante tema mente-cervello vedi l'ottimo articolo di Paul e Patricia Churchland "Il problema mente-cervello".

(*) Per completare il sistema S potremmo aggiungere G come assioma, si otterrebbe però un sistema S' in cui c'è una nuova formula G' indecidibile e così via, senza risolvere il problema.
(**) Qui il punto è proprio quello di voler attribuire alla mente un carattere diverso da quello computazionale (e non tanto la sua eventuale somiglianza ad un computer che è solo un modello interpretativo).
(***) Se la mente segue le leggi della fisica e se queste possono essere trattate computazionalmente, allora cervello e mente sono elementi complementari: la mente (software) è una funzione del cervello (hardware).

(Per chiarimenti su questo post vedi l'ottimo video di Francesco Berto)

I Teoremi di Gödel, l'IA e... un'ipotesi di lavoro! (prima parte)

In questo e nel prossimo post vogliamo mostrare come i due Teoremi di incompletezza di Gödel (del 1931), sebbene non vietino in alcun modo che l'intelligenza artificiale si possa realizzare (nel senso di seguito specificato), ci impongono tuttavia di operare una scelta, o meglio un'ipotesi di lavoro.

Qui con il termine Intelligenza Artificiale (IA) intendiamo un suo aspetto peculiare, secondo cui un sistema meccanico che computa potrebbe pensare in modo umano; in effetti se si suppone che la mente umana non è altro che una macchina computazionale, la tesi dell'IA ne discende direttamente.

Quindi l'IA suppone che pensare è computare e in particolare si pone l'obiettivo di realizzare una macchina che possa pensare umanamente, in modo cioè che "il processo che porta il sistema intelligente a risolvere un problema ricalchi quello umano" (vedi Wkipedia).

Tuttavia prima di introdurre i teoremi di incompletezza, dobbiamo definire cosa si intende con sistema formale (vedi Wkipedia):
"In logica matematica la nozione di sistema formale è utilizzata per fornire una definizione rigorosa del concetto di dimostrazione"; in pratica un sistema formale è un insieme di regole per costruire dimostrazioni.
Nota: si suppone che il sistema sia corretto, cioè se gli assiomi sono veri i teoremi che si deducono con le regole di inferenza sono anch'essi veri.

In breve il problema che Gödel riesce ad esprimere in modo formale nei suoi teoremi è quello dell'autoreferenza (vedi Wikipedia) che si presenta quando una proposizione fa una affermazione su se stessa in modo circolare; un problema già noto agli antichi greci come il Paradosso del mentitore.

Si consideri quindi un sistema formale S, evoluto almeno quanto quello piuttosto semplice dell'aritmetica di Peano; Gödel riesce a formalizzare all'interno del sistema S la seguente frase G che afferma (di se stessa):

(G): G non è dimostrabile in S.

Nota: grazie alla fattorizzazione in numeri primi è possibile assegnare ad una qualsiasi frase formale un numero univoco detto numero di Gödel.

Ora se si suppone che S sia un sistema formale corretto e quindi prova solo cose vere, allora G non è dimostrabile, dunque G è vera: ma allora esiste una verità G che il sistema S non può dimostrare!
Nota: se G fosse dimostrabile allora G (che dice di non essere dimostrabile) sarebbe falsa e quindi S non sarebbe corretto perché dimostra una falsità.

Inoltre se G è vera la sua negazione non-G è falsa (per definizione di negazione), ma allora il sistema S (che prova solo cose vere) non può provare nemmeno non-G: dunque l'enunciato G è indecidibile* in S!

-> Enunciamo quindi il primo teorema di Gödel (vedi Wikipedia):

In ogni teoria matematica S sufficientemente espressiva da contenere l'aritmetica, esiste una formula ${\displaystyle \varphi }$G tale che, se S è coerente**, allora né ${\displaystyle \varphi }$ G né la sua negazione ${\displaystyle \lnot \varphi }$ non-G sono dimostrabili in S.

Nota: con coerente si intende che S non è contraddittorio, d'altra parte se il sistema è corretto è anche coerente (non dimostrando falsità).

Si ricordi che se un sistema è incoerente si può dimostrare una certa proposizione P e la sua negazione non-P ma se così fosse qualsiasi proposizione potrebbe essere dimostrata vera (vedi Wikipedia): sarebbe quindi opportuno riuscire a dimostrare in modo certo la coerenza di S.

Tuttavia Gödel riuscì a mostrare formalmente che il seguente enunciato:

"Se S è coerente allora ciò implica G" 

si può dimostrare in S. Ma allora S non può dimostrare la sua coerenza altrimenti G sarebbe dimostrabile, e ciò è escluso dal primo teorema.

-> Ecco quindi il secondo teorema di Gödel (vedi Wikipedia): 

Sia S una teoria matematica sufficientemente espressiva da contenere l'aritmetica: se S è coerente, non è possibile provare la coerenza di S all'interno di S.

Nota: la coerenza dell'aritmetica fu poi dimostrata nel 1936 in ambito metamatematico da Gerhard Gentzen grazie agli ordinali transfiniti.

Nel prossimo post mostreremo come le argomentazioni espresse nei teoremi di Gödel non siano conclusive sulla possibile realizzazione dell'IA come sopra specificato, ma ci impongano una scelta ben precisa nell'approccio alla comprensione delle nostre capacità cognitive e quindi della nostra mente.

(*) Il risultato è notevole: si dimostra l'indecibilità di una formula G nel sistema S alla quale è tra l'altro collegata la coerenza di S (vedi oltre).

(**) In realtà Gödel richiese la w-coerenza di S che è più forte della sola coerenza, ma poi Rosser dimostrò che non era necessaria.

(Per chiarimenti su questo post vedi l'ottimo video di Francesco Berto)

Stati misti, intrecciati e... il gatto!

Come anticipato nel post "Stati puri, miscele e sovrapposizioni!", ora analizziamo un sistema composto da due elettroni e verifichiamo se si tratta di uno stato di spin puro o misto grazie alla matrice densità prima definita*. 

Consideriamo ad esempio un sistema composto da due elettroni preparati separatamente nei seguenti stati di spin (dove u e d sta per up e down):

|Ψ>=ψ_u|u>+ψ_d|d>   e   |Φ>=φ_u|u>+φ_d|d>.

Lo stato prodotto che descrive il sistema combinato è:

|ΨΦ>=(ψ_u|u>+ψ_d|d>)⊗(φ_u|u>+φ_d|d>)

quindi sviluppando il prodotto tensoriale indicato con ⊗ si ottiene:

|ΨΦ>=ψ_uφ_u|uu>+ψ_uφ_d|ud>+ψ_dφ_u|du>+ψ_dφ_d|dd>.

con le condizioni di normalizzazione:

ψ_uψ_u+ψ_dψ_d=1   e   φ_uφ_u+φ_dφ_d=1.

Nota: ψ_u è il complesso coniugato di ψ_u e lo stesso vale per gli altri valori.

Tuttavia si osservi che in generale un sistema composto da due elettroni è descritto dal seguente stato di spin:

|Ψ>=ψ_uu|uu>+ψ_ud|ud>+ψ_du|du>+ψ_dd|dd>

che non è sempre rappresentabile da uno stato prodotto (vedi sopra) e per il quale vale la condizione di normalizzazione:

ψ_uuψ_uu+ψ_udψ_ud+ψ_duψ_du+ψ_ddψ_dd=1.

Nota: questo stato combinato è detto stato entangled (o intrecciato) quando non può essere fattorizzato in due stati separati.

Ad esempio consideriamo una coppia di elettroni, preparata con spin opposti, il cui stato combinato non fattorizzabile è:

|Ψ>=(1/2)^1/2|ud>+(1/2)^1/2|du>

dove la somma degli stati |ud> e |du> rappresenta due coppie di elettroni con spin opposti in sovrapposizione quantistica, mentre il fattore (1/2)^1/2 indica che la misura di uno dei due stati è equiprobabile poiché:

ψ_udψ_ud=ψ_duψ_du=1/2.

Nota: possiamo ad esempio pensare al caso descritto nell'esperimento EPR (per chiarimenti vedi il post "Un esperimento chiave: EPR").

Calcoliamo quindi la matrice densità, già introdotta nel post "Stati puri, miscele e sovrapposizioni!", che è così definita:

ρ=|Ψ><Ψ|=(1/2)(|ud>+|du>)(<ud|+<du|)

dalla quale svolgendo il prodotto si ottiene:

ρ=(1/2)(|ud><ud|+|ud><du|+|du><ud|+|du><du|).

Premesso che indicheremo i vettori colonna come vettori riga trasposti, scegliamo due vettori di base: |u>=(1,0)^T e |d>=(0,1)^T (dove T indica la matrice trasposta)** e sviluppiamo i prodotti tensoriali:

|ud>=(1,0)^T⊗(0,1)^T=(0,1,0,0)^T   ,   |du>=(0,1)^T⊗(1,0)^T=(0,0,1,0)^T

 <ud|=(1,0)⊗(0,1)=(0,1,0,0)   ,   <du|=(0,1)⊗(1,0)=(0,0,1,0).

Quindi, svolgendo i prodotti sopra definiti, si ottiene la matrice [4x4]:

Infatti come visto nel precedente post, gli elementi di ρ sono i prodotti delle ampiezze di probabilità per i coniugati; in particolare nel nostro caso risulta:

ψ_udψ_ud=ψ_udψ_du=ψ_duψ_ud=ψ_duψ_du=1/2

mentre gli altri elementi di ρ sono tutti nulli (si ricordi che in |Ψ>=ψ_ud|ud>+ψ_du|du> risulta ψ_ud=ψ_du=(1/2)^1/2).

A questo punto possiamo verificare facilmente la relazione ρ=ρ² (basta moltiplicare la matrice ρ per se stessa); ciò significa che siamo in presenza di uno stato puro quindi la conoscenza del sistema combinato è completa***.
Nota: il sistema è stato preparato in uno stato definito di spin perciò è puro, inoltre ciò implica una forte correlazione tra gli spin delle due particelle (poiché se un elettrone è misurato up l'altro è down e viceversa).

Tuttavia la matrice densità ρ riguarda tutto il sistema combinato mentre noi vorremmo descrivere lo stato di ogni singolo elettrone (chiamiamoli A e B) e questo può essere fatto grazie all'introduzione della matrice densità ridotta (che viene giustificata fisicamente dalla teoria della decoerenza).
Nota: come vedremo nella matrice ridotta i termini fuori diagonale, responsabile delle interferenze quantistiche, si annullano.

Introduciamo quindi la matrice densità ridotta che permette di studiare uno dei due sottosistemi (supponiamo A) ed è così definita:

ρ_A=∑<i|ρ|i>=Tr_Bρ

rispetto ad una base di vettori |i> del sistema B.

Nota: Tr_B è l'operatore traccia parziale sulla base di B; in modo equivalente si ha ρ_B=Tr_Aρ. Inoltre se |Ψ> è uno stato prodotto risulta ρ=ρ_A⊗ρ_B.

Perciò nel caso considerato possiamo calcolare la matrice ridotta dello stato di spin dell'elettrone A (oppure di quello B) e risulta:

ρ_A=1/2(|u><u|+|d><d|)=1/2(1,0)^T(1,0)+1/2(0,1)^T(0,1)=(1/2)I₂

dove con I₂ abbiamo indicato la matrice identità [2x2]; da ciò si deduce subito che ρ_A≠ρ_A² cioè siamo in presenza di uno stato composto(!)
Nota: I è una matrice diagonale con tutti gli elementi pari a 1 perciò I²=I.

Ciò significa che gli stati dell'elettrone A (oppure di quello B) non sono in sovrapposizione quantistica, l'incertezza sullo spin è in realtà dovuta alla non completa conoscenza dello stato del sistema totale (poiché stiamo considerando solo il sottositema A) e la probabilità statistica che lo spin di un elettrone sia up oppure down è pari a 1/2.
Nota: a differenza del caso classico però, nemmeno in linea di principio si può stabilire lo stato del sotto sistema A (o B) prima della misura.

È interessante osservare che il famoso Paradosso del gatto di Schrödinger può essere trattato come lo stato entangled che abbiamo ora considerato; ciò significa che anche in questo caso non si ha sovrapposizione di due stati distinti (gatto vivo e gatto morto) poiché il sottosistema "gatto" si trova in uno stato misto di tipo statistico e non è in sovrapposizione quantistica.

(*) Nel precedente post abbiamo definito, per uno stato puro, la matrice densità ρ=|Ψ><Ψ| per la quale risulta ρ²=|Ψ><Ψ|Ψ><Ψ|=ρ; invece per uno stato misto si pone ρ=∑p_i|Ψ_i><Ψ_i| dove p_i è la probabilità che il sistema si trovi nello stato i-esimo e in questo caso risulta ρ≠ρ².
(**) Si osservi che i vettori di base scelti soddisfano correttamente le condizioni di ortonormalità: <u|u>=<d|d>=1 e <u|d>=<d|u>=0.
(***) Le teorie a variabili nascoste affermano invece che la conoscenza quantistica del sistema composto non è completa proprio perché lo stato dei singoli sottositemi non è definito con certezza.

[Si noti che lo stato |Ψ>=(1/2)^1/2|ud>+(1/2)^1/2|du> è lo stato di tripletto con spin S=1 nonostante le particelle abbiano spin opposto (come accade nello stato di singoletto |Ψ>=(1/2)^1/2|ud>-(1/2)^1/2|du> dove però S=0)]

Stati puri, miscele e sovrapposizioni!

Come è noto, definito uno stato |Ψ> di un qualsiasi sistema quantistico la sua evoluzione temporale, una volta fissato lo stato iniziale, è descritta dalla equazione di Schrödinger (vedi Wikipedia), scritta nella notazione di Dirac:

i(h/2π)∂|Ψ>/∂t=H|Ψ>

dove il valore medio dell'operatore hamiltoniano <H> rappresenta il valore di aspettazione dell'energia del sistema.

Si osservi che qui ci limitiamo a trattare il caso di uno spazio finito-dimensionale (cioè definito da n vettori di base |i>) per il quale si ha:

|Ψ>=∑ψ_i|i>

dove ψ_i sono le ampiezze di probabilità relative ai vettori di base |i>.
Nota: per chiarimenti sul vettore di stato di un sistema quantistico vedi i post "I numeri Compessi e la M.Q." e "Le grandezze Osservabili!".

È però possibile dare una descrizione alternativa ma equivalente a quella di Schrödinger definendo il seguente Operatore di densità*:

ρ=|Ψ><Ψ|

dove ρ è rappresentato da una matrice quadrata che si ottiene moltiplicando il vettore colonna |Ψ> per il suo duale vettore riga <Ψ|.
Nota: l'operatore |Ψ><Ψ| è un proiettore poiché applicato ad uno stato |Φ> si ha |Ψ><Ψ|Φ>=k|Ψ> con k=<Ψ|Φ> (cioè proietta |Φ> lungo |Ψ>).

Poiché il vettore di stato è definito rispetto ad una base ortonormale di vettori |i> si ha che gli elementi di matrice ρ_ij sono dati da**

ρ_ij=<i|ρ|j>

da cui segue subito (sostituendo ρ=|Ψ><Ψ| ed essendo ψ_i=<i|Ψ>):

ρ_ij=<i|Ψ><Ψ|j>=ψ_iψ_j

cioè gli elementi di ρ sono i prodotti delle ampiezze di probabilità (associate ai vettori di base dello stato considerato) per i coniugati.
Nota: solo quando i=j il prodotto ψ_iψ_i=|ψ_i|² rappresenta la probabilità che il sistema venga misurato nello stato i-esimo.

Inoltre se abbiamo a che fare con un sistema il cui stato non è ben definito, ma è dato da un ensemble statistico di stati possibili {Ψ_i} si può porre:

ρ=∑p_i|Ψ_i><Ψ_i|

dove p_i è la probabilità statistica che il sistema si trovi nello stato i-esimo: in pratica è la media pesata su tutti gli stati possibili del sistema con ∑p_i=1.
Nota: la probabilità p_i è di tipo statistico poiché è dovuta alla non esatta conoscenza dello stato del sistema (non è una sovrapposizione quantistica).

Si parla quindi di stato puro quando le p_i sono tutte nulle tranne una (e pari a 1), mentre negli altri casi avremo uno stato misto poiché si determina una media pesata su tutti gli stati |Ψ_i> in cui si potrebbe trovare il sistema.

Facciamo subito un esempio di stato puro e consideriamo lo stato di spin di un singolo elettrone (vedi il post "I numeri Complessi e la M.Q.") che può essere descritto in generale (ad es. rispetto alle basi |u> e |d> lungo Z):

|Ψ>=ψ_u|u>+ψ_d|d>.

Ciò significa che, dato uno spin preparato in uno stato qualunque |Ψ> e un apparato di misura orientato lungo l'asse Z, i prodotti ψ_uψ_u e ψ_dψ_d sono le rispettive probabilità che lo spin si trovi nello stato |u> oppure |d>.
Nota: secondo i postulati quantistici, prima della misura lungo l'asse Z gli stati |u> e |d> sono in sovrapposizione quantistica.

Si ricordi infatti che il principale postulato della meccanica quantistica stabilisce che il prodotto ψ_iψ_i (cioè |ψ_i|²) dà la probabilità che il sistema si trovi nello stato i-esimo (ψ_i è il complesso coniugato di ψ_i).

Perciò, come mostrato sopra, le componenti della matrice ρ sono:

ρ_uu=ψ_uψ_u , ρ_ud=ψ_uψ_d , ρ_du=ψ_dψ_u , ρ_dd=ψ_dψ_d

e in particolare risulta: Trρ=∑ρ_ii=∑ψ_iψ_i=1.
Nota: Tr è l'operatore traccia cioè la somma degli elementi ψ_iψ_i della diagonale di ρ quindi si ha Trρ=1 (essendo normalizzata a 1).

Ad esempio se prepariamo (misuriamo) lo stato di spin dell'elettrone nella direzione dell'asse X positivo (detto stato right) allora possiamo scrivere (vedi il post "I numeri Complessi e la M.Q."):

|Ψ_r>=(1/2)^1/2|u>+(1/2)^1/2|d>

quindi risulta per tutti gli elementi della matrice ρ [2x2]:

ρ_uu=ρ_ud=ρ_du=ρ_dd=1/2

dove correttamente si ha Trρ=ρ_uu+ρ_dd=1/2+1/2=1 (cioè la condizione di normalizzazione ψ_uψ_u+ψ_dψ_d=1 è soddisfatta).
Nota: perciò la probabilità che lo spin, misurato lungo l'asse Z, sia up oppure down è pari a 1/2.

Ma ciò che risulta di grande interesse è che per uno stato puro, come quello appena trattato, vale la condizione*** ρ=ρ² e ciò ci permette di distinguere, come vedremo nel prossimo post, uno stato puro da uno stato misto!
Nota: moltiplicando per se stessa la matrice ρ composta da elementi uguali a 1/2 si ottiene di nuovo la matrice ρ.

(*) L'evoluzione di ρ nel tempo è descritta dalla equazione di Von Neumann:
i(h/2π)∂ρ/∂t=[H,ρ] dove [H,ρ]=Hρ-ρH è il commutatore di H e ρ. 
(**) Posto |Ψ>=∑ψ_j|j> si ha <i|Ψ>=∑ψ_j<j|i>=ψ_i poiché solo quando j=i si ha <i|i>=1; inoltre dato un operatore A risulta A_ij=<i|A|j> infatti possiamo scrivere <i|A|Ψ>=∑ψ_j<i|A|j>=∑ψ_jA_ij $<=$$>$ <i|A|Ψ>=<i|∑Φ_j|j>=Φ_i (per j=i), che rappresenta l'equazione A|Ψ>=|Φ> in forma matriciale.
(***) Per uno stato puro si ha ρ=|Ψ><Ψ| quindi ρ²=|Ψ><Ψ|Ψ><Ψ|=ρ essendo per la condizione di normalizzazione <Ψ|Ψ>=1.