Next Article in Journal
Bateman Oscillators: Caldirola-Kanai and Null Lagrangians and Gauge Functions
Previous Article in Journal
Analytical Modeling of the Temporal Evolution of Epidemics Outbreaks Accounting for Vaccinations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spin Current in BF Theory

Department of Mathematical Science, Ritsumeikan University, 6-1-1 Matsugaoka, Otsu, Shiga 520-2102, Japan
Submission received: 10 April 2021 / Revised: 15 May 2021 / Accepted: 26 May 2021 / Published: 7 June 2021
(This article belongs to the Section Classical Physics)

Abstract

:
In this paper, a current that is called spin current and corresponds to the variation of the matter action in BF theory with respect to the spin connection A which takes values in Lie algebra so ( 3 , C ) , in self-dual formalism is introduced. For keeping the 2-form B i constraint (covariant derivation) D B i = 0 satisfied, it is suggested adding a new term to the BF Lagrangian using a new field ψ i , which can be used for calculating the spin current. The equations of motion are derived and the solutions are dicussed. It is shown that the solutions of the equations do not require a specific metric on the 4-manifold M, and one just needs to know the symmetry of the system and the information about the spin current. Finally, the solutions for spherically and cylindrically symmetric systems are found.

1. Introduction

The BF theory on 4-manifold M is a topological theory, which includes constraints when terms turn to gravity theory [1]. The fundamental variables are 2-form B Ω 2 ( M ; so ( 3 , 1 ) ) and spin connection ω , which takes values in Lie algebra so ( 3 , 1 ) , and all derivatives are linear and applied only on ω , which makes it easy for canonical formalism, finding the phase space, Hamiltonian equations, quantization, etc. [1]. This theory does not require a metric to be formulated, as the metric is a derived quantity from the solutions of B. That gives motivation to formulate Einstein’s gravity as a theory of 2-forms rather than the metric tensors, and so no pre-existing geometrical structure is needed to obtain the gravity. Let F ( ω ) Ω 2 ( M ; so ( 3 , 1 ) ) be the curvature of ω . The pure BF theory action is M Tr ( B F ( ω ) ) , which is invariant (symmetric) under local Lorentz transformation (regarded as gauge group) and under arbitrary diffeomorphisms of M, and does not need using a metric. The equations of motion are F ( ω ) = 0 and d ω B = 0 , where d ω stands for covaraint derivative with respect to the connection ω , thus, B defines a twisted de Rham cohomology class [ B ] H D R 2 ( M , so ( 3 , 1 ) ) , and the solution of F ( ω ) = 0 is unique up to gauge and diffeomorphism transformations. There are no local degrees of freedom because the system has so much symmetry that all solutions are locally equivalent under gauge transformation of the group S O ( 3 , 1 ) and under diffeomorphisms of M. Hence, the pure BF theory is a topological theory [2,3].
In constrained BF theory, the Lagrangian includes the constraint term φ I J K L B I J B K L . The traceless matrix φ plays the role of a Lagrangian multiplier that imposes the constraint on the 2-form B I J , so that its solutions are given in terms of 1-forms e I = e μ I d x μ , that is B I J = e I e J , where the capital letters I , J , = 0 , 1 , 2 , 3 are the Lorentz indices and Greek letters μ , ν , = 0 , 1 , 2 , 3 are the space-time tangent indices. The frame fields e μ I d x μ are regarded as gravitational fields, therefore, the constrained BF theory turns to general relativity theory; the reason is that when φ I J K L is not constant (like cosmological constant), the term φ I J K L B I J B K L breaks the diffeomorphisms invariance of BF action, thus, there are non-equivalent local solutions and so local degrees of freedom exist as known in general relativity in the vacuum. Since the field φ I J K L is not a physical variable, the equations of motion of general relativity do not to include it (see Appendix A). The problem with constrained BF theory is that the equation of motion, δ S / δ B = 0 , of the action variation with respect to B contains the non-physical variable φ I J K L , but one can remove it by taking the trace of the equations, but there is also a problem with the trace operation, as it reduces the equations to one equation, which is not enough for getting a solution. For that reason, the solutions of BF theory using the equation δ S / δ ω = 0 are searched for. In general, the equations of motion of constrained BF theory including matter give a relation between the curvature F I J ( ω ) and the frame fields Σ I J = e I e J (the Plebanski 2-form), in matrix notation, that is F = χ Σ + ξ Σ ¯ , where the bar indicates anti-frame field, and χ , ξ are symmetric matrices of scalar fields [4]. Therefore, the problem turns to finding χ and ξ .
Let us start with the definition of the spin current J and discuss its conservation in BF theory including matter (in general, a matter Lagrangian is not specified). The spin current J appears in the equations of motion as a source for d ω B by the equation * d ω B + J = 0 (’*’ is Hodge star operator), and, in order to get d ω B = 0 in this study, a new term is added to BF Lagrangian, like Tr ψ B F ( ω ) , using a new field ψ , which is seen as a redefinition B B + ψ B . One finds that the equation of motion of ψ is the same conservation equation D μ J μ = 0 of the spin current vector field J, where D μ is the covariant derivative. Furthermore, by choosing d ω B = 0 in the equations of motion, the spin current becomes a source for the field ψ instead of B and one gets a new formula (definition) for the spin current using ψ , and since the spin current regards symmetry of the system, the field ψ also regards that symmetry. One can see that the equations of BF theory can be solved only by solving the spin current equation, δ S / δ ω = 0 , J 0 , with d ω B = 0 and without solving the equation δ S / δ B = 0 , which includes the Lagrangian multiplier ϕ I J K L (a non-physical variable), and without using a gravitational metric on M, so that one just needs to use the spin current and know the symmetry of the system. That means that the BF equations can be solved only by using the coupling term M ω μ I J J I J μ , which makes them easy to solve, and makes the theory similar to the gauge theory. Furthermore, since ω I J is 1-form and J I J is a vector field, the term ω μ I J J I J μ is naturally defined on M without needing to use additional structures (like a metric), thus, solving the system equation using only that coupling term gives a topological theory, i.e, the theory turns to finding 1-forms and vector fields, similarly to Chern–Simons theory, which includes the Wilson loops as a source for the gauge field. That makes it easy to solve the equations in different cases of the spin current, e.g, point charge, straight line current, circular current, etc. The lines of the spin current can be described using any coordinates system, e.g, Euclidean coordinates, etc., so the BF theory can be studied in any coordinates system, but in order to avoid an effect of the coordinates on the lines of spin currents, the coordinates are left to be flat (not curved). Furthermore, since the spin current is the source for the field ψ , this field has singularities on the lines of that spin current. One can see that the solution of B I J can always be written as e I e J , so one gets the gravity theory. Finally, an example of explicit solution of the equations in the case of a spherical and cylindrical symmetric systems in static case just by finding the field ψ are given; actually, the field ψ is used for obtaining the spin current.

2. Spin Current in BF Theory

