Physics:Sigma model

Short description: Field theory of a point particle confined to move on a fixed manifold

In physics, a sigma model is a field theory that describes the field as a point particle confined to move on a fixed manifold. This manifold can be taken to be any Riemannian manifold, although it is most commonly taken to be either a Lie group or a symmetric space. The model may or may not be quantized. An example of the non-quantized version is the Skyrme model; it cannot be quantized due to non-linearities of power greater than 4. In general, sigma models admit (classical) topological soliton solutions, for example, the Skyrmion for the Skyrme model. When the sigma field is coupled to a gauge field, the resulting model is described by Ginzburg–Landau theory. This article is primarily devoted to the classical field theory of the sigma model; the corresponding quantized theory is presented in the article titled "non-linear sigma model".

Overview

The name has roots in particle physics, where a sigma model describes the interactions of pions. Unfortunately, the "sigma meson" is not described by the sigma-model, but only a component of it.

The sigma model was introduced by (Gell-Mann Lévy); the name σ-model comes from a field in their model corresponding to a spinless meson called σ, a scalar meson introduced earlier by Julian Schwinger. The model served as the dominant prototype of spontaneous symmetry breaking of O(4) down to O(3): the three axial generators broken are the simplest manifestation of chiral symmetry breaking, the surviving unbroken O(3) representing isospin.

In conventional particle physics settings, the field is generally taken to be SU(N), or the vector subspace of quotient [math]\displaystyle{ (SU(N)_L\times SU(N)_R)/SU(N) }[/math] of the product of left and right chiral fields. In condensed matter theories, the field is taken to be O(N). For the rotation group O(3), the sigma model describes the isotropic ferromagnet; more generally, the O(N) model shows up in the quantum Hall effect, superfluid Helium-3 and spin chains.

In supergravity models, the field is taken to be a symmetric space. Since symmetric spaces are defined in terms of their involution, their tangent space naturally splits into even and odd parity subspaces. This splitting helps propel the dimensional reduction of Kaluza–Klein theories.

In its most basic form, the sigma model can be taken as being purely the kinetic energy of a point particle; as a field, this is just the Dirichlet energy in Euclidean space.

In two spatial dimensions, the O(3) model is completely integrable.

Definition

The Lagrangian density of the sigma model can be written in a variety of different ways, each suitable to a particular type of application. The simplest, most generic definition writes the Lagrangian as the metric trace of the pullback of the metric tensor on a Riemannian manifold. For [math]\displaystyle{ \phi:M\to\Phi }[/math] a field over a spacetime [math]\displaystyle{ M }[/math], this may be written as

[math]\displaystyle{ \mathcal{L} = \frac{1}{2}\sum_{i=1}^n \sum_{j=1}^n g_{ij}(\phi) \; \partial^\mu\phi_i \partial_\mu\phi_j }[/math]

where the [math]\displaystyle{ g_{ij}(\phi) }[/math] is the metric tensor on the field space [math]\displaystyle{ \phi\in\Phi }[/math], and the [math]\displaystyle{ \partial_\mu }[/math] are the derivatives on the underlying spacetime manifold.

This expression can be unpacked a bit. The field space [math]\displaystyle{ \Phi }[/math] can be chosen to be any Riemannian manifold. Historically, this is the "sigma" of the sigma model; the historically-appropriate symbol [math]\displaystyle{ \sigma }[/math] is avoided here to prevent clashes with many other common usages of [math]\displaystyle{ \sigma }[/math] in geometry. Riemannian manifolds always come with a metric tensor [math]\displaystyle{ g }[/math]. Given an atlas of charts on [math]\displaystyle{ \Phi }[/math], the field space can always be locally trivialized, in that given [math]\displaystyle{ U\subset\Phi }[/math] in the atlas, one may write a map [math]\displaystyle{ U\to\mathbb{R}^n }[/math] giving explicit local coordinates [math]\displaystyle{ \phi=(\phi^1, \cdots,\phi^n) }[/math] on that patch. The metric tensor on that patch is a matrix having components [math]\displaystyle{ g_{ij}(\phi). }[/math]

