View source for Geometry of interaction

The ''geometry of interaction'', GoI in short, was defined in the early nineties by Girard as an interpretation of linear logic into operators algebra: formulae were interpreted by Hilbert spaces and proofs by partial isometries.

This was a striking novelty as it was the first time that a mathematical model of logic (lambda-calculus) didn't interpret a proof of <math>A\limp B</math> as a morphism ''from'' <math>A</math> ''to'' <math>B</math><ref>to be precise one should say from ''the space interpreting'' <math>A</math> to the space interpreting'' <math>B</math></ref>, and proof composition (cut rule) as the composition of morphisms. Rather the proof was interpreted as an operator acting ''on'' <math>A\limp B</math>, that is a morphism from <math>A\limp B</math> to <math>A\limp B</math>. For proof composition the problem was then, given an operator on <math>A\limp B</math> and another one on <math>B\limp C</math> to construct a new operator on <math>A\limp C</math>. This problem was solved by the ''execution formula'' that bares some formal analogies with Kleene's formula for recursive functions. For this reason GoI was claimed to be an ''operational semantics'', as opposed to traditionnal [[Semantics|denotational semantics]].

The first instance of the GoI was restricted to the <math>MELL</math> fragment of linear logic (Multiplicative and Exponential fragment) which is enough to encode lambda-calculus. Since then Girard proposed several improvements: firstly the extension to the additive connectives known as ''Geometry of Interaction 3'' and more recently a complete reformulation using Von Neumann algebras that allows to deal with some aspects of [[Light linear logics|implicit complexity]]

The GoI has been a source of inspiration for various authors. Danos and Regnier have reformulated the original model exhibiting its combinatorial nature using a theory of reduction of paths in proof-nets and showing the link with abstract machines; in particular the execution formula appears as the composition of two automata that interact one with the other through their common interface. Also the execution formula has rapidly been understood as expressing the composition of strategies in game semantics. It has been used in the theory of sharing reduction for lambda-calculus in the Abadi-Gonthier-Lévy reformulation and simplification of Lamping's representation of sharing. Finally the original GoI for the <math>MELL</math> fragment has been reformulated in the framework of traced monoidal categories following an idea originally proposed by Joyal.

= The Geometry of Interaction as operators =

The original construction of GoI by Girard follows a general pattern already mentionned in [[coherent semantics]] under the name ''symmetric reducibility''. First set a general space called the ''proof space'' because this is where the interpretations of proofs will live. In the case of GoI, the proof space is the space of bounded operators on <math>\ell^2</math>. Note that the proof space generally contains much more objects than interpretations of proofs; in the GoI case we will see that interpretations of proofs happen to be some very peculiar kind of partial isometries.

Second define a duality on this space that will be denoted as <math>u\perp v</math>. For the GoI, two dualities have proved to work, the first one being nilpotency: two operators <math>u</math> and <math>v</math> are dual if <math>uv</math> is nilpotent, that is, if there is a nonegative integer <math>n</math> such that <math>(uv)^n = 0</math>. We will denote by <math>\bot</math> the set of nilpotent operators so that the duality reads:
: <math>u\perp v</math> iff <math>uv\in\bot</math>.

This duality applies to operators and shouldn't be confused with orthogonality
of vectors. . To enforce this we will reserve the notation <math>\perp</math>
exclusively for the duality of operators and never use it for othogonality of vectors.

Last define a ''type'' as a subset <math>T</math> of the proof space that is equal to its bidual: <math>T = T\biorth</math>. This means that <math>u\in T</math> iff for all operator <math>v</math>, if <math>v\in T\orth</math>, that is if <math>u'v\in\bot</math> for all <math>u'\in T</math>, then <math>uv\in\bot</math>.

It remains now to interpret logical operations, that is associate a type to each formula, an object to each proof and show the ''adequacy lemma'': if <math>u</math> is the interpretation of a proof of the formula <math>A</math> then <math>u</math> belongs to the type associated to <math>A</math>.

== Preliminaries ==

We begin by a brief tour of the operations in Hilbert spaces that will be used in the sequel. In this article <math>H</math> will stand for the Hilbert space <math>\ell^2(\mathbb{N})</math> of sequences <math>(x_n)_{n\in\mathbb{N}}</math> of complex numbers such that the series <math>\sum_{n\in\mathbb{N}}|x_n|^2</math> converges. If <math>x = (x_n)_{n\in\mathbb{N}}</math> and <math>y = (y_n)_{n\in\mathbb{N}}</math> are two vectors of <math>H</math> we denote by <math>\langle x,y\rangle</math> their scalar product:
: <math>\langle x, y\rangle = \sum_{n\in\mathbb{N}} x_n\bar y_n</math>.

