---
title: Extended Quantum Molecular Dynamics (EQMD)
url: https://www.emergentmind.com/topics/extended-quantum-molecular-dynamics-eqmd
type: topic
---

# Extended Quantum Molecular Dynamics (EQMD)

Searching arXiv for recent and foundational EQMD papers relevant to model formulation, applications, and methodological developments.
Extended Quantum Molecular Dynamics (EQMD) is a QMD-like transport framework in which each nucleon is represented by a Gaussian wave packet whose width is itself a dynamical variable rather than a fixed parameter. Across the literature, EQMD is used for stable ground-state nucleus modeling, cluster structure, heavily deformed nuclei, and low- to intermediate-energy nuclear reactions, with particular prominence in studies of giant resonances, clustering phenomena, and structure-sensitive reaction observables [1705.04203] [2405.13652]. Its defining refinements over standard QMD are the use of dynamical wave-packet widths, a phenomenological Pauli potential, frictional cooling for initialization, and subtraction of spurious center-of-mass zero-point motion, which together improve ground-state stability and extend the model’s applicability to collective dynamics and nontrivial nuclear morphologies [2405.13652] [2405.13664].

## 1. Formal structure of the model

In EQMD, the total many-body wave function is taken as a direct product of single-particle Gaussian packets,
\[
\Psi=\prod_i \varphi(\boldsymbol r_i),
\]
with
\[
\varphi(\boldsymbol r_i)=\left(\frac{\nu_i+\nu_i^\ast}{2\pi}\right)^{3/4}
\exp\!\left[-\frac{\nu_i}{2}(\boldsymbol r_i-\boldsymbol R_i)^2+\frac{i}{\hbar}\boldsymbol P_i\cdot \boldsymbol r_i\right],
\]
or, in equivalent notation used in other papers,
\[
v_i=\frac{1}{\lambda_i}+i\delta_i.
\]
Here \(\boldsymbol R_i\) and \(\boldsymbol P_i\) are the coordinate- and momentum-space centroids, while the complex width parameter is dynamical [2206.13921] [1705.04203].

This dynamical width is the main formal distinction from standard QMD emphasized throughout the EQMD literature. Because \(\lambda_i\) and \(\delta_i\) evolve in time, the model can adapt packet localization to the nuclear environment rather than imposing a fixed packet scale. The papers repeatedly connect this feature to improved ground-state stability, better reproduction of static properties, and more realistic collective motion, especially in deformed nuclei and clusterized systems [1705.04203] [2405.13652].

The Hamiltonian expectation value is written as kinetic terms plus an interaction part,
\[
H_{\rm int}=H_{\rm Skyrme}+H_{\rm Coulomb}+H_{\rm Symmetry}+H_{\rm Pauli},
\]
with additional subtraction of spurious center-of-mass zero-point motion. In several implementations, the kinetic energy explicitly contains both the centroid contribution \(\mathbf P_i^2/2m\) and a width-dependent quantum term. This distinction becomes important in applications where internal packet momentum spread affects observables, such as direct hard-photon production [2405.13652] [2008.03514].

## 2. Initialization, frictional cooling, and effective interactions

EQMD ground states are prepared by evolving the variational parameters with damped equations of motion derived from the time-dependent variational principle,
\[
\dot{\mathbf R}_i=\frac{\partial H}{\partial \mathbf P_i}+\mu_{\mathbf R}\frac{\partial H}{\partial \mathbf R_i},\qquad
\dot{\mathbf P}_i=-\frac{\partial H}{\partial \mathbf R_i}+\mu_{\mathbf P}\frac{\partial H}{\partial \mathbf P_i},
\]
\[
\frac{3\hbar}{4}\dot{\lambda}_i=-\frac{\partial H}{\partial \delta_i}+\mu_\lambda\frac{\partial H}{\partial \lambda_i},\qquad
\frac{3\hbar}{4}\dot{\delta}_i=\frac{\partial H}{\partial \lambda_i}+\mu_\delta\frac{\partial H}{\partial \delta_i}.
\]
During cooling, the friction coefficients are negative and drive the system toward a local minimum; during real-time dynamics they are set to zero so that energy is conserved [1705.04203] [2206.13921].

This initialization strategy is central to essentially all EQMD applications. In studies of deformed giant dipole resonances, the cooled states are selected so that the resulting deformation agrees with experiment, enabling axis-resolved resonance analysis in Nd and Sm isotopes [1705.04203]. In broad-nuclide surveys of bubble-like morphologies, frictional cooling is used to obtain relaxed low-energy cluster configurations over the AME2020 chart, with the stopping condition defined by vanishing time derivatives of \(\mathbf R_i\), \(\mathbf P_i\), \(\lambda_i\), and \(\delta_i\) [2605.24676].

The effective interaction is usually decomposed into Skyrme, Coulomb, symmetry, and Pauli terms. In the older EQMD setup used for many GDR calculations, the Skyrme term employs a soft equation of state with
\[
\alpha=-124.3~\text{MeV},\qquad \beta=70.5~\text{MeV},\qquad \gamma=2,
\]
and the symmetry energy coefficient is often taken as \(C_S=25\) MeV or varied to test sensitivity [2206.13921] [1705.04203]. The Pauli potential is introduced phenomenologically to suppress close approach of identical nucleons in phase space; in the deformed-nucleus GDR study it is written as
\[
H_{\rm Pauli}=\frac{c_p}{2}\sum_i (f_i-f_0)^\mu \theta(f_i-f_0),
\]
with \(c_p=15\) MeV, \(f_0=1.0\), and \(\mu=1.3\) [1705.04203].

A recurring practical point is that EQMD’s stability depends not only on the interaction choice but also on how density-dependent terms are propagated. This becomes especially significant once one departs from the original \(\gamma=2\) setup and attempts to use softer, non-integer-exponent equations of state [2405.13664].

## 3. Collective dipole and monopole dynamics

Giant resonances are among the most developed EQMD application domains. In GDR studies, the standard isovector dipole observables are
\[
D_G(t)=\frac{NZ}{A}\big[R_Z(t)-R_N(t)\big],
\qquad
K_G(t)=\frac{NZ}{A\hbar}\left[\frac{P_Z(t)}{Z}-\frac{P_N(t)}{N}\right],
\]
and the \(\gamma\)-ray spectrum is extracted from the Fourier transform of the second derivative of the dipole moment [1705.04203] [2206.13921]. In proton-capture calculations, the compound nucleus is formed dynamically and evolved for about \(700~\text{fm}/c\), with \(15{,}000\) events per energy point and \(b=0\) fm to avoid angular-momentum effects [1706.01012]. In deformed heavy nuclei, the GDR is decomposed into major-axis and minor-axis components, and the splitting satisfies
\[
\frac{\Delta E_m}{\bar E_m}\propto \beta_2,
\]
with fitted relations given separately for Sm and Nd chains [1705.04203].

The same formalism has been used to address deformation, temperature, symmetry energy, angular momentum, and external electromagnetic fields. In deformed Nd and Sm isotopes, EQMD reproduces the transition from a single-hump spectrum in nearly spherical nuclei to double-hump spectra in strongly prolate systems, and for \(^{150}\)Nd the calculation is reported to be perfectly consistent with experiment when \(E_{sym}=32\) MeV [1705.04203]. In proton capture reactions such as \(^{11}\)B\((p,\gamma)^{12}\)C, \(^{27}\)Al\((p,\gamma)^{28}\)Si, \(^{39}\)K\((p,\gamma)^{40}\)Ca, and \(^{67}\)Co\((p,\gamma)^{68}\)Ni, the model yields systematic trends in peak energy, strength, and FWHM, including a temperature-dependent width that is almost constant for \(T<1.5\) MeV, rises sharply for \(1.5\lesssim T\lesssim 3.5\) MeV, and supports saturation above \(4\) MeV [1706.01012].

EQMD has also been used for Coulomb-excited \(^{40}\)Ca in ultra-peripheral \(^{16}\)O+\(^{40}\)Ca collisions with impact parameter \(b=7.13\) fm and beam energies from 50 to 500 MeV/nucleon. In that study, inclusion of the collision-generated magnetic field enhances the GDR peak energy, strength, and width. The magnetic field peaks at about \(2.75\times 10^{15}\) G around \(60~\text{fm}/c\) for 100 MeV/nucleon and remains significant up to roughly \(120~\text{fm}/c\); its effect on GDR broadening is interpreted as a combined consequence of temperature and magnetic-field-enhanced angular momentum rather than heating alone [2206.13921]. The paper further notes that in the low-energy region the observed “GDR” likely contains mixed vibrational and rotational character rather than a pure vibrational dipole mode [2206.13921].

Monopole dynamics provide a complementary benchmark. In the 2024 soft-EOS study, the isoscalar giant monopole resonance of \(^{208}\)Pb is used to test newly introduced Skyrme-based interactions. The original EQMD gives a peak energy around \(20.55\) MeV, far above the RCNP experimental value, which is interpreted as evidence that the original mean field is too stiff. Among the new parameterizations, SkP with \(K_\infty=200\) MeV gives the best agreement for \(^{208}\)Pb, while heavier nuclei are described well overall and lighter nuclei are overestimated by about \(1\) MeV [2405.13652].

## 4. Equation of state, collision terms, and propagation improvements

A major development in recent EQMD work is the relaxation of the model’s original EOS restrictions. One 2024 study replaces the schematic Skyrme-like interaction by a standard Skyrme energy density functional including bulk, gradient, symmetry, and Coulomb terms, and introduces SkP, SkT1, and SKXce parameter sets with incompressibilities \(K_\infty=200\), 236, and 268 MeV, respectively [2405.13652]. Another 2024 study addresses the numerical instability that arises when the density exponent \(\gamma\) is non-integer. It replaces the standard approximate treatment of \(\rho^{\gamma+1}\) by a Monte Carlo integral method and shows that, for \(\gamma=2\), the new method reproduces the original analytical EQMD result, while for non-integer \(\gamma\) it remains stable where the old approximation diverges or leads to non-positive width parameters [2405.13664].

These two developments are complementary. The Skyrme-EDF reformulation provides softer incompressibility with long-time ground-state stability over \(1000~\text{fm}/c\) for \(^{208}\)Pb [2405.13652], while the Monte Carlo propagation method makes it practical to use softer EOS parameterizations in IVGDR calculations without the initialization failures associated with the older approximation [2405.13664]. This suggests a broader transition of EQMD from a historically stiff and limited interaction model toward a more flexible transport framework.

Binary-collision modeling has also been revised. In the original EQMD collision prescription, collisions are handled geometrically with \(d_{\rm coll}=2.0\) fm and
\[
\sigma_{NN}^{\rm EQMD}=\frac{100~{\rm mb}}{1+\epsilon/200~{\rm MeV}}.
\]
A 2025 study of \(^{208}\)Pb GDR replaces this by a stochastic collision treatment based on Gaussian-overlap probabilities and free elastic cross sections from Cugnon et al. [2508.17432]. With the stochastic term, the GDR width becomes strongly dependent on the \(NN\) cross section, and the best agreement with evaluated data is obtained with
\[
c_s=33.2~{\rm MeV},\qquad a=0.57,
\]
reproducing
\[
E_\gamma = 13.38\pm 0.03~{\rm MeV},\qquad \Gamma=4.08\pm 0.09~{\rm MeV}.
\]
The same work concludes that a significant in-medium reduction of free \(NN\) elastic cross sections is needed to reproduce the width [2508.17432].

A different transport refinement appears in direct hard-photon calculations. Because EQMD packets carry intrinsic momentum spread, using only centroid momenta in \(pn\to pn\gamma\) collisions misses a substantial fraction of the available kinetic energy. The 2020 study remedies this by sampling the internal momentum distribution of the wave packet for the first photon-producing \(pn\) collision and modifying Pauli blocking accordingly, improving the yield, inverse slope, and angular distribution for \(^{14}\)N+\(^{12}\)C at \(E/A=20\), 30, and 40 MeV [2008.03514]. The paper’s central point is that the width-dependent quantum kinetic term is not a formal detail but directly affects inelastic observables [2008.03514].

## 5. Clustering, short-range correlations, and nuclear morphology

Beyond giant resonances, EQMD is widely used because it can generate and preserve cluster-sensitive configurations. The 2021 short-range-correlation study modifies EQMD by adding a repulsive short-range term
\[
U_1=\frac{C_{U_1}}{(r+r_0)^{p_1}},
\]
motivated by the absence of a sufficiently strong repulsive core in the original interaction [2112.10348]. In \(^{12}\)C, strengthening the short-range repulsion enhances the high-momentum tail of the nucleon momentum distribution, shifts RMS radii and binding energies, and weakens the emitted proton-pair momentum correlation function calculated with the Lednicky–Lyuboshitz method [2112.10348]. The same work stresses that this remains a simplified SRC picture because explicit tensor forces, spin, and isospin dynamics are absent [2112.10348].

In proton-induced reactions on \(^{12}\)C, EQMD has been coupled to GEMINI to examine how initial structure affects fragment production and event information entropy. For \(p+{}^{12}\)C between 5 and 200 MeV/nucleon, the triangular \(3\alpha\) structure with binding energy \(87.1\) MeV exhibits an extra quasi-elastic \(^{12}\)C branch relative to the spherical structure with binding energy \(104.5\) MeV. This excess is most visible around \(E\approx 40\)–\(80\) MeV and peaks near \(E\approx 50\) MeV; it produces a small dent in both fragment information entropy and multiplicity information entropy, leading the authors to propose event information entropy as a probe for \(\alpha\)-cluster structure [2501.18118].

A more global structural application appears in the 2026 survey of bubble-like nuclei across the AME2020 database. There, EQMD relaxed low-energy cluster configurations are classified by the dimensionless \(BHTU\) scheme: \(B\) counts radial-density inflection points and distinguishes droplet (\(B=0\)), bubble (\(B=1\)), and toroidal bubble (\(B=2\)) nuclei; \(H\) measures central density depletion; \(T\) measures relative surface thickness; and \(U\) measures the relative size of the internal low-density region [2605.24676]. Light nuclei are predominantly droplet-like with \(B=0\), \(H=0\), \(T=1\), \(U=0\), while most medium-mass nuclei have \(B=1\), especially near \(^{40}\)Ca and in neutron-rich regions. Toroidal bubble nuclei emerge for \(Z\approx 25\) and become prevalent in heavy systems, and bubble structures are reported to be widespread in the superheavy region [2605.24676].

These applications indicate the niche EQMD occupies within transport modeling: it is not limited to reaction kinematics, but is repeatedly used as a structure-sensitive dynamical generator for deformed, clustered, and hollow configurations. A plausible implication is that the model’s phenomenological ingredients are being used not only to propagate reactions but also to define effective nuclear morphology classes.

## 6. Limitations, interpretive issues, and related generalizations

Several limitations recur across the EQMD literature. The Pauli principle is enforced through a phenomenological Pauli potential rather than exact antisymmetrization, which is repeatedly presented as a practical approximation rather than a first-principles treatment [1705.04203] [2405.13664]. The original interaction set was effectively limited to hard incompressibility, motivating the introduction of softer Skyrme-based parameterizations [2405.13652]. Momentum-dependent interactions are still absent in the 2024 soft-EOS implementation, and the authors explicitly identify that as a remaining limitation [2405.13652].

Collision modeling is another area of active revision. The original geometric collision term is criticized for possible spurious repeated collisions and for the arbitrariness of \(\sigma_{\rm max}\) and \(d_{\rm coll}\); the stochastic reformulation in the 2025 \(^{208}\)Pb study is presented as more suitable when the cross section is large [2508.17432]. Likewise, the hard-photon work shows that centroid-only collision kinematics can be inadequate when intrinsic wave-packet momentum spread materially contributes to the available collision energy [2008.03514].

Interpretive caution also appears in resonance applications. In the magnetic-field study of \(^{40}\)Ca, the broadening of the fitted GDR spectrum cannot be understood from temperature alone because the field also changes angular momentum, and at low beam energy the resonance likely mixes vibrational and rotational character [2206.13921]. In proton-capture studies, the real process includes strong \(\gamma\) transitions that are difficult to reproduce in full detail, and for \(^{67}\)Co\((p,\gamma)^{68}\)Ni the temperature changes very little, suggesting proton capture may not be an ideal way to excite GDR in nuclei with \(A>68\) [1706.01012].

A related methodological direction is represented by a 2022 wave-packet molecular dynamics framework described as very closely related to EQMD in spirit and in mathematical structure, but generalized to anisotropic packets with full tensorial shape dynamics [2211.08560]. That model allows elongation and rotation in arbitrary directions, uses generalized Ewald summation, and reports a \(15\%\) increase in DC conductivity in dense hydrogen relative to isotropic packets [2211.08560]. Although it is not an EQMD paper in the strict nuclear-transport sense, it indicates how the EQMD philosophy of dynamical packet-shape evolution can be extended when anisotropy itself becomes a leading-order physical degree of freedom.

Taken together, these developments define EQMD as a variational Gaussian-packet transport model whose distinctive contribution lies in the simultaneous treatment of dynamical widths, stabilized initialization, phenomenological fermionic repulsion, and structure-sensitive real-time evolution. The model’s most established achievements are in deformation-resolved GDR systematics, soft-EOS extensions, clustering-sensitive reaction studies, and the identification of exotic density morphologies, while its continuing development is concentrated on collision kernels, EOS flexibility, and more complete microscopic content [1705.04203] [2405.13652] [2605.24676].

Source: https://www.emergentmind.com/topics/extended-quantum-molecular-dynamics-eqmd