het-agi/collideragent/micromegas-calculator

micrOmegas Calculator

Overview

This skill runs micrOmegas for dark matter phenomenology using two Magnus blueprints in sequence:

1. `micromegas-compile` — drops the CalcHEP model into a fresh micrOmegas project, installs the user-supplied main.c, and compiles the executable. The symbolic CalcHEP engine is pre-built in the container, so this step only compiles user code and auto-generated matrix elements.
2. `micromegas-calc` — runs the compiled ./main once, captures its stdout and the structured results.json it writes, and uploads any additional output files.

Both steps execute on the Magnus cloud (see magnus skill).

Prerequisites

CalcHEP model with Z₂-odd marking

micrOmegas identifies dark-matter candidates by the `~` prefix on particle names in the P column of prtcls1.mdl — there is no separate Z₂ quantum number. Every Z₂-odd particle (the DM candidate plus any coannihilation partners) must appear as ~<name> in that column.

Concrete verification step after generate-calchep:

grep -E '\|~[A-Za-z]' path/to/MyModel_CH/prtcls1.mdl

should list every intended Z₂-odd field. Reference examples from shipped micrOmegas projects: SingletDM → ~x1; IDM → ~H3 ~H+ ~X; RDM → ~chi0 ~chi1.

If no ~-prefixed particles appear, the .fr source is still wrong — fix it at authoring time, not by hand-editing .mdl files (the hand-edit gets lost the next time generate-calchep runs). Without the convention, `sortOddParticles` may fail, `darkOmega` may return `NaN`, and every cross section may read zero — yet the blueprint can still see exit code 0. The blueprint detects the known silent-failure markers (Can not compile, Omega=NAN, Omega=nan) and surfaces them as success: false in MAGNUS_RESULT; the caller should still check success before consuming numbers.

Current repo limitation: generate-calchep invokes plain WriteCHOutput[...] and does not pass a dedicated odd-particle list. Therefore the only reliable contract is the generated prtcls1.mdl itself. A FeynRules-side pattern that is verified to work in the current pipeline is to encode the tilde directly in ParticleName (and AntiParticleName for non-self-conjugate fields). Still, treat the exported prtcls1.mdl as the final authority, because export behavior is model-dependent. Likewise, do not assume QuantumNumbers -> {Z2 -> -1} by itself is enough: custom quantum numbers need explicit model-level declarations, and the CalcHEP writer may still not translate them into micrOmegas odd-particle names. Regenerate and re-check prtcls1.mdl until the exported names are correct.

main.c protocol

The user-supplied main.c is responsible for:

1. Setting parameters (either by assignValW / assignValC calls, or by calling slhaRead on the SLHA file that micromegas-calc passes as argv[1])
2. Calling whichever subset of micrOmegas API it needs (sortOddParticles, darkOmega, nucleonAmplitudes, calcSpectrum, …)
3. Writing structured results to `results.json` in the working directory. The blueprint parses this file and surfaces it in MAGNUS_RESULT.

6.3.0 API note: read the dark-matter name from CDM[1] (a char** array). The legacy globals CDM1/CDM2 are still declared in micromegas.h for source compatibility but sortOddParticles no longer fills them — they read as NULL at runtime. Mass is McdmN[1] (still works as before). Same shift applies to nucleonAmplitudes(CDM[1], …).

A minimal skeleton:

#include "micromegas.h"
#include "micromegas_aux.h"
#include <stdio.h>
 
int main(int argc, char** argv) {
    if (argc > 1 && slhaRead(argv[1], 0)) {
        fprintf(stderr, "Failed to read SLHA file: %s\n", argv[1]);
        return 2;
    }
 
    if (sortOddParticles(NULL)) { return 1; }
    if (!CDM[1]) {
        fprintf(stderr, "CDM[1] is null; check ~-prefixed names in prtcls1.mdl\n");
        return 3;
    }
 
    double Xf, Omega;
    Omega = darkOmega(&Xf, 1, 1e-4, NULL);
 
    FILE* out = fopen("results.json", "w");
    if (!out) { return 4; }
    fprintf(out,
            "{\"cdm\":\"%s\",\"cdm_mass_GeV\":%.6e,\"omega_h2\":%.6e,\"xf\":%.3f}\n",
            CDM[1], McdmN[1], Omega, Xf);
    fclose(out);
    return 0;
}

Use bare #include "micromegas.h" as shown above. micromegas-compile exports CPATH=<MICROMEGAS_ROOT>/include to the underlying make invocation, so the same main.c compiles unchanged on the cloud and against a user's local micrOmegas install — no ../include/... prefix needed. Either form works (the legacy #include "../include/micromegas.h" is still resolvable), but bare include is the preferred convention because it stays portable across install layouts.

Also treat parameter names as data, not guesses: assignValW("Name", value) must use the exact Name column from vars1.mdl, and the function returns nonzero on unknown names while leaving defaults in place. For example:

if (assignValW("Mdm1", 500.0)) {
    fprintf(stderr, "Unknown parameter name: Mdm1\n");
    return 5;
}
if (assignValW("laSH", 0.1)) {
    fprintf(stderr, "Unknown parameter name: laSH\n");
    return 6;
}

For anything beyond a protocol smoke test, prefer starting from the closest micrOmegas reference main.c for that model and then trimming it down. Some models need extra setup before heavy calls such as darkOmega or nucleonAmplitudes.

Anything not written to results.json is still preserved in the uploaded output directory, and stdout is tailed into MAGNUS_RESULT for debugging.

Output Paths

All paths are relative to the working directory.

Naming conventions:

• <project_label>: short model name, e.g. singlet_scalar, IDM
• <run_label>: benchmark-point tag, e.g. mS_700GeV, tanb_10, point1

Workflow

Step 1: Compile the project

magnus run micromegas-compile -- \
  --calchep path/to/MyModel_CH \
  --main path/to/main.c \
  --output dm/my_model \
  --project my_model

Parameters:

The CalcHEP directory and main.c are uploaded via FileSecret. On success, the compiled project — including the ./main binary and the auto-generated matrix-element code — is downloaded to --output.

WARNING: If --output points to an existing directory it will be deleted and replaced.

Result (magnus job result <job-id>):

• success (bool)
• project_dir (str): path to the compiled project
• main_binary (str): path to ./main inside the container (informational)
• make_log_tail (str, only on failure): last 3000 chars of make output

Step 2: Run the calculation

magnus run micromegas-calc -- \
  --project dm/my_model \
  --output dm/my_model/run1 \
  --slha path/to/params.slha

Parameters:

Result (magnus job result <job-id>):

• success (bool)
• output_dir (str): path to uploaded run output
• results (object, if results.json present): the parsed JSON written by main.c
• stdout_tail (str): last 4000 chars of ./main stdout
• stderr_tail (str, on failure only)

Parameter scans: compile once, run many

Because the compiled project is returned as a directory and consumed by --project, a typical scan looks like:

# Compile once
magnus run micromegas-compile -- \
  --calchep models/Singlet_CH --main src/main.c \
  --output dm/singlet --project singlet
 
# Scan points by feeding different SLHA cards (or --extra_args if main.c parses argv)
for m in 200 400 600 800 1000 1200; do
  render_slha $m > params_$m.slha
  magnus run micromegas-calc -- \
    --project dm/singlet \
    --slha params_$m.slha \
    --output dm/singlet/mS_${m}GeV
done

Every run reuses the single compiled binary — the expensive symbolic-amplitude step is paid only once.

Examples

Relic density only

main.c:

#include "micromegas.h"
#include "micromegas_aux.h"
#include <stdio.h>
 
int main(int argc, char** argv) {
    if (argc > 1 && slhaRead(argv[1], 0)) return 2;
    if (sortOddParticles(NULL)) return 1;
    if (!CDM[1]) return 3;
    double Xf, Omega = darkOmega(&Xf, 1, 1e-4, NULL);
    FILE* out = fopen("results.json", "w");
    if (!out) return 4;
    fprintf(out,
            "{\"cdm\":\"%s\",\"cdm_mass_GeV\":%.6e,\"omega_h2\":%.6e,\"xf\":%.3f}\n",
            CDM[1], McdmN[1], Omega, Xf);
    fclose(out);
    return 0;
}

Relic + direct detection (SI/SD on p/n)

#include "micromegas.h"
#include "micromegas_aux.h"
#include <stdio.h>
 
int main(int argc, char** argv) {
    if (argc > 1 && slhaRead(argv[1], 0)) return 2;
    if (sortOddParticles(NULL)) return 1;
    if (!CDM[1]) return 3;
 
    double Xf, Omega = darkOmega(&Xf, 1, 1e-4, NULL);
 
    double pA0[2], pA5[2], nA0[2], nA5[2];
    nucleonAmplitudes(CDM[1], pA0, pA5, nA0, nA5);
 
    double mN = 0.939, Mdm = McdmN[1];
    double mu = Mdm * mN / (Mdm + mN);
    double pref = 4.0 / M_PI * mu * mu * 2.568e9;   // GeV^-2 → pb
 
    FILE* out = fopen("results.json", "w");
    if (!out) return 4;
    fprintf(out,
        "{\"cdm\":\"%s\",\"mdm\":%.3f,\"omega_h2\":%.6e,"
        "\"sigma_SI_p\":%.6e,\"sigma_SI_n\":%.6e,"
        "\"sigma_SD_p\":%.6e,\"sigma_SD_n\":%.6e}\n",
        CDM[1], Mdm, Omega,
        pref * pA0[0]*pA0[0], pref * nA0[0]*nA0[0],
        3.0 * pref * pA5[0]*pA5[0], 3.0 * pref * nA5[0]*nA5[0]);
    fclose(out);
    return 0;
}

Key conventions