Two vectors of <math>H</math> are ''othogonal'' if their scalar product is nul. This notion is not to be confused with the orthogonality of operators defined above. The ''norm'' of a vector is the square root of the scalar product with itself:
: <math>\|x\| = \sqrt{\langle x, x\rangle}</math>.

Let us denote by <math>(e_k)_{k\in\mathbb{N}}</math> the canonical hilbertian basis of <math>H</math>: <math>e_k = (\delta_{kn})_{n\in\mathbb{N}}</math> where <math>\delta_{kn}</math> is the Kroenecker symbol. Thus if <math>x=(x_n)_{n\in\mathbb{N}}</math> is a sequence in <math>H</math> we have:
: <math> x = \sum_{n\in\mathbb{N}} x_ne_n</math>.

In this article we call ''operator'' on <math>H</math> a ''continuous'' linear map from <math>H</math> to <math>H</math>. Continuity is equivalent to the fact that operators are ''bounded'', which means that one may define the ''norm'' of an operator <math>u</math> as the sup on the unit ball of the norms of its values:
: <math>\|u\| = \sup_{\{x\in H,\, \|x\| = 1\}}\|u(x)\|</math>.

The set of (bounded) operators is denoted <math>\mathcal{B}(H)</math>. This is our proof space.

The ''range'' or ''codomain'' of the operator <math>u</math> is the set of images of vectors; the ''kernel'' of <math>u</math> is the set of vectors that are anihilated by <math>u</math>; the ''domain'' of <math>u</math> is the set of vectors orthogonal to the kernel:

: <math>\mathrm{Codom}(u) = \{u(x),\, x\in H\}</math>;
: <math>\mathrm{Ker}(u) = \{x\in H,\, u(x) = 0\}</math>;
: <math>\mathrm{Dom}(u) = \{x\in H,\, \forall y\in\mathrm{Ker}(u), \langle x, y\rangle = 0\}</math>.

These three sets are closed subspaces of <math>H</math>.

The ''adjoint'' of an operator <math>u</math> is the operator <math>u^*</math> defined by <math>\langle u(x), y\rangle = \langle x, u^*(y)\rangle</math> for any <math>x,y\in H</math>.

A ''projector'' is an idempotent operator of norm <math>0</math> (the projector
on the null subspace) or <math>1</math>, that is an operator <math>p</math>
such that <math>p^2 = p</math> and <math>\|p\| = 0</math> or <math>1</math>. A projector is auto-adjoint and its domain is equal to its codomain.

A ''partial isometry'' is an operator <math>u</math> satisfying <math>uu^* u = u</math>; as a consequence <math>uu^*</math> is a projector the range of which is the range of <math>u</math>. Similarly <math>u^* u</math> is also a projector the range of which is the domain of <math>u</math>. The restriction of <math>u</math> to its domain is an isometry. Projectors are particular examples of partial isometries.

If <math>u</math> is a partial isometry then <math>u^*</math> is also a partial isometry the domain of which is the codomain of <math>u</math> and the codomain of which is the domain of <math>u</math>.

If the domain of <math>u</math> is <math>H</math> that is if <math>u^* u = 1</math> we say that <math>u</math> has ''full domain'', and similarly for codomain. If <math>u</math> and <math>v</math> are two partial isometries, the equation <math>uu^* + vv^* = 1</math> means that the codomains of <math>u</math> and <math>v</math> are orthogonal and that their direct sum is <math>H</math>.

=== Partial permutations and partial isometries ===

It turns out that most of the operators needed to interpret logical operations are generated by ''partial permutations'' on the basis, which in particular entails that they are partial isometries.

More precisely a partial permutation <math>\varphi</math> on <math>\mathbb{N}</math> is a function defined on a subset <math>D_\varphi</math> of <math>\mathbb{N}</math> which is one-to-one onto a subset <math>C_\varphi</math> of <math>\mathbb{N}</math>. <math>D_\varphi</math> is called the ''domain'' of <math>\varphi</math> and <math>C_\varphi</math> its ''codomain''. Partial permutations may be composed: if <math>\psi</math> is another partial permutation on <math>\mathbb{N}</math> then <math>\varphi\circ\psi</math> is defined by:

