MCPlas, a MATLAB toolbox for reproducible plasma modelling with COMSOL

MCPlas is a transparent and reproducible MATLAB toolbox that automates the generation of fluid-Poisson models for non-thermal plasmas in COMSOL by utilizing structured JSON input data from the LXCat platform, offering advanced electron transport descriptions and validated results for low-pressure argon discharges.

The Gist

Imagine you are trying to bake a very complex, high-tech cake (a plasma discharge) using a super-expensive, professional oven (the software COMSOL). The problem is that the oven's manual is written in a secret code, the recipe is hidden inside a locked box, and every time you want to share your cake with a friend, you have to rewrite the entire recipe from scratch because your friend uses a different brand of oven.

This is the current state of plasma modeling (simulating ionized gas, like in neon signs or fusion reactors). Scientists often struggle to reproduce each other's work because the "recipes" (mathematical equations and chemical reactions) are hidden, messy, or impossible to translate between different computer programs.

Enter MCPlas, a new toolbox described in this paper. Think of MCPlas as a universal translator and automated recipe generator that solves these problems.

Here is a breakdown of what the paper says, using simple analogies:

1. The Problem: The "Black Box" Kitchen

Currently, if a scientist wants to simulate a plasma, they often use commercial software like COMSOL. While powerful, these tools are like "Black Boxes."

• Hidden Recipes: The math behind the simulation is often buried in proprietary files. You can't easily see exactly how the software calculated the result.
• Manual Entry: If you want to add a new chemical reaction (like adding a new spice to the cake), you have to type it in manually, cell by cell. This is slow and prone to typos.
• No Sharing: If you finish your simulation and want to share your data with a colleague using a different software, you can't just send a file. You have to rebuild the whole model from scratch. This violates the modern scientific rule of FAIR (Findable, Accessible, Interoperable, Reusable) data.

2. The Solution: MCPlas (The "Smart Sous-Chef")

MCPlas is a free, open-source tool (a "toolbox") written in a language called MATLAB that talks to COMSOL. It acts as a smart assistant that does three main things:

• The Universal Language (JSON): Instead of writing recipes in secret code, MCPlas uses a standard, open format called JSON. Imagine this as a "Universal Recipe Card." It's a structured list of ingredients and steps that any computer (and any human) can read. It follows the standards of LXCat, a global library of plasma data, ensuring everyone is speaking the same language.
• The Automation (The "Build Bot"): You don't type the equations into COMSOL manually. You feed your "Universal Recipe Card" (the JSON file) into MCPlas. The tool then automatically builds the entire simulation model in COMSOL for you. It sets up the geometry, the math, and the boundaries instantly.
• The Transparency: Because MCPlas writes the code for you, you can see every single equation it used. There are no hidden secrets. If you want to change a boundary condition (like how the plasma hits the wall), you can see exactly where that happens in the code.

3. The "Secret Sauce": Better Physics

The paper highlights that MCPlas isn't just about automation; it also offers better physics than the standard commercial tools.

• The "Drift-Diffusion" Upgrade: Standard tools use a simplified way to calculate how electrons move (like assuming cars on a highway always drive at the same speed). MCPlas offers a more advanced, "novel" way (called DDAn) that accounts for how electrons actually behave, especially near the edges of the plasma.
• The Result: When the authors tested this, they found that the standard tool and MCPlas gave different results near the electrodes (the "walls" of the plasma container). MCPlas's more accurate physics showed that the standard tool was missing some important details about how the plasma behaves in those critical zones.

4. The "Big Cake" Test: Complexity Made Easy

To prove it works, the authors tried two things:

• Simple Cake: They simulated a basic argon plasma (like a neon light). MCPlas built the model perfectly, matching the results of the expensive commercial software.
• Giant Cake: Then, they switched to a massive, complex recipe with 23 different types of particles and 409 different chemical reactions. In the old days, setting this up manually would take weeks and be full of errors. With MCPlas, they just swapped the JSON file, and the tool instantly built the complex model. It handled the complexity effortlessly.

5. The "Cross-Platform" Magic

The most exciting part is Reusability.
The authors took the exact same JSON recipe file they used for MCPlas and fed it into two completely different software programs (PLASIMO and FEDM).

• The Result: All three different software programs produced nearly identical results.
• The Metaphor: It's like writing a recipe in English, sending it to a French chef, a Japanese chef, and an American chef, and having them all bake the exact same cake without needing to translate the instructions. This proves that the data is truly interoperable.

Summary

MCPlas is a game-changer for plasma scientists because it:

1. Opens the Black Box: Makes the math transparent and editable.
2. Saves Time: Automates the boring, error-prone setup of complex models.
3. Improves Accuracy: Uses better physics for electron movement.
4. Connects the World: Allows scientists to share data in a standard format that works across different software, making science more reproducible and collaborative.

In short, it turns plasma modeling from a solitary, error-prone craft into a transparent, automated, and collaborative science.

Technical Summary

1. Problem Statement

The paper addresses critical challenges in the field of low-temperature plasma (LTP) fluid modelling:

