Beam Moment Model Overview

• Beam Moment Model is a family of strategies that use distinct moments—such as bending, velocity, image, covariance, or slice statistics—as primary observables across various disciplines.
• Researchers reconstruct internal moments via equilibrium, constitutive, and kinematic relations in structural mechanics, while employing multibeam or beam-corrected approaches in plasma and radio astronomy.
• The model’s flexibility allows for reduced, interpretable representations that directly link experimental measurements to domain-specific moment hierarchies, enhancing prediction accuracy.

Searching arXiv for the cited papers and related uses of “Beam Moment Model” across domains. arXiv search query: Beam Moment Model (Jirasek et al., 2024, Goldman et al., 2020, Trapman et al., 12 Jun 2025, Serry et al., 2023, Gutman et al., 2021, Ferreira et al., 14 May 2026, Mitchell, 2015) “Beam Moment Model” is used in several non-equivalent ways in recent literature. In structural mechanics it commonly denotes a formulation in which internal bending moments are reconstructed from equilibrium, constitutive, and kinematic relations; in plasma kinetics it denotes multibeam velocity-moment decompositions; in radio astronomy it refers to beam-corrected fitting of moment-zero maps; in paraxial optics it denotes second- and fourth-moment descriptions of beam propagation; and in accelerator physics it denotes slice-moment decompositions of projected emittance (Jirasek et al., 2024, Goldman et al., 2020, Trapman et al., 12 Jun 2025, Garnier et al., 2022, Ferreira et al., 14 May 2026, Mitchell, 2015). The term therefore does not identify a single canonical theory, but a family of moment-based models in which “moment” may mean bending moment, velocity moment, image moment, covariance moment, or slice statistical moment.

1. Terminological scope and disciplinary meanings

The modern usage of the expression is strongly context-dependent. In beam and arch mechanics, the central object is the internal bending moment $M$ and its coupling to curvature, shear, distributed loading, cracking, thermal fields, or boundary feedback. In kinetic and plasma applications, the relevant moments are velocity moments of a distribution $f(\mathbf{v})$. In astronomical image analysis, “moment” refers to moment-zero maps, i.e. velocity-integrated intensity maps. In paraxial optics and accelerator theory, the core objects are second-order covariance moments and their propagation laws.

Domain	Meaning of “beam moment”	Representative paper
Structural mechanics	Internal bending moment, shear, curvature, stability	(Jirasek et al., 2024)
Cracked beams and arches	Moment continuity and slope-jump laws	(Gutman et al., 2021)
Plasma kinetics	Standard and multibeam energy moments	(Goldman et al., 2020)
Protoplanetary disks	Beam-corrected fitting of moment-zero maps	(Trapman et al., 12 Jun 2025)
Paraxial optics	Second- and fourth-order covariance moments	(Garnier et al., 2022, Ferreira et al., 14 May 2026)
Accelerator physics	Slice moments of projected emittance	(Mitchell, 2015)

This multiplicity suggests that “Beam Moment Model” functions less as a standardized formal label than as a shorthand for any framework that elevates moments to the primary state variables or observables.

2. Geometrically exact beam elements and explicit bending-moment reconstruction

A particularly explicit structural interpretation appears in the two-dimensional shear-flexible geometrically exact beam element of “Shear-flexible geometrically exact beam element based on finite differences” (Jirasek et al., 2024). The formulation extends earlier work by including shear distortion and distributed forces and moments acting along the beam. Its governing strategy is flexibility-based: equilibrium is written in integrated resultant form, sectional constitutive equations are inverted into compliance form, and kinematics then map strains back to spatial derivatives of the centerline and sectional inclination.

For centerline coordinates $x_s(\xi), z_s(\xi)$ and sectional inclination $\varphi(\xi)$, the bending moment is common to the Reissner and Ziegler formulations: $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$ with

$$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$

This expression is notable because the distributed-load contribution is accumulated along the current geometry rather than appended as a purely undeformed-span correction. Curvature then follows from the compliance law $\kappa=M/(EI)$, so that moment drives the rotation field through $\varphi'=\kappa$.

The resulting first-order system is discretized by central finite differences on an element-level grid. Because $x_s'$ and $z_s'$ depend only on midpoint rotations, whereas $f(\mathbf{v})$0 depends only on grid-point coordinates, the scheme remains explicit even with central differences. The boundary-value problem is converted into an initial-value problem by a shooting method with unknown left-end generalized forces $f(\mathbf{v})$1. The paper reports high accuracy with quadratic convergence when the spatial discretization is refined, easy modeling of variable stiffness along the element such as rigid joint offsets, and efficient and accurate characterization of buckling and post-buckling behavior (Jirasek et al., 2024).

The same paper uses the moment-based formulation to compare Reissner and Ziegler sectional models. In the Reissner case, stability connections are drawn to Haringx theory; in the Ziegler case, to Engesser theory. This is significant because the moment law is identical at the resultant level, while the distinction enters through the strain measures and compliance mapping.

3. Other structural beam-moment formulations

Several other structural works use “beam moment” in a more classical but still technically distinct sense. In “Equations of motion for cracked beams and shallow arches,” cracks are modeled by massless rotational springs, and the physical bending moment is

$f(\mathbf{v})$2

with moment continuity across cracks and slope jump governed by

$f(\mathbf{v})$3

equivalently $f(\mathbf{v})$4 (Gutman et al., 2021). A specially designed operator $f(\mathbf{v})$5 “absorbs” the crack boundary conditions into the weak form, so the beam moment model becomes an operator-theoretic representation of piecewise Euler–Bernoulli bending with distributional crack effects.

In “An Euler-Bernoulli-Type Beam Model of the Vocal Folds for Describing Curved and Incomplete Glottal Closure Patterns,” the relevant moment is a composite sectional moment,

$f(\mathbf{v})$6

where $f(\mathbf{v})$7 is a composite bending stiffness and $f(\mathbf{v})$8 is a nominal offset moment induced by layer forces (Serry et al., 2023). The paper distinguishes internal moments generated by thyroarytenoid activation from reactive moments generated by the anterior rotational spring. Bowed, concave, and hourglass glottal closure patterns are then interpreted as consequences of the sign and spatial distribution of these competing moments.

In “Some comments on Gao beam model,” a nonlinear Euler–Bernoulli-type formulation is re-examined under vertical and axial loading. The summary states that the cubic slope term implies an effective bending moment

$f(\mathbf{v})$9

and argues that a corrected choice of the integration constant yields the consistent governing equation

$x_s(\xi), z_s(\xi)$0

with $x_s(\xi), z_s(\xi)$1 (Machalová et al., 2019). The central issue is not merely constitutive form, but the consistency of the moment model across pure bending and pure buckling limits.

A thermoelastic generalization appears in “A Beam Theory Consistent with Three-Dimensional Thermo-Elasticity,” where the bending moment in a reduced one-dimensional model is

$x_s(\xi), z_s(\xi)$2

so thermal deviation across the section contributes directly to bending (Favata, 2013). Here the beam moment is derived from a three-dimensional linear thermoelastic parent theory via the Principle of Virtual Powers and cross-sectional reduction, with reactive stresses and hyper-stresses enforcing the internal constraints.

A boundary-control interpretation appears in “Non uniform stability for the Timoshenko beam with tip load,” where the tip rotary inertia is coupled to the boundary bending term $x_s(\xi), z_s(\xi)$3, and the dynamic boundary condition

$x_s(\xi), z_s(\xi)$4

acts as a boundary moment-feedback law (Mercier et al., 2015). Under the equal-speed condition $x_s(\xi), z_s(\xi)$5, the energy decays polynomially rather than exponentially.

Finally, “Image-Based Structural Analysis Using Computer Vision and LLMs: PhotoBeamSolver” operationalizes classical beam moments in an automated solver for idealized drawings. Its mechanical core is entirely conventional: $x_s(\xi), z_s(\xi)$6 with closed-form moment diagrams assembled from detected supports and loads (Fernando, 22 Mar 2026). Here the novelty is not a new moment law, but automatic extraction of the beam model from a drawing.

4. Multibeam energy moments in plasma kinetics

In plasma physics, the phrase is used in a different sense altogether. “Multi-beam Energy Moments of Multibeam Particle Velocity Distributions” defines a multibeam moment framework for a distribution written as a sum of disjoint beams,

$x_s(\xi), z_s(\xi)$7

with beamwise densities $x_s(\xi), z_s(\xi)$8, bulk velocities $x_s(\xi), z_s(\xi)$9, and pressure tensors $\varphi(\xi)$0 (Goldman et al., 2020). The key distinction is between standard moments, computed relative to a single global mean velocity $\varphi(\xi)$1, and multibeam moments, computed beamwise and then summed.

The undecomposed energy density,

$\varphi(\xi)$2

is invariant under either decomposition, but the standard incoherent moments can contain false contributions arising purely from inter-beam drifts. The false thermal energy is defined as

$\varphi(\xi)$3

Likewise, the false pressure tensor is

$\varphi(\xi)$4

The canonical example is a pair of equal and opposite cold beams,

$\varphi(\xi)$5

The global bulk velocity is then $\varphi(\xi)$6, so the standard thermal energy is nonzero even though each beam is individually cold. The multibeam thermal energy vanishes, and the undecomposed energy is entirely coherent in the multibeam sense (Goldman et al., 2020). This directly motivates the multibeam formalism.

The paper also gives a practical workflow for measurements from NASA’s Magnetospheric Multi-Scale Mission and simulations: compute standard moments, identify beams by mixture modeling or clustering, extract $\varphi(\xi)$7 and $\varphi(\xi)$8, compute coherent multibeam moments, subtract false parts from standard incoherent moments, and optionally refine with tri-Maxwellian fits. In this context, a “Beam Moment Model” is therefore a diagnostic framework for reallocating apparent thermal content into coherent multibeam motion.

5. Beam-corrected moment-zero modeling in protoplanetary-disk astronomy

In radio astronomy, the term appears in the image-plane modeling of ALMA data. “The ALMA Survey of Gas Evolution of PROtoplanetary Disks (AGE-PRO): XI. Beam-corrected gas disk sizes from fitting $\varphi(\xi)$9 moment zero maps” uses a beam-convolved forward model to recover intrinsic gas-disk sizes from integrated-intensity maps (Trapman et al., 12 Jun 2025).

The moment-zero map is

$$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$0

constructed from $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$1 cubes with a Keplerian velocity mask and a spatially varying noise map $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$2. The fitting weights are

$$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$3

which is intended to compensate approximately for correlated noise within a synthesized beam. The clean beam is modeled as an elliptical Gaussian, and the forward model is

$$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$4

Three azimuthally symmetric intrinsic profiles are considered: Nuker, Sérsic, and Gaussian, although the Gaussian is generally inadequate for CO. Most sources are fitted with a Nuker profile, while two require Sérsic fits. The beam-corrected gas size $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$5 is measured not on the convolved image, but on the unconvolved best-fit model through a curve-of-growth condition,

$$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$6

This operationally defines the beam moment model in this paper: fit in the beam-convolved image plane, measure size in the unconvolved model plane.

The reported gas-to-dust size ratio

$$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$7

spans approximately $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$8 to $$M(\xi)=-M_{ab}+X_{ab}\bigl(z_s(\xi)-z_a\bigr)-Z_{ab}\bigl(x_s(\xi)-x_a\bigr)+M_p(\xi),$$9, with a sample median $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$0. The younger Lupus disks have $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$1, while the older Upper Sco disks have $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$2, and the paper reports no clear increase of $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$3 with age (Trapman et al., 12 Jun 2025). Cloud contamination, azimuthal asymmetry, elevated emitting surfaces, and the locality of the Levenberg–Marquardt fitter are identified as limitations.

6. Optical beam moments: covariance, fourth moments, and second-order equivalence