: <math>n\in D_{\varphi\circ\psi}</math> iff <math>n\in D_\psi</math> and <math>\psi(n)\in D_\varphi</math>;
: if <math>n\in D_{\varphi\circ\psi}</math> then <math>\varphi\circ\psi(n) = \varphi(\psi(n))</math>;
: the codomain of <math>\varphi\circ\psi</math> is the image of the domain.

Partial permutations are well known to form a structure of ''inverse monoid'' that we detail now.

A ''partial identitie'' is a partial permutation <math>1_D</math> whose domain and codomain are both equal to a subset <math>D</math> on which <math>1_D</math> is the identity function. Partial identities are idempotent for composition.

Among partial identities one finds the identity on the empty subset, that is the empty map, that we will denote as <math>0</math> and the identity on <math>\mathbb{N}</math> that we will denote by <math>1</math>. This latter permutation is the neutral for composition.

If <math>\varphi</math> is a partial permutation there is an inverse partial permutation <math>\varphi^{-1}</math> whose domain is <math>D_{\varphi^{-1}} = C_{\varphi}</math> and who satisfies:

: <math>\varphi^{-1}\circ\varphi = 1_{D_\varphi}</math>
: <math>\varphi\circ\varphi^{-1} = 1_{C_\varphi}</math>

Given a partial permutation <math>\varphi</math> one defines a partial isometry <math>u_\varphi</math> by:
: <math>u_\varphi(e_n) = 
   \begin{cases}
     e_{\varphi(n)} & \text{ if }n\in D_\varphi,\\
     0              & \text{ otherwise.}
   \end{cases}
  </math>
In other terms if <math>x=(x_n)_{n\in\mathbb{N}}</math> is a sequence in <math>\ell^2</math> then <math>u_\varphi(x)</math> is the sequence <math>(y_n)_{n\in\mathbb{N}}</math> defined by:
: <math>y_n = x_{\varphi^{-1}(n)}</math> if <math>n\in C_\varphi</math>, <math>0</math> otherwise.

We will (not so abusively) write <math>e_{\varphi(n)} = 0</math> when <math>\varphi(n)</math> is undefined.

The domain of <math>u_\varphi</math> is the subspace spaned by the family <math>(e_n)_{n\in D_\varphi}</math> and the codomain of <math>u_\varphi</math> is the subspace spaned by <math>(e_n)_{n\in C_\varphi}</math>. As a particular case if <math>\varphi</math> is <math>1_D</math> the partial identity on <math>D</math> then <math>u_\varphi</math> is the projector on the subspace spaned by <math>(e_n)_{n\in D}</math>.

If <math>\psi</math> is another partial permutation then we have:
: <math>u_\varphi u_\psi = u_{\varphi\circ\psi}</math>.

If <math>\varphi</math> is a partial permutation then the adjoint of <math>u_\varphi</math> is:
: <math>u_\varphi^* = u_{\varphi^{-1}}</math>.

In particular the projector on the domain of <math>u_{\varphi}</math> is given by:
: <math>u^*_\varphi u_\varphi = u_{1_{D_\varphi}}</math>.

and similarly the projector on the codomain of <math>u_\varphi</math> is:
: <math>u_\varphi u_\varphi^* = u_{1_{C_\varphi}}</math>.

{{Proposition|
Let <math>u_\varphi</math> and <math>u_\psi</math> be two partial isometries generated by partial permutations. Then we have:
: <math>u_\varphi + u_\psi = 0</math> iff <math>u_\varphi = u_\psi = 0</math>,
that is iff <math>\varphi</math> and <math>\psi</math> are the nowhere defined partial permutation.
}} Indeed suppose <math>u_\varphi + u_\psi = 0</math> then for any <math>n</math> we have <math>u_\varphi(e_n) + u_\psi(e_n) = e_{\varphi(n)} + e_{\psi(n)} = 0</math> which is possible only if <math>\varphi(n)</math> and <math>\psi(n)</math> are undefined.

=== From operators to matrices: internalization/externalization ===

It will be convenient to view operators on <math>H</math> as acting on <math>H\oplus H</math>, and conversely. For this purpose we define an isomorphism <math>H\oplus H \cong H</math> by <math>x\oplus y\rightsquigarrow p(x)+q(y)</math> where <math>p:H\mapsto H</math> and <math>q:H\mapsto H</math> are partial isometries given by:

