Physics for Maple 2023

All Products Maple MapleSim

Home : Support : Online Help : Physics : Updates : Physics for Maple 2023

Physics

Maple provides a state-of-the-art environment for algebraic computations in Physics, with emphasis on ensuring that the computational experience is as natural as possible. The theme of the Physics project for Maple 2023 has been the use of the package in Education, new commands, the consolidation of the functionality introduced in previous releases, also related to the StandardModel package, and several enhancements necessary for producing the new Courseware-Support, Mechanics material.

As part of its commitment to providing the best possible computational environment in Physics, Maplesoft launched a Maple Physics: Research and Development website in 2014, which enabled users to download research versions of the package, ask questions, and provide feedback. The results from this accelerated exchange have been incorporated into the Physics package in Maple 2023. The presentation below illustrates both the novelties and the kind of mathematical formulations that can now be performed.

Courseware support: Mechanics

New command Lagrange Equations

New command Substitute

New command StandardModel:-Lagrangian

New Physics:-Library commands: FactorSum and RemoveCommonFactors

New coordinates labels

Documentation advanced examples

See Also

Courseware support: Mechanics

The Physics package, with its various subpackages, is now a resourceful, mature project. For Maple 2023, part of the focus has been in its use in Education in Physics, starting with Mechanics. The resulting material is now part of Maple's help system and can be opened to interact with it as a Maple document. The presentation is organized as a hyperlinked syllabus covering typical topics appearing in Mechanics undergrad courses of the first years. For the Maple 2023 release, Part I covers:

1	Position, velocity and acceleration in Cartesian, cylindrical and spherical coordinates

a	The position $\vec{r} (t)$ as a function of time

b	The velocity $\vec{v} (t)$

c	The acceleration $\vec{a} (t)$

d	Deriving these formulas

e	Velocity and acceleration in the case of 2-dimensional motion on the x, y plane

1	The equations of motion

a	A single particle

i	The equations of motion - vectorial form

ii	The case of constant acceleration

iii	Motion under gravitational force close to the Earth's surface

iv	Motion under gravitational force not close to the Earth's surface

A	Circular motion

B	Escape velocity

i	Different acceleration in different regions

ii	The equations of motion using tensor notation

A	Cartesian coordinates

B	Curvilinear coordinates

a	Many-particle systems

i	Center of mass

ii	The equations of motion

iii	Static: reactions of planes and tensions on cables

a	Lagrange equations

i	Motion of a pendulum

1	Conservation laws

Work

b	Conservation of the total energy of a closed system or a system in a constant external field

c	Conservation of the total momentum of a closed system

d	Conservation of angular momentum

e	Cyclic coordinates

1	Integration of the equations of motion

a	Motion in one dimension

b	Reduced mass

i	The two-body problem

ii	A many-body problem

a	Motion in a central field

b	Kepler's problem

1	Small Oscillations

a	Free oscillations in one dimension

b	Forced oscillations

c	Oscillations of systems with many degrees of freedom

1	Rigid-body motion

a	Angular velocity

b	Inertia tensor

c	Angular momentum of a rigid body

d	The equations of motion of a rigid body

1	Non-inertial coordinate systems

a	Coriolis force and centripetal force

New command Lagrange Equations

There is now a new Physics command, LagrangeEquations, taking advantage of the functional differentiation capabilities. LagrangeEquations receives an expression representing a Lagrangian and returns a sequence of Lagrange equations with as many equations as coordinates are indicated. The formula in the traditional case where the Lagrangian depends on 1st order derivatives of the coordinates and there is only one parameter, $t$ , is

$\frac{&DifferentialD;}{&DifferentialD; t} \frac{\partial}{\partial {\vec{v}}_{i}} L = \frac{\partial}{\partial {\vec{r}}_{i}} L$

where $\frac{\partial}{\partial {\vec{r}}_{i}} L$ formally represents the derivative with respect to the coordinates of the $i^{th}$ particle, equal to the Gradient when working in Cartesian coordinates; $\frac{\partial}{\partial {\vec{v}}_{i}} L$ represents the equivalent operation, replacing each coordinate by the corresponding generalized velocity and $\frac{&DifferentialD;}{&DifferentialD; t}$ represents the total derivative with respect to $t$ , the parameter parametrizing the coordinates. In more general cases the number of parameters can be many. For example, in electrodynamics, the "coordinate" is a tensor field $A_{μ} (x, y, z, t)$ , there are then four coordinates, one for each of the values of the index $μ$ , and there are four parameters $(x, y, z, t)$ . LagrangeEquations can handle tensors and vectors of the Physics package as well as derivatives using vectorial differential operators (see d_ and Nabla), works by performing functional differentiation (see Fundiff), and handles 1st, and higher order derivatives of the coordinates in the Lagrangian automatically.

Examples

$with (Physics) &colon;$

$Setup (mathematicalnotation = true, coordinates = cartesian)$

$[coordinatesystems = \{X\}, mathematicalnotation = true]$

(1)

The Lagrangian of a one-dimensional oscillator - small oscillations

$L ≔ \frac{1}{2} {diff (x (t), t)}^{2} - \frac{1}{2} k {x (t)}^{2}$

$L ≔ \frac{{\dot{x} (t)}^{2}}{2} - \frac{k {x (t)}^{2}}{2}$

(2)

The corresponding Lagrange equation gives Newton's second law, a 2nd order linear ODE for $x (t)$

$LagrangeEquations (L, x)$

$\overset{..}{x} (t) + x (t) k = 0$

(3)

The Lagrangian of a pendulum of mass $m$ and length $l$ where the suspension point moves uniformly over a vertical circumference centered at the origin, with a constant frequency $ω$

$CompactDisplay (φ (t))$

$φ (t) will now be displayed as φ$

(4)

$L ≔ \frac{1}{2} m (- 2 diff (φ (t), t) a l ω \sin (ω t - φ (t)) + {diff (φ (t), t)}^{2} l^{2} + 2 \cos (φ (t)) g l)$

$L ≔ \frac{m (- 2 \dot{φ} a l ω \sin (ω t - φ) + {\dot{φ}}^{}^{2} l^{2} + 2 \cos (φ) g l)}{2}$

(5)

The Lagrange equations

$LagrangeEquations (L, φ)$

$- m l (- a ω^{2} \cos (ω t - φ) + \overset{..}{φ} l + \sin (φ) g) = 0$

(6)

The Maxwell equations can be derived as Lagrange equations as follows. For simplicity, consider Maxwell equations in vacuum. Define first a tensor representing the 4D electromagnetic field potential

