Benutzer:CoTangent/Butcher-Gruppe
Im mathematischen Teilgebiet der Numerik bezeichnet die Butcher-Gruppe eine bestimmte Gruppe, die beim Studium von Runge-Kutter-Verfahren zur Lösung von gewöhnlichen Differenzialgleichungen auftritt. Sie ist nach dem neuseeländischen Mathematiker John C. Butcher benannt.
Connes und Kreimer haben 1999 herausgefunden, dass diese Gruppe außerdem die Charaktergruppe einer bestimmten Hopf-Algebra über gewurzelten Bäumen auftritt und dass diese bereits vollkommen unabhängig von der Verbindung zur Numerik in der Renormierung von Quantenfeldtheorien in der Physik verwendet wurde.
Motivation: Anfangswertproblem
Viele Vorgänge in Natur und Technik lassen sich mit Hilfe von gewöhnlichen Differentialgleichungen formulieren. Eine solche Gleichung lässt sich immer in die folgende Form bringen:
Hierbei ist ein (zeitunabhängiges) Vektorfeld auf einer offenen Menge . Wenn wie hier der Wert ist vorgegeben ist, spricht man auch von einem Anfangswertproblem.
Wenn das Vektorfeld beliebig oft differenzierbar ist, dann folgt aus dem Satz von Picard-Lindelöf, dass dieses Anfangswertproblem eine eindeutige Lösung besitzt und diese Lösungskurve selbst auch beliebig oft differenzierbar ist. Allerdings ist es meistens nicht möglich, die Lösung in einer geschlossenen Form anzugeben.
Eine Möglichkeit, die Lösung zu untersuchen ist ein Potenzreihenansatz um den Punkt : Dafür benötigt man die höheren Ableitungen der gesuchten Funktion, also . Hierzu kann man grundsätzlich die gegebene Differentialgleichung einfach mehrfach ableiten, aber da jedes Mal mehr Terme hinzukommen, wird die Sache schnell unübersichtlich. Im Jahre 1857 fand Arthur Cayley eine kompkate Formel, um die höheren Ableitungen einer Lösung eines Anfangswertproblems zu berechnen, hierfür benötigen wir aber als kombinatorische Hilfsmittel gewurzelte Bäume:
Differentiale und gewurzelte Bäume
Ein gewurzelter Baum ist ein Graph mit einem ausgezeichneten Knoten, genannt Wurzel, in dem jeder andere Knoten mit der Wurzel durch einen eindeutigen Weg verbunden ist. Wenn die Wurzel eines gewurzelten Baumes (mit allen Kanten, die an sie angrenzen) entfernt wird, dann zerfällt der Originalbaum in gewurzelte Bäume , wobei als Wurzeln der neuen Bäume immer die Knoten genommen werden, die mit der Orginalwurzel verbunden waren. Umgekehrt kann man immer aus verschiedenen gewurzelten Bäumen einen neuen Baum konstruieren, indem man eine neue Wurzel hinzufügt und diese mit allen bisherigen Wurzeln verbindet. Diesen neuen Baum bezeichnet man auch mit eckigen Klammern: .
Die Anzahl der Knoten eines Baumes bezeichnet man mit . Der Baum mit nur einem Knoten (die isolierte Wurzel) wird mit bezeichnet.
Eine Heap-Ordnung auf einem Baum ist eine Verteilung der Zahlen auf die Knoten des Graphen, sodass die Zahlen immer in aufsteigender Reihenfolge erscheinen, wenn man sich auf einem Weg von der Wurzel wegbewegt. Zwei Heap-Ordnungen heißen äquivalent, wenn einen Automorphismus von gewurzelten Bäumen gibt, der die eine Heap-Ordnung auf die andere überträgt. Die Anzahl der Äquivalenzklassen von Heap-Ordnungen wird mit bezeichnet und kann mit Butchers Formel berechnet werden:
wobei die Symmetriegruppe von ist und die Fakultät eines Baumes rekursiv definiert wird durch
und die Fakultät einer isolierten Wurzel auf gesetzt wird:
Als nächstes definieren wir die sogenannten elementaren Differentiale rekursiv:
Unter Verwendung dieser Notation gibt es nun eine kompakte Darstellung der -ten Ableitung der Lösungskurve:
Im einfachsten Fall N = 1 (dann sind x und f reellwertige Funktionen einer reellen Variable) ergibt die obige Formel:
wobei die vier Terme genau den vier gewurzelten Bäumen auf vier Knoten entsprechen (siehe Diagramm). Im Falle N = 1 entspricht diese Formel der bereits früher bekannten Formel_von_Faà_di_Bruno von 1855, in der die Baumstruktur keine Rolle spielt. Wenn aber , so muss man die Formel vorsichtiger in der folgenden Form schreiben: wobei die Baumstruktur entscheidend wird.
Damit erhalten wir eine geschlossene Formel für die Taylor-Reihe der Lösungskurve mit Entwicklungspunkt 0:
--
Butcher series and Runge–Kutta method
Das oben schon erwähnte Anfangswertproblem
lässt sich näherungsweise mit dem Runge-Kutta-Verfahren lösen. Bei diesem nimmt man eine -Matrix
und einen Vektor
mit m Komponenten.
Bei diesem Verfahren werden iterativ Näherungslösungen gefunden, indem zuerst eine LösungX1, ... , Xm of
and then setting
Vorlage:Harvtxt showed that the solution of the corresponding ordinary differential equations
has the power series expansion
where φj and φ are determined recursively by
and
The power series above are called B-series or Butcher series.[1][2] The corresponding assignment φ is an element of the Butcher group. The homomorphism corresponding to the actual flow has
Butcher showed that the Runge-Kutta method gives an nth order approximation of the actual flow provided that φ and Φ agree on all trees with n nodes or less. Moreover, Vorlage:Harvtxt showed that the homomorphisms defined by the Runge-Kutta method form a dense subgroup of the Butcher group: in fact he showed that, given a homomorphism φ', there is a Runge-Kutta homomorphism φ agreeing with φ' to order n; and that if given homomorphims φ and φ' corresponding to Runge-Kutta data (A, b) and (A' , b' ), the product homomorphism corresponds to the data
Vorlage:Harvtxt proved that the Butcher group acts naturally on the functions f. Indeed, setting
they proved that
As an example when N = 1, so that x and f are real-valued functions of a single real variable, the formula yields
where the four terms correspond to the four rooted trees from left to right in Figure 3 above.
In a single variable this formula is the same as Faà di Bruno's formula of 1855; however in several variables it has to be written more carefully in the form
where the tree structure is crucial.
Definition using Hopf algebra of rooted trees
The Hopf algebra H of rooted trees was defined by Vorlage:Harvtxt in connection with Kreimer's previous work on renormalization in quantum field theory. It was later discovered that the Hopf algebra was the dual of a Hopf algebra defined earlier by Vorlage:Harvtxt in a different context. The characters of H, i.e. the homomorphisms of the underlying commutative algebra into R, form a group, called the Butcher group. It corresponds to the formal group structure discovered in numerical analysis by Vorlage:Harvtxt.
The Hopf algebra of rooted trees H is defined to be the polynomial ring in the variables t, where t runs through rooted trees.
- Its comultiplication is defined by
where the sum is over all proper rooted subtrees s of t; is the monomial given by the product the variables ti formed by the rooted trees that arise on erasing all the nodes of s and connected links from t. The number of such trees is denoted by n(t\s).
- Its counit is the homomorphism ε of H into R sending each variable t to zero.
- Its antipode S can be defined recursively by the formula
The Butcher group is defined to be the set of algebra homomorphisms φ of H into R with group structure
The inverse in the Butcher group is given by
and the identity by the counit ε.
Using complex coefficients in the construction of the Hopf algebra of rooted trees one obtains the complex Hopf algebra of rooted trees. Its C-valued characters form a group, called the complex Butcher group GC. The complex Butcher group GC is an infinite-dimensional complex Lie group[3] which appears as a toy model in the Vorlage:Section link of quantum field theories.
Butcher series and Runge–Kutta method
The non-linear ordinary differential equation
can be solved approximately by the Runge-Kutta method. This iterative scheme requires an m x m matrix
and a vector
with m components.
The scheme defines vectors xn by first finding a solution X1, ... , Xm of
and then setting
Vorlage:Harvtxt showed that the solution of the corresponding ordinary differential equations
has the power series expansion
where φj and φ are determined recursively by
and
The power series above are called B-series or Butcher series.[1][4] The corresponding assignment φ is an element of the Butcher group. The homomorphism corresponding to the actual flow has
Butcher showed that the Runge-Kutta method gives an nth order approximation of the actual flow provided that φ and Φ agree on all trees with n nodes or less. Moreover, Vorlage:Harvtxt showed that the homomorphisms defined by the Runge-Kutta method form a dense subgroup of the Butcher group: in fact he showed that, given a homomorphism φ', there is a Runge-Kutta homomorphism φ agreeing with φ' to order n; and that if given homomorphims φ and φ' corresponding to Runge-Kutta data (A, b) and (A' , b' ), the product homomorphism corresponds to the data
Vorlage:Harvtxt proved that the Butcher group acts naturally on the functions f. Indeed, setting
they proved that
Lie algebra
Vorlage:Harvtxt showed that associated with the Butcher group G is an infinite-dimensional Lie algebra. The existence of this Lie algebra is predicted by a theorem of Vorlage:Harvtxt: the commutativity and natural grading on H implies that the graded dual H* can be identified with the universal enveloping algebra of a Lie algebra . Connes and Kreimer explicitly identify with a space of derivations θ of H into R, i.e. linear maps such that
the formal tangent space of G at the identity ε. This forms a Lie algebra with Lie bracket
is generated by the derivations θt defined by
for each rooted tree t.
The infinite-dimensional Lie algebra from Vorlage:Harvtxt and the Lie algebra L(G) of the Butcher group as an infinite-dimensional Lie group are not the same. The Lie algebra L(G) can be identified with the Lie algebra of all derivations in the dual of H (i.e. the space of all linear maps from H to R), whereas is obtained from the graded dual. Hence turns out to be a (strictly smaller) Lie subalgebra of L(G).[3]
Renormalization
Vorlage:Harvtxt provided a general context for using Hopf algebraic methods to give a simple mathematical formulation of renormalization in quantum field theory. Renormalization was interpreted as Birkhoff factorization of loops in the character group of the associated Hopf algebra. The models considered by Vorlage:Harvtxt had Hopf algebra H and character group G, the Butcher group. Vorlage:Harvtxt has given an account of this renormalization process in terms of Runge-Kutta data.
In this simplified setting, a renormalizable model has two pieces of input data:[5]
- a set of Feynman rules given by an algebra homomorphism Φ of H into the algebra V of Laurent series in z with poles of finite order;
- a renormalization scheme given by a linear operator R on V such that R satisfies the Rota-Baxter identity
- and the image of R – id lies in the algebra V+ of power series in z.
Note that R satisfies the Rota-Baxter identity if and only if id – R does. An important example is the minimal subtraction scheme
In addition there is a projection P of H onto the augmentation ideal ker ε given by
To define the renormalized Feynman rules, note that the antipode S satisfies
so that
The renormalized Feynman rules are given by a homomorphism of H into V obtained by twisting the homomorphism Φ • S. The homomorphism is uniquely specified by
Because of the precise form of Δ, this gives a recursive formula for .
For the minimal subtraction scheme, this process can be interpreted in terms of Birkhoff factorization in the complex Butcher group. Φ can be regarded as a map γ of the unit circle into the complexification GC of G (maps into C instead of R). As such it has a Birkhoff factorization
where γ+ is holomorphic on the interior of the closed unit disk and γ– is holomorphic on its complement in the Riemann sphere C with γ–(∞) = 1. The loop γ+ corresponds to the renormalized homomorphism. The evaluation at z = 0 of γ+ or the renormalized homomorphism gives the dimensionally regularized values for each rooted tree.
In example, the Feynman rules depend on additional parameter μ, a "unit of mass". Vorlage:Harvtxt showed that
so that γμ– is independent of μ.
The complex Butcher group comes with a natural one-parameter group λw of automorphisms, dual to that on H
for w ≠ 0 in C.
The loops γμ and λw · γμ have the same negative part and, for t real,
defines a one-parameter subgroup of the complex Butcher group GC called the renormalization group flow (RG).
Its infinitesimal generator β is an element of the Lie algebra of GC and is defined by
It is called the beta-function of the model.
In any given model, there is usually a finite-dimensional space of complex coupling constants. The complex Butcher group acts by diffeomorphims on this space. In particular the renormalization group defines a flow on the space of coupling constants, with the beta function giving the corresponding vector field.
More general models in quantum field theory require rooted trees to be replaced by Feynman diagrams with vertices decorated by symbols from a finite index set. Connes and Kreimer have also defined Hopf algebras in this setting and have shown how they can be used to systematize standard computations in renormalization theory.
Example
Vorlage:Harvtxt has given a "toy model" involving dimensional regularization for H and the algebra V. If c is a positive integer and qμ = q / μ is a dimensionless constant, Feynman rules can be defined recursively by
where z = 1 – D/2 is the regularization parameter. These integrals can be computed explicitly in terms of the Gamma function using the formula
In particular
Taking the renormalization scheme R of minimal subtraction, the renormalized quantities are polynomials in when evaluated at z = 0.
Notes
References
- Bergbauer, Christoph; Kreimer, Dirk (2005), “The Hopf Algebra of Rooted Trees in Epstein-Glaser Renormalization”, in Annales Henri Poincaré, volume 6, issue 2, Bibcode: 2005AnHP....6..343B, DOI: , pages 343–367
- Boutet de Monvel, Louis (2003), “Algèbre de Hopf des diagrammes de Feynman, renormalisation et factorisation de Wiener-Hopf (d'après A. Connes et D. Kreimer). [Hopf algebra of Feynman diagrams, renormalization and Wiener-Hopf factorization (following A. Connes and D. Kreimer)]”, in Astérisque (Séminaire Bourbaki)[1], volume 290, pages 149–165
- Brouder, Christian (2000), “Runge–Kutta methods and renormalization”, in Eur.Phys.J., volume C12, Bibcode: 2000EPJC...12..521B, DOI: , pages 521–534
- Bogfjellmo, G.; Schmeding, A. (2015), “The Lie group structure of the Butcher group”, in Foundations of Computational Mathematics[2], DOI:
- Brouder, Christian (2004), “Trees, Renormalization and Differential Equations”, in BIT Numerical Mathematics[3], volume 44, issue 3, DOI: , pages 425–438
- Butcher, J.C (1963), “Coefficients for the study of Runge-Kutta integration processes”, in J. Austral. Math. Soc., volume 3, issue 2, DOI: , pages 185–201
- Butcher, J.C (1972), “An algebraic theory of integration methods”, in Math. Comput., volume 26, issue 117, DOI: , JSTOR 2004720, pages 79–106
- Numerical methods for ordinary differential equations, 2nd edition, John Wiley & Sons Ltd., 2008, ISBN 978-0-470-72335-7
- Butcher, J.C (2009), “Trees and numerical methods for ordinary differential equations”, in Numerical Algorithms[4], Springer online
- Cayley, Arthur (1857), “On the theory of analytic forms called trees”, in Philosophical Magazine[5], volume XIII, pages 172–176 (also in Volume 3 of the Collected Works of Cayley, pages 242–246)
- Connes, Alain; Kreimer, Dirk (1998), “Hopf Algebras, Renormalization and Noncommutative Geometry”, in Communications in Mathematical Physics[6], volume 199, Bibcode: 1998CMaPh.199..203C, DOI: , pages 203–242
- Connes, Alain; Kreimer, Dirk (1999), “Lessons from quantum field theory: Hopf algebras and spacetime geometries”, in Letters in Mathematical Physics, volume 48, DOI: , pages 85–96
- Connes, Alain; Kreimer, Dirk (2000), “Renormalization in quantum field theory and the Riemann-Hilbert problem. I. The Hopf algebra structure of graphs and the main theorem”, in Comm. Math. Phys.[7], volume 210, Bibcode: 2000CMaPh.210..249C, DOI: , pages 249–273
- Connes, Alain; Kreimer, Dirk (2001), “Renormalization in quantum field theory and the Riemann-Hilbert problem. II. The β-function, diffeomorphisms and the renormalization group”, in Comm. In Math. Phys.[8], volume 216, Bibcode: 2001CMaPh.216..215C, DOI: , pages 215–241
- Gracia-Bondía, José; Várilly, Joseph C.; Figueroa, Héctor (2000) Elements of noncommutative geometry, Birkhäuser, ISBN 0-8176-4124-6, Chapter 14.
- Grossman, R.; Larson, R. (1989), “Hopf algebraic structures of families of trees”, in Journal Algebra[9], volume 26, archived from the original on 2008-08-20, pages 184–210
- Hairer, E.; Wanner, G. (1974), “On the Butcher group and general multi-value methods”, in Computing[10], volume 13, DOI: , pages 1–15
- Kreimer, Dirk (1998), “On the Hopf algebra structure of perturbative quantum field theories”, in Adv. Theor. Math. Phys., volume 2, Bibcode: 1997q.alg.....7029K, pages 303–334
- Kreimer, Dirk (1999), “Chen's iterated integral represents the operator product expansion”, in Adv. Theor. Math. Phys., volume 3, Bibcode: 1999hep.th....1099K, pages 627–670
- Kreimer, Dirk (2007) Factorization in Quantum Field Theory: An Exercise in Hopf Algebras and Local Singularities (Frontiers in Number Theory, Physics, and Geometry II), Springer, Bibcode: 2003hep.th....6020K, pages 715–736
- “On the structure of Hopf algebras”, in Annals of Mathematics (Second Series), volume 81, issue 2, 1965, DOI: , JSTOR 1970615, pages 211–264
[[Category:Combinatorics]] [[Category:Numerical analysis]] [[Category:Quantum field theory]] [[Category:Renormalization group]] [[Category:Hopf algebras]]
- ↑ a b Referenzfehler: Ungültiges
<ref>
-Tag; kein Text angegeben für Einzelnachweis mit dem Namen Butcher2008. - ↑ Jackson, K. R.; Kværnø, A.; Nørsett, S.P. (1994), “The use of Butcher series in the analysis of Newton-like iterations in Runge-Kutta formulas”, in Applied Numerical Mathematics, volume 15, issue 3, DOI: , pages 341–356 (Special issue to honor professor J. C. Butcher on his sixtieth birthday)
- ↑ a b Referenzfehler: Ungültiges
<ref>
-Tag; kein Text angegeben für Einzelnachweis mit dem Namen :0. - ↑ Jackson, K. R.; Kværnø, A.; Nørsett, S.P. (1994), “The use of Butcher series in the analysis of Newton-like iterations in Runge-Kutta formulas”, in Applied Numerical Mathematics, volume 15, issue 3, DOI: , pages 341–356 (Special issue to honor professor J. C. Butcher on his sixtieth birthday)
- ↑ Vorlage:Harvnb