Two optical works use beam moments in the statistical-optics sense. “Partially Coherent Electromagnetic Beam Propagation in Random Media” studies a white-noise paraxial regime in which the complex wave amplitude satisfies an Itô–Schrödinger equation. Because the noise is Gaussian and Markov in propagation distance $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$4, moment equations close at all orders, and the paper analyzes the second and fourth moments explicitly (Garnier et al., 2022). In the scintillation regime, the intensity covariance is written in terms of the fundamental second-order functions $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$5: $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$6 In the strongly scattering smooth case, the beam radius obeys

$$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$7

with asymptotic deep-propagation laws

$$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$8

Here the beam moment model is a hierarchy of field moments used to characterize spread, coherence loss, and scintillation.

A different but related second-order framework is developed in “Second-order moment equivalence of twisted Gaussian Schell model beams and orbital angular momentum eigenmodes” (Ferreira et al., 14 May 2026). For the phase-space vector $$M_p(\xi)=-\int_0^\xi m(s)\,ds+\int_0^\xi P_x(s)\,z_s'(s)\,ds-\int_0^\xi P_z(s)\,x_s'(s)\,ds.$$9, the covariance matrix

$\kappa=M/(EI)$0

of any cylindrically symmetric coherent OAM eigenmode with azimuthal index $\kappa=M/(EI)$1 takes the universal form

$\kappa=M/(EI)$2

This has the same pattern as a twisted Gaussian Schell-model beam at the zero-curvature plane, yielding the term-by-term mapping

$\kappa=M/(EI)$3

Under any ABCD system represented by a symplectic matrix $\kappa=M/(EI)$4, the moments propagate as

$\kappa=M/(EI)$5

Consequently, any matched OAM and TGSM beams are second-order indistinguishable at every propagation plane and share identical beam-width evolution, far-field divergence, and $\kappa=M/(EI)$6 beam-quality factor (Ferreira et al., 14 May 2026). This is a strict covariance-moment notion of beam equivalence.

7. Slice moments and projected emittance in accelerator physics

In accelerator beam dynamics, “A General Slice Moment Decomposition of RMS Beam Emittance” treats the beam distribution $\kappa=M/(EI)$7 on six-dimensional phase space $\kappa=M/(EI)$8 as a probability density and decomposes projected emittance into slice contributions (Mitchell, 2015). For the horizontal plane, the projected rms emittance is

$\kappa=M/(EI)$9

where $\varphi'=\kappa$0 is the $\varphi'=\kappa$1 horizontal covariance block. With longitudinal slices indexed by $\varphi'=\kappa$2, one defines the slice covariance $\varphi'=\kappa$3, slice centroid $\varphi'=\kappa$4, and slice emittance

$\varphi'=\kappa$5

The law of total covariance gives

$\varphi'=\kappa$6

where $\varphi'=\kappa$7 is the covariance matrix of slice centroids. The central result is the four-term nonnegative decomposition

$\varphi'=\kappa$8

These terms are

$\varphi'=\kappa$9

$x_s'$0

$x_s'$1

$x_s'$2

The interpretation is geometric. $x_s'$3 is the mean intrinsic slice emittance. $x_s'$4 is a mismatch term arising from slice-to-slice variation of Twiss functions or, equivalently, from variation of $x_s'$5. $x_s'$6 is a linear centroid-misalignment term, and $x_s'$7 is a nonlinear centroid-misalignment term. The paper emphasizes that these terms can be reconstructed from slice statistics, which makes the decomposition useful for diagnosing projected emittance growth due to wakefields, chromaticity, time-dependent focusing, or transverse–longitudinal coupling (Mitchell, 2015).

Taken together, these works show that a beam moment model is best understood as a modeling strategy rather than a single formula. Its defining feature is the choice of moment variables—mechanical resultants, velocity moments, intensity moments, covariance matrices, or slice statistics—as the primary carriers of structure. The diversity of current usage also implies a recurring methodological advantage: once the correct moment hierarchy is identified, it often yields reduced models with strong interpretability, explicit propagation laws, and direct access to experimentally observable quantities.