Dear All,

I would like to simulate the permeation of a solvent through a membrane. To do this, I first modeled the polymer on top of a graphene surface, which acts as a wall to prevent the solvent molecules from crossing to the other side. However, after relaxation, the graphene surface is no longer flat — it appears curved. When I tried to fix the graphene surface using the Fix Bonds option, I encountered an error:

---- Dynamics parameters ----

Ensemble                        : NPT
Temperature                     : 300.00 K
 Control method                : Nose
 Q ratio                       : 0.0100000
Pressure                        : 0 GPa
 Control method                : Berendsen
 Decay constant                : 0.1000000 ps
Timestep                        : 5.00 fs
Number of steps                 : 2000
Duration                        : 10 ps
Initial velocities              : Random
Random number seed              : 1761653801

---- Energy parameters ----

Forcefield                      : COMPASSIII (Version 1.2)
Charges                         : Forcefield assigned
Electrostatic terms:
 Summation method              : Group based
 Cutoff distance               : 12.5 A
 Spline width                  : 1 A
 Buffer width                  : 2 A

van der Waals terms:
 Summation method              : Group based
 Cutoff distance               : 12.5 A
 Spline width                  : 1 A
 Long range correction         : YES
 Buffer width                  : 2 A

---- Thermostat ----

Nose mass parameter, Q mass: 52.9808 kcal/mol*ps^2

---- Constraints ----

Distance constraints            : 6140

Dynamics calculation - Constraint solver failed to converge

Termination status : Error