"Conserved En." is increasing linearly with simulation time - Is this normal?

GROMACS version:2020.3
GROMACS modification: No

My aim is to compute the self-diffusivity of methanol in zeolite Y.
To achieve this, the first task is to equilibrate the system. To assess whether the system is equilibrated or not, I plot various quantities as a function of simulation time. These quantities are:
Potential energy, Kinetic energy, Total energy, Conserved energy, Temperature, and Pressure.
I get these quantities using the `gmx energy` command.

In spite of its name, the “conserved energy” is increasing linearly with simulation time (see attached plot). Is this a normal behavior?

I must admit I am not familiar with the quantity “conserved energy”. What is it describing?

Dear Dr. Lemkul,

Unfortunately, I am not able to understand the information in the link. Is there not any review article or text book that explains what this quantity is and/or its expected behavior?

Thanks

Maybe I should add that I am doing an NVT simulation, and I am using the Nose-Hoover thermostat. I also attach the mdp file in case it contains relevant information.

Regards,

This answer on the mailing list from 2013 is incorrect.

The conserved energy quantity should be conserved analytically. However, discretization as well as approximations lead to small integration errors causing a linear drift in most cases. How large this drift can be be and what one can tolerate depends on the system and is often difficult to say.

I don’t know the size of your system, so I can’t say what reasonable for your case.

My system contains 4800 atoms.
The length of the simulation is 40 ns.
The total drift of the conserved energy is 1.11767e+06 kJ/mol

Normalizing, the drift of the conserved energy per ns and atoms becomes:
5.82 kJ mol-1 ns-1 atom-1

Does this look like a reasonable drift of the conserved energy?

A late answer: yes, that is a reasonable order of magnitude.