$Define (A [μ])$

$\{A_{μ}, γ_{μ}, σ_{μ}, \partial_{μ}, g_{μ, ν}, ε_{α, β, μ, ν}, X_{μ}\}$

(7)

The electromagnetic field tensor

$F [μ, ν] ≔ d_[μ] (A [ν] (X)) - d_[ν] (A [μ] (X))$

$F_{μ, ν} ≔ \partial_{μ} (A_{ν} (X)) - \partial_{ν} (A_{μ} (X))$

(8)

The Lagrangian

$L ≔ {F [μ, ν]}^{2}$

$L ≔ (\partial_{μ} (A_{ν} (X)) - \partial_{ν} (A_{μ} (X))) (\partial_{}^{μ} (A_{}^{ν} (X)) - \partial_{}^{ν} (A_{}^{μ} (X)))$

(9)

Maxwell equations in 4D tensorial notation

$LagrangeEquations (L, A)$

$- 4 □ (A_{}^{α} (X)) + 4 \partial_{μ} (\partial_{}^{α} (A_{}^{μ} (X))) = 0$

(10)

New command Substitute

The command for doing syntactical exact-match substitutions is subs. Frequently, however, what we intend to accomplish with a substitution is more of a mathematical substitution. For different kinds of mathematical substitutions Maple has the commands algsubs, Physics:-SubstituteTensor, Physics:-SubstituteTensorIndices, Physics:-Library:-SubstituteOperator and Physics:-Library:-SubstituteMatrix, resulting in a myriad of commands, useful, however difficult to remember. A new single command, Physics:-Substitute, unifies all of those so that one does not need to remember which one to use.

Examples

$restart &semi; with (Physics) &colon;$

•	Substitutions of sub-products and sub-sums, possibly involving noncommutative objects

$eq ≔ F + A B (x) = E [j]$

$eq ≔ F + A B (x) = E_{j}$

(11)

$ee ≔ A B (x) C$

(12)

$Substitute (eq, ee)$

$- C (F - E_{j})$

(13)

•	Substitution of matrices and vectors taking mathematical equality into account

$eq ≔ alpha + Matrix (2, symbol = a) = β + Matrix (2, symbol = b)$

$eq ≔ α + [\begin{array}{c} a_{1, 1} & a_{1, 2} \\ a_{2, 1} & a_{2, 2} \end{array}] = β + [\begin{array}{c} b_{1, 1} & b_{1, 2} \\ b_{2, 1} & b_{2, 2} \end{array}]$

(14)

$ee ≔ Matrix (2, symbol = a)$

$ee ≔ [\begin{array}{c} a_{1, 1} & a_{1, 2} \\ a_{2, 1} & a_{2, 2} \end{array}]$

(15)

$Substitute (eq, ee)$

$- α + β + [\begin{array}{c} b_{1, 1} & b_{1, 2} \\ b_{2, 1} & b_{2, 2} \end{array}]$

(16)

•	Set up some non-commutative operands for substitution

$Setup (quantumoperators = \{A, B, C, E, F\}) &colon;$

Some typical examples not handled by subs, eval or algsubs

$eq ≔ F + A B = E$

(17)

$ee ≔ \exp (λ A B C) B \exp (- λ A B C)$

$ee ≔ {&ExponentialE;}^{λ A B C} B {&ExponentialE;}^{- λ A B C}$

(18)

$Substitute (eq, ee)$

${&ExponentialE;}^{λ (E - F) C} B {&ExponentialE;}^{- λ (E - F) C}$

(19)

•	Substituting tensorial sub-expressions

$Define (A, B, C, F, G) &colon;$

$Defined objects with tensor properties$

(20)

$eq ≔ A [μ] = G [ν, α] A [α] F [μ, ν]$

$eq ≔ A_{μ} = A_{α} F_{μ, ν} G_{}^{ν, α}$

(21)

$ee ≔ A [ν] A [ν]$

$ee ≔ A_{ν} A_{}^{ν}$

(22)

$Substitute (eq, ee)$

$A_{α} F_{ν, β} G_{}^{β, α} A_{κ} F_{λ}^{ν} G_{}^{λ, κ}$

(23)

$Setup (spinorindices = lowercaselatin_is)$

$[spinorindices = lowercaselatin_is]$

(24)

$Define (T [μ, i, j])$

$\{A, B, C, F, G, γ_{μ}, σ_{μ}, T_{μ, i, j}, \partial_{μ}, g_{μ, ν}, ε_{α, β, μ, ν}\}$

(25)

$Substitute (Dgamma [μ] [i, j] = T [μ, i, j], Dgamma [~alpha] [~k, m], disregardfreeindices)$

$T_{m}^{α, k}$

(26)

•	Substitutions of tensor indices

$g_[α, μ] A [~mu] g_[~alpha, ~nu] B [ν, σ, ~rho]$

$g_{α, μ} g_{}^{α, ν} A_{}^{μ} B_{ν, σ}^{ρ}$

(27)

$Check (, all)$

$[\{α, μ, ν\}], \{σ, ~rho\}$

(28)

$Substitute (α = β, mu = delta, sigma = rho,)$

$g_{β, δ} g_{}^{β, ν} A_{}^{δ} B_{ν, ρ}^{ρ}$

(29)

New command StandardModel:-Lagrangian

One of the distinctive aspects of the Standard Model is the complexity of its Lagrangian. In this context, the new StandardModel:-Lagrangian returns the Lagrangian of the model after symmetry breaking, optionally restricted to only the interaction terms, or only one of its QED, QCD or electroweak sectors, or only one of the different sub-terms involved in the electroweak part; all of that with the covariant derivatives, and sums over leptons and quarks optionally expanded.

The algebraic expressions returned by Lagrangian are fully computable, so you can use them as starting point to construct other Lagrangians (add or subtract terms), or the Action and related field equations (see d_, D_ for covariant derivatives, diff and Fundiff for functional differentiation), or to compute scattering amplitudes (see FeynmanDiagrams and FeynmanIntegral).

Examples

$restart &semi; with (Physics) &colon; with (StandardModel)$

$_______________________________________________________$

$Setting lowercaselatin_is letters to represent Dirac spinor indices$

$Setting lowercaselatin_ah letters to represent SU(3) adjoint representation, (1..8) indices$

$Setting uppercaselatin_ah letters to represent SU(3) fundamental representation, (1..3) indices$

$Setting uppercaselatin_is letters to represent SU(2) adjoint representation, (1..3) indices$

$Setting uppercasegreek letters to represent SU(2) fundamental representation, (1..2) indices$

$_______________________________________________________$

$Defined as the electron, muon and tau leptons and corresponding neutrinos: e_{j}, μ_{j}, τ_{j}, {ν^{(e)}}_{j}, {ν^{(μ)}}_{j}, {ν^{(τ)}}_{j}$

$Defined as the up, charm, top, down, strange and bottom quarks: u_{A, j}, c_{A, j}, t_{A, j}, d_{A, j}, s_{A, j}, b_{A, j}$

$Defined as gauge tensors: B_{μ}, 𝔹_{μ, ν}, A_{μ}, 𝔽_{μ, ν}, W_{μ, J}, 𝕎_{μ, ν, J}, {W^{+}}_{μ}, {𝕎^{+}}_{μ, ν}, {W^{-}}_{μ}, {𝕎^{-}}_{μ, ν}, Z_{μ}, ℤ_{μ, ν}, G_{μ, a}, 𝔾_{μ, ν, a}$

$Defined as Gell-Mann (Glambda), Pauli (Psigma) and Dirac (Dgamma) matrices: λ_{a}, σ_{J}, γ_{μ}$

$Defined as the electric, weak and strong coupling constants: g__e, g__w, g__s$

$Defined as the charge in units of | g__e | for 1) the electron, muon and tauon, 2) the up, charm and top, and 3) the down, strange and bottom: q__e = −1, q__u = \frac{2}{3}, q__d = - \frac{1}{3}$

$Defined as the weak isospin for 1) the electron, muon and tauon, 2) the up, charm and top, 3) the down, strange and bottom, and 4) all the neutrinos: I__e = - \frac{1}{2}, I__u = \frac{1}{2}, I__d = - \frac{1}{2}, I__n = \frac{1}{2}$

$You can use the active form without the % prefix, or the 'value' command to give the corresponding value to any of the inert representations q__e, q__u, q__d, I__e, I__u, I__d, I__n$

$_______________________________________________________$

$Default differentiation variables for d_, D_ and dAlembertian are: \{X = (x, y, z, t)\}$

$Minkowski spacetime with signatre (- - - +)$

$_______________________________________________________$

$[I__d, I__e, I__n, I__u, q__d, q__e, q__u, BField, BFieldStrength, Bottom, CKM, Charm, Down, ElectromagneticField, ElectromagneticFieldStrength, Electron, ElectronNeutrino, FSU3, Glambda, GluonField, GluonFieldStrength, HiggsBoson, Lagrangian, Muon, MuonNeutrino, Strange, Tauon, TauonNeutrino, Top, Up, WField, WFieldStrength, WMinusField, WMinusFieldStrength, WPlusField, WPlusFieldStrength, WeinbergAngle, ZField, ZFieldStrength, g__e, g__s, g__w]$

(30)

The massless fields of the model are the electromagnetic and gluon fields and the three neutrinos

$Setup (massless)$

$[masslessfields = \{G, ν^{(μ)}, ν^{(τ)}, A, ν^{(e)}\}]$

(31)

The Leptons and Quarks of the model are

$StandardModel :- Leptons$

$[e, μ, τ, ν^{(e)}, ν^{(μ)}, ν^{(τ)}]$

(32)

$StandardModel :- Quarks$

$[u, c, t, d, s, b]$

(33)

The Gauge fields, and their related field strengths displayed with Open Face type fonts

$StandardModel :- GaugeFields$

$[A, 𝔽, B, 𝔹, W, 𝕎, G, 𝔾, W^{-}, 𝕎^{-}, W^{+}, 𝕎^{+}, Z, ℤ]$

(34)

To represent the interaction Lagrangians for the QCD and electroweak sectors as sums over leptons and quarks, all of them fermions, it is useful to introduce four anticommutative prefixes, used below as summation indices in the formulas

$Setup (anticommutativeprefix = \{f__D, f__L, f__Q, f__U\})$

$[anticommutativeprefix = \{f__D, f__L, f__Q, f__U\}]$

(35)

For readability, omit from the display of formulas the functionality of all the fields entering the Standard Model (see CompactDisplay) and use the lowercase i instead of the uppercase I to represent the imaginary unit

$CompactDisplay ((StandardModel :- Leptons, StandardModel :- Quarks, StandardModel :- GaugeFields, HiggsBoson, f__L, f__Q, f__U, f__D) (X), quiet) &colon; interface (imaginaryunit = i) &colon;$

The Lagrangian of the whole Standard Model after symmetry breaking, in its most compact form:

$Lagrangian ()$

$\sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} (ⅈ {(γ_{μ})}_{j, k} ▿_{μ} - m_{f__Q} δ_{j, k}) {f__Q}_{A, k} - \frac{𝔾_{μ, ν, a}^{2}}{4} - \frac{𝔽_{μ, ν}^{2}}{4} - \frac{{𝕎^{+}}_{μ, ν} {𝕎^{-}}_{μ, ν}}{2} + m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} - \frac{ℤ_{μ, ν}^{2}}{4} + \frac{m_{Z}^{2} Z_{μ}^{2}}{2} + \frac{{\partial_{μ} (Φ)}^{2}}{2} - \frac{m_{Φ}^{2} Φ^{2}}{2} + \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__L}_{k}) - m_{f__L} {f__L}_{j}) + \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} ⅈ {(γ_{μ})}_{j, k} \overset{&conjugate0;}{{f__L}_{j}} \partial_{μ} ({f__L}_{k}) + \sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__Q}_{A, k}) - m_{f__Q} {f__Q}_{A, j}) + g__e {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k}) A_{μ} + \frac{g__w ({(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) (I__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__n \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l} + I__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l}) - {\sin (θ_{w})}^{2} {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k})) Z_{μ}}{\cos (θ_{w})} - \frac{g__w \sqrt{2} {(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) ((\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} 𝕄_{f__U, f__D} \overset{&conjugate0;}{{f__U}_{A, j}} {f__D}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{1})}_{j} (X)} {({f__L}_{2})}_{l} (X)) {W^{+}}_{μ} + (\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} \overset{&conjugate0;}{𝕄_{f__U, f__D}} \overset{&conjugate0;}{{f__D}_{A, j}} {f__U}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{2})}_{j} (X)} {({f__L}_{1})}_{l} (X)) {W^{-}}_{μ})}{2} - \frac{g__w m_{Φ}^{2} (Φ^{3} + \frac{Φ^{4}}{8 m_{W}})}{4 m_{W}} + (\frac{g__w Φ}{m_{W}} + \frac{{g__w}^{2} Φ^{2}}{4 m_{W}^{2}}) (m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} + \frac{{Z_{μ}}^{2} m_{Z}^{2}}{2}) - ⅈ g__w (({𝕎^{+}}_{μ, ν} {W^{-}}_{μ} - {W^{+}}_{μ} {𝕎^{-}}_{μ, ν}) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})) + {W^{-}}_{ν} {W^{+}}_{μ} (𝔽_{μ, ν} \sin (θ_{w}) - ℤ_{μ, ν} \cos (θ_{w}))) - \frac{{g__w}^{2} ((2 {W^{+}}_{μ} {W^{-}}_{μ} + {(A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w}))}^{2}) (2 {W^{+}}_{ν} {W^{-}}_{ν} + {(A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w}))}^{2}) + {({W^{+}}_{μ} {W^{-}}_{ν} + {W^{+}}_{ν} {W^{-}}_{μ} + (A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w})) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})))}^{2})}{4} - \frac{g__w (\sum_{f__L = [e, μ, τ, ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} m_{f__L} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{j} + \sum_{f__Q = [u, c, t, d, s, b]} m_{f__Q} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, j}) Φ}{2 m_{W}}$