Let M be a connected oriented smooth 4-manifold and P M be an S O ( 3 , 1 ) -principal bundle with a spin connection ω , which is locally a 1-form with values in so ( 3 , 1 ) and F Ω 2 ( M ; so ( 3 , 1 ) P ) is its curvature. The BF theory action is invariant under global and local Lorentz transformation, which gives a conserved current, and it is called here a spin current. Before discussing the conservation of the spin current, let us introduce the self-dual formalism.
Definition 1.
The self-dual projection is a homomorphism
so ( 3 , 1 ) P = P × S O ( 3 , 1 ) so ( 3 , 1 ) so ( 3 , C ) P = P × S O ( 3 , 1 ) so ( 3 , C )
defined by
P I J i : so ( 3 , 1 ) so ( 3 , 1 ) C = so ( 3 , C ) so ( 3 , C ) so ( 3 , C ) ,
for P × S O ( 3 , 1 ) so ( 3 , 1 ) = P × so ( 3 , 1 ) P × so ( 3 , 1 ) S O ( 3 , 1 ) S O ( 3 , 1 ) , where Latin letters i , j , = 1 , 2 , 3 , with using the matrices [5]
P I J i = 1 2 ε i j k , for I = i , J = j , and P 0 j i = P j 0 i = i 2 δ j i , f o r I = 0 , J = j 0 ,
where ε i j k is the totally anti-symmetric Levi-Civita tensor, and δ j i is the Kronecker delta.
This relates to the fact that the complexified Lie algebra of S O ( 3 , 1 ) has the decomposition so ( 3 , 1 ) C = so ( 3 , C ) so ( 3 , C ) [6]. The new connection is locally an so ( 3 , C ) -valued 1-form A on M whose components are
A μ i = P I J i ω μ I J = 1 2 ε i j k ω μ j k i ω μ 0 i ,
and its curvature is
F i ( A ) = P I J i F I J ( ω ) = d A i + ε i j k A j A k ,
where d is the exterior derivative. The two form B I J is mapped to B i = P I J i B I J . The covariant derivative D μ = μ + ω μ I J acting on sections of T M so ( 3 , 1 , C ) P becomes D μ = μ + A μ , with A μ i j = ε i j k A μ k , where μ is the affine connection on the tangent (T) space T M .
Using the new variables one can write the Lagrangian of matter (without specifying matter fields), L matter ( e I , ω I J ) , as L matter ( B i , A i , B ¯ i , A ¯ i ) , where B ¯ i and A ¯ i (anti-self-dual representation) are the complex conjugation of B i and A i . The Urbantke formula (Equation (A4), Appendix B) writes the metric g μ ν using only the constrained B i without using the constrained B ¯ i . Furthermore, the self-dual connection A i is compatible with B i via d A B i = 0 , while the anti-self-dual connection A ¯ i is compatible with B ¯ i via d A ¯ B ¯ i = 0 . Here d A is the exterior covariant derivative with respect to the connection A. By that one may suppose
δ δ B ¯ i L matter = 0 , δ δ A ¯ i L matter = 0 ,
or just writing L matter ( B i , A i ) .
Definition 2.
Let A be the self-dual connection on the so ( 3 , C ) -bundle so ( 3 , C ) P M . Let L matter be the Lagrangian of matter fields on M. Then the spin current J i μ is defined as
J i μ = δ δ A μ i L matter .
The matter action, S matter , is required to be invariant under any infinitesimal local Lorentz transformation ω μ I J ω μ I J + D μ Λ I J for infinitesimal transformation parameter Λ I J Ω 0 ( M ; so ( 3 , 1 ) P ) . Now, let us assume S matter has this property. Then one gets the following.
Lemma 1.
The spin current J i μ given by J i μ = δ δ A μ i L matter in gravity theory is conserved [7].
Proof. 
Since S matter is invariant under infinitesimal gauge transformation Λ I J , it is invariant under (one may suggest the condition (4))
A μ i A μ i + D μ Λ i for Λ i = P I J i Λ I J Ω 0 ( M ; so ( 3 , C ) P ) .
The variation
S matter A i + D Λ i S matter A i = M d 4 x D μ Λ i δ δ A μ i L matter = M d 4 x Λ i D μ δ δ A μ i L matter + M d 3 x μ Λ i δ δ A μ i L matter
vanishes for arbitrary Λ i only when D μ δ δ A μ i L matter = 0 , where Λ i is considered to vanish on the boundary M . Thus, the current J i μ = δ δ A μ i L matter is conserved. Actually, the previous calculation based on the idea that A i and A ¯ i transform independently under infinitesimal local Lorentz transformation ω μ I J ω μ I J + D μ Λ I J , therefore, there is another current that associates with the connection A ¯ i when the matter Lagrangian depends also on A ¯ i . □
One finds the same for the general relativity (GR) action; by using the variables ( Σ i , A i ) , one obtains the equation
M d 4 x D μ Λ i δ δ A μ i S GR = 0 M d 4 x Λ i D μ δ δ A μ i S GR = 0 .
In 3 + 1 decomposition of the space-time manifold M = Σ × R , let Σ t be space-like slice of constant time t, with the coordinates x a , a = 1 , 2 , 3 , let 0 be the time index. In the Hamilton–Jacobi system, by using the variables ( E i a , A a i , A 0 i ) on the slice of constant time Σ t , the equation becomes
Σ × R d 4 x Λ i D a δ δ A a i S GR Σ × R d 4 x Λ i D 0 δ δ A 0 i S GR = const . × Σ × R d 4 x Λ i D a E i a + D 0 ( D a E i a ) = 0 ,
which is satisfied when D a E i a = 0 , where E i a is conjugate momentum to A a i , and the relations E i a = const . × δ δ A a i S GR and D a E i a = const . × δ δ A 0 i S GR [5,8,9] are used. In this paper, one fixes D a E i a = 0 .
Remark 1.
To note is that the current J i μ is similar to the currents in Yang–Mills theory of the gauge fields, and one can see this clearly when regards the connection A μ i as a gauge field, by that the current J i μ relates to the local Lorentz invariance (local symmetry). The metric g μ ν = e μ I e ν J η I J ( η I J is Lorentz metric) is invariant under arbitrary local Lorentz transformations, like e μ I ( x ) U I J ( x ) e μ J ( x ) , for U ( x ) S O ( 3 , 1 ) , therefore, the local Lorentz symmetry is an internal degree of freedom.
Definition 3.
The action of BF theory including matter (without cosmological constant) on S O ( 3 , 1 ) -principal bundle P M is defined to be [10]
S = S topological + S constraints + S matter ,
with
S topological = M B i F i ( A ) , a n d S constraints = 1 2 M φ i j B i B j ,
where φ Γ M ; End ( so ( 3 , C ) P ) is a traceless matrix of scalar fields φ i j . Actually, it is not required to be symmetric since a new term to be added to BF Lagrangian (see the discussion below Equation (23)). The connection A on the Lie algebra bundle so ( 3 , C ) P , which is locally a 1-form with values in so ( 3 , C ) and its curvature F ( A ) Ω 2 M ; so ( 3 , C ) P are defined in the Equations (2) and (3). The index contraction is done by using δ i j , the Killing form on so ( 3 , C ) .
Hence,
S = M B i F i ( A ) + 1 2 φ i j B i B j + S matter .
Since the matrix φ is traceless, one can write φ i j = m i j ( m 11 + m 22 + m 33 ) δ i j / 3 , for some not traceless matrix ( m i j ) . The variation of the action with respect to m i j produces a quadratic equation in B i whose solution turns the theory into general relativity. These are
B i B j = 1 3 δ i j B i B j .
The solutions to this are all of the following form B i = P I J i e I e J , in which the gravitational fields e μ I are considered as frame fields [11]. Using the self-dual formula (Equation (1)), the constrained 2-form B i is written as
B i = 1 2 ε i j k e j e j i e 0 e i = Σ i ,
this is B i constrained = Σ i , using the notation Σ i = P I J i e I e J .
The equation of motion with respect to B i is
F i ( A ) + φ i j B j + δ S matter δ B i = 0 ,
or
F i ( A ) = φ i j B j δ S matter δ B i .
Since F i ( A ) Ω 2 M ; so ( 3 , C ) P is 2-form with values in so ( 3 , C ) , the δ δ B i S matter is also 2-form with values in so ( 3 , C ) .
Lemma 2.
In constrained B i , the variation δ δ B i S matter Ω 2 M ; so ( 3 , C ) P has the form
δ S matter δ B i constraint = T i j Σ j + ξ i j Σ ¯ j ,
for some matrices T i j , ξ i j Γ M ; End ( so ( 3 , C ) P ) , with T i j = T j i and ξ i j = ξ ¯ j i (see Appendix B, for more details).
Therefore, in the vacuum, T i j = 0 is set. Using this formula in Equation (7) implies
F i ( A ) = φ i j Σ j T i j Σ j ξ i j Σ ¯ j ,
or
F i ( A ) = ψ i j Σ j ξ i j Σ ¯ j ,
for some matrix ψ = φ T Γ M ; End ( so ( 3 , C ) P ) [12].
Since Tr ( φ ) = 0 , so Tr ( ψ ) = Tr ( T ) , Equation (8) yields
Σ i μ ν F μ ν i = Tr ( T ) for Σ μ ν i Σ j μ ν = δ j i , and Σ μ ν i Σ ¯ j μ ν = 0 .
Thus, in the vacuum, T i j = 0 , one has: Σ i μ ν F μ ν i = 0 . Σ i μ ν F μ ν i is called the Tr F here. Equation (9) does not contain the non-physical variable φ , but the problem with it is that the trace process decreasing the number of equations. Therefore, Σ i μ ν F μ ν i = Tr ( T ) is a condition on the solutions. The (0,2) tensor Σ j μ ν is inverse of the 2-form Σ μ ν i (Appendix C).
The equation of motion with respect to the connection A i is
D B i + δ δ A i S matter = 0 ,
or
ε μ ν ρ σ D ν B ρ σ i + J μ i = 0 .
where D B i = d B i + ε i j k A j B k .
One can see that ε μ ν ρ σ D ν B ρ σ i = 0 cannot be chosen when J μ i 0 , but the condition ε μ ν ρ σ D ν B ρ σ i = 0 leads to the constraint D a E a i = 0 which is satisfied in the Hamilton–Jacobi system, Equation (5). One gets D a E a i = 0 from ε μ ν ρ σ D ν B ρ σ i = 0 by setting μ = 0 , so ε 0 a b c D a B b c i = 0 , then ε 0 a b c ε a b c is used to get ε a b c D a B b c i = 2 D a E a i = 0 , where E a i = ε a b c B b c i / 2 is conjugated to the connection A a i on space-like slice of constant time on which the coordinates x a are used. Furthermore, the condition D B i = 0 is necessary when the connection A i is flat, by that the 2-form B i belongs to the twisted de Rham cohomology classes H D R 2 ( M , so ( 3 , C ) P ) , and this is necessary for getting a topological theory. One can solve that problem by adding new terms to the BF action (6) with which there are many possibilities for controlling Equation (10) for J μ i 0 with choosing D B i = 0 . Only some simple possibilities are chosen below in order to get simple results.
By acting by D μ on Equation (10), one gets:
ε μ ν ρ σ D μ D ν B ρ σ i + D μ J μ i = 0 ,
and using D μ J μ i = 0 , one obtains
ε μ ν ρ σ D μ D ν B ρ σ i = ε μ ν ρ σ [ D μ , D ν ] B ρ σ i / 2 = 0 ,
but [ D μ , D ν ] = F μ ν ( A ) , therefore,
ε μ ν ρ σ ( F μ ν ( A ) ) i j B ρ σ j = 0 .
Then using ( F μ ν ( A ) ) i j = ε i j k F μ ν k ( A ) , implies
ε μ ν ρ σ ε i j k F μ ν j ( A ) B ρ σ k = 0 .
As it is shown below, one can regard Equation (11) as an equation of motion with respect to a new field ψ i , with the possibility of choosing D B i = 0 with J μ i 0 .
In order to include the constraints D μ J μ i = 0 and D B i = 0 in BF theory, the following action is suggested.
Definition 4.
A new term is added to the BF action (6) to get
S = M B i F i ( A ) + 1 2 φ i j B i B j + M ε i j k ψ i B j F k ( A ) + S matter ,
in which M ε i j k ψ i B j F k ( A ) is added, for some vector field ψ i Γ M ; so ( 3 , C ) P . One can relate the new term to a redefinition like B i B i + ε i j k ψ j B k in pure BF Lagrangian B i F i ( A ) .
The equation of motion of this action with respect to the field ψ i is
ε i j k ε μ ν ρ σ B μ ν j F ρ σ k ( A ) = 0 ,
which is the same as Equation (11). By using ( F μ ν ( A ) ) i j = ε i j k F μ ν k ( A ) , one gets
ε μ ν ρ σ B μ ν j ( F ρ σ ( A ) ) i j = 0 ,
but [ D μ , D ν ] = F μ ν ( A ) . Therefore,
ε μ ν ρ σ [ D ρ , D σ ] B μ ν i = 0 ,
but ε μ ν ρ σ D ρ B μ ν i = 0 ( D B i = 0 ) is chosen as suggested before, thus,
δ δ ψ i S matter = 0
is satisfied.
The equation of motion of this action with respect to the connection A k is
D B k + D ε k i j ψ i B j + δ δ A k S matter = 0 ,
or
ε μ ν ρ σ D ν B ρ σ k + ε μ ν ρ σ D ν ε k i j ψ i B ρ σ j + δ δ A μ k S matter = 0 .
In this equation, one can choose the condition ε μ ν ρ σ D ν B ρ σ k = 0 , which is equivalent to D μ Σ i μ ν = 0 in constrained B i , since B i constraint = Σ i = P I J i e I e J and [13,14]
1 2 ! e 1 ε μ ν ρ σ Σ ρ σ i = Σ i μ ν = i Σ i μ ν = i Σ μ ν i , e = det ( e μ I ) ,
where the Hodge duality theory between the forms and the tensor fields is used; here, Σ μ ν i is 2-form and Σ i μ ν is (0,2)-tensor field. One can see that Σ i μ ν is inverse of Σ μ ν i , so that Σ j μ ν Σ μ ν i = δ j i (see Appendix C for more details).
With that, the term ε μ ν ρ σ D ν B ρ σ k in constrained B i becomes i D ν Σ μ ν i , and so one can choose the condition D ν Σ μ ν i = 0 , which is locally equivalent to D B i = 0 in constrained B i .
The remaining equation of (14) in constrained B i is
D ε k i j ψ i Σ j + δ δ A k S matter = 0 ,
or
ε μ ν ρ σ ε k i j D ν ψ i Σ ρ σ j + δ δ A μ k S matter = 0 ,
hence,
2 i ε k i j D ν ψ i Σ μ ν j + J μ k = 0 ,
in which the spin current J i μ = δ S matter / δ A μ i is used and the condition D ν Σ μ ν j = 0 is imposed. Below, the condition D μ J μ k = 0 is discussed. Here, both Σ μ ν j and J μ k are tensor fields.
Remark 2.
To note is that Equation (16) is similar to the current J μ n = μ φ i T i j n φ j in scalar field theory with symmetry and generators T i j n , so one has J 0 n = 0 φ i T i j n φ j = π i T i j n φ j , where π i is conjugate momentum to φ i . Similarly, Equation (16) gives J 0 i = 2 i ( D μ ψ k ) ε k i j Σ j μ 0 (for μ = a = 1 , 2 , 3 ), so here D a ψ i is conjugate momentum to Σ i a 0 , noting that the indices raising in Σ j μ ν is done by using a metric g μ ν .
The equation of motion of the action (12) with respect to B i in constrained BF (like deriving Equation (8)) is
F i ( A ) + φ i j Σ j + ε i j k ψ j F k ( A ) + T i j Σ j + ξ i j Σ ¯ j = 0 .
Multiplying by Σ i μ ν , summing over the indices and using Σ i μ ν Σ ¯ μ ν j = 0 , one obtains:
Σ i μ ν F μ ν i + φ i j Σ i μ ν Σ μ ν j + ε i j k ψ j Σ i μ ν F μ ν k + T i j Σ i μ ν Σ μ ν j = 0 .
Then, using Σ i μ ν Σ μ ν j = δ i j to obtain
Σ i μ ν F μ ν i + Tr ( φ ) + ε i j k ψ j Σ i μ ν F μ ν k + Tr ( T ) = 0 .
Since Tr ( φ ) = 0 and ε i j k Σ i μ ν F μ ν k = 0 (see Equations (13) and (15)), one finds:
Σ i μ ν F μ ν i + Tr ( T ) = Tr ( F ) + Tr ( T ) = 0 .
Equation (17) allows us to write F i ( A ) in terms of Σ i and Σ ¯ i , and since ε i j k Σ i μ ν F μ ν k = 0 , one can write
F i ( A ) = χ i j Σ j + χ i j Σ ¯ j ,
for some symmetric matrix ( χ i j ) and skew-hermitian matrix ( χ i j ) . Using this equation in Equation (19), one obtains:
Tr ( χ i j ) + Tr ( T i j ) = 0 .
In addition to this relation, there is another relation between the vector field ψ i and the symmetric matrix χ i j when T i j 0 and J μ i 0 , from the conservation of the current (16), D ν J ν i = 0 , one has (for D μ Σ μ ν i = 0 ):
i 2 D ν J ν i = D ν D μ ψ k ε k i j Σ μ ν j = 1 2 D ν , D μ ψ k ε k i j Σ μ ν j = 1 2 F μ ν ( A ) ε k m ψ m ε k i j Σ μ ν j = 1 2 F μ ν ( A ) ε k m ψ m ε k i j Σ j μ ν = 1 2 F μ ν ( A ) ε k m ψ m ε k i j Σ j μ ν = 1 2 F μ ν ( A ) ψ m δ i δ m j δ m i δ j Σ j μ ν = 1 2 F μ ν i ( A ) ψ j Σ j μ ν 1 2 F μ ν j ( A ) ψ i Σ j μ ν ,
and using Equation (20), one gets:
i 2 D ν J ν i = 1 2 χ i m Σ μ ν m ψ j Σ j μ ν 1 2 χ j m Σ μ ν m ψ i Σ j μ ν = 1 2 χ i m ψ j δ j m 1 2 χ j m ψ i δ j m = 1 2 χ i j ψ j 1 2 ψ i t r χ i j = 0 , for J μ i 0 .
This is another relation between the vector field ψ i and the symmetric matrix χ i j in existence of matter T i j 0 with J μ i 0 . In this case, the matrix χ i j has to satisfy det χ i j Tr ( χ i j ) = 0 in order to get ψ i 0 ; of course, this condition is not needed in the vacuum T i j = 0 , J μ i = 0 .
Using Equation (20) in (17), one obtains:
χ i j Σ j + φ i j Σ j + ε i j k ψ j χ k Σ + T i j Σ j + ( ) i j Σ ¯ j = 0 .
That yields
χ i + φ i + ε i j k ψ j χ k + T i = 0 .
This equation relates to the equation of motion δ S / δ B = 0 , and it includes the Lagrangian multiplier φ i j , which is a non-physical variable that makes (23) difficult to solve. Therefore, one needs to find χ and ψ using the other equations of motion obtained above. One can see that φ is not required to be symmetric matrix since the third term in (23) is not symmetric in general. The symmetric matrix T i j is assumed to be given using the matter Lagrangian (Appendix B), thus, the total unknown variables are 3 + 5 + 8 = 16 of the vector ψ , the symmetric matrix χ (with (21)) and the traceless matrix φ . Equation (23) gives 9 equations, therefore, there are 16 9 = 7 unknown variables, but when J μ i 0 , they reduce to 6 unknown variables (regarding Equation (22)). However, if one chooses a solution for which the symmetric matrix χ i j becomes diagonal, like
χ = K i δ j i = diag ( K 1 , K 2 , K 3 ) ,
for some scalar functions K 1 , K 2 and K 3 on M. Thus, the unknown variables reduce to 4 variables and to 3 variables when J μ i 0 .
Remark 3.
The field ψ i is a solution of D μ D μ ψ i = 0 (see Appendix C), so if D μ v i = 0 , then ψ i + v i is another solution, and that makes the components ψ 1 , ψ 2 and ψ 3 of the vector field ψ i independent variables, therefore, one can regard them as the degrees of freedom of the system and solve the equations of motions in terms of them. Note that ψ i ψ i + v i ( D v i = 0 ) does not change the current J ν i = 2 i D μ ψ k ε k i j Σ μ ν j .
The Bianchi identity D F i = 0 implies D χ i j Σ j = 0 (for D Σ i = 0 ), hence ( d K i ) δ i j Σ j = 0 , where D δ i j = 0 is used along with the covariant derivative D v i = d v i + ε i j k A j v k . Therefore, one obtains:
ε μ ν ρ σ ( ν K i ) Σ ρ σ i = 0 ,
where ν d / d x ν . In 3 + 1 decomposition of the space-time manifold M = Σ × R , let Σ t be the space-like slice of constant time t with the coordinates x a with a = 1 , 2 , 3 , and 0 is the time index. The equation ϵ μ ν ρ σ D μ Σ ν ρ i = 0 ( D B i = 0 ) decomposes into two equations,
D a E a i = 0 and ϵ a b c D b B c i = 0 ,
in which the vector field E i and the 1-form B i :
E a i = ϵ 0 a b c Σ b c i / 2 = ϵ a b c Σ b c i / 2 , B c i = Σ 0 c i ,
are introduced on the space-like slice Σ t (the field E a i is conjugate to the connection A a i ).
Equation (25) decomposes into (for 0 K i = 0 )
ε 0 a b c ( a K i ) Σ b c i = ( a K i ) ε a b c Σ b c i = 2 ( a K i ) E a i = 0 , ε a b c 0 ( b K i ) Σ c 0 i = ε a b c ( b K i ) B c i = 0 .
One can solve them by writing (for non-zero curvature F i ( A ) )
E a i = 1 2 ε a b c ( b K i ) r c i , B c i = ( c K i ) u i ,
for some r i Ω 1 M ; so ( 3 , C ) P and u i Γ M ; so ( 3 , C ) P . The functions K i are scalars, the indices are just for distinguishing each from the others. Thus, one gets the solutions
Σ a b i = ( [ a K i ) r b ] i , and Σ 0 a i = ( a K i ) u i .
Equation (26) implies D r i = 0 and D u i = 0 .
In the static case J i a = 0 (zero current) with J i 0 0 (non-zero charge), the spin current formula J k ν = 2 i ε k i j D μ ψ i Σ μ ν j decomposes into two equations:
J k b = 2 i ε k i j D a ψ i Σ a b j = 0 , and J k 0 = 2 i ε k i j D a ψ i Σ a 0 j 0 .
One can solve the first equation in terms of ψ i by writing Σ a b i = f v [ a D b ] ψ i , for some vector v Γ M ; T Σ that satisfies v a D a ψ i = 0 , and f is scalar function on M. Including f in v, one can just write Σ a b i = v [ a D b ] ψ i . Let us note that D a ψ i = g a b D b ψ i without a need the used metric g a b to be specified. Regarding the second equation of (30), when J 0 i = 0 , one gets the solution, Σ 0 a i = Σ a 0 i = f D a ψ i . Furthermore, when J 0 i 0 , we let Σ 0 a i = f D a ξ i for a vector field ξ i ψ i . To note is that no specific metric g a b is required for raising and lowering the indices a , b , on Σ t , so let it be the metric coming from pulling back of the Lorentz metric, where Σ t is kept to be immersed in R 4 .
Lemma 3.
By comparing the solutions Σ a b i = v [ a D b ] ψ i and Σ 0 a i = f D a ξ i of Equations (30) with the solutions (29), and in order to get a correspondence between that solutions, one finds that
ψ i = K i b i , r c i = v c b i , ξ i = K i u i , D b i = 0 , D u i = 0 , d v = 0 , f = 1 ,
for some vector fields, b i , u i Γ M ; so ( 3 , C ) P .
By that, one obtains the solutions,
E i = 1 2 ϵ a b c Σ a b i c = 1 2 ϵ a b c v a b K i b i c Γ M ; T Σ so ( 3 , C ) P , B i = Σ 0 a i d x a = d K i u i Γ M ; T * Σ so ( 3 , C ) P ,
without needing to use a specific metric.
Remark 4.
Regarding the solutions of Equations (29) and (30), let us note that for every two solutions of Σ a b i and Σ a b i , the metric g a b satisfies Σ a b i = g a a g b b Σ a b i . Furthermore, the metric used in Σ a b i = g a a g b b Σ a b i is not necessarily the same metric used in D a ψ i = g a b D b ψ i for getting the solutions of (30). It is convenient to start from a solution of Σ a b i , and by using a metric g a b , to obtain the corresponding solution of Σ a b i .
Remark 5.
In solution (32), one can see that Σ I J can be written as e I e J , as required in constraint BF theory to get gravity theory, that is, according to self-dual projection, there are vector fields b I and K I satisfying K i b i = P I J i b I ( K J b J ) , therefore, Σ a b I J = b I v [ a D b ] K J b J , then one can write e a I = v a b I and e b J = b K J b J . Furthermore, from Σ 0 a I J = u I a K J u J , one gets: e a J = a K J u J and e 0 I = v 0 u I , for v 0 = 1 . A more general case is to find three vector fields b 1 I , b 2 I and K I satisfying K i b i = P I J i b 1 I ( K J b 2 J ) , therefore, Σ a b I J = b 1 I v [ a D b ] K J b 2 J , then one can write e a I = v a b 1 I and e b J = b K J b 2 J . Furthermore, from Σ 0 a I J = u 1 I a K J u 2 J , one gets e a J = a K J u 2 J and e 0 I = v 0 u 1 I , for v 0 = 1 . By that, (32) can be written as Σ i = P I J i Σ I J for Σ I J = e I e J . However, to note is that solution (32) is a general solution and one has to find a special solution, like to let b i and u i be constant fields.
Using the solution of Σ 0 a i in Equation (30), one obtains:
J k 0 = Q k = 2 i ϵ k i j b i a K i u j a K j = i ϵ k i j b i a K i u j a K j b j a K j u i a K i = i ϵ k i j a K i a K j ( b i u j b j u i ) .
One can see that J k 0 0 takes place only when b u . Therefore, in the vacuum it must be b = u . If J k a 0 , it must be ψ i = g i b i with g i K i .
If the charges J 0 i 0 are given as functions on M, then letting b i Γ M ; so ( 3 , C ) P be constant field on M, one can determine the scalar functions K i using Equation (33), and obtaining the vector v Γ M ; T Σ using v a a K i = 0 . However, to satisfy D b i = d b i + ε j k i A j b k = 0 for a constant vector field b i , the connection A μ i must be written as A μ i = A μ b i . Furthermore, u i = b i + f ( x ) b i + a i is chosen, for a constant a i Γ M ; so ( 3 , C ) P satisfying a i b i = 0 , the function f is needed for satisfying D u i = 0 . Examples of determining b i and u i in spherical and cylindrical symmetries are given below. Then, one obtains B a i and E a i using Equation (32), and obtain the matrix χ using Equation (24), thus, obtaining the curvature F = χ Σ . Note that v a E a i = 0 , v a B a i = 0 and v a a K i = 0 depend on the symmetry of the system, for example, spherical symmetry, cylindrical symmetry, and so on. Thus, one sees that the equations of motion of BF theory can be solved without needing to use a gravitational metric on the manifold M.

