liquid sodium MD, r_c = 1.971 bohr (temperature derivatives from r_c = 1.971), N = 2000, 30 ps per run after 6 ps equilibration
  left out (still crystalline, never melted): na_1.971_1.0_346_1.pkl
  left out (still crystalline, never melted): na_1.971_1.0_346_1.pkl
  density x0.98, target 371 K: mean T  370.5 K  P = -0.1798 +- 0.0001 GPa  E = -220.6161 mHa
  density x1.00, target 371 K: mean T  373.0 K  P = -0.0584 +- 0.0002 GPa  E = -220.6871 mHa
  density x1.00, target 396 K: mean T  397.1 K  P = -0.0321 +- 0.0001 GPa  E = -220.4204 mHa
  density x1.02, target 371 K: mean T  375.2 K  P = +0.0688 +- 0.0002 GPa  E = -220.7471 mHa
  cut-off (impulsive) pressure added at density x0.98: +0.0338 GPa
  cut-off (impulsive) pressure added at density x1.00: +0.0400 GPa
  cut-off (impulsive) pressure added at density x1.00: +0.0400 GPa
  cut-off (impulsive) pressure added at density x1.02: +0.0463 GPa
  temperatures used: 373.0, 397.1 K;  (dP/dT)_V = 1.0925 +- 0.0275 MPa/K;  c_v = 3.498 +- 0.083 k_B per ion (measured, Fink & Leibowitz: 3.49)
  isothermal B_T = 6.090 +- 0.010 GPa;  adiabatic B_S = 6.470 GPa;  Cp/Cv = 1.0623
  thermodynamic route: speed of sound c = 2641.9 m/s at T = 373.0 K, pressure -0.058 GPa
  structural route: S(k) at k = 0.076, 0.108, 0.132, 0.153 bohr^-1: 0.0115, 0.0123, 0.0132, 0.0112; S(0) = 0.0121
    -> isothermal B from the structure alone = 10.36 GPa (thermodynamic 6.09); a pulse carried by the pair forces alone would run at about 3445 m/s
  measured: S(0) = n kT beta_T = 0.0235 (beta_T from Fink & Leibowitz-type data, 1.886e-10 /Pa)
BEST: c = 2641.9 +- 4.4 m/s