(36)

In the output above we see, among other things, the $γ_{5}$ Dirac matrix, and the Cabibbo - Kobayashi - Maskawa matrix $𝕄$ , and the tensor indices of different kinds all explicit. See StandardModel for the notational conventions used, which are standard in the literature but for a few things, like a sign in the definition of $γ_{5}$ , that depend on the reference.

The Quantum Electrodynamics (QED) Lagrangian

The simplest sector of this Lagrangian (8) is the QED one

$Lagrangian (QED)$

$\overset{&conjugate0;}{e_{j}} (ⅈ {(γ_{μ})}_{j, k} ▿_{}^{μ} - m_{e} δ_{j, k}) e_{k} - \frac{𝔽_{μ, ν} 𝔽_{}^{μ, ν}}{4}$

(37)

The applied form can be obtained using the Library command ApplyProductsOfDifferentialOperators over the output (9) or passing the optional argument applied

$Lagrangian (QED, applied)$

$\overset{&conjugate0;}{e_{j}} (ⅈ ▿_{μ} (e_{k}) {(γ_{}^{μ})}_{j, k} - m_{e} δ_{j, k} e_{k}) - \frac{𝔽_{μ, ν} 𝔽_{}^{μ, ν}}{4}$

(38)

Only the interaction part of this Lagrangian is relevant when computing scattering amplitudes. To get that part, you can either expand the covariant derivative operator or or pass the optional keyword expanded, in which case also the trace of $&Fopf;__μ,ν$ gets expanded

$Lagrangian (QED, expanded)$

$ⅈ \overset{&conjugate0;}{e_{j}} \partial_{μ} (e_{k}) {(γ_{}^{μ})}_{j, k} - g__e \overset{&conjugate0;}{e_{j}} e_{k} A_{μ} {(γ_{}^{μ})}_{j, k} - m_{e} \overset{&conjugate0;}{e_{j}} e_{j} - \frac{(\partial_{μ} (A_{ν}) - \partial_{ν} (A_{μ})) (\partial_{}^{μ} (A_{}^{ν}) - \partial_{}^{ν} (A_{}^{μ}))}{4}$

(39)

then discard the non-interaction terms

$remove (has,, [d_, m])$

$- {(γ_{}^{μ})}_{j, k} g__e \overset{&conjugate0;}{e_{j}} e_{k} A_{μ}$

(40)

or simpler: pass the keyword interaction

$Lagrangian (QED, interaction)$

$- {(γ_{}^{μ})}_{j, k} g__e \overset{&conjugate0;}{e_{j}} e_{k} A_{μ}$

(41)

All the algebraic expressions returned by Lagrangian are fully computable in that further calculations can proceed starting from them. For example (see FeynmanDiagrams), this is the self-energy of the electron

$FeynmanDiagrams (, incoming = [Electron], outgoing = [Electron], numberofloops = 1, diagrams)$

$- \int - \frac{{(u_{e})}_{l} ({\vec{P}}_{1}) \overset{&conjugate0;}{{(u_{e})}_{m} ({\vec{P}}_{2})} {g__e}^{2} {(γ_{}^{α})}_{m, n} {(γ_{}^{ν})}_{p, l} (({P__1}_{β} + {p__2}_{β}) {(γ_{}^{β})}_{n, p} + m_{e} δ_{n, p}) g_{α, ν} δ (- P__2 + P__1)}{8 π^{3} ({(P__1 + p__2)}^{2} - m_{e}^{2} + ⅈ ε) ({p__2}^{2} + ⅈ ε)} &DifferentialD; {p__2}^{4}$

(42)

The Quantum Chromodynamics (QCD) Lagrangian

Next in complexity is the QCD Lagrangian

$Lagrangian (QCD)$

$\sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} (ⅈ {(γ_{μ})}_{j, k} ▿_{}^{μ} - m_{f__Q} δ_{j, k}) {f__Q}_{A, k} - \frac{𝔾_{μ, ν, a} 𝔾_{a}^{μ, ν}}{4}$

(43)

$Lagrangian (QCD, expanded)$