3. Solutions for Spherically Symmetric System

It was shown above that one can solve the equations of motion in BF theory by using a complex vector field ψ i = K i b i Γ M ; so ( 3 , C ) P , which allows us to obtain v , E i , B i and J 0 i , according to Equations (31)–(33). Here, the solutions to be found for spherically symmetric system in the vacuum ( u = b ) and then apply it for matter located at a point. As it was seen above, the solution of the system regards the symmetry of that system, since one searches for a vector v Γ M ; T Σ that satisfies v a D a ψ i = 0 , v a E a i = 0 and v a B a i = 0 . For example, in spherical symmetry, the spherical coordinates ( r , θ , φ ) to be used on the space-like slice Σ t = Σ = R 3 . Letting the vector field ψ i to depend only on the radius r, one gets (for D b i = 0 ):
D 2 ψ i = D 2 ( K i b i ) = b i 2 K i = b i 1 r 2 r r 2 r K i = 0 K i = c i r ,
so ψ i = c i b i / r , for some constants c i R . Actually, one can include c i in b i and just write ψ i = b i / r . Therefore,
D ψ i = b i ( d K i ) = b i d r r + d θ θ + d φ φ 1 r = b i r 2 d r ,
thus, the 1-form v ( d v = 0 , g a b v a D b ψ i = 0 ) is
v = a 1 d θ + a 2 d φ , a 1 , a 2 R ,
where, in the spherical symmetry, a 1 and a 2 kept to not depend on the coordinates θ and φ . The values of the constants a 1 and a 2 are not significant since a i = g i j a j is the Killing vector, thus, set a 1 = a 2 = 1 . The used metric g a b here is the standard metric in the spherical coordinates, because no any other metric is defined. In what follows in this Section, the indices r, θ and φ denote the spherical components.
Using Equation (32), one gets the solutions of the 1-form B i and the vector field E i ,
B i = ( d K i ) u i = b i r 2 d r , f o r u = b ,
E i = 1 2 ε a b c Σ b c i a = 1 2 ε a b c v b c K i b i a = ε φ θ r b i 2 r 2 φ ε θ φ r b i 2 r 2 θ = ε φ r θ b i 2 r 2 φ + ε θ r φ b i 2 r 2 θ .
By that, one gets:
Σ 0 r i = b i r 2 , Σ r θ i = b i 2 r 2 , Σ r φ i = b i 2 r 2 ,
while the other components like Σ 0 θ i , Σ 0 φ i , are zeros.
One obtains the matrix χ using Equation (24) with the solution (34),
χ = K j δ j i = 1 r diag ( c 1 , c 2 , c 3 ) ,
where the constants c i have to be determined in order to satisfy the condition Tr χ = 0 (in the vacuum), so i = 1 3 c i = 0 . Thus, one gets the curvature F = χ Σ + χ Σ ¯ (with setting χ = 0 in the vacuum [15]),
F 0 r i = χ i j Σ 0 r j = c i b i r 3 , F r θ i = χ i j Σ r θ j = c i b i 2 r 3 , F r φ i = χ i j Σ r φ j = c i b i 2 r 3 .
Now, let us calculate the connection A i and the field b i , which satisfies D b i = 0 . Using F i = d A i + ε i j k A j A k , one obtains:
F 0 r i = 1 2 0 A r i r A 0 i + ε i j k A 0 j A r k , F θ r i = 1 2 θ A r i r A θ i + ε i j k A θ j A r k , F φ r i = 1 2 φ A r i r A φ i + ε i j k A φ j A r k .
Since a spherically symmetric system is under consideration, the connection A i are considered depending on r only. If the gauge A r i = 0 is chosen, then:
F 0 r i = 1 2 r A 0 i , F θ r i = 1 2 r A θ i , F φ r i = 1 2 r A φ i ,
and, therefore, using the solution (36), one obtains:
1 2 r A 0 i = c i b i r 3 , 1 2 r A θ i = c i b i 2 r 3 , 1 2 r A φ i = c i b i 2 r 3 .
However, D b i = d b i + ε j k i A j b k = 0 and μ b i = 0 for μ r , therefore,
r b i + ε i j k A r j b k = r b i = 0 , 0 b i + ε i j k A 0 j b k = ε i j k A 0 j b k = 0 , θ b i + ε i j k A θ j b k = ε i j k A θ j b k = 0 , φ b i + ε i j k A φ j b k = ε i j k A φ j b k = 0 .
Therefore, the field b i Γ M ; so ( 3 , C ) P is constant, and one gets the solution,
A 0 i = c i b i r 2 , A θ i = c i b i 2 r 2 , A φ i = c i b i 2 r 2 ,
where ε i j k b j b k = 0 is sued. Thus, in this solution the field b i Γ M ; so ( 3 , C ) P is constant on M = Σ × R . Next is to find b i in the case of matter located at a point.

