[Mon Jul 03 17:53:27 CST 2023] [MD] [warn] 'Starting MedeA Core 3.6.0' Forcefield set to: pcff+ File: Local::/home/medea/MD/Jobs/dir0/222/pcff+.frc md5 sum: 271FB553B535FE0449E31A6F8A0898D3 MedeA version 3.6.0 #------------------------------------------------------------------------------- # This calculation has 2 stages #------------------------------------------------------------------------------- Stage 1: Variables tstep = 1 fs T = 298.2 K P = 1 atm Stage 2: This LAMMPS calculation has 6 stages: Stage 2.1: Initialize LAMMPS for a Class 2 type of forcefield 'pcff+' (from Local::/home/medea/MD/Jobs/dir0/222/pcff+.frc) Current system formula: K94Na8Li30Mg616S93Cl1178O8675H16606 (K94Na8Li30Mg616S93Cl1178O8675H16606) The nonbond terms are handled with a simple cutoff at 9.5 Angstrom. The long-range van der Waals interactions are included via tail corrections. Stage 2.2: Minimize using conjugate gradients and linesearch fast Relative energy convergence: 0.0 Force convergence: 1.0 Maximum atomic motion: 0.05 Maximum number of iterations: 1000 Maximum number of force evaluations: 10000 Stage 2.3: Set the velocities for 298.2 K Stage 2.4: NPT integration for 2 ns with a timestep of 1 fs P is 1 atm with cell constrained isotropically T is 298.2 K Stage 2.5: NVT integration for 5 ns with a timestep of 1 fs T is 298.2 K Stage 2.6: Surface tension calculation using NVT integration for 2 ns with a timestep of 1 fs, T is 298.2 K #------------------------------------------------------------------------------- # Running the calculation #------------------------------------------------------------------------------- Stage 1: Variables tstep = 1 fs T = 298.2 K P = 1 atm Stage 2: This LAMMPS calculation has 6 stages: ***** LAMMPS error in LAMMPS.out ***** ERROR: Non-numeric pressure - simulation unstable (../fix_nh.cpp:1061) ************************************** Stage 2.1: Initialize LAMMPS for a Class 2 type of forcefield 'pcff+' (from Local::/home/medea/MD/Jobs/dir0/222/pcff+.frc) Current system formula: K94Na8Li30Mg616S93Cl1178O8675H16606 (K94Na8Li30Mg616S93Cl1178O8675H16606) The nonbond terms are handled with a simple cutoff at 9.5 Angstrom. The long-range van der Waals interactions are included via tail corrections. Stage 2.2: Minimize using conjugate gradients and linesearch fast Relative energy convergence: 0.0 Force convergence: 1.0 Maximum atomic motion: 0.05 Maximum number of iterations: 1000 Maximum number of force evaluations: 10000 Note: 2.2_Minimization.out contains warning messages from LAMMPS - examine this file for additional information. The minimization finished after 41 steps ending because the linesearch alpha is zero Property Value Units ------------------ ---------- ----------- nSteps: 41 Initial Fmax: 33854.8 kJ/mol/Ang Fmax: 337478.1 kJ/mol/Ang Initial Frms: 1644.3 kJ/mol/Ang Frms: 3387.7 kJ/mol/Ang P: 15335.3 atm V: 278556.9 Ang^3 rho: 1.3074 g/mL Sxx: -16354.6 atm Syy: -15467.4 atm Szz: -14183.8 atm Syz: 1953.3 atm Sxz: -1142.3 atm Sxy: -438.8 atm Initial Epot: 1571063.3 kJ/mol Epot: -1683939.9 kJ/mol a: 79.1497 Ang b: 76.7349 Ang c: 45.8640 Ang alpha: 90.00 degree beta: 90.00 degree gamma: 90.00 degree sum_fx: -0.0 kJ/mol/Ang sum_fy: 0.0 kJ/mol/Ang sum_fz: -0.0 kJ/mol/Ang Stage 2.3: Set the velocities for 298.2 K Stage 2.4: NPT integration for 2 ns with a timestep of 1 fs P is 1 atm with cell constrained isotropically T is 298.2 K Property Value +/- Uncertainty Units After Steps % Run ------------------ ---------- ----------- ------------- ------------ ------- t: fs Error with T: The run is not long enough to analyze. Sorry! Error with P: The run is not long enough to analyze. Sorry! Error with V: The run is not long enough to analyze. Sorry! Error with rho: The run is not long enough to analyze. Sorry! Error with Etotal: The run is not long enough to analyze. Sorry! Error with a: The run is not long enough to analyze. Sorry! Error with b: The run is not long enough to analyze. Sorry! Error with c: The run is not long enough to analyze. Sorry! Error with alpha: The run is not long enough to analyze. Sorry! Error with beta: The run is not long enough to analyze. Sorry! Error with gamma: The run is not long enough to analyze. Sorry! Error with Epot: The run is not long enough to analyze. Sorry! Error with Ekin: The run is not long enough to analyze. Sorry! Error with Evdw: The run is not long enough to analyze. Sorry! Error with Ecoul: The run is not long enough to analyze. Sorry! Error with Sxx: The run is not long enough to analyze. Sorry! Error with Syy: The run is not long enough to analyze. Sorry! Error with Szz: The run is not long enough to analyze. Sorry! Error with Syz: The run is not long enough to analyze. Sorry! Error with Sxz: The run is not long enough to analyze. Sorry! Error with Sxy: The run is not long enough to analyze. Sorry! Stage 2.5: NVT integration for 5 ns with a timestep of 1 fs T is 298.2 K Stage 2.6: Surface tension calculation using NVT integration for 2 ns with a timestep of 1 fs, T is 298.2 K ********************************************* WARNING: this LAMMPS stage completed on Mon 03 July 2023 at 17:53:55 CST after 14 s (0:00:14) with the following error: ERROR: Non-numeric pressure - simulation unstable (../fix_nh.cpp:1061) ********************************************* Entire job completed on Mon 03 July 2023 at 17:53:55 CST after 14 s (0:00:14) and running 1 tasks.