$\overset{&conjugate0;}{u_{A, j}} (ⅈ (\partial_{μ} (u_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} u_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{u} δ_{j, k} u_{A, k}) + \overset{&conjugate0;}{c_{A, j}} (ⅈ (\partial_{μ} (c_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} c_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{c} δ_{j, k} c_{A, k}) + \overset{&conjugate0;}{t_{A, j}} (ⅈ (\partial_{μ} (t_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} t_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{t} δ_{j, k} t_{A, k}) + \overset{&conjugate0;}{d_{A, j}} (ⅈ (\partial_{μ} (d_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} d_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{d} δ_{j, k} d_{A, k}) + \overset{&conjugate0;}{s_{A, j}} (ⅈ (\partial_{μ} (s_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} s_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{s} δ_{j, k} s_{A, k}) + \overset{&conjugate0;}{b_{A, j}} (ⅈ (\partial_{μ} (b_{A, k}) - \frac{ⅈ g__s {(λ_{a})}_{A, B} b_{B, k} G_{μ, a}}{2}) {(γ_{}^{μ})}_{j, k} - m_{b} δ_{j, k} b_{A, k}) - \frac{(\partial_{μ} (G_{ν, a}) - \partial_{ν} (G_{μ, a}) + g__s {f_{su3}}_{a, b, c} G_{μ, b} G_{ν, c}) (\partial_{}^{μ} (G_{a}^{ν}) - \partial_{}^{ν} (G_{a}^{μ}) + g__s {f_{su3}}_{a, d, e} G_{d}^{μ} G_{e}^{ν})}{4}$

(44)

For computing scattering amplitudes, only the interaction part of this Lagrangian is relevant

$Lagrangian (QCD, interaction)$

$\frac{g__s {(λ_{a})}_{A, B} \sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{B, k} G_{μ, a} {(γ_{}^{μ})}_{j, k}}{2} - g__s {f_{su3}}_{a, b, c} (\partial_{μ} (G_{ν, a}) G_{b}^{μ} G_{c}^{ν} - \frac{g__s {f_{su3}}_{c, d, e} G_{μ, a} G_{α, b} G_{e}^{μ} G_{d}^{α}}{4})$

(45)

$Lagrangian (QCD, interaction, expanded)$

$\frac{g__s {(λ_{a})}_{A, B} (\overset{&conjugate0;}{u_{A, j}} u_{B, k} + \overset{&conjugate0;}{c_{A, j}} c_{B, k} + \overset{&conjugate0;}{t_{A, j}} t_{B, k} + \overset{&conjugate0;}{d_{A, j}} d_{B, k} + \overset{&conjugate0;}{s_{A, j}} s_{B, k} + \overset{&conjugate0;}{b_{A, j}} b_{B, k}) G_{μ, a} {(γ_{}^{μ})}_{j, k}}{2} - g__s {f_{su3}}_{a, b, c} (\partial_{μ} (G_{ν, a}) G_{b}^{μ} G_{c}^{ν} - \frac{g__s {f_{su3}}_{c, d, e} G_{μ, a} G_{α, b} G_{e}^{μ} G_{d}^{α}}{4})$

(46)

The amplitude at tree level for the process with two incoming and two outgoing Up quarks (particle and antiparticle) exchanging a gluon

$FeynmanDiagrams (, incomingparticles = [Up, conjugate (Up)], outgoingparticles = [Up, conjugate (Up)], numberofloops = 0, diagrams)$

$- \frac{ⅈ {(u_{u})}_{C, l} ({\vec{P}}_{1}) \overset{&conjugate0;}{{(v_{u})}_{E, m} ({\vec{P}}_{2})} \overset{&conjugate0;}{{(u_{u})}_{F, n} ({\vec{P}}_{3})} {(v_{u})}_{G, p} ({\vec{P}}_{4}) {g__s}^{2} {(λ_{g})}_{F, G} {(γ_{}^{κ})}_{n, p} {(λ_{f})}_{E, C} {(γ_{}^{β})}_{m, l} g_{β, κ} δ_{f, g} δ (- {P__3}_{}^{λ} - {P__4}_{}^{λ} + {P__1}_{}^{λ} + {P__2}_{}^{λ})}{16 π^{2} (({P__1}_{σ} + {P__2}_{σ}) ({P__1}_{}^{σ} + {P__2}_{}^{σ}) + ⅈ ε)} + \frac{ⅈ {(u_{u})}_{C, l} ({\vec{P}}_{1}) \overset{&conjugate0;}{{(v_{u})}_{E, m} ({\vec{P}}_{2})} \overset{&conjugate0;}{{(u_{u})}_{F, n} ({\vec{P}}_{3})} {(v_{u})}_{G, p} ({\vec{P}}_{4}) {g__s}^{2} {(λ_{g})}_{E, G} {(γ_{}^{κ})}_{m, p} {(λ_{f})}_{F, C} {(γ_{}^{β})}_{n, l} g_{β, κ} δ_{f, g} δ (- {P__3}_{}^{λ} - {P__4}_{}^{λ} + {P__1}_{}^{λ} + {P__2}_{}^{λ})}{16 π^{2} (({P__1}_{σ} - {P__3}_{σ}) ({P__1}_{}^{σ} - {P__3}_{}^{σ}) + ⅈ ε)}$

(47)

The Electro-Weak Lagrangian

The electroweak sector of the Standard Model Lagrangian is significantly more complicated

$Lagrangian (electroweak)$

$- \frac{𝔽_{μ, ν}^{2}}{4} - \frac{{𝕎^{+}}_{μ, ν} {𝕎^{-}}_{μ, ν}}{2} + m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} - \frac{ℤ_{μ, ν}^{2}}{4} + \frac{m_{Z}^{2} Z_{μ}^{2}}{2} + \frac{{\partial_{μ} (Φ)}^{2}}{2} - \frac{m_{Φ}^{2} Φ^{2}}{2} + \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__L}_{k}) - m_{f__L} {f__L}_{j}) + \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} ⅈ {(γ_{μ})}_{j, k} \overset{&conjugate0;}{{f__L}_{j}} \partial_{μ} ({f__L}_{k}) + \sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__Q}_{A, k}) - m_{f__Q} {f__Q}_{A, j}) + g__e {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k}) A_{μ} + \frac{g__w ({(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) (I__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__n \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l} + I__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l}) - {\sin (θ_{w})}^{2} {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k})) Z_{μ}}{\cos (θ_{w})} - \frac{g__w \sqrt{2} {(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) ((\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} 𝕄_{f__U, f__D} \overset{&conjugate0;}{{f__U}_{A, j}} {f__D}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{1})}_{j} (X)} {({f__L}_{2})}_{l} (X)) {W^{+}}_{μ} + (\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} \overset{&conjugate0;}{𝕄_{f__U, f__D}} \overset{&conjugate0;}{{f__D}_{A, j}} {f__U}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{2})}_{j} (X)} {({f__L}_{1})}_{l} (X)) {W^{-}}_{μ})}{2} - \frac{g__w m_{Φ}^{2} (Φ^{3} + \frac{Φ^{4}}{8 m_{W}})}{4 m_{W}} + (\frac{g__w Φ}{m_{W}} + \frac{{g__w}^{2} Φ^{2}}{4 m_{W}^{2}}) (m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} + \frac{{Z_{μ}}^{2} m_{Z}^{2}}{2}) - ⅈ g__w (({𝕎^{+}}_{μ, ν} {W^{-}}_{μ} - {W^{+}}_{μ} {𝕎^{-}}_{μ, ν}) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})) + {W^{-}}_{ν} {W^{+}}_{μ} (𝔽_{μ, ν} \sin (θ_{w}) - ℤ_{μ, ν} \cos (θ_{w}))) - \frac{{g__w}^{2} ((2 {W^{+}}_{μ} {W^{-}}_{μ} + {(A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w}))}^{2}) (2 {W^{+}}_{ν} {W^{-}}_{ν} + {(A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w}))}^{2}) + {({W^{+}}_{μ} {W^{-}}_{ν} + {W^{+}}_{ν} {W^{-}}_{μ} + (A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w})) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})))}^{2})}{4} - \frac{g__w (\sum_{f__L = [e, μ, τ, ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} m_{f__L} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{j} + \sum_{f__Q = [u, c, t, d, s, b]} m_{f__Q} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, j}) Φ}{2 m_{W}}$