The base manifold [math]\displaystyle{ M }[/math] must be a differentiable manifold; by convention, it is either Minkowski space in particle physics applications, flat two-dimensional Euclidean space for condensed matter applications, or a Riemann surface, the worldsheet in string theory. The [math]\displaystyle{ \partial_\mu \phi = \partial\phi/\partial x^\mu }[/math] is just the plain-old covariant derivative on the base spacetime manifold [math]\displaystyle{ M. }[/math] When [math]\displaystyle{ M }[/math] is flat, [math]\displaystyle{ \partial_\mu \phi = \nabla\phi }[/math] is just the ordinary gradient of a scalar function (as [math]\displaystyle{ \phi }[/math] is a scalar field, from the point of view of [math]\displaystyle{ M }[/math] itself.) In more precise language, [math]\displaystyle{ \partial_\mu\phi }[/math] is a section of the jet bundle of [math]\displaystyle{ M\times\Phi }[/math].

Example: O(N) non-linear sigma model

Taking [math]\displaystyle{ g_{ij}=\delta_{ij} }[/math] the Kronecker delta, i.e. the scalar dot product in Euclidean space, one gets the [math]\displaystyle{ O(n) }[/math] non-linear sigma model. That is, write [math]\displaystyle{ \phi =\hat{u} }[/math] to be the unit vector in [math]\displaystyle{ \mathbb{R}^{n} }[/math], so that [math]\displaystyle{ \hat{u}\cdot\hat{u}=1 }[/math], with [math]\displaystyle{ \cdot }[/math] the ordinary Euclidean dot product. Then [math]\displaystyle{ \hat{u}\in S^{n-1} }[/math] the [math]\displaystyle{ (n-1) }[/math]-sphere, the isometries of which are the rotation group [math]\displaystyle{ O(n) }[/math]. The Lagrangian can then be written as

[math]\displaystyle{ \mathcal{L} = \frac{1}{2} \nabla_\mu\hat{u} \cdot \nabla_\mu\hat{u} }[/math]

For [math]\displaystyle{ n=3 }[/math], this is the continuum limit of the isotropic ferromagnet on a lattice, i.e. of the classical Heisenberg model. For [math]\displaystyle{ n=2 }[/math], this is the continuum limit of the classical XY model. See also the n-vector model and the Potts model for reviews of the lattice model equivalents. The continuum limit is taken by writing

[math]\displaystyle{ \delta_h[\hat{u}](i,j)=\frac{\hat{u}_i - \hat{u}_j}{h} }[/math]

as the finite difference on neighboring lattice locations [math]\displaystyle{ i,j. }[/math] Then [math]\displaystyle{ \delta_h[\hat{u}]\to\partial_\mu\hat{u} }[/math] in the limit [math]\displaystyle{ h\to 0 }[/math], and [math]\displaystyle{ \hat{u}_i\cdot \hat{u}_j\to \partial_\mu\hat{u} \cdot \partial_\mu\hat{u} }[/math] after dropping the constant terms [math]\displaystyle{ \hat{u}_i\cdot\hat{u}_i=1 }[/math] (the "bulk magnetization").

In geometric notation

The sigma model can also be written in a more fully geometric notation, as a fiber bundle with fibers [math]\displaystyle{ \Phi }[/math] over a differentiable manifold [math]\displaystyle{ M }[/math]. Given a section [math]\displaystyle{ \phi:M\to\Phi }[/math], fix a point [math]\displaystyle{ x\in M. }[/math] The pushforward at [math]\displaystyle{ x }[/math] is a map of tangent bundles

[math]\displaystyle{ \mathrm{d}_x\phi:T_xM\to T_{\phi(x)}\Phi\quad }[/math] taking [math]\displaystyle{ \quad \partial_\mu \mapsto \frac{\partial \phi^i}{\partial x^\mu} \partial_i }[/math]

where [math]\displaystyle{ \partial_\mu=\partial/\partial x^\mu }[/math] is taken to be an orthonormal vector space basis on [math]\displaystyle{ TM }[/math] and [math]\displaystyle{ \partial_i=\partial/\partial q^i }[/math] the vector space basis on [math]\displaystyle{ T\Phi }[/math]. The [math]\displaystyle{ \mathrm{d}\phi }[/math] is a differential form. The sigma model action is then just the conventional inner product on vector-valued k-forms

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \int_M \mathrm{d}\phi \wedge {*\mathrm{d}\phi} }[/math]