Solutions for Matter Located at a Point

If there is matter located at a point in Σ t = Σ = R 3 , one, thus, has a spherically symmetric system in a static case J i 0 0 , J i a = 0 . Let the origin ( 0 , 0 , 0 ) R 3 to be that point, therefore, the charge (33) is given by Q i ( x ) = Q 0 i δ 3 ( x ) , so R 3 Q 0 i δ 3 ( x ) = Q 0 i = const . (conservation of the charges). In order to get the same solution as in Equations (35) and (36), the field b i is kept to be a constant, and in Equation (33),
J k 0 ( x ) = Q k ( x ) = i ϵ k i j a K i a K j ( b i u j b j u i ) ,
Q i ( x ) = Q 0 i δ 3 ( x ) is used.
For spherical symmetry, the functions K i are given by K i = c i / r (Equation (34)), therefore,
Q k ( x ) = 2 i g r r r K i r K j ϵ k i j b i u j = i 2 c i c j r 4 ϵ k i j b i u j .
Therefore, in order to get Q i = Q 0 i δ 3 ( x ) , one replaces 1 / r 4 with 1 / ( r 4 + ϵ 4 ) , for some infinitesimal parameter ϵ 0 + , and looks for a solution for the field u j like
u j = b j + ϵ f b j + ϵ a j / ( 2 i π 2 2 ) ,
for some function f on M that is needed for satisfying D u i = 0 and a constant vector field a i Γ M ; so ( 3 , C ) P . With that, one obtains (for i , j k ):
Q k ( r ) = 2 i c i c j r 4 + ϵ 4 ϵ k i j b i b j + ϵ f b j + ϵ 2 i π 2 2 a j = 1 π 2 2 ϵ r 4 + ϵ 4 ϵ k i j ( c i b i ) ( c j a j ) .
Comparing with Q i ( x ) = Q 0 i δ 3 ( x ) , one finds ϵ k i j ( c i b i ) ( c j a j ) = Q 0 k = const . , and, by imposing ( c i a i ) ( c i a i ) = 1 with ( c i b i ) ( c i a i ) = 0 :
( c b ) 2 = ( c i b i ) ( c i b i ) = Q 0 k Q 0 k ,
thus, c i b i = e i Q 0 k Q 0 k are chosen for e i = 1 , so that the constant field b i is determined by Q 0 k Q 0 k with free S O ( 3 , C ) rotation.
By that (for r > 0 , ϵ 0 + ), one obtains the same solutions as in Equations (35) and (36), but with c i b i = e i Q 0 k Q 0 k for e i = 1 and i = 1 3 c i = 0 . Since c i b i is a finite value, it is not sufficient to let the constants c i take arbitrary values, so chosen them to be ( c i ) = ( 1 , 1 , 2 ) .
By that, an example for the possibility of solving the equations of motion in BF theory without needing to use a gravitational metric on M is given, so that one needs just to use a vector field ψ i Γ M ; so ( 3 , C ) P , which is defined in the spin current J k μ = 2 i ε k i j D ν ψ i Σ μ ν j of matter, Equation (16). Furthermore, it is shown that the solutions depend on the symmetry of the system, since one needs some vector v that satisfies v a a K i = 0 , v a E a i = 0 , v a B a i = 0 for obtaining the solutions.