(48)

To decipher this result it is useful to see the structure of physically recognizable terms. For that purpose, you can use the keyword showterms

$Lagrangian (electroweak, showterms)$

$L__K + L__N + L__C + L__H + L__HV + L__WWV + L__WWVV + L__Y where [L__K = - \frac{𝔽_{μ, ν}^{2}}{4} - \frac{{𝕎^{+}}_{μ, ν} {𝕎^{-}}_{μ, ν}}{2} + m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} - \frac{ℤ_{μ, ν}^{2}}{4} + \frac{m_{Z}^{2} Z_{μ}^{2}}{2} + \frac{{\partial_{μ} (Φ)}^{2}}{2} - \frac{m_{Φ}^{2} Φ^{2}}{2} + \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__L}_{k}) - m_{f__L} {f__L}_{j}) + \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} ⅈ {(γ_{μ})}_{j, k} \overset{&conjugate0;}{{f__L}_{j}} \partial_{μ} ({f__L}_{k}) + \sum_{f__Q = [u, c, t, d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} (ⅈ {(γ_{μ})}_{j, k} \partial_{μ} ({f__Q}_{A, k}) - m_{f__Q} {f__Q}_{A, j}), L__N = g__e {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k}) A_{μ} + \frac{g__w ({(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) (I__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__n \sum_{f__L = [ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{l} + I__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l} + I__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, l}) - {\sin (θ_{w})}^{2} {(γ_{μ})}_{j, k} (q__e \sum_{f__L = [e, μ, τ]} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{k} + q__u \sum_{f__Q = [u, c, t]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k} + q__d \sum_{f__Q = [d, s, b]} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, k})) Z_{μ}}{\cos (θ_{w})}, L__C = - \frac{g__w \sqrt{2} {(γ_{μ})}_{j, k} (δ_{k, l} + {(γ_{5})}_{k, l}) ((\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} 𝕄_{f__U, f__D} \overset{&conjugate0;}{{f__U}_{A, j}} {f__D}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{1})}_{j} (X)} {({f__L}_{2})}_{l} (X)) {W^{+}}_{μ} + (\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} \overset{&conjugate0;}{𝕄_{f__U, f__D}} \overset{&conjugate0;}{{f__D}_{A, j}} {f__U}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{2})}_{j} (X)} {({f__L}_{1})}_{l} (X)) {W^{-}}_{μ})}{2}, L__H = - \frac{g__w m_{Φ}^{2} (Φ^{3} + \frac{Φ^{4}}{8 m_{W}})}{4 m_{W}}, L__HV = (\frac{g__w Φ}{m_{W}} + \frac{{g__w}^{2} Φ^{2}}{4 m_{W}^{2}}) (m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{μ} + \frac{{Z_{μ}}^{2} m_{Z}^{2}}{2}), L__WWV = −ⅈ g__w (({𝕎^{+}}_{μ, ν} {W^{-}}_{μ} - {W^{+}}_{μ} {𝕎^{-}}_{μ, ν}) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})) + {W^{-}}_{ν} {W^{+}}_{μ} (𝔽_{μ, ν} \sin (θ_{w}) - ℤ_{μ, ν} \cos (θ_{w}))), L__WWVV = - \frac{{g__w}^{2} ((2 {W^{+}}_{μ} {W^{-}}_{μ} + {(A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w}))}^{2}) (2 {W^{+}}_{ν} {W^{-}}_{ν} + {(A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w}))}^{2}) + {({W^{+}}_{μ} {W^{-}}_{ν} + {W^{+}}_{ν} {W^{-}}_{μ} + (A_{μ} \sin (θ_{w}) - Z_{μ} \cos (θ_{w})) (A_{ν} \sin (θ_{w}) - Z_{ν} \cos (θ_{w})))}^{2})}{4}, L__Y = - \frac{g__w (\sum_{f__L = [e, μ, τ, ν^{(e)}, ν^{(μ)}, ν^{(τ)}]} m_{f__L} \overset{&conjugate0;}{{f__L}_{j}} {f__L}_{j} + \sum_{f__Q = [u, c, t, d, s, b]} m_{f__Q} \overset{&conjugate0;}{{f__Q}_{A, j}} {f__Q}_{A, j}) Φ}{2 m_{W}}]$

(49)

In this result we see a sum of $L_{terms}$ , and after 'where' there is a list of equations with the formulas represented by each $L_{term}$ . Take from the above, for instance, only the charged current $L_{C}$ term that involves interaction between the leptons and the corresponding neutrinos: you can do that with the mouse, copy and paste, or using the $term = .. &period;$ option of Lagrangian

$Lagrangian (electroweak, term = L [C])$

$L__C = - \frac{g__w \sqrt{2} (δ_{k, l} + {(γ_{5})}_{k, l}) ((\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} 𝕄_{f__U, f__D} \overset{&conjugate0;}{{f__U}_{A, j}} {f__D}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{1})}_{j} (X)} {({f__L}_{2})}_{l} (X)) {W^{+}}_{μ} + (\sum_{f__D = [d, s, b]} \sum_{f__U = [u, c, t]} \overset{&conjugate0;}{𝕄_{f__U, f__D}} \overset{&conjugate0;}{{f__D}_{A, j}} {f__U}_{A, l} + \sum_{f__L = [[ν^{(e)}, e], [ν^{(μ)}, μ], [ν^{(τ)}, τ]]} \overset{&conjugate0;}{{({f__L}_{2})}_{j} (X)} {({f__L}_{1})}_{l} (X)) {W^{-}}_{μ}) {(γ_{}^{μ})}_{j, k}}{2}$

(50)

The same term in expanded form, useful for computing scattering amplitudes

$Lagrangian (electroweak, term = L [C], expanded)$