where the [math]\displaystyle{ \wedge }[/math] is the wedge product, and the [math]\displaystyle{ * }[/math] is the Hodge star. This is an inner product in two different ways. In the first way, given any two differentiable forms [math]\displaystyle{ \alpha,\beta }[/math] in [math]\displaystyle{ M }[/math], the Hodge dual defines an invariant inner product on the space of differential forms, commonly written as

[math]\displaystyle{ \langle\!\langle\alpha,\beta\rangle\!\rangle \ =\ \int_M \alpha\wedge {*\beta} }[/math]

The above is an inner product on the space of square-integrable forms, conventionally taken to be the Sobolev space [math]\displaystyle{ L^2. }[/math] In this way, one may write

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \langle\!\langle \mathrm{d}\phi, \mathrm{d}\phi\rangle\!\rangle }[/math]

This makes it explicit and plainly evident that the sigma model is just the kinetic energy of a point particle. From the point of view of the manifold [math]\displaystyle{ M }[/math], the field [math]\displaystyle{ \phi }[/math] is a scalar, and so [math]\displaystyle{ \mathrm{d}\phi }[/math] can be recognized just the ordinary gradient of a scalar function. The Hodge star is merely a fancy device for keeping track of the volume form when integrating on curved spacetime. In the case that [math]\displaystyle{ M }[/math] is flat, it can be completely ignored, and so the action is

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \int_M \Vert \nabla\phi\Vert^2 d^mx }[/math]

which is the Dirichlet energy of [math]\displaystyle{ \phi }[/math]. Classical extrema of the action (the solutions to the Lagrange equations) are then those field configurations that minimize the Dirichlet energy of [math]\displaystyle{ \phi }[/math]. Another way to convert this expression into a more easily-recognizable form is to observe that, for a scalar function [math]\displaystyle{ f:M\to\mathbb{R} }[/math] one has [math]\displaystyle{ \mathrm{d}*f=0 }[/math] and so one may also write

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \langle\!\langle \phi, \Delta\phi\rangle\!\rangle }[/math]

where [math]\displaystyle{ \Delta }[/math] is the Laplace–Beltrami operator, i.e. the ordinary Laplacian when [math]\displaystyle{ M }[/math] is flat.

That there is another, second inner product in play simply requires not forgetting that [math]\displaystyle{ \mathrm{d}\phi }[/math] is a vector from the point of view of [math]\displaystyle{ \Phi }[/math] itself. That is, given any two vectors [math]\displaystyle{ v,w\in T\Phi }[/math], the Riemannian metric [math]\displaystyle{ g_{ij} }[/math] defines an inner product

[math]\displaystyle{ \langle v,w\rangle = g_{ij}v^i w^j }[/math]

Since [math]\displaystyle{ \mathrm{d}\phi }[/math] is vector-valued [math]\displaystyle{ \mathrm{d}\phi = (\mathrm{d}\phi^1,\cdots, \mathrm{d}\phi^n) }[/math] on local charts, one also takes the inner product there as well. More verbosely,

[math]\displaystyle{ \mathcal{S} = \frac{1}{2}\int_M g_{ij}(\phi)\; \mathrm{d}\phi^i \wedge {*\mathrm{d}\phi^j} }[/math]

The tension between these two inner products can be made even more explicit by noting that

[math]\displaystyle{ B_{\mu\nu}(\phi)= g_{ij} \partial_\mu\phi^i \partial_\nu\phi^j }[/math]

is a bilinear form; it is a pullback of the Riemann metric [math]\displaystyle{ g_{ij} }[/math]. The individual [math]\displaystyle{ \partial_\mu\phi^i }[/math] can be taken as vielbeins. The Lagrangian density of the sigma model is then

[math]\displaystyle{ \mathcal{L} = \frac{1}{2} g^{\mu\nu}B_{\mu\nu} }[/math]

for [math]\displaystyle{ g_{\mu\nu} }[/math] the metric on [math]\displaystyle{ M. }[/math] Given this gluing-together, the [math]\displaystyle{ \mathrm{d}\phi }[/math] can be interpreted as a solder form; this is articulated more fully, below.

Motivations and basic interpretations

