Physics[dAlembertian] - d'Alembertian differential operator with respect to the spacetime coordinates
any mathematical expression or relation between expressions
(optional) list of differentiation variables with respect to which the d'Alembertian is to be computed
The dAlembertian command is a computational representation for the d'Alembert differential operator, displayed as a square box □. The definition used is
where the differential operator ∂__mu, represented in Maple by d_, is defined by
and where xμ represents the contravariant spacetime vector of a coordinate system defined using Coordinates or Setup.
In a galilean system (Euclidean or Minkowski), both xμ and the differential of the coordinates ∂xμ are vectors (tensors with 1 index), and so ∂μ is a vector representing
where xμ=gμ,ν⁢xν is the covariant spacetime vector, and so ∂μ⁡∂μ⁡A=∂μ⁡∂μ⁡A.
Unlike the galilean case, in a curvilinear system of coordinates, neither x__mu nor ∂μ are vectors and the formula above for ∂μ loses its meaning; instead, the convention used in the Physics package (it becomes the one above only in the galilean case) is
Hence in a non-galilean spacetime ∂μ⁡∂μ⁡A≠∂μ⁡∂μ⁡A, because the metric g⁢μ,ν depends on the coordinates and so the two factors in the right-hand-side above do not commute, and the formulas relating the second application of d_[mu] to the dAlembertian are
reflecting the nontensorial character in curvilinear coordinates of ∂μ and of the dAlembertian.
The %dAlembertian command is the inert form of dAlembertian; that is, it represents the same mathematical operation while displaying the operation unevaluated. To evaluate the operation, use the value command. To obtain the form of dAlembertian or its inert representation as a sum of diff constructions, use convert/diff.
As in the case of d_[mu], when only one argument is given to dAlembertian, the differentiation variables are the ones set by the Setup command, typically x1,x2,x3,x4 (x0 is automatically assigned to x4), represented by X (see Coordinates). To change these default differentiation variables, see Setup.
Regardless of the existence of default differentiation variables, you can always call dAlembertian with two arguments, where the second argument is a list with the differentiation variables you want; in this case, the list should have as many symbols as the spacetime dimension, which by default is 4 but can be set to any value by the Setup command.
Some automatic simplifications are carried out each time dAlembertian(A) is called, as follows:
- If A does not depend on the differentiation variables, then 0 is returned.
- If A is an unknown function (the rule for its derivative is unknown), a Dirac delta function, or a derivative, then the result is returned unevaluated, as dAlembertian(A).
- If A is of the form d_μ⁡B, then d_[mu](dAlembertian(B)) is returned.
- Otherwise, dAlembertian(A) is computed calling d_ as in ∂μ⁡∂μ⁡A.
In general, to accomplish differentiation, dAlembertian calls d_, which in turn makes calls to the Physics/diff command, which in turn uses the standard Maple diff command with appropriate arguments. In this way, any user-defined differentiation rule in the library or that you created, such as for a function foo of the form `diff/foo`, is automatically taken into account by dAlembertian.
Setup(mathematicalnotation = true);
Set the default differentiation variables for dAlembertian and d_, define a spacetime tensor function A, and use the enhanced display scheme of the differential equation packages.
Setup(differentiationvariables = X);
⁢Default differentiation variables for d_, D_ and dAlembertian are:⁢X=x1,x2,x3,x4
⁢Systems of spacetime coordinates are:⁢X=x1,x2,x3,x4
Defined objects with tensor properties
PDEtools[declare]((f, A, B, Phi)(X));
f⁡X⁢will now be displayed as⁢f
A⁡X⁢will now be displayed as⁢A
B⁡X⁢will now be displayed as⁢B
Phi⁡X⁢will now be displayed as⁢Φ
The dAlembertian operator is the double application of d_[mu], with the index contracted.
In the default 4 = 3+1 dimensional Minkowski spacetime, the form of dAlembertian as a sum of diff constructions is
The dAlembertian deals normally with derivatives expressed in terms of any of the Maple differential operators, including D, diff, and d_. It also distributes over sums and products.
d_[mu](f(X)*A[mu](X) + B[mu](X));
Because dAlembertian is a second order differential operator (dAlembertian(f) = d_[mu](d_[mu](f))), the differentiation rule, when applied to a product as in above, requires decomposing the dAlembertian operation into the sequence of d_[mu] operations. Note also the introduction of a dummy index alpha, which could be any spacetime index; these indices are chosen after checking that they are not assigned and not already present in the given expression at the time of introducing them.
Regardless of having set the default differentiation variables to X, you can call dAlembertian or d_ with other coordinates as differentiation variables. For example, set one more coordinate system:
⁢Systems of spacetime coordinates are:⁢X=x1,x2,x3,x4,Y=y1,y2,y3,y4
Note that the output above displays the differentiation variables y1,y2,y3,y4. This is so because they are not the "default differentiation variables;" if you set them to be y1,y2,y3,y4, then the omitted differentiation variables are y1,y2,y3,y4, and the other differentiation variables (namely, x1,x2,x3,x4) are displayed.
Setup(differentiationvariables = Y);
⁢Default differentiation variables for d_, D_ and dAlembertian are:⁢Y=y1,y2,y3,y4
The dAlembertian enters the classical field equations in various models; this is the electromagnetic field tensor.
F[mu, nu] := d_[mu](A[nu](X)) - d_[nu](A[mu](X));
Maxwell equations result from taking the functional derivative of the Action. Use delay evaluation quotes to display the operation unevaluated, then perform the operation on the next line.
'Fundiff'( Intc(F[mu,nu]^2, X), A[rho](Y));
subs(Y = X, (20));
To simplify the contracted spacetime indices, use the Simplify command.
The Lagrangean of the lambdaPhi^4 model, the corresponding Action, and the field equations:
L := 1/2 * d_[mu](Phi(X)) * d_[mu](Phi(X)) - m^2/2 * Phi(X)^2 + lambda/4*Phi(X)^4;
subs(Y = X, (24));
Coordinates, d_, Define, diff, Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Setup
Download Help Document
What kind of issue would you like to report? (Optional)