: <math>p(e_n) = e_{2n}</math>,
: <math>q(e_n) = e_{2n+1}</math>.

From the definition <math>p</math> and <math>q</math> have full domain, that is satisfy <math>p^* p = q^* q = 1</math>. On the other hand their codomains are orthogonal, thus we have <math>p^* q = q^* p = 0</math>. Note that we also have <math>pp^* + qq^* = 1</math>.

The choice of <math>p</math> and <math>q</math> is actually arbitrary, any two partial isometries with full domain and orthogonal codomains would do the job.

Let <math>U</math> be an operator on <math>H\oplus H</math>. We can write <math>U</math> as a matrix:
: <math>U = \begin{pmatrix}
  u_{11} & u_{12}\\
  u_{21} & u_{22}
  \end{pmatrix}</math>
where each <math>u_{ij}</math> operates on <math>H</math>.

Now through the isomorphism <math>H\oplus H\cong H</math> we may transform <math>U</math> into the operator <math>u</math> on <math>H</math> defined by:

: <math>u = pu_{11}p^* + pu_{12}q^* + qu_{21}p^* + qu_{22}q^*</math>.

We call <math>u</math> the ''internalization'' of <math>U</math>. Internalization is compatible with composition (functorial so to speak): if <math>V</math> is another operator on <math>H\oplus</math> then the internalization of the matrix product <math>UV</math> is the product <math>uv</math>.

Conversely given an operator <math>u</math> on <math>H</math> we may externalize it obtaining an operator <math>U</math> on <math>H\oplus H</math>:
: <math>u_{11} = p^*up</math>;
: <math>u_{12} = p^*uq</math>;
: <math>u_{21} = q^*up</math>;
: <math>u_{22} = q^*uq</math>.

== Interpreting the multiplicative connectives ==

Recall that when <math>u</math> and <math>v</math> are operators we denote by <math>u\perp v</math> the fact that <math>uv</math> is nilpotent, and that <math>\bot</math> denotes the set of nilpotent operators so that <math>u\perp v</math> iff <math>uv\in\bot</math>.

If <math>X</math> is  set of operators also recall that <math>X\orth</math> denotes the set of dual operators:
: <math>X\orth = \{v\in \mathcal{B}(H) \text{ such that }\forall u\in X, uv \in\bot\}</math>.

There are a few properties of this duality that we will use without mention in the sequel; let <math>X</math> and <math>Y</math> be sets of operators:
: <math>X\subset X\biorth</math>;
: <math>X\orth = X\triorth</math>.
: if <math>X\subset Y</math> then <math>Y\orth\subset X\orth</math>;

In particular <math>X\orth</math> is always a type (equal to its biorthogonal). We say that <math>X</math> ''generates'' the type <math>X\biorth</math>.

=== The tensor and the linear application ===

Given two types <math>A</math> and <math>B</math> two types, we define their tensor by:

: <math>A\tens B = \{pup^* + qvq^*, u\in A, v\in B\}\biorth</math>

Note the closure by biorthogonal to make sure that we obtain a type. From what precedes we see that <math>A\tens B</math> is generated by the internalizations of operators on <math>H\oplus H</math> of the form:
: <math>\begin{pmatrix}
   u & 0\\
   0 & v
  \end{pmatrix}</math>

This is an abuse of notations as this operation is more like a direct sum than a tensor. We will stick to this notation though because it defines the interpretation of the tensor connective of linear logic.

The linear implication is derived from the tensor by duality: given two types <math>A</math> and <math>B</math> the type <math>A\limp B</math> is defined by:
: <math>A\limp B = (A\tens B\orth)\orth</math>.

Unfolding this definition we see that we have:
: <math>A\limp B = \{u\in\mathcal{B}(H)\text{ such that } \forall v\in A, \forall w\in B\orth,\, u.(pvp^* + qwq^*) \in\bot\}</math>.

=== The idendity ===

The interpretation of the identity is an example of the internalization/externalization procedure. Given a type <math>A</math> we are to find an operator <math>\iota</math> in type <math>A\limp A</math>, thus satisfying:
: <math>\forall u\in A, v\in A\orth,\, \iota(pup^* + qvq^*)\in\bot</math>.

An easy solution is to take <math>\iota = pq^* + qp^*</math>. In this way we get <math>\iota(pup^* + qvq^*) = qup^* + pvq^*</math>. Therefore <math>(\iota(pup^* + qvq^*))^2 = quvq^* + pvup^*</math>, from which one deduces that this operator is nilpotent iff <math>uv</math> is nilpotent. It is the case since <math>u</math> is in <math>A</math> and <math>v</math> in <math>A\orth</math>.

It is interesting to note that the <math>\iota</math> thus defined is actually the internalization of the operator on <math>H\oplus H</math> given by the matrix:
: <math>\begin{pmatrix}0 & 1\\1 & 0\end{pmatrix}</math>.

We will see once the composition is defined that the <math>\iota</math> operator is the interpretation of the identity proof, as expected.

=== Interpreting the cut rule: the execution formula ===

Let <math>A</math> and <math>B</math> be two types and <math>u</math> an operator in <math>A\limp B</math>. By definition this means that given <math>v</math> in <math>A</math> and <math>w</math> in <math>B\orth</math> the operator <math>u.(pvp^* + qwq^*)</math> is nilpotent.

Let us define <math>u_{11}</math> to <math>u_{22}</math> by externalization as above. If we compute <math>(u.(pvp^* + qwq^*))^n</math> we see that this is a finite sum of operators of the form:
: <math>q(u_{22}w)^{k_0}u_{21}v(u_{11}v)^{k_1}u_{12}w\dots u_{12}w(u_{22}w)^{k_{p+1}}q^*</math>,
: <math>p(u_{11}v)^{k_1}u_{12}w\dots u_{12}w(u_{22}w)^{k_{p+1}}q^*</math>,
: <math>q(u_{22}w)^{k_0}u_{21}v(u_{11}v)^{k_1}u_{12}w\dots (u_{11}v)^{k_p}p^*</math> or
: <math>p(u_{11}v)^{k_1}u_{12}w\dots (u_{11}v)^{k_p}p^*</math>
where each of these monimials has exactly <math>n</math> factors of the form <math>u_{ij}(v\text{ or }w)</math>.

Let us suppose that <math>u</math>, <math>v</math> and <math>w</math> are partial isometries generated by partial permutations of the basis. Then this is also the case of all these monomials. For <math>n</math> big enough, we know that <math>(u.(pvp^* + qwq^*))^n = 0</math>, thus it is also the case that all these monomials are null for <math>n</math> big enough. Now this is equivalent to the conjonction of two facts:
: <math>u_{11}v</math> is nilpotent and
: <math>\bigl(u_{22} + u_{21}v\sum(u_{11}v)^ku_{12}\bigr).w</math> is nilpotent.

=== Interpreting the tensor rule ===

Let now <math>A, A', B</math> and <math>B'</math> be types and consider two operators <math>u</math> and <math>u'</math> respectively in <math>A\limp B</math> and <math>A\limp B'</math>. We define an operator denoted <math>u\tens u'</math> by:
: <math>\begin{align}
    u\tens u' &= ppp^*upp^*p^* + qpq^*upp^*p^* + ppp^*uqp^*q^* + qpq^*uqp^*q^*\\
              &+ pqp^*vpq^*p^* + qqq^*vpq^*p^* + pqp^*vqq^*q^* + qqq^*vqq^*q^*
  \end{align}</math>

To understand this formula it is convenient to think <math>u</math> and <math>u'</math> as the internalizations of the matrices:
: <math>U  = \begin{pmatrix}u_{11}   & u_{12}\\
                            u_{21}   & u_{22}
             \end{pmatrix}
  </math> and <math>U' = \begin{pmatrix}u'_{11} & u'_{12}\\
                            u'_{21} & u'_{22}
                         \end{pmatrix}</math>
where the <math>u_{ij}</math>'s and the <math>u'_{ij}</math>'s are defined by the formula above, eg <math>u_{11} = p^*up</math>.

Then <math>u\tens u'</math> is actually the internalization of the matrix <math>U\tens U'</math> given by:

: <math>
    U\tens U' =
      \begin{pmatrix}
        u_{11} & 0       & u_{12}  & 0       \\
        0      & u'_{11} & 0       & u'_{12} \\
        u_{21} & 0       & u_{22}  & 0       \\
        0      & u'_{21} & 0       & u'_{22} \\
      \end{pmatrix}
  </math>

It remains to show that, given that <math>u</math>and <math>u'</math> are in types <math>A\limp B</math> and <math>A'\limp B'</math>, then <math>u\tens u'</math> is in <math>A\tens A'\limp B\tens B'</math>. We postpone this for after the definition of the execution.

= The Geometry of Interaction as an abstract machine =