Several interpretational and foundational remarks can be made about the classical (non-quantized) sigma model. The first of these is that the classical sigma model can be interpreted as a model of non-interacting quantum mechanics. The second concerns the interpretation of energy.

Interpretation as quantum-mechanics

This follows directly from the expression

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \langle\!\langle \phi, \Delta\phi\rangle\!\rangle }[/math]

given above. Taking [math]\displaystyle{ \Phi=\mathbb{C} }[/math], the function [math]\displaystyle{ \phi:M\to\mathbb{C} }[/math] can be interpreted as a wave function, and its Laplacian the kinetic energy of that wave function. The [math]\displaystyle{ \langle\!\langle \cdot, \cdot\rangle\!\rangle }[/math] is just some geometric machinery reminding one to integrate over all space. The corresponding quantum mechanical notation is [math]\displaystyle{ \phi=|\psi\rangle. }[/math] In flat space, the Laplacian is conventionally written as [math]\displaystyle{ \Delta=\nabla^2 }[/math]. Assembling all these pieces together, the sigma model action is equivalent to

[math]\displaystyle{ \mathcal{S} = \frac{1}{2} \int_M \langle \psi| \nabla^2 |\psi\rangle dx^m = \frac{1}{2} \int_M \psi^\dagger(x) \nabla^2 \psi(x) dx^m }[/math]

which is just the grand-total kinetic energy of the wave-function [math]\displaystyle{ \psi(x) }[/math], up to a factor of [math]\displaystyle{ \hbar/m }[/math]. To conclude, the classical sigma model on [math]\displaystyle{ \mathbb{C} }[/math] can be interpreted as the quantum mechanics of a free, non-interacting quantum particle. Obviously, adding a term of [math]\displaystyle{ V(\phi) }[/math] to the Lagrangian results in the quantum mechanics of a wave-function in a potential. Taking [math]\displaystyle{ \Phi=\mathbb{C}^n }[/math] is not enough to describe the [math]\displaystyle{ n }[/math]-particle system, in that [math]\displaystyle{ n }[/math] particles require [math]\displaystyle{ n }[/math] distinct coordinates, which are not provided by the base manifold. This can be solved by taking [math]\displaystyle{ n }[/math] copies of the base manifold.

The solder form

It is very well-known that the geodesic structure of a Riemannian manifold is described by the Hamilton–Jacobi equations. In thumbnail form, the construction is as follows. Both [math]\displaystyle{ M }[/math] and [math]\displaystyle{ \Phi }[/math] are Riemannian manifolds; the below is written for [math]\displaystyle{ \Phi }[/math], the same can be done for [math]\displaystyle{ M }[/math]. The cotangent bundle [math]\displaystyle{ T^*\Phi }[/math], supplied with coordinate charts, can always be locally trivialized, i.e.

[math]\displaystyle{ \left. T^*\Phi\right|_U \cong U \times \mathbb {R}^n }[/math]

The trivialization supplies canonical coordinates [math]\displaystyle{ (q^1,\cdots,q^n,p_1,\cdots,p_n) }[/math] on the cotangent bundle. Given the metric tensor [math]\displaystyle{ g_{ij} }[/math] on [math]\displaystyle{ \Phi }[/math], define the Hamiltonian function

[math]\displaystyle{ H(q,p)=\frac{1}{2}g^{ij}(q)p_ip_j }[/math]

where, as always, one is careful to note that the inverse of the metric is used in this definition: [math]\displaystyle{ g^{ij}g_{jk}=\delta^i_k. }[/math] Famously, the geodesic flow on [math]\displaystyle{ \Phi }[/math] is given by the Hamilton–Jacobi equations

[math]\displaystyle{ \dot{q}^i=\frac{\partial H}{\partial p_i}\quad }[/math] and [math]\displaystyle{ \quad\dot{p}_i=-\frac{\partial H}{\partial q^i} }[/math]

The geodesic flow is the Hamiltonian flow; the solutions to the above are the geodesics of the manifold. Note, incidentally, that [math]\displaystyle{ dH/dt=0 }[/math] along geodesics; the time parameter [math]\displaystyle{ t }[/math] is the distance along the geodesic.

