GMX-PME_ERROR(1)

✨ GMX-PME_ERROR ✨

okay bestie, so this is a page from the GROMACS manual β€” it's gmx-pme_error(1), and it is literally the girl who tells you how wrong your PME is being. πŸ’…


πŸ’– NAME πŸ’–

gmx-pme_error β€” estimate the error of using PME with a given input file. that's it, that's the job, no cap.


πŸ¦‹ SYNOPSIS πŸ¦‹

gmx pme_error [-s [<.tpr>]] [-o [<.out>]] [-so [<.tpr>]] [-beta <real>]
              [-[no]tune] [-self <real>] [-seed <int>] [-[no]v]

🌸 DESCRIPTION 🌸

so gmx pme_error estimates the error of the electrostatic forces if you're using the sPME algorithm. it's giving "vibe check your electrostatics" and honestly? iconic.

the flag -tune will determine the splitting parameter such that the error is equally distributed over the real and reciprocal space part. balanced. symmetrical. main character energy. ✨

now here's the tea β˜•: the part of the error that stems from self interaction of the particles is computationally demanding. like, it's a lot. BUT β€” a good approximation is to just use a fraction of the particles for this term, and you can indicate that fraction with the flag -self. slay, honestly. work smarter not harder. πŸŽ€


πŸ‘‘ OPTIONS πŸ‘‘

πŸ“₯ Options to specify input files

  • -s [<.tpr>] (topol.tpr) β€” Portable xdr run input file. it's the input, bestie, we love a portable xdr file here πŸ’–

πŸ“€ Options to specify output files

  • -o [<.out>] (error.out) β€” Generic output file. generic on the outside, holding all your error data on the inside 🫢
  • -so [<.tpr>] (tuned.tpr) β€” Portable xdr run input file. (Optional) β€” yes she's optional, she's the tuned-up version, no pressure ✨

🌈 Other options

  • -beta <real> (-1) β€” if it's positive, it overwrites ewald_beta from the .tpr file with this value. default is -1, which means hands off πŸ”₯
  • -[no]tune (no) β€” tune the splitting parameter such that the error is equally distributed between real and reciprocal space. default is no because we keep things as-is unless you say otherwise πŸ’…
  • -self <real> (1) β€” if it's between 0.0 and 1.0, determine the self interaction error from just this fraction of the charged particles. default is 1 (all of 'em) πŸ“
  • -seed <int> (0) β€” random number seed used for the Monte Carlo algorithm when -self is set to a value between 0.0 and 1.0. default seed is 0, and yes it's giving Monte Carlo so the seed matters fr fr 🎲
  • -[no]v (no) β€” be loud and noisy. default is no, so she's quiet unless you ask her to spill πŸ—£οΈ
⚠️ heads up: -self only does its fraction-of-the-particles thing when it's between 0.0 and 1.0, and -seed only matters when -self is in that range. and -beta only overwrites when it's positive. iykyk 🧚

πŸ”— SEE ALSO πŸ”—

more information about GROMACS is available at <http://www.gromacs.org/>. go read, bestie πŸ“š


2026, GROMACS development team


and that's the whole page, angel 🫢 this lovely lil manual is from GROMACS 2026.4, dated Oct 02, 2026. go tune those splitting parameters, you absolute icon πŸ’…βœ¨

πŸ“¦ where this page comes from