4. Solutions for Cylindrically Symmetric System

In a cylindrically symmetric system, the matter is considered to be homogeneously located along the Z-axis. Similar to the above-considered spherical symmetry, one searches for the field ψ i = K i b i Γ M ; so ( 3 , C ) P , and then calculates v , E i , B i and J 0 i , according to Equations (31)–(33). The vector v Γ M ; T Σ satisfies v a D a ψ i = 0 , v a E a i = 0 and v a B a i = 0 , thus, it is the Killing vector. Here, the solution in the vacuum ( u = b ) to be found and, then, to be to be applied to a matter located homogeneously along the Z-axis. The needed information for solving the equations of motion is only the spin charge Q i ( x ) , Equation (33). As it was mentioned before, there is no need to use a gravitational metric, a standard metric to be used instead. In cylindrical symmetry, the cylindrical coordinates ( ρ , φ , z ) on the space-like slice Σ t = Σ = R 3 of constant time t to be used. Letting the vector field ψ i to depend only on the radius ρ , one gets (for D 2 ψ i = 0 and D b i = 0 ):
D 2 ψ i = D 2 ( K i b i ) = b i 2 K i = b i 1 ρ ρ ρ ρ K i = 0 K i = c i log ( ρ ) ,
so ψ i = c i b i log ( ρ ) , for some constants c i R . Therefore,
D ψ i = c i b i ( d K i ) = c i b i d ρ ρ + d φ φ + d z z log ( ρ ) = c i b i ρ d ρ ,
thus, the 1-form v ( d v = 0 , g a b v a D b ψ i = 0 ) is
v = a 1 d φ + a 2 d z , a 1 , a 2 R ,
where, in the cylindrical symmetry, a 1 and a 2 are kept not depending on the coordinates z and φ . The values of the constants a 1 and a 2 are not significant since a i = g i j a j is Killing vector, thus, are set as a 1 = a 2 = 1 . The used metric g a b here is the standard metric in the cylindrical coordinates because no any other metric is defined. In what follows in this Section, the indices ρ , θ and φ denote the cylindrical components.
Using Equation (32), one gets the solutions of the 1-form B i and the vector field E i , namely,
B i = c i u i ρ d ρ = c i b i ρ d ρ , for u = b ,
E i = 1 2 ε a b c Σ b c i a = 1 2 ε a b c v b c K i b i a = ε φ z ρ c i b i 2 ρ φ + ε z φ ρ c i b i 2 ρ z .
By that, one gets:
Σ 0 ρ i = Σ ρ 0 i = c i b i ρ , Σ ρ z i = Σ z ρ i = c i b i 2 ρ , Σ ρ φ i = Σ φ ρ i = c i b i 2 ρ ,
while the other components like Σ 0 z i , Σ 0 φ i , are zeros.
The matrix χ is obtained using Equation (24) with the solution (38),
χ = K j δ j i = log ( ρ ) diag ( c 1 , c 2 , c 3 ) ,
where the constants c i have to be determined in order to satisfy the condition Tr χ = 0 (in the vacuum), so i = 1 3 c i = 0 . Thus, one gets the curvature F = χ Σ + χ Σ ¯ (with setting χ = 0 in the vacuum [15]),
F 0 ρ i = χ i j Σ 0 ρ j = log ( ρ ) ρ c i c i b i , F ρ z i = χ i j Σ ρ z j = log ( ρ ) 2 ρ c i c i b i , F ρ φ i = χ i j Σ ρ φ j = log ( ρ ) 2 ρ c i c i b i .
Using the gauge A ρ i = 0 , with letting A i depend only on ρ , one obtains
A 0 i = ( log ( ρ ) ) 2 ( c i ) 2 b i , A z i = 1 2 ( log ( ρ ) ) 2 ( c i ) 2 b i , A φ i = 1 2 ( log ( ρ ) ) 2 ( c i ) 2 b i ,
where the field b i Γ M ; so ( 3 , C ) P is constant on M = Σ × R .
As in the spherical symmetry case, one finds b i by using the spin charge Q i ( x ) , which is given by Equation (33). Since the system is static and the matter homogeneously located along the Z-axis, the spin charge Q i ( ρ , φ , z ) is given by Q i ( ρ , φ , z ) = Q 0 i δ ( ρ ) / 2 π ρ , which yields 0 0 2 π ρ d φ d ρ Q 0 i δ ( ρ ) / 2 π ρ = Q 0 i for each point of Z. Here, Q 0 i is the point charge located at each point of the Z-axis.
In order to get the same solution, as in Equations (39) and (40), the field b i is kept to be constant, and in Equation (33):
J k 0 ( x ) = Q k ( x ) = i ϵ k i j a K i a K j ( b i u j b j u i ) ,
Q i ( x ) = Q 0 i δ ( ρ ) / 2 π ρ is used.
In cylindrical symmetry, the functions K i are given by K i = c i log ( ρ ) , Equation (38), therefore,
Q k ( x ) = 2 i g ρ ρ ρ K i ρ K j ϵ k i j b i u j = 2 i c i c j ρ 2 ϵ k i j b i u j .
Therefore, in order to get Q i = Q 0 i δ ( ρ ) / 2 π ρ , one replaces 1 / ρ 2 with 1 / ρ ρ 1 ϵ , for some infinitesimal parameter ϵ 0 + , and chooses a solution for the field u j like
u j = b j + ϵ f b j + ϵ a j / ( 4 i π ) ,
for some scalar function f on M that is needed for satisfying D u i = 0 , with a constant vector field a i Γ M ; so ( 3 , C ) P . By that, one obtains (for i , j k )
Q k ( r ) = 2 i c i c j ρ ρ 1 ϵ ϵ k i j b i b j + ϵ f b j + ϵ 4 i π a j = 1 2 π ρ ϵ ρ 1 ϵ ϵ k i j ( c i b i ) ( c j a j ) .
Comparing with Q i = Q 0 i δ ( ρ ) / 2 π ρ , one finds: ϵ k i j ( c i b i ) ( c j a j ) = Q 0 k = const . , and, by imposing ( c i a i ) ( c i a i ) = 1 with ( c i b i ) ( c i a i ) = 0 :
( c b ) 2 = ( c i b i ) ( c i b i ) = Q 0 k Q 0 k ,
thus, choosing c i b i = e i Q 0 k Q 0 k for e i = 1 , so the constant field c i b i is determined by Q 0 k Q 0 k with free S O ( 3 , C ) rotation.
By that (for ρ > 0 , ϵ 0 + ), one obtains the same solutions as in Equations (39) and (40), but with c i b i = e i Q 0 k Q 0 k for e i = 1 and i = 1 3 c i = 0 . Since c i b i is a finite value, it is not sufficient to let the constants c i take arbitrary values, so chosen to be ( c i ) = ( 1 , 1 , 2 ) .