The sigma model takes the momenta in the two manifolds [math]\displaystyle{ T^*\Phi }[/math] and [math]\displaystyle{ T^*M }[/math] and solders them together, in that [math]\displaystyle{ \mathrm{d}\phi }[/math] is a solder form. In this sense, the interpretation of the sigma model as an energy functional is not surprising; it is in fact the gluing together of two energy functionals. Caution: the precise definition of a solder form requires it to be an isomorphism; this can only happen if [math]\displaystyle{ M }[/math] and [math]\displaystyle{ \Phi }[/math] have the same real dimension. Furthermore, the conventional definition of a solder form takes [math]\displaystyle{ \Phi }[/math] to be a Lie group. Both conditions are satisfied in various applications.

Results on various spaces

The space [math]\displaystyle{ \Phi }[/math] is often taken to be a Lie group, usually SU(N), in the conventional particle physics models, O(N) in condensed matter theories, or as a symmetric space in supergravity models. Since symmetric spaces are defined in terms of their involution, their tangent space (i.e. the place where [math]\displaystyle{ \mathrm{d}\phi }[/math] lives) naturally splits into even and odd parity subspaces. This splitting helps propel the dimensional reduction of Kaluza–Klein theories.

On Lie groups

For the special case of [math]\displaystyle{ \Phi }[/math] being a Lie group, the [math]\displaystyle{ g_{ij} }[/math] is the metric tensor on the Lie group, formally called the Cartan tensor or the Killing form. The Lagrangian can then be written as the pullback of the Killing form. Note that the Killing form can be written as a trace over two matrices from the corresponding Lie algebra; thus, the Lagrangian can also be written in a form involving the trace. With slight re-arrangements, it can also be written as the pullback of the Maurer–Cartan form.

On symmetric spaces

A common variation of the sigma model is to present it on a symmetric space. The prototypical example is the chiral model, which takes the product

[math]\displaystyle{ G=SU(N)\times SU(N) }[/math]

of the "left" and "right" chiral fields, and then constructs the sigma model on the "diagonal"

[math]\displaystyle{ \Phi=\frac{SU(N)\times SU(N)}{SU(N)} }[/math]

Such a quotient space is a symmetric space, and so one can generically take [math]\displaystyle{ \Phi=G/H }[/math] where [math]\displaystyle{ H\subset G }[/math] is the maximal subgroup of [math]\displaystyle{ G }[/math] that is invariant under the Cartan involution. The Lagrangian is still written exactly as the above, either in terms of the pullback of the metric on [math]\displaystyle{ G }[/math] to a metric on [math]\displaystyle{ G/H }[/math] or as a pullback of the Maurer–Cartan form.

Trace notation

In physics, the most common and conventional statement of the sigma model begins with the definition

[math]\displaystyle{ L_\mu=\pi_\mathfrak{m} \circ \left(g^{-1}\partial_\mu g \right) }[/math]

Here, the [math]\displaystyle{ g^{-1}\partial_\mu g }[/math] is the pullback of the Maurer–Cartan form, for [math]\displaystyle{ g\in G }[/math], onto the spacetime manifold. The [math]\displaystyle{ \pi_\mathfrak{m} }[/math] is a projection onto the odd-parity piece of the Cartan involution. That is, given the Lie algebra [math]\displaystyle{ \mathfrak{g} }[/math] of [math]\displaystyle{ G }[/math], the involution decomposes the space into odd and even parity components [math]\displaystyle{ \mathfrak{g}=\mathfrak{m}\oplus\mathfrak{h} }[/math] corresponding to the two eigenstates of the involution. The sigma model Lagrangian can then be written as

[math]\displaystyle{ \mathcal{L}=\frac{1}{2}\mathrm{tr}\left(L_\mu L^\mu\right) }[/math]

This is instantly recognizable as the first term of the Skyrme model.

Metric form

The equivalent metric form of this is to write a group element [math]\displaystyle{ g\in G }[/math] as the geodesic [math]\displaystyle{ g=\exp(\theta^i T_i) }[/math] of an element [math]\displaystyle{ \theta^i T_i\in \mathfrak{g} }[/math] of the Lie algebra [math]\displaystyle{ \mathfrak{g} }[/math]. The [math]\displaystyle{ [T_i,T_j]={f_{ij}}^k T_k }[/math] are the basis elements for the Lie algebra; the [math]\displaystyle{ {f_{ij}}^k }[/math] are the structure constants of [math]\displaystyle{ \mathfrak{g} }[/math].

