Physics
Maple provides a state-of-the-art environment for algebraic computations in physics, with emphasis on ensuring the computational experience is as natural as possible. The theme of the Physics project for Maple 18 has been the consolidation and integration of the Physics package with the rest of the Maple library, making it even easier to combine standard Maple commands and techniques with Physics-specific computations. With more than 500 enhancements throughout the entire package to increase robustness and versatility, an extension of its typesetting capabilities to support even more standard notation, as well 17 new Physics:-Library commands to support further explorations and extensions, Maple 18 extends the range of physics-related algebraic formulations that can be done in a natural way inside Maple. The impact of these changes is across the board, from vector analysis to quantum mechanics, relativity and field theory.
As part of its commitment to providing the best possible environment for algebraic computations in physics, Maplesoft has launched a Maple Physics: Research and Development web site, where users can download research versions, ask questions, and provide feedback. The results from this accelerated exchange with people around the world have been incorporated into the Physics package in Maple 18.
Simplify
4-Vectors, Substituting Tensors
Functional Differentiation
More Metrics in the Database of Solutions to Einstein's Equations
Commutators, AntiCommutators
Expand and Combine
New Enhanced Modes in Physics Setup
Dagger
Vectors Package
Library
Miscellaneous
Simplification is perhaps the most common operation performed in a computer algebra system. In Physics, this typically entails simplifying tensorial expressions, or expressions involving noncommutative operators that satisfy certain commutator/anticommutator rules, or sums and integrals involving quantum operators and Dirac delta functions in the summands and integrands. Relevant enhancements were introduced in Maple 18 for all these cases.
Examples
restart; withPhysics: Setupmathematicalnotation = true;
mathematicalnotation=true
Simplification of sums when the summand is linear in KroneckerDeltas:
Sumsqrtn+1 KroneckerDeltam,n+1+sqrtn KroneckerDeltam,n−1,n = 0 .. infinity
∑n=0∞⁡n+1⁢δm,n+1+n⁢δm,n−1
m+m+1
Simplification of tensorial expressions. To facilitate typing, set the spacetime indices to be lowercaselatin:
Setupspacetimeindices = lowercaselatin
spacetimeindices=lowercaselatin
Define a tensor Fa:
DefineFa
Defined objects with tensor properties
γμ,Fa,σμ,∂μ,gμ,ν,δμ,ν,εα,β,μ,ν
The following tensorial expression,
−1010921950548⁢Ff2⁢Fd2⁢Fj2⁢Fa⁢Fb⁢Fc−2023081095744⁢Ff2⁢Fd2⁢Fa⁢Fb⁢Fc+1701⁢Fc⁢Fe⁢Ff⁢ge,f−81⁢Fe⁢Ff2⁢δc,e−1350⁢Fe⁢δc,f⁢δe,f⁢81⁢Fd⁢Fi⁢Fj⁢gi,j−162⁢Fj⁢δd,i⁢δi,j+108⁢Fj⁢δd,j⁢2511⁢Fa⁢Fh⁢gd,h+27⁢Fd⁢Fh⁢ga,h+324⁢Fd⁢Fh⁢δa,h⁢432⁢Fb⁢Fn2+27⁢Fb⁢gl,n2+6156⁢Fn⁢δb,n−81⁢Fa⁢Fd⁢Fk2+81⁢Fd⁢Fk⁢δa,j⁢δj,k+8100⁢Fj⁢Fk⁢ga,k⁢δd,j⁢81⁢Fb⁢Fd⁢Fi2−5832⁢Fb⁢Fd⁢δh,i2+135⁢Fi2⁢δb,d⁢−642978⁢Fe2⁢Ff⁢gc,f+8748⁢Ff⁢gc,f+128755884390192⁢Fe2⁢Fh2⁢Fa⁢Fb⁢Fc−25470904248⁢Fh2⁢Fn2⁢Fi2⁢Fa⁢Fb⁢Fc−148090286178⁢Fe2⁢Fn2⁢Fh2⁢Fi2⁢Fa⁢Fb⁢Fc−−5022⁢Fc⁢Fi⁢gh,i⁢δd,h−1458⁢Fd⁢Fh⁢gc,i⁢δh,i+81⁢Fh⁢Fi⁢δc,i⁢δd,h⁢−5994⁢Fe⁢Fg⁢gb,e⁢δa,g−2700⁢Fa⁢Fg⁢δb,g−9396⁢Fg2⁢ga,b⁢−7695⁢Fd⁢Fk⁢Fm⁢δk,m−30780⁢Fm⁢δd,k⁢δk,m+500455863936⁢Fd2⁢Fa⁢Fb⁢Fc−1850647623120⁢ga,b⁢Fh2⁢Fg2⁢Fc−18366600960⁢Fe2⁢Fn2⁢Fi2⁢Fd2⁢Fa⁢Fb⁢Fc−355087618560⁢Fe2⁢Fd2⁢Fa⁢Fb⁢Fc−5179618847700⁢Fh2⁢Fa⁢Fb⁢Fc−290804515200⁢Fe2⁢Fd2⁢Fj2⁢Fa⁢Fb⁢Fc−462661905780⁢ga,b⁢Fh2⁢Fg2⁢Fk2⁢Fc−1378008798300⁢Fh2⁢Fn2⁢Fa⁢Fb⁢Fc
−1010921950548⁢Fd⁢F⁢d⁢d⁢Ff⁢F⁢f⁢f⁢Fj⁢F⁢j⁢j⁢Fa⁢Fb⁢Fc−2023081095744⁢Fd⁢F⁢d⁢d⁢Ff⁢F⁢f⁢f⁢Fa⁢Fb⁢Fc+1701⁢ge,f⁢Fc⁢F⁢e⁢e⁢F⁢f⁢f−81⁢Ff⁢F⁢f⁢f⁢Fe⁢δcece−1350⁢Fe⁢δc,f⁢δ⁢e,f⁢e,f⁢81⁢gi,j⁢Fd⁢F⁢i⁢i⁢F⁢j⁢j−162⁢Fj⁢δd,i⁢δ⁢i,j⁢i,j+108⁢Fj⁢δdjdj⁢2511⁢Fa⁢F⁢h⁢h⁢g⁢dh⁢dh+324⁢Fh⁢F⁢d⁢d⁢δahah+27⁢F⁢d⁢d⁢F⁢h⁢h⁢ga,h⁢432⁢Fn⁢F⁢n⁢n⁢Fb+27⁢gl,n⁢g⁢l,n⁢l,n⁢Fb+6156⁢Fn⁢δbnbn−81⁢Fk⁢F⁢k⁢k⁢Fa⁢Fd+81⁢Fd⁢Fk⁢δa,j⁢δ⁢j,k⁢j,k+8100⁢ga,k⁢Fj⁢F⁢k⁢k⁢δdjdj⁢81⁢Fb⁢Fi⁢F⁢d⁢d⁢F⁢i⁢i−5832⁢Fb⁢F⁢d⁢d⁢δh,i⁢δ⁢h,i⁢h,i+135⁢Fi⁢F⁢i⁢i⁢δbdbd⁢−642978⁢Fe⁢F⁢e⁢e⁢gc,f⁢F⁢f⁢f+8748⁢gc,f⁢F⁢f⁢f+128755884390192⁢Fe⁢F⁢e⁢e⁢Fh⁢F⁢h⁢h⁢Fa⁢Fb⁢Fc−25470904248⁢Fh⁢F⁢h⁢h⁢Fi⁢F⁢i⁢i⁢Fn⁢F⁢n⁢n⁢Fa⁢Fb⁢Fc−148090286178⁢Fe⁢F⁢e⁢e⁢Fh⁢F⁢h⁢h⁢Fi⁢F⁢i⁢i⁢Fn⁢F⁢n⁢n⁢Fa⁢Fb⁢Fc−−5022⁢gh,i⁢Fc⁢F⁢i⁢i⁢δdhdh−1458⁢gc,i⁢Fd⁢Fh⁢δ⁢h,i⁢h,i+81⁢Fh⁢Fi⁢δcici⁢δdhdh⁢−5994⁢gb,e⁢F⁢e⁢e⁢Fg⁢δagag−2700⁢Fa⁢Fg⁢δbgbg−9396⁢Fg⁢F⁢g⁢g⁢ga,b⁢−7695⁢Fk⁢Fm⁢F⁢d⁢d⁢δ⁢k,m⁢k,m−30780⁢Fm⁢δ⁢dk⁢dk⁢δ⁢k,m⁢k,m+500455863936⁢Fd⁢F⁢d⁢d⁢Fa⁢Fb⁢Fc−1850647623120⁢Fg⁢F⁢g⁢g⁢Fh⁢F⁢h⁢h⁢ga,b⁢Fc−18366600960⁢Fd⁢F⁢d⁢d⁢Fe⁢F⁢e⁢e⁢Fi⁢F⁢i⁢i⁢Fn⁢F⁢n⁢n⁢Fa⁢Fb⁢Fc−355087618560⁢Fd⁢F⁢d⁢d⁢Fe⁢F⁢e⁢e⁢Fa⁢Fb⁢Fc−5179618847700⁢Fh⁢F⁢h⁢h⁢Fa⁢Fb⁢Fc−290804515200⁢Fd⁢F⁢d⁢d⁢Fe⁢F⁢e⁢e⁢Fj⁢F⁢j⁢j⁢Fa⁢Fb⁢Fc−462661905780⁢Fg⁢F⁢g⁢g⁢Fh⁢F⁢h⁢h⁢Fk⁢F⁢k⁢k⁢ga,b⁢Fc−1378008798300⁢Fh⁢F⁢h⁢h⁢Fn⁢F⁢n⁢n⁢Fa⁢Fb⁢Fc
has various terms with contracted indices. In each term, {a,b,c} are free indices:
Check, all
The repeated indices per term are: ...,...,...; the free indices are: ...
d,f,j,d,f,d,e,f,h,i,j,l,n,d,e,f,h,i,j,k,e,h,h,i,n,e,h,i,n,d,e,g,h,i,k,m,d,g,h,d,e,i,n,d,e,h,d,e,j,g,h,k,h,n,a,b,c
Taking into account Einstein's sum rule for contracted (repeated) indices, the symmetry properties of gi,j and δi,j, this tensorial expression is equal to zero:
0
The simplification of integrals and sums involving quantum operators that satisfy algebra rules is now more powerful, both in the continuous and discrete case. Consider a field, ψ, and its expansion in terms in a basis of functions, φ using operators, a and a†, that satisfy:an,a⁡p†−=δ⁡n−p
Setupop = a, psi, algebrarules = %Commutatoran, %Daggerap = Diracn−p, %Commutatoran, ap = 0
* Partial match of 'op' against keyword 'quantumoperators'
algebrarules=a⁡n,a⁡p−=0,a⁡n,a⁡p†−=δ⁡n−p,quantumoperators=a,ψ
The expansion of terms ψ and ψ†is given by:
ψr=∫−∞∞φn,r⁢anⅆn
ψ⁡r=∫−∞∞φ⁡n,r⁢a⁡nⅆn
subsn=p,r=s,Dagger
ψ⁡s†=∫−∞∞φ⁡p,s&conjugate0;⁢a⁡p†ⅆp
The commutator ψr,ψ†s is equal to:
Commutator,
ψ⁡r,ψ⁡s†−=∫−∞∞φ⁡n,r⁢a⁡nⅆn,∫−∞∞φ⁡p,s&conjugate0;⁢a⁡p†ⅆp−
expand
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∫−∞∞φ⁡n,r⁢a⁡nⅆn⁢∫−∞∞φ⁡p,s&conjugate0;⁢a⁡p†ⅆp−∫−∞∞φ⁡p,s&conjugate0;⁢a⁡p†ⅆp⁢∫−∞∞φ⁡n,r⁢a⁡nⅆn
The products of integrals on the right-hand side can both be combined into double integrals, then recombined into a single integral and simplified taking into account the algebra rule stated: an,a⁡p†−= δn−p.
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∫−∞∞φ⁡p,s&conjugate0;⁢φ⁡p,rⅆp
The step involving only the combination of the integrals can now also be performed separately:
combine
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∫−∞∞∫−∞∞−φ⁡p,s&conjugate0;⁢φ⁡n,r⁢a⁡p†⁢a⁡n+φ⁡n,r⁢φ⁡p,s&conjugate0;⁢a⁡n⁢a⁡p†ⅆnⅆp
The extended capabilities in Simplify regarding integration also work in the discrete case, over sums. Redo the algebra rule now considering the same relations but in the discrete case.
Setupredo, quantumoperators = a, psi, algebrarules=%Commutatoran,%Daggerap=KroneckerDeltan,p,%Commutatoran,ap=0, spacetimeindices = greek
algebrarules=an,ap−=0,an,ap†−=δn,p,quantumoperators=a,ψ,spacetimeindices=greek
The following sum can now be simplified by combining the sums and taking into account the new (discrete) algebra rules, or just performing the combination step:
ψ⁡r=Sum⁡φn⁡r⁢an,n=−∞..∞
ψ⁡r=∑n=−∞∞⁡φn⁡r⁢an
ψ⁡s†=∑p=−∞∞⁡φp⁡s&conjugate0;⁢ap†
expandCommutator,
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∑n=−∞∞⁡φn⁡r⁢an⁢∑p=−∞∞⁡φp⁡s&conjugate0;⁢ap†−∑p=−∞∞⁡φp⁡s&conjugate0;⁢ap†⁢∑n=−∞∞⁡φn⁡r⁢an
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∑p=−∞∞⁡φp⁡s&conjugate0;⁢φp⁡r
ψ⁡r⁢ψ⁡s†−ψ⁡s†⁢ψ⁡r=∑p=−∞∞⁡∑n=−∞∞⁡−φp⁡s&conjugate0;⁢φn⁡r⁢ap†⁢an+φn⁡r⁢φp⁡s&conjugate0;⁢an⁢ap†
Improvements in the simplification of annihilation and the creation of fermionic operators, as well as the related occupation number operator:
Setupanticommutativeprefix = psi
anticommutativeprefix=_λ,ψ
am ≔ Annihilationpsi, notation = explicit
am≔a−ψ1
ap ≔ Creationpsi, notation = explicit
ap≔a+ψ1
The related occupation number operator:
N ≔ ap . am
N≔a+ψ1⁢a−ψ1
Consider the application of these fermionic operators to a related state vector:
Ketpsi, 1
ψ1
am .
ψ0
Increasing the occupation number,
ap .
In other words, powers of annihilation and creation fermionic operators are equal to zero:
am2
a−ψ12
Simplifyam2
Simplifyap2
The occupation number operator is also idempotent:
N N N = N
a+ψ1⁢a−ψ1⁢a+ψ1⁢a−ψ1⁢a+ψ1⁢a−ψ1=a+ψ1⁢a−ψ1
These expressions can now be simplified:
Simplify%
a+ψ1⁢a−ψ1=a+ψ1⁢a−ψ1
The simplification of vectorial expressions was also enhanced. For example:
withVectors :
B0→·v→&xB0→+B1→
B0→·v→×B0→+B1→
B0→·v→×B1→
In Maple 17, it is possible to define a tensor with a tensorial equation, where the tensor being defined is on the left-hand side. Then, on the right-hand side, you write either a tensorial expression with free and repeated indices, or a Matrix or Array with the components themselves. In Maple 18, you can also define a 4-Vector with a tensorial equation, where you indicate the vector's components on the right-hand side as a list.
One new Library routine specialized for tensor substitutions was added to the Maple library: SubstituteTensor, which substitutes the equation(s) Eqs into an expression, taking care of the free and repeated indices, such that: 1) equations in Eqs are interpreted as mappings having the free indices as parameters, and 2) repeated indices in Eqs do not clash with repeated indices in the expression. This new routine can also substitute algebraic sub-expressions of type product or sum within the expression, generalizing and unifying the functionality of the subs and algsubs commands for algebraic tensor expressions.
Define a contravariant 4-vector with components p__x, p__y, p__z,p__t:
Definep~mu = p__x, p__y, p__z,p__t
γμ,σμ,∂μ,gμ,ν,p⁢μ⁢μ,δμ,ν,εα,β,μ,ν
You can retrieve the components in different ways:
p~mu = Library:-TensorComponentsp~mu
p⁢μ⁢μ=p__x,p__y,p__z,p__t
pmu = Library:-TensorComponentspmu
pμ=−p__x,−p__y,−p__z,p__t
Or, indexing p with a contravariant or covariant integer value of the index:
p~1≠p1
p__x≠−p__x
You can compute with p⁢μ⁢μ algebraically; p[~mu] will return its components only when the index assumes integer values 0≤μ ≤4.
pmu pnu ep_mu,nu,alpha,beta
εα,β,μ,ν⁢p⁢μ⁢μ⁢p⁢ν⁢ν
pmu2
pμ⁢p⁢μ⁢μ
SumOverRepeatedIndices
p__t2−p__x2−p__y2−p__z2
Define some tensors for experimentation with the new Library:-SubstituteTensor command:
DefineA,B,F,G
A,B,F,G,γμ,σμ,∂μ,gμ,ν,p⁢μ⁢μ,δμ,ν,εα,β,μ,ν
A substitution equation:
Aμ=Gν,α⁢Aα⁢Fμ,ν
Aμ=Gν,α⁢A⁢α⁢α⁢Fμνμν
Substitute into Aν⁢A⁢ν⁢ν: the free indices of (48) are taken as parameters, repeated indices in the substitution equation do not repeat more than once in the result:
Library:-SubstituteTensor, Aν⁢A⁢ν⁢ν
Gβ,α⁢A⁢α⁢α⁢Fνβνβ⁢Gλ,κ⁢A⁢κ⁢κ⁢F⁢ν,λ⁢ν,λ
When the left-hand side of the substitution equation is a tensor function, the functionality is also taken as a parameter,
Library:-SubstituteTensorAμ⁡X=Bμ⁡X,Aν⁡Y
Bν⁡Y
SubstituteTensor can also substitute sub-expressions of type product or sum, similar to what algsubs does, so for example substitute:
Amu Bnu = Gmu,nu
Aμ⁢Bν=Gμ,ν
into,
Aalpha Fmu,nu Bbeta Arho Brho Gmu,nu
Aα⁢Fμ,ν⁢G⁢μ,ν⁢μ,ν⁢Bβ⁢Aρ⁢B⁢ρ⁢ρ
Library:-SubstituteTensor,
Gα,β⁢Gρρρρ⁢Fμ,ν⁢G⁢μ,ν⁢μ,ν
Check the free and repeated indices of this result and verify that the free indices of (52) are the same:
Check,all
μ,ν,ρ,α,β
Check,free
The free indices are: ...
α,β
The Physics:-Fundiff command for functional differentiation has been extended to handle all the complex components (abs, argument, conjugate, Im, Re, signum) and vectorial differential operators in order to compute field equations using variational principles when the field function enters the Lagrangian together with its conjugate. For an example illustrating the use of the new capabilities in the context of a more general problem, see the MaplePrimes™ post Quantum Mechanics using Computer Algebra.
restart; withPhysics:Setupmathematicalnotation = true
A function and its conjugate are considered independent from each other regarding functional differentiation:
%Fundiff = Fundiffconjugatefx, fy
δδ⁢f⁡y⁡f⁡x&conjugate0;=0
Fundiff can now compute functional derivatives of expressions involving vectorial differential operators and the corresponding conjugate functions.
withVectors:
%Fundiff = Fundiff%Gradientfx, fy
δδ⁢f⁡y⁡∇f⁡x=δ′⁡x−y⁢i∧
Set a system of coordinates for functional differentiation with many variables:
CoordinatesQ = X, Y, Z, T
Default differentiation variables for d_, D_ and dAlembertian are: Q=X,Y,Z,T
Systems of spacetime Coordinates are: Q=X,Y,Z,T
Q
The Action for a system:
S ≔ IntcNorm%GradientPhix,y,z,t2, x,y,z,t
S≔∫−∞∞∫−∞∞∫−∞∞∫−∞∞∇Φ⁡x,y,z,t2ⅆxⅆyⅆzⅆt
The equations of motion through functional differentiation:
%Fundiff = FundiffS, PhiX,Y,Z,T
δδ⁢Φ⁡Q⁡∫−∞∞∫−∞∞∫−∞∞∫−∞∞∇Φ⁡x,y,z,t2ⅆxⅆyⅆzⅆt=−ⅆ2ⅆZ2Φ⁡Q&conjugate0;−ⅆ2ⅆY2Φ⁡Q&conjugate0;−ⅆ2ⅆX2Φ⁡Q&conjugate0;
The Action for a complex scalar field with a Lagrangian quadratic in the derivatives: note that abs now automatically maps into the Norm of the vector.
S≔Intc−12absGradientΦx,y,z,t2, x,y,z,t
S≔∫−∞∞∫−∞∞∫−∞∞∫−∞∞−ⅆⅆxΦ⁡x,y,z,t&conjugate0;⁢∂∂xΦ⁡x,y,z,t2−ⅆⅆyΦ⁡x,y,z,t&conjugate0;⁢∂∂yΦ⁡x,y,z,t2−ⅆⅆzΦ⁡x,y,z,t&conjugate0;⁢∂∂zΦ⁡x,y,z,t2ⅆxⅆyⅆzⅆt
The corresponding field equation:
FundiffS, Φ⁡X,Y,Z,T = 0
ⅆ2ⅆZ2Φ⁡Q&conjugate0;2+ⅆ2ⅆY2Φ⁡Q&conjugate0;2+ⅆ2ⅆX2Φ⁡Q&conjugate0;2=0
A database of solutions to Einstein's equations was added to the Maple library in Maple 15 with a selection of metrics from "Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C.; and Herlt, E., Exact Solutions to Einstein's Field Equations" and "Hawking, Stephen; and Ellis, G. F. R., The Large Scale Structure of Space-Time". More metrics from these two books were added for Maple 16 and Maple 17. These metrics can be searched using the command DifferentialGeometry:-Library:-MetricSearch, or directly using g_ (the Physics command representing the spacetime metric that also sets the metric to your choice).
For Maple 18, fifty more metrics were added to the database from Chapter 28 of the aforementioned book entitled "Exact Solutions to Einstein's Field Equations".
It is now possible to list all the metrics of a chapter by indexing the metric command with the chapter's number, for example, entering g_["28"].
By default, the metric is a Minkowski type:
g_
You can query about metrics directly from the metric command g_
g_Bajer
____________________________________________________________
28,58.2,1=Authors=Bajer, Kowalezynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=All tensor components given with respect to the anholonomic frame,The coordinates xi and xi1 are complex conjugates
28,58.3,1=Authors=Bajer, Kowalezynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Stationary,Comments=Case 1 of 2,The coordinates xi and xi1 are complex conjugates, the parameters _Q0 and _Q01 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,58.3,2=Authors=Bajer, Kowalezynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 1 of 2,The coordinates xi and xi1 are complex conjugates, the parameters _Q0 and _Q01 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,58.4,1=Authors=Bajer, Kowalezynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
Warning, found more than one match for the keyword 'Bajer', as seen above. Please refine your 'keyword' or re-enter the metric 'g_[...]' with the list of three numbers identifying the metric, for example as in g_[[28, 58.2, 1]] or Setup(metric = [28, 58.2, 1])
When you identified the metric, you can set it directly from g_ (alternatively, you can do that using Setup):
g_28, 58.2, 1
Systems of spacetime Coordinates are: X=r,ξ,ξ1,u
Default differentiation variables for d_, D_ and dAlembertian are: X=r,ξ,ξ1,u
The Bajer, Kowalezynski (1985) metric in coordinates r,ξ,ξ1,u
Parameters: κ0,Q0,Q01,m0,b,Q
Comments: All tⅇnsor componⅇnts gⅈvⅇn wⅈth rⅇspⅇct to thⅇ anholonomⅈc framⅇ
New in Maple 18, you can now also list all the metrics of a chapter. For example, for the metrics of Chapter 28,
g_28
28,16,1=Authors=Robinson-Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,17,1=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence
28,21,1=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi
28,21,2=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi
28,21,3=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Static,Comments=The coordinate zeta is changed to xi
28,21,4=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi,This is _a special case of Kasner spacetime Stephani [13, 51,1], [13, 53,1]
28,21,5=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Static,Comments=The coordinate zeta is changed to xi
28,21,6=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi
28,21,7=Authors=Robinson (1975),Foster, Newman (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Static,Comments=The coordinate zeta is changed to xi
28,24,1=Authors=Collinson, French (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=Stephani claims this metric is static which is false since the orbit type is generically Riemannian,The assumption _b >0 and _c >0 is made so that the given base point is in the domain
28,25,1=Authors=Collinson, French (1967),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,1=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1],The metric is type D at points where r = 3*_m/(xi1+xi2) and type II on either side of this hypersurface. For convenience, it is assumed that 3*_m - r*(xi1 + xi2) > 0,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,2=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1].,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,3=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1].,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,41,1=Authors= Bartrum (1967),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=PureRadiation,RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,43,1=Authors=Robinson, Trautman (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=PureRadiation,RobinsonTrautman,Comments=h1(u) is the conjugate of h(u)
28,44,1=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 1 of 6. K = 1, Riemmannian Orbits,_Q1 is the conjugate of Q,The metric is defined for all r > 0. The restriction 2*r^2 - 4*_m*r +_kappa0*_Q*_Q1 < 0 gives Riemannian orbits for the isometry group
28,44,2=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Static,Comments=Case 2 of 6 K = 1, PseudoRiemmannian Orbits,_Q1 is the conjugate of Q,The metric is defined for all r > 0. The restriction 2*r^2 - 4*_m*r +_kappa0*_Q*_Q1 > 0 gives PseudoRiemannian orbits for the isometry group
28,44,3=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 3 of 6 K = 0, Riemmannian Orbits,_Q1 is the conjugate of Q,The metric is defined for all r > 0. The restriction _Q*_Q1*_kappa0-4*_m*r < 0 gives Riemannian orbits for the isometry group
28,44,4=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Static,Comments=Case 4 of 6 K = 0, Pseudo-Riemmannian Orbits,_Q1 is the conjugate of Q,The metric is defined for all r > 0. The restriction _Q*_Q1*_kappa0-4*_m*r > 0 gives Pseudo-Riemannian orbits for the isometry group
28,44,5=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case of 6. K = 1, Riemmannian Orbits,_Q1 is the conjugate of _Q,The metric is defined for all r > 0. The restriction _Q*_Q1*_kappa0 - 4*m*r- 2*r^2 < 0 gives Riemannian orbits for the isometry group
28,44,6=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Static,Comments=Case 6 of 6. K = 1, Pseudo-Riemmannian Orbits,_Q1 is the conjugate of _Q,The metric is defined for all r > 0. The restriction 2*r^2 - 4*_m*r +_kappa0*_Q*_Q1 > 0 gives Pseudo-Riemannian orbits for the isometry group
28,45,1=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Static,Comments=Case 1 of 2,Note that the metric only depends on the square of the function P. If the funtion P^(-2) is negative, then the solution is static.,The parameter _q determines _a duality rotation on the electromagnetic field.
28,45,2=Authors=Leroy (1976),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=,Comments=Case 2 of 2,Note that the metric only depends on the square of the function P. If the funtion P^(-2) is negative, then the solution is static.,The parameter _q determines _a duality rotation on the electromagnetic field.
28,46,1=PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 1 of 2
28,46,2=PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Stationary,Comments=Case 2 of 2
28,53,1=Authors=Bartrum (1967),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,PureRadiation,Stationary,Comments=Case 1 of 2. We take _m > 0 for simplicity. AlternativeNullTetrad1 is the standard Robinson Trautman null tetrad. The orbit type is pseudo-Riemannian at the given base point if _f1(0)^2*_f2(0)^2 > 2_m. If alpha = constant then the electromagnetic field is inheriting.
28,53,2=Authors=Bartrum (1967),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,PureRadiation,Comments=Case 2 of 2. We take _m > 0 for simplicity. AlternativeNullTetrad1 is the standard Robinson Trautman null tetrad. The orbit type is pseudo-Riemannian at the given base point if _f1(0)^2*_f2(0)^2 < 2_m. If alpha = constant then the electromagnetic field is inheriting.
28,55,1=PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Stationary,Comments=Case 1 of 2. We choose _A > 0 and u > 0 for simplicity. AlternativeNullTetrad1 is the standard Robinson Trautman null tetrad.
28,55,2=PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 2 of 2. We choose _A > 0 and u > 0 for simplicity. AlternativeNullTetrad1 is the standard Robinson Trautman null tetrad.
28,60,1=Authors=Kowalczynski (1978),Kowalczynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates, the parameters _Q0 and _Q01 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,61,1=Authors=Kowalczynski (1978),Kowalczynski (1985),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates, the parameters _Q0 and _Q01 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,64,1=Authors=Herlt and Stephani (1984),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,66,1=Authors=Herlt and Stephani (1984),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,67,1=Authors=Herlt and Stephani (1984),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,68,1=Authors=Herlt and Stephani (1984),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,72,1=PrimaryDescription=PureRadiation,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,73,1=Authors=Frolov and Khlebnikov (1975),PrimaryDescription=PureRadiation,SecondaryDescription=RobinsonTrautman,Comments=The coordinates xi and xi1 are complex conjugates
28,74,1=Authors=Frolov and Khlebnikov (1975),PrimaryDescription=PureRadiation,SecondaryDescription=RobinsonTrautman,Comments=With _m(u) = constant, the metric is Ricci flat and becomes 28.24 in Stephani.
28,56.1,1=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Stationary,Comments=Case 1 of 3, epsilon = 0, stationary,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.2,2=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 2 of 3, epsilon = 0,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.2,3=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=Case 3 of 3, epsilon = 1,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.3,1=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.4,1=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.5,1=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,56.6,1=Authors=Leroy, 1976,PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
Warning, found more than one match for the keyword '28', as seen above. Please refine your 'keyword' or re-enter the metric 'g_[...]' with the list of three numbers identifying the metric, for example as in g_[[28, 16, 1]] or Setup(metric = [28, 16, 1])
When computing with products of noncommutative operators, the results depend on the algebra of commutators and anticommutators that you previously set. Besides that, in Physics, various mathematical objects themselves satisfy specific commutation rules. You can query about these rules using the Library commands Commute and Anticommute. Previously existing functionality and enhancements in this area were refined and implemented in Maple 18. Among them:
Both Commutator and AntiCommutator now accept matrices as arguments.
The AntiCommutator of products of fermionic operators - for instance annihilation and creation operators - is now derived automatically from the intrinsic anticommutation rules they satisfy.
Commutators and Anticommutators of vectorial quantum operators A→,B→, are now implemented and expressed using the dot (scalar) product, as in A→,B→−=A→·B→−B→·A→
If two noncommutative operators a and S satisfy a†,S−=0 , then the commutator a,S†− is automatically taken equal to 0; if in addition S is Hermitian, then a,S−is also automatically taken equal to zero.
Commutator and AntiCommutator now also operate on matrices:
M__1≔ Matrix2,2,a,b,c,d
M__1≔abcd
M__2≔ Matrix2,2,alpha,beta,gamma,delta
M__2≔αβγδ
%Commutator = CommutatorM__1, M__2
abcd,αβγδ−=b⁢γ−c⁢βa⁢β−α⁢b+b⁢δ−β⁢d−γ⁢a+c⁢α−δ⁢c+d⁢γ−b⁢γ+c⁢β
Commutators of vectorial operators were implemented, expressed using the scalar (dot) product:
withVectors: Setupop = A_, B_
quantumoperators=A→,B→
Commutator = expand@CommutatorA_, B_
A→,B→−=A→·B→−B→·A→
Define 4 pairs of annihilation/creation operators for fermionic particles:
for j to 4 do apj ≔ Creationpsi, j, notation = explicit; amj ≔ Annihilationpsi, j, notation = explicit end do;
ap1≔a+ψ1
am1≔a−ψ1
ap2≔a+ψ2
am2≔a−ψ2
ap3≔a+ψ3
am3≔a−ψ3
ap4≔a+ψ4
am4≔a−ψ4
For these operators, the system knows about the anticommutator rules satisfied between any two of them, the algebra is set on the fly when you define them.
Setupalgebra
* Partial match of 'algebra' against keyword 'algebrarules'
algebrarules=a−ψ1,a+ψ1+=1,a−ψ2,a+ψ2+=1,a−ψ3,a+ψ3+=1,a−ψ4,a+ψ4+=1
Using that information, the system now also knows about the anticommutator of products of these fermionic operators. For example, on the left-hand side: inert, on the right-hand side: computed.
%AntiCommutator = AntiCommutatorap1,am1⁢am2⁢ap2
a+ψ1,a−ψ1⁢a−ψ2⁢a+ψ2+=a−ψ2⁢a+ψ2
This new functionality is automatically used when sorting products of non-commutative operators according to a specified ordering using the new Library:-SortProducts routine; let P be a product:
P ≔ ap2 am2 am1 ap1
P≔a+ψ2⁢a−ψ2⁢a−ψ1⁢a+ψ1
Rewrite this product using the following ordering: with the creation operators to the left, using the anticommutator relations between them.
ApAm ≔ seqapj, j = 1 .. 2, seqamj, j = 1 .. 2
ApAm≔a+ψ1,a+ψ2,a−ψ1,a−ψ2
Library:-SortProductsP, ApAm, useanticommutator
a+ψ1⁢a+ψ2⁢a−ψ1⁢a−ψ2+a+ψ2⁢a−ψ2
Changes in design: