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;$

Maple

Maple Add-Ons

MapleSim

MapleSim Add-Ons

システムズエンジニアリング

コンサルティングサービス

オンライン教育製品

教育利用

産業分野

自動車及び航空宇宙

ロボティクス

マシンデザイン
& 産業オートメーション

その他

アプリケーション

製品価格

購入

教育機関向け
学生ライセンス

Maplesoft Elite Maintenance (EMP)

サポート

製品のトレーニング

オンライン製品ヘルプ

Webセミナー
& イベント

出版物

コンテンツ Hub

事例とアプリケーション

コミュニティ

Maplesoftについて

メディアセンター

ユーザコミュニティ

お問合せ

Online Help

All Products Maple MapleSim

Maple

シンプルな操作性の高度数学ソフトウェア

Maple Add-Ons

MapleSim

先進的なシステムレベルモデリング

MapleSim Add-Ons

システムズエンジニアリング

コンサルティングサービス

オンライン教育製品

教育利用

産業分野

自動車及び航空宇宙

ロボティクス

マシンデザイン & 産業オートメーション

その他

アプリケーション

製品価格

購入

教育機関向け 学生ライセンス

Maplesoft Elite Maintenance (EMP)

サポート

製品のトレーニング

オンライン製品ヘルプ

Webセミナー & イベント

出版物

コンテンツ Hub

事例とアプリケーション

コミュニティ

Maplesoftについて

メディアセンター

ユーザコミュニティ

お問合せ

Online Help

All Products Maple MapleSim

マシンデザイン
& 産業オートメーション

教育機関向け
学生ライセンス

Webセミナー
& イベント