Plugging this directly into the above and applying the infinitesimal form of the Baker–Campbell–Hausdorff formula promptly leads to the equivalent expression

[math]\displaystyle{ \mathcal{L} = \frac{1}{2} g_{ij}(\phi) \; \mathrm{d}\phi_i \wedge {*\mathrm{d}\phi_j} = \frac{1}{2} \; {W_i}^m {W^n}_j \;\; \mathrm{d}\phi_i \wedge {*\mathrm{d}\phi_j} \;\;\mathrm{tr} (T_m T_n) }[/math]

where [math]\displaystyle{ \mathrm{tr} (T_m T_n) }[/math] is now obviously (proportional to) the Killing form, and the [math]\displaystyle{ {W_i}^m }[/math] are the vielbeins that express the "curved" metric [math]\displaystyle{ g_{ij} }[/math] in terms of the "flat" metric [math]\displaystyle{ \mathrm{tr} (T_m T_n) }[/math]. The article on the Baker–Campbell–Hausdorff formula provides an explicit expression for the vielbeins. They can be written as

[math]\displaystyle{ W=\sum_{n=0}^\infty \frac{(-1)^nM^n}{(n+1)!} = (I-e^{-M})M^{-1} }[/math]

where [math]\displaystyle{ M }[/math] is a matrix whose matrix elements are [math]\displaystyle{ {M_j}^k = \theta^i{f_{ij}}^k }[/math].

For the sigma model on a symmetric space, as opposed to a Lie group, the [math]\displaystyle{ T_i }[/math] are limited to span the subspace [math]\displaystyle{ \mathfrak{m} }[/math] instead of all of [math]\displaystyle{ \mathfrak{g}=\mathfrak{m}\oplus\mathfrak{h} }[/math]. The Lie commutator on [math]\displaystyle{ \mathfrak{m} }[/math] will not be within [math]\displaystyle{ \mathfrak{m} }[/math]; indeed, one has [math]\displaystyle{ [ \mathfrak{m},\mathfrak{m}]\subset \mathfrak{h} }[/math] and so a projection is still needed.

Extensions

The model can be extended in a variety of ways. Besides the aforementioned Skyrme model, which introduces quartic terms, the model may be augmented by a torsion term to yield the Wess–Zumino–Witten model.

Another possibility is frequently seen in supergravity models. Here, one notes that the Maurer–Cartan form [math]\displaystyle{ g^{-1}dg }[/math] looks like "pure gauge". In the construction above for symmetric spaces, one can also consider the other projection

[math]\displaystyle{ A_\mu = \pi_\mathfrak{h}\circ\left(g^{-1}\partial_\mu g\right) }[/math]

where, as before, the symmetric space corresponded to the split [math]\displaystyle{ \mathfrak{g}= \mathfrak{m}\oplus \mathfrak{h} }[/math]. This extra term can be interpreted as a connection on the fiber bundle [math]\displaystyle{ M\times H }[/math] (it transforms as a gauge field). It is what is "left over" from the connection on [math]\displaystyle{ G }[/math]. It can be endowed with its own dynamics, by writing

[math]\displaystyle{ \mathcal{L}=g_{ij}F^i\wedge *F^j }[/math]

with [math]\displaystyle{ F^i=dA^i }[/math]. Note that the differential here is just "d", and not a covariant derivative; this is not the Yang–Mills stress-energy tensor. This term is not gauge invariant by itself; it must be taken together with the part of the connection that embeds into [math]\displaystyle{ L_\mu }[/math], so that taken together, the [math]\displaystyle{ L_\mu }[/math], now with the connection as a part of it, together with this term, forms a complete gauge invariant Lagrangian (which does have the Yang–Mills terms in it, when expanded out).

References

• Gell-Mann, M.; Lévy, M. (1960), "The axial vector current in beta decay", Il Nuovo Cimento 16 (4): 705–726, doi:10.1007/BF02859738, Bibcode: 1960NCim...16..705G 

• Ketov, Sergei (2009). "Nonlinear Sigma model". Scholarpedia 4 (1): 8508. doi:10.4249/scholarpedia.8508. Bibcode: 2009SchpJ...4.8508K.