$L__C = - \frac{g__w \sqrt{2} (δ_{k, l} + {(γ_{5})}_{k, l}) ((𝕄_{u, d} \overset{&conjugate0;}{u_{A, j}} d_{A, l} + 𝕄_{c, d} \overset{&conjugate0;}{c_{A, j}} d_{A, l} + 𝕄_{t, d} \overset{&conjugate0;}{t_{A, j}} d_{A, l} + 𝕄_{u, s} \overset{&conjugate0;}{u_{A, j}} s_{A, l} + 𝕄_{c, s} \overset{&conjugate0;}{c_{A, j}} s_{A, l} + 𝕄_{t, s} \overset{&conjugate0;}{t_{A, j}} s_{A, l} + 𝕄_{u, b} \overset{&conjugate0;}{u_{A, j}} b_{A, l} + 𝕄_{c, b} \overset{&conjugate0;}{c_{A, j}} b_{A, l} + 𝕄_{t, b} \overset{&conjugate0;}{t_{A, j}} b_{A, l} + \overset{&conjugate0;}{{ν^{(e)}}_{j}} e_{l} + \overset{&conjugate0;}{{ν^{(μ)}}_{j}} μ_{l} + \overset{&conjugate0;}{{ν^{(τ)}}_{j}} τ_{l}) {W^{+}}_{μ} + (\overset{&conjugate0;}{𝕄_{u, d}} \overset{&conjugate0;}{d_{A, j}} u_{A, l} + \overset{&conjugate0;}{𝕄_{c, d}} \overset{&conjugate0;}{d_{A, j}} c_{A, l} + \overset{&conjugate0;}{𝕄_{t, d}} \overset{&conjugate0;}{d_{A, j}} t_{A, l} + \overset{&conjugate0;}{𝕄_{u, s}} \overset{&conjugate0;}{s_{A, j}} u_{A, l} + \overset{&conjugate0;}{𝕄_{c, s}} \overset{&conjugate0;}{s_{A, j}} c_{A, l} + \overset{&conjugate0;}{𝕄_{t, s}} \overset{&conjugate0;}{s_{A, j}} t_{A, l} + \overset{&conjugate0;}{𝕄_{u, b}} \overset{&conjugate0;}{b_{A, j}} u_{A, l} + \overset{&conjugate0;}{𝕄_{c, b}} \overset{&conjugate0;}{b_{A, j}} c_{A, l} + \overset{&conjugate0;}{𝕄_{t, b}} \overset{&conjugate0;}{b_{A, j}} t_{A, l} + \overset{&conjugate0;}{e_{j}} {ν^{(e)}}_{l} + \overset{&conjugate0;}{μ_{j}} {ν^{(μ)}}_{l} + \overset{&conjugate0;}{τ_{j}} {ν^{(τ)}}_{l}) {W^{-}}_{μ}) {(γ_{}^{μ})}_{j, k}}{2}$

(51)

A process at tree level with a positron and electronic neutrino incoming and the antiparticle of the muon and a muon neutrino outgoing after exchanging a W boson

$FeynmanDiagrams (rhs (), incoming = [conjugate (Electron), ElectronNeutrino], outgoing = [conjugate (Muon), MuonNeutrino], numberofloops = 0, diagrams)$

$\frac{- \frac{ⅈ}{4} \overset{&conjugate0;}{{(v_{e})}_{m} ({\vec{P}}_{1})} {(u_{ν^{(e)}})}_{n} ({\vec{P}}_{2}) {(v_{μ})}_{p} ({\vec{P}}_{3}) \overset{&conjugate0;}{{(u_{ν^{(μ)}})}_{q} ({\vec{P}}_{4})} (- \frac{ⅈ \sqrt{2} g__w {(γ_{5})}_{r, p} {(γ_{}^{α})}_{q, r}}{2} - \frac{ⅈ \sqrt{2} g__w {(γ_{}^{α})}_{q, p}}{2}) (- \frac{ⅈ \sqrt{2} g__w {(γ_{5})}_{s, n} {(γ_{}^{ν})}_{m, s}}{2} - \frac{ⅈ \sqrt{2} g__w {(γ_{}^{ν})}_{m, n}}{2}) (- g_{α, ν} + \frac{({P__1}_{ν} + {P__2}_{ν}) ({P__1}_{α} + {P__2}_{α})}{m_{W^{-}}^{2}}) δ (- {P__3}_{}^{β} - {P__4}_{}^{β} + {P__1}_{}^{β} + {P__2}_{}^{β})}{π^{2} (({P__1}_{κ} + {P__2}_{κ}) ({P__1}_{}^{κ} + {P__2}_{}^{κ}) - m_{W^{-}}^{2} + ⅈ ε)}$

(52)

The term $L_{HV}$ of the electroweak Lagrangian contains the interaction between the Higgs and the Z and W bosons

$Lagrangian (electroweak, term = L [HV], interaction)$

$L__HV = (\frac{g__w Φ}{m_{W}} + \frac{{g__w}^{2} Φ^{2}}{4 m_{W}^{2}}) (m_{W}^{2} {W^{+}}_{μ} {W^{-}}_{}^{μ} + \frac{m_{Z}^{2} Z_{μ} Z_{}^{μ}}{2})$

(53)

The probability density at one loop for a process with three Higgs incoming and outgoing: to omit the large - two pages - algebraic result containing Feynman integrals, end the input line with ":", or remove that ending to see and manipulate the integrals. Note this result involves Feynman diagrams with two, three, four, five and six vertices. New in Maple 2023 are the diagrams with more than four vertices.

$FeynmanDiagrams (rhs (), incoming = [HiggsBoson, HiggsBoson, HiggsBoson], outgoing = [HiggsBoson, HiggsBoson, HiggsBoson], numberofloops = 1, diagrams, output = probabilitydensity) &colon;$

New Physics:-Library commands: FactorSum and RemoveCommonFactors

It is sometimes the case that one needs to factor out something without changing or complicating the algebraic structure of each factor. For this purpose, there is a new Physics:-Library command, FactorSum.

Examples

$restart &semi;$

$f ≔ A \sin (x) θ (x) k - \frac{A \sin (x) \cdot θ (x) m ω^{2}}{B + C}$

$f ≔ A \sin (x) θ (x) k - \frac{A \sin (x) θ (x) m ω^{2}}{B + C}$

(54)

We want to extract the factors $A \sin (x) θ (x)$ . The standard approach is to use factor or simplify

$factor (f)$

$\frac{A \sin (x) θ (x) (- m ω^{2} + k B + k C)}{B + C}$

(55)

$simplify (f, size)$

$\frac{A \sin (x) θ (x) (k (B + C) - m ω^{2})}{B + C}$