• Lack of Reproducibility and Transparency: Commercial software (e.g., COMSOL Plasma Module) and in-house tools often hide governing equations, boundary conditions, and reaction kinetics within proprietary files or ad-hoc scripts. This makes it difficult to reconstruct, verify, or modify published models.
• Manual Implementation Errors: Defining complex Reaction Kinetic Models (RKMs) with numerous species and reactions is time-consuming and prone to human error when entered manually.
• Interoperability Issues: There is no machine-readable community standard for plasma chemistry data. Consequently, the same RKM cannot be reliably reused across different modelling platforms, hindering cross-validation and conflicting with FAIR (Findable, Accessible, Interoperable, Reusable) data principles.
• Limited Access to Advanced Physics: Commercial packages often restrict users from implementing advanced electron transport descriptions (beyond standard drift-diffusion) or custom boundary conditions.

2. Methodology

The authors developed MCPlas, an open-source MATLAB toolbox that automates the generation of equation-based fluid-Poisson models for non-thermal plasmas within the COMSOL Multiphysics environment.

• Core Architecture:

  • Input Standardization: All model-defining data (species properties, RKM, transport/rate coefficients, geometry, and boundary conditions) are stored in structured JSON documents.
  • Schema Validation: These JSON files are validated against schemas aligned with the LXCat platform (for plasma chemistry) and the Plasma-MDS (Plasma Metadata Schema). This ensures data integrity and interoperability.
  • Automation: The toolbox uses the LiveLink™ for MATLAB® module to automatically compile a fully functional COMSOL model from the validated JSON inputs. This eliminates manual setup in the GUI.
  • Code Transparency: The generated model provides explicit access to all implemented equations, boundary conditions, and source terms via the MATLAB source code.

• Advanced Physics Implementation:

  • Electron Transport: MCPlas supports three drift-diffusion approximations:
    1. Conventional (DDAc).
    2. Simplified (DDA53, commonly used in commercial codes).
    3. Novel (DDAn): A higher-order approximation derived from the expansion of the Electron Velocity Distribution Function (EVDF) in Legendre polynomials. This offers improved accuracy at low and atmospheric pressures and is unique to MCPlas.
  • Boundary Conditions: Implements advanced boundary treatments for electron density and energy flux (based on Hagelaar et al.) and surface charge accumulation on dielectrics.
  • Stabilization: Includes techniques like logarithmic transformation of densities and source term stabilization to ensure numerical stability in regions with steep gradients or low concentrations.

• Supported Geometries: 1D (Cartesian), 1p5D (Polar), 2D (Cartesian), and 2p5D (Cylindrical).

3. Key Contributions

• FAIR-Compliant Workflow: MCPlas establishes a transparent, reproducible workflow where model definitions are externalized into machine-readable JSON files, decoupling the physics definition from the simulation software.
• Cross-Platform Interoperability: By adhering to the LXCat JSON schema, the same RKM input files can be used across different software platforms without modification.
• Advanced Physics Accessibility: It democratizes access to advanced electron transport models (DDAn) and custom boundary conditions that are typically locked behind proprietary interfaces in commercial software.
• Automation of Complexity: The toolbox can automatically generate models for highly complex chemistries (e.g., 23 species, 409 reactions) that would be tedious and error-prone to set up manually.

4. Results

The toolbox was validated through two DC and two RF argon glow discharge test cases, comparing MCPlas against the COMSOL Plasma Module (CPM), PLASIMO, and FEDM.

• Verification against CPM:

  • When configured with identical physics (standard drift-diffusion and CPM boundary conditions), MCPlas results matched CPM results perfectly, verifying the correctness of the generated equations.
  • Impact of Advanced Physics: When MCPlas utilized the novel DDAn transport and advanced boundary conditions, significant differences emerged compared to CPM. Specifically, MCPlas predicted lower particle densities and broader cathode-sheath regions, highlighting the physical significance of these advanced formulations.

• Handling Complex RKMs:

  • The toolbox successfully simulated a 23-species argon RKM (including individual 1s/2p states and molecular species) with 409 processes.
  • Comparisons showed that simplified 4-species models failed to capture the detailed dynamics of excited states and molecular ions, demonstrating the necessity of complex chemistries for accurate plasma composition analysis.

• Cross-Platform Reusability:

  • The same JSON input files (for both 4-species and 23-species models) were used in PLASIMO (Finite Volume Method) and FEDM (Finite Element Method).
  • All three tools (MCPlas, PLASIMO, FEDM) produced mutually consistent results (differences < few percent), proving that the JSON schema successfully enables interoperable, cross-platform plasma modelling.

5. Significance

• Reproducibility: MCPlas solves the "black box" problem in plasma modelling by making the entire model definition explicit and editable.
• Scientific Rigor: It facilitates the verification of published models and allows researchers to test the sensitivity of results to specific transport formulations or boundary conditions.
• Community Standardization: By adopting the LXCat JSON format, MCPlas promotes the adoption of FAIR data principles in the LTP community, enabling the sharing and reuse of complex reaction kinetics data across different research groups and software ecosystems.
• Future-Proofing: The separation of data (JSON) from the solver (COMSOL/PLASIMO/FEDM) ensures that models can be easily migrated to new software or updated as physics descriptions improve.