parse_msd_output
Compute self- or collective-diffusion coefficients by parsing RASPA2 MSD output files from MD simulations.
Instructions
Parse RASPA2 MSD files and compute self- or collective-diffusion coefficients (2-1/3-1/3-2).
Requires a completed NVT-MD or NPT-MD simulation with 'ComputeMSD yes' set. RASPA2 writes MSDSelf_.dat (self-diffusion) and MSDCollective_.dat.
The Einstein relation D = MSD(t) / 6t is fitted to the linear regime (latter 50% of trajectory). Result is given in A²/ps and m²/s.
Typical self-diffusivities in MOFs:
Fast gas (H2, He): 10⁻⁸ – 10⁻⁷ m²/s
CO2, CH4 in wide pores: 10⁻⁹ – 10⁻⁸ m²/s
Slow diffusers (large MOF pores, tight channels): < 10⁻¹¹ m²/s
Args: output_dir: Path to RASPA2 Output directory. molecule: Optional filter by molecule name (e.g. "CO2"). diffusion_type: "self" (MSDSelf) or "collective" (MSDCollective).
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| molecule | No | ||
| output_dir | Yes | ||
| diffusion_type | No | self |