5. Conclusions

In this paper, the BF theory has been studied including matter by redefinition of the 2-form B i as B i + ε i j k ψ j B k , so that one can get D B i = 0 , with D being a covariant derivative and Latin letters i , j , = 1 , 2 , 3 , in the case of non-zero spin current of matter fields. The new field ψ i is defined using the spin current vector, J ν i = ( D μ ψ k ) ε k i j Σ μ ν j , Σ μ ν j being the frame field and Greek letters μ , ν , = 0 , 1 , 2 , 3 , are the space-time tangent indices. It is shown that one can solve the BF equations by using only the spin current of matter, that is, it is enough to solve the equations δ S / δ A i = 0 (for the action variation with the spin connection A i ), D B i = 0 and J ν i = ( D μ ψ k ) ε k i j Σ μ ν j without using a gravitational metric on the M 4-manifold and without a need to solve the equation δ S / δ B i = 0 , which includes the Lagrangian multiplier ϕ i j (a non-physical variable), so that one gets ϕ i j by using the solutions in δ S / δ B i = 0 . It is found that to obtain the solutions of BF theory, it is enough to use (find) the field ψ i and the Killing vector v (satisfies v a D a ψ i = 0 ) in Euclidean coordinates, where it is convenient to describe the spin currents and their lines in Euclidean coordinates with no need to describe them in curved coordinates. Furthermore, it is possible to obtain the solutions of BF theory using only the charges J 0 i 0 when they are given as functions on M in the static case (discussion below Equation (33)). It is shown that the singularities appear in solution of ψ i , that is related to the idea that the spin current J μ i is the source for ψ i , therefore, ψ i has singularities on the line of that spin current, and the singularities appear by that and not by using a gravitational metric. It is found that the solutions of BF theory equations depend on the symmetry of the system and every two solutions of ( Σ a b i , Σ 0 a i ) and ( Σ a b i , Σ 0 a i ) determine a metric (Remark 4), and those solutions are able to be written as e I e J using 1-forms e I , where I , J , = 0 , 1 , 2 , 3 are the Lorentz indices. Finally, the solutions of BF theory are applied to spherically and cylindrically symmetric systems in static case of matter.

Funding

This research received no external funding.

Acknowledgments

I am grateful to Yoshihiro Fukumoto at Ritsumeikan University for many useful notes and advice on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Verifying the General Relativity Equations in the Vacuum Using BF Equations

Satisfying of the equation of motion of classical general relativity (GR) in the vacuum is tested using the solutions in BF theory; the equations do not include the non-physical variable ϕ . The GR equations in the vacuum are [5]:
D a E i a = 0 , E i a F a b i = 0 , C = ε i j k E a i E b j F a b k = 0 .
Here, the gauge e a 0 = 0 is used with the metric g a b = δ i j e a i e b j . Note that D a E i a = 0 is satisfied when D B i = 0 is satisfied (see the discussion below Equation (10)).
Using Equation (20), in the vacuum χ i j = 0 , it reduces to F i ( A ) = χ i j Σ j , on the spacelike surface of constant time Σ t , it becomes F i ( A ) a b = χ i j Σ a b j . Using E a i = ε a b c Σ b c i , one obtains:
F i ( A ) a b = χ i j Σ a b j = χ i j ε a b c E j c .
Multiplying it by E i a and summing over the contracted indices, one gets:
E i a F i ( A ) a b = E i a χ i j ε a b c E j c = χ i j ε a b c E a i E j c = 0 ,
where the fact that the matrix χ i j is symmetric is used. Therefore, the second constraint (A1) is satisfied. Using F i ( A ) a b = χ i j ε a b c E j c in C, yields
C = ε i j k E a i E b j χ k i ε a b c E i c = χ k i ε i j k ε a b c E i a E j b E i c .
Then using E i a = e e i a , where e i a is the inverse of the gravitational field e a i and e = det ( e a i ) , one obtains:
C = e 3 χ k i ε i j k ε a b c e i a e j b e i c = e 3 χ k i ε i j k e 1 ε i j i = 2 e 2 Tr χ .
However, in the vacuum Tr χ = 0 (Equation (21)), thus, C = ε i j k E a i E b j F a b k = 0 is satisfied. By that one finds that the general relativity constraints are satisfied in the vacuum using the equations of motion of BF theory. In existence of matter, the first two equations are still satisfied, the third equation becomes C = ε i j k E a i E b j F a b k = 2 e 2 Tr χ = 2 e 2 Tr ( T i j ) .

Appendix B. Calculating δ S matter δ Σ i