(56)

Both performed the operation requested, but, in addition, performed additional undesired operations or normalizations. An alternative to accomplish this without changing the algebraic structure of each factor, is

$Physics :- Library :- FactorSum (f)$

$(k - \frac{m ω^{2}}{B + C}) A \sin (x) θ (x)$

(57)

The relevance of this is in that more often than not, textbooks write this expression as in (57), not (56), not (55), and when studying a book or trying to follow its computations getting exact match to what we see on paper, $Physics :- Library :- FactorSum$ is of help in a way that the other commands are not.

Similarly, when we have equations with both sides having a common factor and therefore we want to just cancel out that factor without changing anything else, the routine to use is the new $Physics :- Library :- RemoveCommonFactors$ . For example, in

$= A \sin (x) θ (x)$

$(k - \frac{m ω^{2}}{B + C}) A \sin (x) θ (x) = A \sin (x) θ (x)$

(58)

we see, on the left-hand and right-hand sides, the common factors $A \sin (x) θ (x)$ . One simple way of performing this cancellation of common factors on the left and right hand sides, sufficiently identical to what we do with paper and pencil, is

$Physics :- Library :- RemoveCommonFactors ()$

$k - \frac{m ω^{2}}{B + C} = 1$

(59)

Note that, as shown above, you don't need to load the Physics to use its Library commands

New coordinates labels

Notation is one of the most important things to communicate with others in science. In the context of the Physics package, strong emphasis is put on using textbook notation as much as possible regarding input and output. Still, for historical reasons, when using the Physics package, the labels used to refer to a coordinate system had been a single Capital Letter, as in X, Y, ...It was not possible to use, e.g. $X'$ , or $x$ .

That has changed in Maple 2023: any symbol can be used as a coordinate system label.

Examples

$with (Physics) &colon;$

A lowercase letter can be used as a coordinate's label

$Setup (coordinates = x)$

$[coordinatesystems = \{x\}]$

(60)

In these cases where the coordinate's system label is not a single Capital Letter, say X, the default coordinates automatically associated to the label, instead of $[x1, x2, .. &period;]$ , are now literal subscripts $[x__1, x__2, .. &period;]$ , that is: symbols with a root, in this example x, followed by the double underscore __, then a number. The automatic mapping $x0 \to x4$ is in these cases $x__0 \to x__4$ , also using literal subscripts x__0 and x__4.

Hence

$x__0$

$x__4$

(61)

Everything else works as with the traditional Capital Letter coordinate-system labels. So $x_{}^{μ}$ is a tensor because it is the 4-vector related to the system of coordinates $x = (x__1, x__2, x__3, x__4)$

$type (x [mu], Library :- PhysicsType :- Tensor)$

$true$

(62)

Consequently, $x {[mu]}^{2}$ also receives the automatic reformatting of repeated indices making one of then contravariant

$x {[mu]}^{2}$

$x_{μ} x_{}^{μ}$

(63)

$SumOverRepeatedIndices ()$

$- {x__1}^{2} - {x__2}^{2} - {x__3}^{2} + {x__4}^{2}$

(64)

A typical textbook case example: the label has a prime using the new feature of Maple 2023 that allows for primed variables. For example,

$Coordinates (X')$

$\{X', x\}$

(65)

Note that the coordinates automatically associated to $X'$ also have a prime; that is also as in textbooks. One can as well use $x$ and $x'$ at the same time to describe transformations of coordinates between two reference systems

$Coordinates (x')$

$\{x, x'\}$

(66)

(in the above, the previous label X' is automatically removed to reuse the coordinates $(x'__1, x'__2, x'__3, x'__4)$ with the new label x' )

The case of a coordinate system where the label has a prime is also special. For example, suppose you set a coordinate system without a prime, as in

$Coordinates (X = cartesian)$

$\{X, x, x'\}$

(67)

If you now set the label $X'$ , because the label $X$ is already set, the coordinates associated to $X'$ are those of $X$ with a prime

$Coordinates (X')$

$\{X, X', x, x'\}$

(68)

Finally, it is now also possible to use literalsubscripts with any label. For example, for a Capital Letter label, say instead of $Y = (y1, y2, y3, y4)$ , you can set Y with literal subscripts:

$Coordinates (Y, literalsubscripts)$

$\{X, X', Y, x, x'\}$

(69)

Summarizing:

1. The label can now be any symbol; particularly useful is the case of a lowercase letter as a label: we can now represent position vectors in tensorial notation using a more natural textbook-notation (lowercase).

2. We can now use the traditional x and x', or X and X' notations, at the same time, as frequently seen in textbooks when presenting transformations of coordinates

3. The coordinates automatically associated to these new coordinate labels use literal subscripted letters as in $x_{1}, x__2, .. &period;$ , a textbook-like typesetting, instead of the old, not so good looking, $x1, x2, .. &period;$

4. It is now possible (optional) to use literal subscripts also with Capital Letter coordinate-system labels

5. Coordinate-system labels with a prime are treated in a special way when the same label without a prime already exists, by just concatenating a prime to the coordinates of the label without prime.

Documentation advanced examples

One of the most important parts of the Physics project is its documentation; the illustration of the use of the package in different scenarios. The three relevant help pages for that are

•	The Physics,Examples

•	The Physics,Tensors

•	The Physics,Updates

For Maple 2021, the first of these pages got extended with four sections covering new material in Vector Analysis, Mechanics and Classical Field Theory.
1. Vectors in Spherical Coordinates using Tensor Notation
2. The equations of motion in curvilinear coordinates, tensor notation and Coriolis force

3. The EnergyMomentum tensor for the Proca Lagrangian

4. The Gross-Pitaevskii field equations for a quantum system of identical particles"

For Maple 2022, the first of these pages got extended with three new sections:
1. Parametrization of Curves, Surfaces and Volumes
2. Integral Vector Calculus
3. The StandardModel in Particle Physics

For Maple 2023, the first of these pages got extended with three new sections:

1	. Moving charges: the retarded and Liénard-Wiechert potentials, electric and magnetic fields $\vec{E}$ and $\vec{H}$

2	. Deriving the mathematical form of 4D relativistic Lorentz transformations

3	. Einstein's principle of relativity

	See Also
	Index of New Maple 2023 Features, Physics , Computer Algebra for Theoretical Physics, The Physics project, The Physics Updates

Download Help Document

Maple

Maple Add-Ons

MapleSim

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim

Maple

Powerful math software that is easy to use

Maple Add-Ons

MapleSim

Advanced System Level Modeling

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim