pith. sign in

arxiv: 2510.27446 · v2 · pith:UPCJA7QMnew · submitted 2025-10-31 · ⚛️ physics.plasm-ph · quant-ph

Modeling partially-ionized dense plasma using wavepacket molecular dynamics

Pith reviewed 2026-05-21 20:06 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph quant-ph
keywords dense plasmawave packet molecular dynamicspartially ionized plasmahydrogenfree energy minimizationpath integral Monte Carlocharge state distributionplasma structure
0
0 comments X

The pith

Wave-packet molecular dynamics with explicit bound states produces charge distributions in dense hydrogen plasma that can be compared directly to path-integral Monte Carlo results.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a wave packet molecular dynamics framework for partially-ionized dense plasmas that rests on a chemical model including bound-state wavefunctions explicitly. For hydrogen it obtains self-consistent charge-state distributions by minimizing the free energy following an earlier approach. This construction permits direct benchmarking of static equilibrium properties against path-integral Monte Carlo data. The central aim is to test the model's approximations while tracing the link between ionization and plasma structure. A reader would value a route to these quantities that stays computationally lighter than full quantum simulations.

Core claim

We develop a wave packet molecular dynamics framework for modeling the structural properties of partially-ionized dense plasmas, based on a chemical model that explicitly includes bound state wavefunctions. Using hydrogen as a representative system, we compute self-consistent charge state distributions through free energy minimization, following the approach of Plummer et al. This enables a direct comparison of static equilibrium properties with path integral Monte Carlo data, facilitating an evaluation of the model's underlying approximations and its ability to capture the complex interplay between ionization and structure in dense plasma environments.

What carries the argument

Wave-packet molecular dynamics framework built on a chemical model that includes bound-state wavefunctions explicitly, together with free-energy minimization to determine charge-state distributions.

If this is right

  • The model's approximations can be evaluated by direct comparison of static properties to path-integral Monte Carlo data.
  • The framework reveals how ionization levels couple to the spatial structure of the plasma.
  • Self-consistent charge states are obtained without fixing ionization fractions in advance.
  • The method is shown for hydrogen as a test case representative of partially ionized dense plasma conditions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same construction could be applied to time-dependent or transport properties once the equilibrium version is validated.
  • If the comparisons succeed, the method offers a way to scan wide ranges of density and temperature where full path-integral Monte Carlo remains expensive.
  • Explicit wavefunctions for bound states appear necessary for accurate equilibrium modeling in this regime.
  • The approach might be adapted to multi-species plasmas to examine how impurities alter ionization-structure relations.

Load-bearing premise

The chemical model that explicitly includes bound-state wavefunctions, combined with free-energy minimization following Plummer et al. (2025), accurately represents the equilibrium charge-state distribution and structural properties of partially ionized dense hydrogen.

What would settle it

Marked disagreement between the wave-packet model's computed radial distribution functions or ionization fractions and the corresponding path-integral Monte Carlo results at matching densities and temperatures would show the central claim does not hold.

Figures

Figures reproduced from arXiv: 2510.27446 by Daniel Plummer, Gianluca Gregori, Patrick Hollebon, Pontus Svensson, Sam M. Vinko, Wiktor Jasniak.

Figure 1
Figure 1. Figure 1: FIG. 1. Radial distribution function (RDF) dependence on confinement parameter [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Each step of the ionization calculation for [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Therefore, we proceed with an in-depth analy [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Minimized ionization state ¯z [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Radial distribution functions (RDFs) from the bound-WPMD model for different values of the confinement parameter [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Radial distribution functions (RDFs) at [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Spin-resolved neutral-neutral radial distribution func [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Decomposition of the ion-neutral kernel used inter [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Ensemble average of Coulomb energy at zero cou [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
read the original abstract

We develop a wave packet molecular dynamics framework for modeling the structural properties of partially-ionized dense plasmas, based on a chemical model that explicitly includes bound state wavefunctions. Using hydrogen as a representative system, we compute self-consistent charge state distributions through free energy minimization, following the approach of Plummer et al. [Phys. Rev. E 111, 015204 (2025)]. This enables a direct comparison of static equilibrium properties with path integral Monte Carlo data, facilitating an evaluation of the model's underlying approximations and its ability to capture the complex interplay between ionization and structure in dense plasma environments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a wave-packet molecular dynamics (WPMD) framework for partially ionized dense plasmas. It employs a chemical model that includes bound-state wavefunctions explicitly for hydrogen, obtains self-consistent charge-state distributions via free-energy minimization following the approach of Plummer et al. (Phys. Rev. E 111, 015204, 2025), and uses the resulting static equilibrium properties for direct comparison against path-integral Monte Carlo (PIMC) data to assess the model's approximations and its treatment of the ionization-structure interplay.

Significance. If the central consistency between the static free-energy minimization and the dynamic WPMD propagation holds, the work would provide a valuable bridge between chemical models and explicit wavefunction dynamics for dense-plasma modeling. The explicit inclusion of bound-state wavefunctions and the direct PIMC comparison are strengths that could help evaluate approximations in regimes where ionization and structure are coupled. The extension of the Plummer et al. framework is noted as a technical contribution.

major comments (2)
  1. [§2.2 and §3.1] §2.2 (Free-energy minimization) and §3.1 (Wave-packet propagation): The charge-state distribution is obtained from static minimization of the Plummer free-energy functional that treats bound states via explicit wavefunctions. However, the time-dependent Schrödinger evolution of those same wave packets during MD can produce spreading, overlap, and effective ionization not folded back into the functional. The manuscript does not describe a feedback mechanism or constraint that enforces the minimized distribution throughout the dynamics; without this, the structural observables compared to PIMC are no longer guaranteed to correspond to the claimed self-consistent charge states.
  2. [§4] §4 (Comparison to PIMC): The direct comparison of pair-correlation functions and other structural quantities is presented as validation of the model's ability to capture ionization-structure coupling. Because the charge-state input is fixed from the static minimization, any discrepancy or agreement with PIMC could arise from the missing dynamic feedback rather than from the underlying approximations; the paper should quantify how sensitive the reported structural properties are to small changes in the charge-state distribution.
minor comments (2)
  1. [§2.1] Notation for the wave-packet width parameter is introduced in §2.1 but used without redefinition in later equations; a brief reminder or table of symbols would improve readability.
  2. [Introduction] The abstract states that the framework 'enables a direct comparison' with PIMC, but the manuscript should clarify in the introduction whether the comparison is performed at fixed density and temperature or whether additional thermodynamic constraints are applied.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the major comments point by point below, indicating where revisions will be made to clarify the hybrid static-dynamic nature of the model.

read point-by-point responses
  1. Referee: [§2.2 and §3.1] §2.2 (Free-energy minimization) and §3.1 (Wave-packet propagation): The charge-state distribution is obtained from static minimization of the Plummer free-energy functional that treats bound states via explicit wavefunctions. However, the time-dependent Schrödinger evolution of those same wave packets during MD can produce spreading, overlap, and effective ionization not folded back into the functional. The manuscript does not describe a feedback mechanism or constraint that enforces the minimized distribution throughout the dynamics; without this, the structural observables compared to PIMC are no longer guaranteed to correspond to the claimed self-consistent charge states.

    Authors: We thank the referee for identifying this key aspect of the hybrid approach. The free-energy minimization determines the equilibrium charge-state distribution, which sets the initial wave-packet populations and widths for the WPMD trajectories. The subsequent propagation generates configurations for computing static structural observables such as pair correlations. We agree that the manuscript does not explicitly describe a feedback mechanism to readjust populations during dynamics, and this constitutes an approximation of the chemical model. In the revised manuscript we will add a dedicated paragraph in §3.1 clarifying that the charge states are held fixed at their minimized values for the purpose of the PIMC comparison, together with a brief discussion of the expected validity of this approximation in the density-temperature regime examined. revision: yes

  2. Referee: [§4] §4 (Comparison to PIMC): The direct comparison of pair-correlation functions and other structural quantities is presented as validation of the model's ability to capture ionization-structure coupling. Because the charge-state input is fixed from the static minimization, any discrepancy or agreement with PIMC could arise from the missing dynamic feedback rather than from the underlying approximations; the paper should quantify how sensitive the reported structural properties are to small changes in the charge-state distribution.

    Authors: The referee is correct that the reported structural quantities are computed with a fixed charge-state distribution obtained from the static minimization. To address the concern, we will add a sensitivity study in the revised §4 in which the hydrogen charge fractions are varied by a few percent around the minimized values while keeping all other parameters fixed, and the resulting changes in the pair-correlation functions are quantified. These additional results will demonstrate the robustness of the comparison with PIMC data. revision: yes

Circularity Check

1 steps flagged

Self-citation to prior free-energy minimization provides charge-state input while wave-packet dynamics add independent structural content

specific steps
  1. self citation load bearing [Abstract]
    "we compute self-consistent charge state distributions through free energy minimization, following the approach of Plummer et al. [Phys. Rev. E 111, 015204 (2025)]. This enables a direct comparison of static equilibrium properties with path integral Monte Carlo data"

    The charge-state distribution required for the claimed direct PIMC comparison is obtained solely by adopting the free-energy minimization procedure from a prior publication by the same lead author, without an independent derivation or verification step within this manuscript's wave-packet dynamics.

full rationale

The paper's central claim rests on computing self-consistent charge-state distributions via free-energy minimization from a 2025 paper by the same lead author, then using the new wave-packet MD to model structure and compare to PIMC. This self-citation is load-bearing for the equilibrium properties but does not reduce the entire derivation to a tautology, as the wave-packet propagation and structural observables introduce new elements. No fitted-input-as-prediction or self-definitional loops appear in the provided text. The approach is therefore partially dependent on prior work but retains independent content, warranting a moderate score.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the ledger is therefore minimal and provisional. The model implicitly assumes that a chemical picture with explicit bound states remains valid at the densities considered and that free-energy minimization yields the correct equilibrium charge states.

axioms (1)
  • domain assumption A chemical model that treats bound electrons via explicit wavefunctions plus free-energy minimization accurately describes equilibrium in partially ionized dense hydrogen.
    Invoked when the authors state they follow the Plummer et al. (2025) approach to obtain self-consistent charge-state distributions.

pith-pipeline@v0.9.0 · 5639 in / 1276 out tokens · 33059 ms · 2026-05-21T20:06:22.032417+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  1. [1]

    PREXTREME

    for example. Molecular dynamics simulations were performed at dif- ferent values ofλto sample the potential energy function. For a given ionization state ¯zand value of the coupling parameterλ, 12 independent runs were performed for N= 1024 total electrons. After an initial 25 fs of velocity scaling at full coupling (λ= 1), the coupling parameter was sequ...

  2. [2]

    Once radial distances are sampled, angular co- ordinates are generated from a uniform distribution over the unit sphere to findz

    This enables a fast sampling routine, given that the Gamma distri- bution can be directly sampled with common numerical libraries. Once radial distances are sampled, angular co- ordinates are generated from a uniform distribution over the unit sphere to findz. These are then linearly trans- lated to find samples fory=z+r ik. Appendix E: Self energy Coulom...

  3. [3]

    (E1) holds within statistical error

    As can be seen, Eq. (E1) holds within statistical error

  4. [4]

    Bonitz, T

    M. Bonitz, T. Dornheim, Z. A. Moldabekov, S. Zhang, P. Hamann, H. K¨ ahlert, A. Filinov, K. Ramakrishna, and J. Vorberger, Ab initio simulation of warm dense matter, Physics of Plasmas27, 042710 (2020)

  5. [5]

    Bonitz, J

    M. Bonitz, J. Vorberger, M. Bethkenhagen, M. P. B¨ ohme, D. M. Ceperley, A. Filinov, T. Gawne, F. Graziani, G. Gregori, P. Hamann, S. B. Hansen, M. Holzmann, S. X. Hu, H. K¨ ahlert, V. V. Karasiev, U. Kleinschmidt, L. Kordts, C. Makait, B. Militzer, Z. A. Moldabekov, C. Pierleoni, M. Preising, K. Ramakrishna, R. Redmer, S. Schwalbe, P. Svensson, and T. Do...

  6. [6]

    J. T. Su and W. A. Goddard, Excited electron dynamics modeling of warm dense matter, Phys. Rev. Lett.99, 185003 (2007)

  7. [7]

    Lavrinenko, P

    Y. Lavrinenko, P. R. Levashov, D. V. Minakov, I. V. Mo- rozov, and I. A. Valuev, Equilibrium properties of warm dense deuterium calculated by the wave packet molecular dynamics and density functional theory method, Phys. Rev. E104, 045304 (2021)

  8. [8]

    Jakob, P.-G

    B. Jakob, P.-G. Reinhard, C. Toepffer, and G. Zwick- nagel, Wave packet simulation of dense hydrogen, Phys. Rev. E76, 036406 (2007)

  9. [9]

    Svensson, Y

    P. Svensson, Y. Aziz, T. Dornheim, S. Azadi, P. Holle- bon, A. Skelt, S. M. Vinko, and G. Gregori, Modeling of warm dense hydrogen via explicit real-time electron dy- namics: Dynamic structure factors, Phys. Rev. E110, 055205 (2024)

  10. [10]

    Zwicknagel and T

    G. Zwicknagel and T. Pschiwul, Wpmd simulations of a two-component plasma, Journal of Physics A: Mathe- matical and General39, 4359 (2006)

  11. [11]

    R. A. Davis, W. A. Angermeier, R. K. T. Hermsmeier, and T. G. White, Ion modes in dense ionized plasmas through nonadiabatic molecular dynamics, Phys. Rev. Res.2, 043139 (2020)

  12. [12]

    Klakow, C

    D. Klakow, C. Toepffer, and P. Reinhard, Semiclassi- cal molecular dynamics for strongly coupled Coulomb systems, The Journal of Chemical Physics101, 10766 (1994)

  13. [13]

    Svensson, P

    P. Svensson, P. Hollebon, D. Plummer, S. M. Vinko, and G. Gregori, Modeling of warm dense hydrogen via explicit real-time electron dynamics: Electron transport proper- ties, Phys. Rev. E111, 045208 (2025)

  14. [14]

    I. V. Morozov and I. A. Valuev, Localization constraints in gaussian wave packet molecular dynamics of nonideal plasmas, Journal of Physics A: Mathematical and Theo- retical42, 214044 (2009)

  15. [15]

    W. A. Angermeier, B. S. Scheiner, N. R. Shaffer, and T. G. White, Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen, Philosophical Transactions of the Royal Soci- ety A: Mathematical, Physical and Engineering Sciences 381, 20230034 (2023)

  16. [16]

    Y. Yao, Q. Zeng, K. Chen, D. Kang, Y. Hou, Q. Ma, and J. Dai, Reduced ionic diffusion by the dynamic elec- tron–ion collisions in warm dense hydrogen, Physics of Plasmas28, 012704 (2021)

  17. [17]

    P. E. Grabowski, A review of wave packet molecular dynamics, inFrontiers and Challenges in Warm Dense Matter, edited by F. Graziani, M. P. Desjarlais, R. Red- mer, and S. B. Trickey (Springer International Publish- ing, Cham, 2014) pp. 265–282

  18. [18]

    Ebeling and B

    W. Ebeling and B. Militzer, Quantum molecular dynam- ics of partially ionized plasmas, Physics Letters A226, 298 (1997)

  19. [19]

    Plummer, P

    D. Plummer, P. Svensson, D. O. Gericke, P. Hollebon, S. M. Vinko, and G. Gregori, Ionization calculations using classical molecular dynamics, Phys. Rev. E111, 015204 (2025)

  20. [20]

    Dornheim, S

    T. Dornheim, S. Schwalbe, M. P. B¨ ohme, Z. A. Mold- abekov, J. Vorberger, and P. Tolias, Ab initio path integral monte carlo simulations of warm dense two- 16 component systems without fixed nodes: Structural properties, The Journal of Chemical Physics160, 164111 (2024)

  21. [21]

    Marx and J

    D. Marx and J. Hutter,Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods(Cambridge Uni- versity Press, 2009)

  22. [22]

    Feldmeier and J

    H. Feldmeier and J. Schnack, Molecular dynamics for fermions, Rev. Mod. Phys.72, 655 (2000)

  23. [23]

    Runge and E

    E. Runge and E. K. U. Gross, Density-functional the- ory for time-dependent systems, Phys. Rev. Lett.52, 997 (1984)

  24. [24]

    Kononov, C.-W

    A. Kononov, C.-W. Lee, T. P. dos Santos, B. Robin- son, Y. Yao, Y. Yao, X. Andrade, A. D. Baczewski, E. Constantinescu, A. A. Correa, Y. Kanai, N. Modine, and A. Schleife, Electron dynamics in extended systems within real-time time-dependent density-functional the- ory, MRS Communications12, 1002 (2022)

  25. [25]

    Svensson, T

    P. Svensson, T. Campbell, F. Graziani, Z. Moldabekov, N. Lyu, V. S. Batista, S. Richardson, S. M. Vinko, and G. Gregori, Development of a new quantum trajectory molecular dynamics framework, Philosophical Transac- tions of the Royal Society A381(2023)

  26. [26]

    Ebeling, V

    W. Ebeling, V. E. Fortov, and V. Filinov,Quan- tum Statistics of Dense Gases and Nonideal Plasmas (Springer International Publishing, 2017)

  27. [27]

    W. A. Angermeier and T. G. White, An investigation into the approximations used in wave packet molecular dynamics for the study of warm dense matter, Plasma4, 294 (2021)

  28. [28]

    Klakow, C

    D. Klakow, C. Toepffer, and P.-G. Reinhard, Hydrogen under extreme conditions, Physics Letters A192, 55 (1994)

  29. [29]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Comp. Phys. Comm.271,...

  30. [30]

    Ditchfield, W

    R. Ditchfield, W. J. Hehre, and J. A. Pople, Self- consistent molecular orbital methods. vi. energy opti- mized gaussian atomic orbitals, The Journal of Chemical Physics52, 5001 (1970)

  31. [31]

    D. A. McQuarrie,Quantum chemistry, 2nd ed. (Univer- sity Science, Sausalito, Calif, 2007)

  32. [32]

    Ebeling, A

    W. Ebeling, A. Filinov, M. Bonitz, V. Filinov, and T. Pohl, The method of effective potentials in the quantum-statistical theory of plasmas, Journal of Physics A: Mathematical and General39, 4309 (2006)

  33. [33]

    Knaup, P.-G

    M. Knaup, P.-G. Reinhard, and C. Toepffer, Wave packet molecular dynamics simulations of hydrogen near the transition to a metallic fluid, Contributions to Plasma Physics39, 57 (1999)

  34. [34]

    A. V. Filinov and M. Bonitz, Equation of state of par- tially ionized hydrogen and deuterium plasma revisited, Physical Review E108, 055212 (2023)

  35. [35]

    Z. A. Moldabekov, T. Dornheim, and M. Bonitz, Screen- ing of a test charge in a free-electron gas at warm dense matter and dense non-ideal plasma conditions, Contribu- tions to Plasma Physics62, e202000176 (2022)

  36. [36]

    Frenkel and B

    D. Frenkel and B. Smit, Chapter 7 - free energy calcu- lations, inUnderstanding Molecular Simulation (Second Edition), edited by D. Frenkel and B. Smit (Academic Press, San Diego, 2002) second edition ed., pp. 167–200

  37. [37]

    H. M. Bellenbaum, M. P. B¨ ohme, M. Bonitz, T. D¨ oppner, L. B. Fletcher, T. Gawne, D. Kraus, Z. A. Moldabekov, S. Schwalbe, J. Vorberger, and T. Dornheim, Estimating ionization states and continuum lowering from ab initio path integral monte carlo simulations for warm dense hy- drogen, Phys. Rev. Res.7, 033016 (2025)

  38. [38]

    Bonitz and L

    M. Bonitz and L. Kordts, Ionization potential depres- sion and fermi barrier in warm dense matter: A first- principles approach, Contributions to Plasma Physics , e70001 (2025)

  39. [39]

    Dornheim, T

    T. Dornheim, T. D¨ oppner, P. Tolias, M. P. B¨ ohme, L. B. Fletcher, T. Gawne, F. R. Graziani, D. Kraus, M. J. Mac- Donald, Z. A. Moldabekov, S. Schwalbe, D. O. Gericke, and J. Vorberger, Unraveling electronic correlations in warm dense quantum plasmas, Nature Communications 16, 5103 (2025)

  40. [40]

    K. P. Huber and G. H. Herzberg, Constants of diatomic molecules, inNIST Chemistry WebBook, NIST Standard Reference Database Number 69, edited by P. J. Linstrom and W. G. Mallard (National Institute of Standards and Technology, Gaithersburg, MD, 20899, 2025) data pre- pared by Jean W. Gallagher and Russell D. Johnson, III

  41. [41]

    H. Xiao, A. Jaramillo-Botero, P. L. Theofanis, and W. A. Goddard, Non-adiabatic dynamics modeling framework for materials in extreme conditions, Mechanics of Mate- rials90, 243 (2015), proceedings of the IUTAM Sympo- sium on Micromechanics of Defects in Solids

  42. [42]

    Campbell, P

    T. Campbell, P. Svensson, B. Larder, D. Plummer, S. M. Vinko, and G. Gregori, Molecular dynamics frame- work coupled with smoothed particle hydrodynamics for quantum plasma simulations, Phys. Rev. Res.7, 023286 (2025)

  43. [43]

    Ono and H

    A. Ono and H. Horiuchi, Antisymmetrized molecular dy- namics for heavy ion collisions, Progress in Particle and Nuclear Physics53, 501 (2004)

  44. [44]

    Ono and K

    J. Ono and K. Ando, Semiquantal molecular dynam- ics simulations of hydrogen-bond dynamics in liquid wa- ter using multi-dimensional gaussian wave packets, The Journal of Chemical Physics137, 174503 (2012)

  45. [45]

    H. V. Henderson and S. R. Searle, On deriving the inverse of a sum of matrices, SIAM Review23, 53 (1981)

  46. [46]

    Svensson,Novel Modelling Techniques for Charged Many-Body Systems with Quantum and Relativistic Ef- fects, Ph.D

    P. Svensson,Novel Modelling Techniques for Charged Many-Body Systems with Quantum and Relativistic Ef- fects, Ph.D. thesis, Worcester College, University of Ox- ford, Oxford (2024)

  47. [47]

    W. L. Slattery, G. D. Doolen, and H. E. DeWitt, Im- proved equation of state for the classical one-component plasma, Phys. Rev. A21, 2087 (1980)