LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] corner atoms fly away for core-shell model
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] corner atoms fly away for core-shell model

From: Axel Kohlmeyer <akohlmey@...24...>
Date: Fri, 18 Aug 2017 11:53:27 -0400

On Fri, Aug 18, 2017 at 9:31 AM, bo He <bohe@...6615...> wrote:
Dear Lammps users,

I am using core-shell model for Barium titanat oxide, i would like to build the model same as figure 6 in It's a cantilever beam, one side is fixed, the another side is free and after relaxation a downward force will apply to the free end. However during relaxation the corner atoms at the free end fly away and blowing up the temperature.

​are you certain this is an effect​ of the core-shell model?
have you tried the same geometry with a non-polarizable setup?


Details of my model:

1.Core-shell model for
Barium titanat oxide, the core-shell requires periodic boundary for all direction. In order to study the size effect, sufficient vacuum space created around the model.

2. temperature is applied through nvt. I also tried with temp/rescale, same thing happens.

3. The system is not charge neutral due to free surfaces.

4. Potential used from DOI: 10.1063/1.4827475

5. code of relaxation:

# compute temperature only for the mobile part (exclude the fixed end)
compute CSequ mobile temp/cs mobile_cores mobile_shells

# velocity bias option

variable T_start equal 415

velocity mobile create ${T_start} 134 dist gaussian mom yes rot yes bias yes temp CSequ

velocity mobile scale ${T_start} temp CSequ

velocity fixed set 0 0 0

fix setforce_fixed fixed setforce 0.0 0.0 0.0

fix nve all nve

fix temp_mobile mobile temp/rescale 100 ${T_start} ${T_start} 1 0.5

fix mom_mobile mobile momentum 10 linear 0 0 0 angular

fix_modify temp_mobile temp CSequ

thermo_modify temp CSequ  

timestep 0.004

run 1000

Any suggestion is appreciated.

Best regards


Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list

Dr. Axel Kohlmeyer  akohlmey@...12...24...
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.