<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><span style="font-size:12.8px">Dear DFTB+ community,</span><br style="font-size:12.8px"><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I am a user of DFTB+1.2, and currently I have encountered some problems with the barostat in DFTB+ when I try to simulate a gas phase system.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I am trying to simulate the association and crystallization behavior of atomic borons and nitrogens at a given temperature. So I constructed a cubic supercell system with 343 atoms, with about equal numbers of B and N atoms. Initial distances between atoms are 4nm, therefore they don't interact. The system is subject to thermostat and barostat. Below is my input file for the test.</div><div style="font-size:12.8px">*******************</div><div><div style="font-size:12.8px">Geometry = GenFormat {</div><div style="font-size:12.8px"> <<< "bnnt.gen"</div><div style="font-size:12.8px"> }</div><div style="font-size:12.8px"> Driver = VelocityVerlet{</div><div style="font-size:12.8px">  MovedAtoms = 1:-1</div><div style="font-size:12.8px">  TimeStep [Femtosecond] = 1.0</div><div style="font-size:12.8px">  Thermostat = NoseHoover {</div><div style="font-size:12.8px">    Temperature [K] = TemperatureProfile {</div><div style="font-size:12.8px">        constant      1 8000</div><div style="font-size:12.8px">        linear    40000 2000</div><div style="font-size:12.8px">        constant  10000000 2000</div><div style="font-size:12.8px">        }</div><div style="font-size:12.8px">    CouplingStrength [cm^-1] = 3000</div><div style="font-size:12.8px">  }</div><div style="font-size:12.8px">  Barostat = {</div><div style="font-size:12.8px">    Pressure [Pa]=500000.0</div><div style="font-size:12.8px">    Timescale [ps]=0.01</div><div style="font-size:12.8px">    }</div><div style="font-size:12.8px">  MDRestartFrequency = 500</div><div style="font-size:12.8px"> OutputPrefix = "geo_stp"</div><div style="font-size:12.8px"> }</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"> Hamiltonian = DFTB{</div><div style="font-size:12.8px">    Charge = 0</div><div style="font-size:12.8px">    Filling = Fermi{</div><div style="font-size:12.8px">        Temperature [K] = 3000</div><div style="font-size:12.8px">    }</div><div style="font-size:12.8px">    KPointsAndWeights = SupercellFolding{</div><div style="font-size:12.8px">        1       0       0</div><div style="font-size:12.8px">        0       1       0</div><div style="font-size:12.8px">        0       0       1</div><div style="font-size:12.8px">        0.0     0.0     0.0</div><div style="font-size:12.8px">    }</div><div style="font-size:12.8px">    MaxAngularMomentum = {</div><div style="font-size:12.8px">       B = "p"</div><div style="font-size:12.8px">       N = "p"</div><div style="font-size:12.8px">  }</div><div style="font-size:12.8px">  SlaterKosterFiles = Type2FileNames {</div><div style="font-size:12.8px">    Prefix = "/gpfs/home/lonhan/software/ex<wbr>ecutable/SK-files/matsci-0-3/"</div><div style="font-size:12.8px">    Separator = "-"</div><div style="font-size:12.8px">    Suffix = ".skf"</div><div style="font-size:12.8px">    LowerCaseTypeName = No</div><div style="font-size:12.8px">  }</div><div style="font-size:12.8px"> }</div><div><div><span style="font-size:12.8px">Options={</span></div><div><span style="font-size:12.8px">  ShowFoldedCoords=Yes</span></div><div><span style="font-size:12.8px">}</span></div></div></div><div style="font-size:12.8px">****************</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">The pressure reported in the MD output is orders of magnitude higher than it should be. And after some thousands of steps, the run will fail with the following error:</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><div><i>ERROR!</i></div><div><i>-> Failure in diagonalisation routine dsygvd, non-positive definite overlap! Minor 485107 responsible.</i></div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I tried to increase the coupling strength up to 10000 (cm^-1), but I still get this error. </div><div style="font-size:12.8px">The input geometry and output files are attached.<br></div><div style="font-size:12.8px">Could anyone​​​​​​ help me identify if this is a problem with my simulation approach or with the barostat?</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Best regards,</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Longtao Han</div><div style="font-size:12.8px">PhD student</div><div style="font-size:12.8px">IACS, Stony Brook University</div><div style="font-size:12.8px">Email: <a href="mailto:longtao.han@stonybrook.edu" target="_blank">longtao.han@stonybrook.<wbr>edu</a></div></div>
</div><br></div>