Energy Decomposition Method Overview

• Energy Decomposition Method is a technique that partitions energy functionals into additive or projected components to enhance numerical stability and physical interpretability.
• It is applied across diverse fields such as low-beta MHD, dual-energy CT, quantum chemistry, and continuum mechanics, demonstrating tailored formulations for each domain.
• The method involves trade-offs like relaxed total-energy conservation or dependence on projection operators, balancing stability with precise energy-channel analysis.

Searching arXiv for recent and foundational uses of “Energy Decomposition Method” across domains. In arXiv usage, the expression “Energy Decomposition Method” denotes a family of domain-specific procedures that partition an energy functional, an interaction energy, a transport quantity, or a measurement model into components chosen for numerical stability, physical interpretability, inverse reconstruction, or scalable approximation. The phrase appears in time-evolving coronal magnetohydrodynamics, dual-energy CT, energy decomposition analysis in electronic structure theory, morphoelasticity and chemo-/poro-mechanics, energy-preserving empirical mode decomposition, exact nonlinear MHD-wave analysis, reaction-coordinate analysis, energy-based fragmentation of the Born–Oppenheimer potential, and decomposition-based distributed energy-resource coordination (Wang et al., 28 Aug 2025, Zhao et al., 2015, Tahmasbi et al., 24 Sep 2025, Chenchiah et al., 2012, Chua et al., 2023, Singh et al., 2015, Raboonik et al., 2024, Li, 2022, Barker et al., 2024, Andrianesis et al., 2021).

1. Scope and core mathematical pattern

Across these literatures, the common operation is to replace a single undifferentiated quantity by additive or projected components that are easier to evolve, interpret, or optimize. In low-$\beta$ coronal MHD, the total energy is written as

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$

and the model advances $E_1$ rather than $E$ to avoid catastrophic cancellation in pressure recovery (Wang et al., 28 Aug 2025). In dual-energy CT, attenuation is represented in a basis form,

$$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$

so that spatial dependence and energy dependence can be separated and later recombined into task-weighted projections (Zhao et al., 2015).

In quantum chemistry, the same phrase usually refers to decomposition of interaction or binding energy into physically interpretable channels. ALMO-EDA uses

$E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$

while periodic EDA and ETS-EDA use

$$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$

In continuum mechanics, the decomposition may act directly at the energy level, as in

$W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$

or as an additive split of Helmholtz free energy,

$$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$

rather than a kinematic multiplicative split (Chenchiah et al., 2012, Chua et al., 2023).

This range of meanings makes the topic intrinsically context-dependent. In some papers the method is primarily a stabilization device, in others a physically motivated analysis of bonding, and in others a reconstruction or approximation framework. The shared structure is decomposition, but the decomposed object, admissible projections, and correctness criteria are determined by the governing equations of each field.

2. Low-$\beta$ MHD and time-evolving coronal simulation

In the coronal-modeling literature, the method is introduced to address the ill-conditioning of pressure recovery when the plasma beta is extremely small, $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$0. In a conservative ideal-MHD scheme that advances total energy density, thermal pressure is recovered by subtracting magnetic and kinetic energies from $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$1. When $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$2 is orders of magnitude smaller than $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$3, discretization error in $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$4 can make the subtraction

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$5

suffer catastrophic cancellation and even yield negative pressures. The COCONUT formulation therefore advances the non-magnetic energy $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$6 and recovers pressure from

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$7

which removes the subtraction of $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$8 in low-$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$9 regions (Wang et al., 28 Aug 2025).

The resulting decomposed-energy equation is not a purely conservative total-energy equation; it includes explicit source-like terms representing magnetic work,

$E_1$0

or equivalently

$E_1$1

For $E_1$2, the right-hand side reduces to the mechanical work done by magnetic stress. In COCONUT, the induction equation remains unchanged, hyperbolic GLM cleaning is retained, and only the advanced energy variable and its source terms are modified (Wang et al., 28 Aug 2025).

The numerical implementation uses a finite volume Godunov method on unstructured, geodesic pyramidal meshes, second-order BDF2 in time with Newton iterations, HLL inviscid fluxes, piecewise linear reconstruction, Venkatakrishnan limiting on density, pressure, and energy, HGLM cleaning with $E_1$3, and positivity-preserving safeguards for $E_1$4 and $E_1$5. Stability in the decomposed-energy equation is further reinforced by blending a small LLF-like dissipation into the energy flux only where $E_1$6 is extremely small (Wang et al., 28 Aug 2025).

Validation in COCONUT is explicit. During solar minimum, decomposed-energy and full-energy solutions are nearly identical: for CR 2073 the average relative differences over the domain are $E_1$7 for density, $E_1$8 for speed, and $E_1$9 for $E$0. For the rising phase CR 2248, the differences are larger but remain modest: $E$1 for density, $E$2 for speed, and $E$3 for $E$4. For the solar-maximum time-evolving run of CR 2296, full-energy COCONUT crashed because of low-$E$5 instabilities, whereas the decomposed-energy version remained robust with near-surface fields exceeding $E$6 G and reproduced streamer and coronal-hole structure seen in EUV and $E$7 observations (Wang et al., 28 Aug 2025).

The main trade-off is conservation. Because $E$8 is evolved with explicit magnetic-work terms, and because HGLM cleaning and low-$E$9 flux dissipation are retained, exact algebraic total-energy conservation is relaxed for stability. The paper presents this as acceptable for thermodynamic coronal models that already include gravity, conduction, radiative losses, and empirical heating, but it also states that advancing $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$0 may remain preferable when strict total-energy conservation is required for fine energy-budget analysis across shocks or resistive reconnection (Wang et al., 28 Aug 2025).

3. Dual-energy CT, material decomposition, and spectral reweighting

In dual-energy CT, “energy decomposition” usually refers to the representation of attenuation in a basis-material model and the subsequent recovery of material-selective images or spectra. One formulation uses dual-energy acquisition to estimate water and iodine basis images, forward-projects them polychromatically to recover the outgoing spectrum,

$$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$1

and then reweights the extracted spectrum to improve soft-tissue contrast. The reweighted projection is

$$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$2

with weights

$$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$3

In the reported phantom study, the extracted spectrum achieved RMSE $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$4 and mean energy difference $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$5 keV relative to an energy-resolved photon-counting detector spectrum; dose-normalized CNR improved from $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$6 to $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$7 for the $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$8 mg/mL iodine insert and from $$\mu(E,\mathbf r)=\sum_m \alpha_m(\mathbf r) f_m(E),$$9 to $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$0 for the $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$1 mg/mL insert (Zhao et al., 2015).

A second line of work replaces ill-conditioned image-domain inversion by a learned nonlinear map $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$2 from low- and high-energy reconstructions to basis images. The dual interactive Wasserstein GAN framework uses two interactive U-Net generators, two WGAN-GP critics, and a selector that alternates training according to the hybrid losses of the two generators. Its generator objective is

$E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$3

with $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$4. On XCAT phantoms, DIWGAN produced the best soft-tissue metrics among the reported methods: for lung soft tissue, RMSE $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$5, PSNR $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$6 dB, and SSIM $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$7; for head soft tissue, RMSE $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$8, PSNR $E_{\rm int}=E_{\rm frz}+E_{\rm pol}+E_{\rm deloc},$9 dB, and SSIM $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$0. The reported ROI-based noise standard deviation reductions averaged $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$1 versus direct inversion, $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$2 versus iterative decomposition, $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$3 versus FCN, and $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$4 versus Butterfly-Net (Shi et al., 2020).

The provided literature also describes a one-step formulation for dual-energy propagation-based phase-contrast CT in which phase retrieval, reconstruction, and material decomposition are performed directly from intensity data using a Fresnel-diffraction forward model, rather than as a pre-reconstruction or post-reconstruction two-step procedure. In that setting, material decomposition is embedded in a joint optimization over basis coefficients $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$5 through

$$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$6

with the propagated intensity fit iteratively to the measurements (Liao et al., 2023).

These CT uses share a specific inverse-problem logic. The decomposition does not merely relabel already reconstructed energy; it makes the energy dependence of attenuation explicit so that spectrum estimation, task-based reweighting, or direct material recovery becomes possible. The limitations are likewise explicit in the cited works: basis sufficiency, spectral calibration, beam hardening, detector nonidealities, and propagation of decomposition error into spectrum estimation or reconstruction all remain central constraints (Zhao et al., 2015, Shi et al., 2020).

4. Quantum chemistry and electronic-structure energy decomposition

In quantum chemistry, the phrase is closely associated with partitioning interaction energy into chemically interpretable terms. In periodic ALMO-EDA for aqueous systems, the interaction energy is decomposed as

$$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$7

Here $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$8 is the frozen-density interaction, $$\Delta E_{\rm bond}=\Delta E_{\rm prep}+\Delta E_{\rm int}+\Delta E_{\rm disp},\qquad \Delta E_{\rm int}=\Delta E_{\rm elstat}+\Delta E_{\rm Pauli}+\Delta E_{\rm orb}.$$9 is polarization under the ALMO constraint, and $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$0 is the additional stabilization from interfragment occupied–virtual mixing. The cited work then learns the dominant donor charge-transfer contributions with an oxygen-centered SOAP descriptor and an MLP $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$1, predicting the first and second strongest donor CT interactions. The single-output model achieved RMSE $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$2 mHartree on a $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$3-sample test set, and the two-output model RMSE $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$4 mHartree, while retaining the locality-based scaling advantage over periodic ALMO-EDA/DFT (Tahmasbi et al., 24 Sep 2025).

Periodic EDA for extended systems adopts a different but related partition. For fragments $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$5 and $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$6 forming $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$7,

$W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$8

The pEDA implementation handles restricted and unrestricted fragments for $W_g(F)=\sum_{k=1}^K W_k(FG_k^{-1}),$9D to $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$0D systems, with or without fragment occupations, and includes reciprocal-space sampling when fragment occupations are not fixed. It was validated on molecular donor–acceptor systems and then applied to surface-adsorbate bonding scenarios including CO on MgO(001), $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$1 on Pd(001) and Cu(001), and CO and $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$2 on Si(001)c(4x2), with reported convergence using TZ2P, suitable integration accuracy, and k-space sampling (Raupach et al., 2015).

ETS-EDA within KS-DFT uses the same electrostatic/Pauli/orbital pattern plus an explicit dispersion term,

$$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$3

and the cited functional-sensitivity study reports that different density functionals produce significant variation in EDA terms, with dispersion correction terms showing the highest variability. At the same time, a machine-learning analysis found the functional label to be the least important feature for bonding-group classification, whereas $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$4 and $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$5 were the most informative. Hydrogen-bonded, covalent, and donor–acceptor classes remained distinguishable despite cross-functional variation (Oestereich et al., 2023).

Two further quantum-chemical uses widen the meaning of decomposition. One is localized operator partitioning for electronic energy transfer, where fragment Hamiltonians are defined by symmetric projector partitioning,

$$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$6

so that subsystem energies and energy currents can be defined consistently under antisymmetrization, with the construction reducing to Förster and Dexter limits in their respective regimes (Khan et al., 2012). The other is energy-based fragmentation of the Born–Oppenheimer potential, where Möbius inversion on posets yields contribution terms

$$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$7

and ML-SUPANOVA extends the decomposition across fragment, electronic-method, and basis-set hierarchies using product-poset Möbius coefficients (Barker et al., 2024).

In this research tradition, energy decomposition is primarily an interpretive and approximation framework. Its central questions are not pressure positivity or image noise, but which channels stabilize a bond, how subsystem energy should be assigned in overlapping quantum systems, and how many-body or multilevel contributions can be combined without double counting.

5. Continuum mechanics: growth, poro-mechanics, and operator-level splits

A mechanically distinct use of the term appears in morphoelasticity. There, the classical multiplicative decomposition $$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$8 is shown to be fundamentally incompatible with shear-resistance except in dilational growth for the crystalline case discussed. The proposed replacement is an energy–deformation decomposition,

$$\psi=\psi_{\rm solid}+\psi_{\rm poro}+\psi_{\rm chem}+\psi_{\rm grad}+\psi_{\rm couple},$$9

or equivalently

$\beta$0

The cited square-lattice theorem states that if a single $\beta$1 exists such that $\beta$2 for all $\beta$3, then necessarily the growth is isotropic dilation, $\beta$4 and $\beta$5. EDD therefore shifts the decomposition from kinematics to mechanical energy and preserves shear-resistance by assigning different growth descriptors to different energetic channels (Chenchiah et al., 2012).

In chemo- and poro-mechanics, the additive split is stated directly at the free-energy level. For a single fluid phase under the affine-deformation assumption,

$\beta$6

This yields

$\beta$7

and

$\beta$8

The argument is that additive energy decomposition is more transparent than a multiplicative inelastic split in poro- and chemo-mechanics because the transported phase has its own energy and stress. The same structure is also presented as advantageous for coupling to phase-field fracture, where only the solid contribution needs to be degraded (Chua et al., 2023).

A third mechanical use is purely algebraic. The multiplying decomposition of stress and strain employs the projection tensors

$\beta$9

or their Voigt-matrix counterparts, to decompose not only $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$00 and $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$01 but also constitutive and compliance operators. With $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$02, $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$03, $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$04, and $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$05, the strain energy splits into volumetric and deviatoric parts in a form directly usable in FEM and BEM. This operator-level decomposition is presented as particularly useful for selective reduced integration and for separating volumetric and shear contributions during assembly (Lee et al., 2012).

These continuum-mechanical variants share a strong formal preference for additive or projected energy structure over a single undifferentiated constitutive object. The decomposed quantities are not measurement signals but stored energy, free energy, or operator action, and the principal motivations are shear-resistance, thermodynamic transparency, and direct compatibility with numerical formulation.

6. Dynamical, signal, and control-theoretic variants

In exact nonlinear ideal-MHD wave analysis, the decomposition is applied to the instantaneous time variation of total energy density. Using the local eigensystem of the ideal-MHD flux Jacobians,

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$06

the total energy-rate is decomposed into contributions from slow, Alfvén, and fast waves plus entropy and field-divergence pseudo-eigenmodes,

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$07

The cited simulations show that the method uniquely identifies component wave modes of a composite wavefield, including mode conversion, and that entropy and divergence pseudo-modes serve as diagnostics of non-adiabaticity and $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$08 control in numerical data (Raboonik et al., 2024).

In empirical mode decomposition, the objective is exact energy preservation rather than physical-wave separation. The energy-preserving EMD algorithms produce linearly independent, non-orthogonal yet energy-preserving components satisfying a Parseval-type identity. For EPIMFs, the decomposition

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$09

is constructed so that $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$10, which guarantees

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$11

A second variant uses reverse-order Gram–Schmidt orthogonalization to obtain orthogonal IMFs while preserving IMF-like properties and exact energy additivity, enabling cleaner Hilbert spectra and more reliable instantaneous frequency estimates (Singh et al., 2015).

Along a reaction coordinate, energy decomposition is defined through reactive-flux weighting. The flux-weighted potential energy profile is

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$12

and the per-coordinate contribution density is

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$13

In plateau regions of reactive flux, the derivative of the flux-weighted energy equals the sum of these per-coordinate terms. The same framework introduces an energy-weighted reactive current and a directional derivative in collective-variable space that is proposed as useful for identifying the reaction coordinate and the minimum free energy path (Li, 2022).

A control-theoretic use appears in distribution-system operation. There the “decomposition method” couples a centralized AC OPF with parallel DER self-dispatch through distribution locational marginal costs. The planner problem is written as

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$14

and solved by alternating a fixed-injection network OPF with proximal DER updates,

$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$15

The reported method preserves feasibility at every iteration, incorporates voltage and ampacity congestion and dynamic transformer degradation into DLMCs, and reaches near-optimal coordination in the case studies (Andrianesis et al., 2021).

These variants extend the expression beyond static energy partitioning. In all four cases, decomposition is used to isolate a dynamically meaningful channel: eigenmodes of an MHD solution, IMF components of a signal, coordinate-level energetic driving along reactive flux, or system-versus-agent structure in AC OPF.

7. Cross-cutting structure, misconceptions, and limitations

A recurrent misconception is that an energy decomposition method is a single transferable algorithm. The cited works show otherwise. In COCONUT, the decomposed variable is a non-magnetic energy introduced to avoid subtractive loss of significance (Wang et al., 28 Aug 2025). In CT, the decomposition is a basis-material inversion coupled to spectrum extraction or learned reconstruction (Zhao et al., 2015, Shi et al., 2020). In ALMO-EDA and pEDA, it is an interaction-energy analysis that separates electrostatics, Pauli repulsion, polarization, charge transfer, orbital relaxation, and dispersion (Tahmasbi et al., 24 Sep 2025, Raupach et al., 2015). In morphoelasticity and poro-mechanics, it is a constitutive reformulation at the free-energy level (Chenchiah et al., 2012, Chua et al., 2023). In signal and wave analysis, it is a decomposition of energy content or energy rate into exactly additive modal channels (Singh et al., 2015, Raboonik et al., 2024).

Another misconception is that decomposition automatically improves all desiderata simultaneously. The record is more specific. COCONUT improves numerical stability in low-$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$16 regions but relaxes strict total-energy conservation (Wang et al., 28 Aug 2025). DIWGAN suppresses noise and beam-hardening artifacts but depends on representative paired training data and on how accurately simulated labels reflect scanner physics (Shi et al., 2020). The ML surrogate for ALMO-EDA scales linearly with the number of molecules but predicts only the dominant donor CT components unless complementary models are added (Tahmasbi et al., 24 Sep 2025). In multilevel fragmentation, combination-consistency depends on using meet-subsemilattice structure; the paper explicitly shows that connected induced subgraphs on cyclic graphs can miscount contributions, motivating convex subgraphs and Möbius-based correction (Barker et al., 2024).

A further cross-cutting point is that decomposition frequently introduces explicit source terms, coupling coefficients, or projection operators that must be handled with care. In the coronal model, magnetic work appears on the right-hand side of the $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$17 equation (Wang et al., 28 Aug 2025). In reactive-flux decomposition, coordinate-wise contributions require stable estimates of $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$18, $$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$19, and current-density weighting (Li, 2022). In localized operator partitioning for electronic energy transfer, subsystem Hamiltonians must be symmetrically partitioned so that fragment energies remain real and additive under antisymmetrization (Khan et al., 2012). In multiplying decomposition for constitutive operators, the projection matrices are singular and therefore suitable for splitting but not for inversion (Lee et al., 2012).

The topic is therefore best understood as a methodological motif rather than a unitary theory. The motif is the explicit partition of an energetically relevant object into channels that align with the numerical, physical, or inferential structure of a problem. What counts as a correct decomposition is domain-specific: pressure positivity in low-$$E=\frac{1}{2}\mathbf{B}^2+E_1,\qquad E_1=\frac{p}{\gamma-1}+\frac{1}{2}\rho \mathbf{v}^2,$$20 MHD, contrast-to-noise improvement and basis sufficiency in CT, chemically interpretable interaction channels in EDA, shear-resistance or thermodynamic transparency in mechanics, exact additivity in signal analysis, or feasibility-preserving coordination in networked optimization.