Starting from
δ S matter δ B σ 1 ρ 1 constraint = δ S matter δ Σ σ 1 ρ 1 = 2 g δ S matter δ g μ 1 ν 1 g 2 δ g μ 1 ν 1 δ Σ σ 1 ρ 1 = T μ 1 ν 1 g 2 δ g μ 1 ν 1 δ Σ σ 1 ρ 1 ,
where T μ ν is energy-momentum tensor. Using
δ g μ 1 ν 1 = 1 g δ g g μ 1 ν 1 1 g g μ 1 ν 1 δ g ,
one obtains:
g δ g μ 1 ν 1 δ Σ σ 1 ρ 1 = δ g g μ 1 ν 1 δ Σ σ 1 ρ 1 g μ 1 ν 1 δ g δ Σ σ 1 ρ 1 = I 1 I 2 .
Using the Urbantke formula [12],
g g μ 1 ν 1 = ε μ ν ρ σ ε i j k Σ μ 1 μ i Σ ν ρ j Σ σ ν 1 k ,
one gets:
δ g g μ 1 ν 1 = δ ε μ ν ρ σ ε i j k Σ μ 1 μ i Σ ν ρ j Σ σ ν 1 k = ε μ ν ρ σ ε i j k δ Σ μ 1 μ i Σ ν ρ j Σ σ ν 1 k + ε μ ν ρ σ ε i j k Σ μ 1 μ i δ Σ ν ρ j Σ σ ν 1 k + ε μ ν ρ σ ε i j k Σ μ 1 μ i Σ ν ρ j δ Σ σ ν 1 k .
Calculating I 1 in (A3):
I 1 = δ δ Σ σ 1 ρ 1 g g μ 1 ν 1 = ε μ ν ρ σ ε i j k Σ ν ρ j Σ σ ν 1 k δ i δ μ 1 σ 1 δ μ ρ 1 + ε μ ν ρ σ ε i j k Σ μ 1 μ i Σ σ ν 1 k δ j δ ν σ 1 δ ρ ρ 1 + ε μ ν ρ σ ε i j k Σ μ 1 μ i Σ ν ρ j δ k δ σ σ 1 δ ν 1 ρ 1 = ε ρ 1 ν ρ σ ε j k Σ ν ρ j Σ σ ν 1 k δ μ 1 σ 1 + ε μ σ 1 ρ 1 σ ε i k Σ μ 1 μ i Σ σ ν 1 k + ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 = ε ρ 1 ν ρ σ ε i j Σ ν ρ i Σ σ ν 1 j δ μ 1 σ 1 + ε μ σ 1 ρ 1 σ ε i j Σ μ 1 μ i Σ σ ν 1 j + ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 = ε ρ 1 ν ρ σ ε i j Σ ν ρ i Σ σ ν 1 j δ μ 1 σ 1 ε μ σ 1 ρ 1 σ ε i j Σ μ 1 μ i Σ σ ν 1 j + ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 .
To calculate I 2 , one uses
g = i 6 ε μ ν ρ σ δ i j Σ μ ν i Σ ρ σ j ,
hence,
I 2 = g μ 1 ν 1 δ δ Σ σ 1 ρ 1 g = i 6 g μ 1 ν 1 ε μ ν ρ σ δ i j δ i δ μ σ 1 δ ν ρ 1 Σ ρ σ j + i 6 g μ 1 ν 1 ε μ ν ρ σ δ i j Σ μ ν i δ j δ ρ σ 1 δ σ ρ 1 = i 6 g μ 1 ν 1 ε σ 1 ρ 1 ρ σ δ i Σ ρ σ i + i 6 g μ 1 ν 1 ε μ ν σ 1 ρ 1 δ i Σ μ ν i = i 3 g μ 1 ν 1 ε μ ν σ 1 ρ 1 δ i Σ μ ν i .
Therefore,
g δ g μ 1 ν 1 δ Σ σ 1 ρ 1 = I 1 I 2 = ε ρ 1 ν ρ σ ε i j Σ ν ρ i Σ σ ν 1 j δ μ 1 σ 1 ε μ σ 1 ρ 1 σ ε i j Σ μ 1 μ i Σ σ ν 1 j + ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 i 3 g μ 1 ν 1 ε μ ν σ 1 ρ 1 δ i Σ μ ν i .
By this, Equation (A2) reads:
δ S matter δ B σ 1 ρ 1 constraint = 1 2 T μ 1 ν 1 I 1 I 2 = 1 2 T μ 1 ν 1 ε μ σ 1 ρ 1 σ ε i j Σ μ 1 μ i Σ σ ν 1 j 1 2 T μ 1 ν 1 ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 1 2 T μ 1 ν 1 ε ρ 1 ν ρ σ ε i j Σ ν ρ i Σ σ ν 1 j δ μ 1 σ 1 + i 6 T μ 1 ν 1 g μ 1 ν 1 ε μ ν σ 1 ρ 1 δ i Σ μ ν i = I 1 + I 2 + I 3 + I 4 .
Using T μ 1 ν 1 = T I J e I μ 1 e J ν 1 , the first term of (A6) reads:
2 I 1 = T μ 1 ν 1 ε μ ν σ 1 ρ 1 ε i j Σ μ 1 μ i Σ ν ν 1 j = T I J e I μ 1 e J ν 1 ε μ ν σ 1 ρ 1 ε i j P K L i P K 1 L 1 j e μ 1 K e μ L e ν K 1 e ν 1 L 1 = T I J δ I K δ J L 1 ε μ ν σ 1 ρ 1 ε i j P K L i P K 1 L 1 j e μ L e ν K 1 = T I J ε μ ν σ 1 ρ 1 ε i j P I L i P K 1 J j e μ L e ν K 1 = T I J ε μ ν σ 1 ρ 1 ε i j P I L i P K 1 J j P n L K 1 Σ μ ν n + P ¯ n L K 1 Σ ¯ μ ν n ,
where
e [ μ L e ν ] K 1 = P n L K 1 Σ μ ν n + P ¯ n L K 1 Σ ¯ μ ν n
is used.
One uses the following property of the self-dual projection:
1 2 ε i j P I J i P K L j = 1 4 η J K P I L η J L P I K 1 4 I J ,
which can be easily checked when I = 0 and J , K , L are spatial indices, and when I = K = 0 and J , L are spatial indices, so the Lorentz invariance asserts that this property is also satisfied when I , J , K , L are all spatial indices. By using this property, one obtains:
2 I 1 = T I J ε μ ν σ 1 ρ 1 1 2 η L K 1 P I J η L J P I K 1 η I K 1 P L J + η I J P L K 1 e μ L e ν K 1 = T I J ε μ ν σ 1 ρ 1 1 2 P I J g μ ν η L J P I K 1 e μ L e ν K 1 η I K 1 P L J e μ L e ν K 1 + η I J Σ μ ν = T I J ε μ ν σ 1 ρ 1 1 2 η L J P I K 1 e μ L e ν K 1 η I K 1 P L J e μ L e ν K 1 + η I J Σ μ ν = 1 2 T I J ε μ ν σ 1 ρ 1 η L J P I K 1 + η I K 1 P L J e [ μ L e ν ] K 1 1 2 η I J T I J ε μ ν σ 1 ρ 1 Σ μ ν ,
hence,
2 I 1 = 1 2 T I J ε μ ν σ 1 ρ 1 η L J P I K 1 + η I K 1 P L J P i L K 1 Σ μ ν i + P ¯ i L K 1 Σ ¯ μ ν i 1 2 η I J T I J ε μ ν σ 1 ρ 1 Σ μ ν .
Finally,
I 1 = 1 4 T I J ε μ ν σ 1 ρ 1 η L J P I K 1 + η I K 1 P L J P i L K 1 Σ μ ν i + P ¯ i L K 1 Σ ¯ μ ν i 1 4 η I J T I J ε μ ν σ 1 ρ 1 Σ μ ν .
Furthermore, the second term of (A6) becomes
2 I 2 = T μ 1 ν 1 ε μ ν ρ σ 1 ε i j Σ μ 1 μ i Σ ν ρ j δ ν 1 ρ 1 = T I J e I μ 1 e J ρ 1 ε μ ν ρ σ 1 ε i j P K L i P K 1 L 1 j e μ 1 K e μ L e ν K 1 e ρ L 1 = T I J δ I K e J ρ 1 ε μ ν ρ σ 1 ε i j P K L i P K 1 L 1 j e μ L e ν K 1 e ρ L 1 = T I J e J ρ 1 ε μ ν ρ σ 1 ε i j P I L i P K 1 L 1 j e μ L e ν K 1 e ρ L 1 = T I J ε μ ν ρ σ 1 1 2 η L K 1 P I L 1 η L L 1 P I K 1 η I K 1 P L L 1 + η I L 1 P L K 1 e J ρ 1 e μ L e ν K 1 e ρ L 1 = 1 2 T I J ε μ ν ρ σ 1 P I L 1 g μ ν e ρ L 1 P I K 1 g μ ρ e ν K 1 η I K 1 P L L 1 e μ L e ν K 1 e ρ L 1 + η I L 1 P L K 1 e μ L e ν K 1 e ρ L 1 e J ρ 1 = 1 2 T I J ε μ ν ρ σ 1 η I K 1 P L L 1 + η I L 1 P L K 1 e J ρ 1 e μ L e ν K 1 e ρ L 1 .
Using ε μ ν ρ σ 1 e μ L e ν K 1 e ρ L 1 = e ε L K 1 L 1 M e M σ 1 , where e is the determinant of ( e μ I ) , one obtains:
2 I 2 = 1 2 e T I J η I K 1 P L L 1 + η I L 1 P L K 1 e J ρ 1 ε L K 1 L 1 M e M σ 1 = 1 2 e T I J η I K 1 P L L 1 + η I L 1 P L K 1 ε L K 1 L 1 M e J [ ρ 1 e M σ 1 ] = 1 2 e T I J η I K 1 P L L 1 + η I L 1 P L K 1 ε L K 1 L 1 M e [ J ρ 1 e M ] σ 1 = 1 2 e T I J η I K 1 P L L 1 + η I L 1 P L K 1 ε L K 1 L 1 M P i J M Σ i ρ 1 σ 1 + P ¯ i J M Σ ¯ i ρ 1 σ 1 .
One uses the self-dual property,
P I J i ε I J K L = 2 i P i K L ,
to get
2 I 2 = i e T I J η I K 1 P K 1 M + η I L 1 P L 1 M P i J M Σ i ρ 1 σ 1 + P ¯ i J M Σ ¯ i ρ 1 σ 1 .
Furthermore, using
1 2 ! e 1 ε μ ν ρ σ Σ ρ σ i = Σ i μ ν = i Σ i μ ν = i Σ μ ν i , e = g ,
one obtains:
I 2 = 1 4 T I J η I K 1 P K 1 M + η I L 1 P L 1 M P i J M ε ρ 1 σ 1 ρ σ Σ ρ σ i + P ¯ i J M ε ρ 1 σ 1 ρ σ Σ ¯ ρ σ i .
One gets the third term of (A6) by the replacing σ 1 ρ 1 in I 2 , and, with reversing its sign:
I 3 = + 1 4 T I J η I K 1 P K 1 M + η I L 1 P L 1 M P i J M ε σ 1 ρ 1 ρ σ Σ ρ σ i + P ¯ i J M ε σ 1 ρ 1 ρ σ Σ ¯ ρ σ i = I 2 ,
therefore,
I 2 + I 3 = + 1 2 T I J η I K 1 P K 1 M + η I L 1 P L 1 M P i J M ε σ 1 ρ 1 ρ σ Σ ρ σ i + P ¯ i J M ε σ 1 ρ 1 ρ σ Σ ¯ ρ σ i .
The fourth term of (A6) is
I 4 = i 6 T μ 1 ν 1 g μ 1 ν 1 ε μ ν σ 1 ρ 1 Σ μ ν .
Using T μ ν g μ ν = T I J η I J = T , Equation (A6) reads:
δ S matter δ B σ 1 ρ 1 constraint = I 1 + I 2 + I 3 + I 4 = 1 4 T I J η L J P I K + η I K P L J P i L K ε μ ν σ 1 ρ 1 Σ μ ν i + P ¯ i L K ε μ ν σ 1 ρ 1 Σ ¯ μ ν i , 1 4 T ε μ ν σ 1 ρ 1 Σ μ ν + 1 2 T I J η I K P K M + η I L P L M P i J M ε σ 1 ρ 1 ρ σ Σ ρ σ i + P ¯ i J M ε σ 1 ρ 1 ρ σ Σ ¯ ρ σ i + i 6 T ε μ ν σ 1 ρ 1 Σ μ ν .
For brevity, this form is written as
δ S matter δ B σ 1 ρ 1 constraint = ε μ ν σ 1 ρ 1 T i Σ μ ν i + ε μ ν σ 1 ρ 1 ξ i Σ ¯ μ ν i ,
or
δ S matter δ B i constraint = δ S m δ Σ i = T i j Σ j + ξ i j Σ ¯ j ,
with the complex matrices T and ξ given by
T i j = 1 4 T I J η L J P i I K + η I K P i L J P j L K 1 2 T I J η I K P i K M + η I L P i L M P j J M + T i 6 1 4 δ i j = 1 2 P i I K P j J K T I J + i 6 1 4 T δ i j ,
and
ξ i j = 1 4 T I J η L J P i I K + η I K P i L J P ¯ j L K + 1 2 T I J η I K P i K M + η I L P i L M P ¯ j J M = 3 2 P i I K P ¯ j J K T I J .
The self-dual projection matrices P I J i are given in Equation (1),
P I J i = 1 2 ε i j k , for I = i , J = j , and P 0 j i = P j 0 i = i 2 δ j i , for I = 0 , J = j 0 .
and P ¯ I J i are their complex conjugate. One finds:
P i I K P j J K T I J = P i I 0 P j J 0 + P i I P j J T I J = P i 0 P j m 0 T m + P i 0 P j 0 T 00 + P i m P j 0 T m 0 + P i 0 P j m T 0 m + P i n P j m T n m = P i 0 P j m 0 T m + P i 0 P j 0 T 00 + P i m P j 0 T m 0 + P i 0 P j m T 0 m + P i n P j m T n m = i 2 δ i i 2 δ j m T m + i 2 δ i i 2 δ j T 00 + 1 2 ε i m i 2 δ j T m 0 + i 2 δ i 1 2 ε j m T 0 m + 1 2 ε i n 1 2 ε j m T n m = 1 4 T i j 1 4 δ i j T 00 + i 4 ε i m j T m 0 + i 4 ε j m i T 0 m + 1 4 δ i j δ n m δ i m δ n j T n m = 1 4 T i j 1 4 δ i j T 00 + i 4 ε i m j T 0 m + i 4 ε j m i T 0 m + 1 4 δ i j δ n m T n m 1 4 T j i = 1 4 δ i j T 00 + 1 4 δ i j δ n m T n m = 1 4 δ i j η 00 T 00 + 1 4 δ i j η n m T n m = 1 4 δ i j η I J T I J = 1 4 δ i j T .
Using this in Equation (A9), one gets:
T i j = 1 8 δ i j T + i 6 1 4 T δ i j = i 6 3 8 T δ i j .

