# Temperature

Hi all
I have a question about how imposing the temperature of the environment in my simulation.
If I want to simulate a functional CNT at 600K (the temperature of the environment) and a random velocity is adopted to reach the room temperature and the simulation commences by relaxing the CNT at room temperature for 5 ps,
Is it true to consider 600K as room temperature?
Indeed, I use this approach but interestingly, I see the strange behavior in my simulation. Some functional group is dissociated from CNT’s wall! (I use AIREBO potential and CH2=C as functional group)

Thanks.

Hi all
I have a question about how imposing the temperature of the environment in
my simulation.
If I want to simulate a functional CNT at 600K (the temperature of the
environment) and a random velocity is adopted to reach the room temperature
and the simulation commences by relaxing the CNT at room temperature for 5
ps,

​not sure what you are saying here. it would be better, if you provide your
input and explain what​ you are referring to with which statement.

Is it true to consider 600K as room temperature?

Indeed, I use this approach but interestingly, I see the strange behavior
in my simulation. Some functional group is dissociated from CNT's wall! (I
use AIREBO potential and CH2=C as functional group)

​doesn't seem very interesting to me. you could just have an unrelaxed or
unstable initial geometry. also, at elevated temperatures, bonds will break
more easily., especially, when they are already strained in some way.

axel.

This is my input file:

# ---------- Initialize Simulation ---------------------

units metal
dimension 3
atom_style atomic
boundary s s s
atom_modify map array
timestep 1e-3

# ---------- Define Interatomic Potential --------------

pair_style airebo 3.0 1 1
pair_coeff * * CH.airebo C H

# ---------- Define Settings ----------------------------

#Initial velocities------------------------------------------------

velocity all create 300 511124 rot yes mom yes

# Relaxing the Imported structure in the desired temperature------------------

fix relax all nvt temp 300 300 0.01
compute mystress middle_atoms stress/atom myTemp
compute reducedstress middle_atoms reduce sum c_mystress[3]

thermo_style custom c_reducedstress
thermo 1000
dump relax_position all atom 5000 dump2.atom #For having video!
run 5000
undump relax_position
unfix relax

# Zeroing Stress----------------------------------------------------

thermo_style custom fmax fnorm
min_style cg
minimize 0 0 1000 1000000

# Defining fixes-------------------------------------------------------

variable stepnum equal 1000
variable Nstep loop \${stepnum}
variable disp_upper_atoms equal 0.01*v_Nstep

fix fixing_Left_side lower_atoms setforce 0 0 0
fix fixing_Right_side upper_atoms setforce 0 0 0

label looplable

displace_atoms upper_atoms move 0 0 -0.01
fix fix1 middle_atoms nvt temp 300 300 0.01
thermo 5000
thermo_style custom step temp pe v_sigma
run 5000

next Nstep
jump in.CNT looplable

Is it true to change my input file to:

#Initial velocities------------------------------------------------

velocity all create 300 511124 rot yes mom yes =====> velocity all create 600 511124 rot yes mom yes

# Relaxing the Imported structure in the desired temperature------------------

fix relax all nvt temp 300 300 0.01 ======> fix relax all nvt temp 600 600 0.01

# Defining fixes-------------------------------------------------------

label looplable

displace_atoms upper_atoms move 0 0 -0.01
fix fix1 middle_atoms nvt temp 300 300 0.01 ======> fix fix1 middle_atoms nvt temp 600 600 0.01
thermo 5000
thermo_style custom step temp pe v_sigma
run 5000

next Nstep
jump in.CNT looplable

in order to simulate my functional CNT in the thermal environment? (600K is the temperature of the environment. Indeed, I want to simulate the functional CNT in thermal environment)
I hope I can explain my purpose truly.
Thanks.

Is it true to change my input file to:

#Initial velocities------------------------------------------------

velocity all create 300 511124 rot yes mom yes =====>
velocity all create 600 511124 rot yes mom yes

# Relaxing the Imported structure in the desired
temperature------------------

fix relax all nvt temp 300 300 0.01 ======> fix relax
all nvt temp 600 600 0.01

...

# Defining fixes-------------------------------------------------------
....

label looplable

displace_atoms upper_atoms move 0 0 -0.01
fix fix1 middle_atoms nvt temp 300 300 0.01 ======>
fix fix1 middle_atoms nvt temp 600 600 0.01
thermo 5000
thermo_style custom step temp pe v_sigma
run 5000

next Nstep
jump in.CNT looplable

in order to simulate my functional CNT in the thermal environment? (600K
is the temperature of the environment. Indeed, I want to simulate the
functional CNT in thermal environment)
I hope I can explain my purpose truly.

​whether this does what you want or not, is not something that you should