Appendix C. Calculating ψi

Here, it is verifired that the (0,2)-tensor field Σ μ ν i defined in
i Σ μ ν i = 1 2 ! e 1 ε μ ν ρ σ Σ ρ σ i , e = det ( e μ I ) ,
is inversion of the 2-form Σ μ ν i , that is Σ μ ν i Σ j μ ν = δ j i . Multiplying by Σ μ ν j and summing up over contracted indices, one gets:
i Σ μ ν j Σ μ ν i = 1 2 e 1 ε μ ν ρ σ Σ μ ν j Σ ρ σ i .
Then, using Σ μ ν i = P I J i e μ I e ν J , one finds:
i Σ μ ν j Σ μ ν i = 1 2 e 1 ε μ ν ρ σ P j I J P K L i e μ I e ν J e ρ K e σ L .
Using ε μ ν ρ σ e μ I e ν J e ρ K e σ L = e ε I J K L , one gets:
i Σ μ ν j Σ μ ν i = 1 2 e 1 P j I J P K L i e ε I J K L = 1 2 P j I J P K L i ε I J K L .
Then, one uses the self-dual projection property,
P I J i ε I J K L = 2 i P i K L ,
to obtain
i Σ μ ν j Σ μ ν i = 1 2 P j I J ( 2 i P i I J ) = P j I J ( i P i I J ) .
Therefore,
Σ μ ν j Σ μ ν i = P j I J P i I J = δ j i ,
where the self-dual projection property, P j I J P i I J = δ j i is used. The sum is over the contracted indices.
Now, let us calculate the vector field ψ i , which is given in J k μ = 2 i ε k i j D ν ψ i Σ μ ν j , therefore, it is related to the current J μ i , or, in other words, the current is source for ψ i . If J k μ = 2 i ε k i j D ν ψ i Σ μ ν j is multiplied by Σ μ ρ k and summed up over contracted indices, one obtains:
Σ μ ρ k J k μ = 2 i ε k i j D ν ψ i Σ μ ρ k Σ μ ν j = 2 i ε k i j D ν ψ i P I J k e μ I e ρ J P K L j e μ K e ν L = 2 i ε k i j P I J k P K L j D ν ψ i η I K e ρ J e ν L = 2 i η I K ε i k j P I J k P K L j D ν ψ i e ρ J e ν L .
Using Equation (A7):
ε i k j P I J k P K L j = 1 2 η J K P i I L η J L P i I K + 1 2 η I K P i J L η I L P i J K .
Multiplying by η I K and summing over I , K , one gets:
η I K ε i k j P I J k P K L j = 1 2 P i J L 0 + 1 2 4 P i J L P i J L = P i J L .
Using this formula in (A12), one finds:
Σ μ ρ k J k μ = P I J k e μ I e ρ J J k μ = 2 i P i J L D ν ψ i e ρ J e ν L .
Multiplying by e K ρ (inverse of e ρ I ), one gets P I K k e μ I J k μ = 2 i P i K L D ν ψ i e ν L . From this equation, one can see that J μ i = 2 i g μ ν D ν ψ i + v μ i , for some vector v μ i satisfying v μ i e μ I = 0 . However, D μ Σ μ ν i = 0 , so D ν e μ I = 0 , therefore, e μ I D ν v μ i = 0 , and thus, one can choose D ν v μ i = 0 which allows us to calculate v μ i and then, to calculate ψ i in terms of J μ i and A μ i . The equation D μ J μ i = 0 implies D μ D μ ψ i = D 2 ψ i = 0 . Let us note that a specific metric g μ ν defined by e μ I = g μ ν e ν I is used, but, according to the Remark 4, there is a metric for every two solutions of ( Σ a b i , Σ 0 a i ) and ( Σ a b i , Σ 0 a i ) , and also the equation D 2 ψ i = 0 is invariant under any coordinates transformation x x , ψ i ( x ) = ψ i ( x ) , so that the same equation D 2 ψ i = 0 is valid for other metric g μ ν .

References

  1. Freidel, L.; Speziale, S. On the relations between gravity and BF theories. SIGMA 2012, 8, 032. [Google Scholar] [CrossRef] [Green Version]
  2. Baez, J.C. An introduction to spin foam models of quantum gravity and BFtheory. In Geometry and Quantum Physics; Gausterer, H., Pittner, L., Grosse, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2000; pp. 25–93. [Google Scholar] [CrossRef]
  3. Cattaneo, A.S.; Cotta-Ramusino, P.; Froehlich, J.; Martellini, M. Topological BF theories in 3 and 4 dimensions. J. Math. Phys. 1995, 36, 6137–6160. [Google Scholar] [CrossRef] [Green Version]
  4. Celada, M.; Gonzáez, D.; Montesinos, M. Plebanski-like action for general relativity and anti-self-dual gravity. Phys. Rev. D 2016, 93, 104058. [Google Scholar] [CrossRef] [Green Version]
  5. Rovelli, C. Quantum Gravity; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar] [CrossRef]
  6. Herfray, Y. Pure connection formulation. Twistors and the shase for a twistor action for general relativity. J. Math. Phys. 2017, 58, 112505. [Google Scholar] [CrossRef]
  7. Forger, M.; Römer, H. Currents and the energy-momentum tensor in classical field theory: A fresh look at an old problem. Ann. Phys. 2004, 309, 306–389. [Google Scholar] [CrossRef] [Green Version]
  8. Yoneda, G.; Shinkai, H. Constraints and reality conditions in the Ashtekar formulation of general relativity. Class. Quant. Grav. 1996, 13, 783–790. [Google Scholar] [CrossRef] [Green Version]
  9. Rovelli, C. Ashtekar formulation of general relativity and loop-space nonperturbative quantum gravity: A report. Class. Quant. Grav. 1991, 8, 1613. [Google Scholar] [CrossRef]
  10. Smolin, L. An invitation to loop quantum gravity. In Quantum Theory Symmetries; Argyres, P.C., Hodges, T.J., Mansouri, F., Scanio, J.J., Suranyi, P., Wijewardhana, L.C.R., Eds.; World Scientific: Singapore, 2004; pp. 655–682. [Google Scholar] [CrossRef] [Green Version]
  11. Freidel, L.; Krasnov, K.; Puzio, R. BF description of higher-dimensional gravity theories. Adv. Theor. Math. Phys. 1999, 3, 1289–1324. [Google Scholar] [CrossRef] [Green Version]
  12. Bengtsson, I. 2-form geometry and the ’t Hooft-Plebanski action. Class. Quant. Grav. 1995, 12, 1581. [Google Scholar] [CrossRef] [Green Version]
  13. Tennie, F.; Wohlfarth, M.N.R. Consistent matter couplings for Plebanski gravity. Phys. Rev. D 2010, 28, 104052. [Google Scholar] [CrossRef] [Green Version]
  14. Bennett, D.L.; Laperashvili, L.V.; Nielsen, H.B.; Tureanu, A. Gravity and mirror gravity in Plebanski formulation. Int. J. Mod. Phys. A 2013, 28, 1350035. [Google Scholar] [CrossRef] [Green Version]
  15. Krasnov, K. Plebański formulation of general relativity: A practical introduction. Gen. Relativ. Grav. 2011, 43, 1–15. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almatwi, M. Spin Current in BF Theory. Physics 2021, 3, 427-448. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020029

AMA Style

Almatwi M. Spin Current in BF Theory. Physics. 2021; 3(2):427-448. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020029

Chicago/Turabian Style

Almatwi, Malik. 2021. "Spin Current in BF Theory" Physics 3, no. 2: 427-448. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020029

Article Metrics

Back to TopTop