Homec4science
Diffusion lammps (lammm-devel) Browse

lammps/examples/USER/atc/elastic/in.electron_density
1ca6fb0f8160lammm-devel

Last Change
Raw File

in.electron_density
View Options

# needs description
echo both
units metal
# PARAMETERS-----------------------------
variable s equal 1
variable L equal 10
variable e equal 4
variable E equal 0.0001
variable V equal $E*$L
# END -----------------------------------
atom_style atomic
lattice diamond 1.0
boundary f p p
region box block -$L $L 0 1 0 1
create_box 1 box
group box region box
atom_modify sort 0 1
timestep 0.0
mass * 12.01
# coupling
### NOTE ATC: material cnt does not provide all interfaces for charge_density physics and will be treated as null
fix AtC box atc electrostatic-equilibrium CNT.mat
fix_modify AtC internal_quadrature off
#fix_modify AtC atom_weight constant internal 1.0 NOTE penultimate is a group
fix_modify AtC atom_weight constant box 1.0
fix_modify AtC omit atomic_charge
fix_modify AtC mesh create $e 1 1 box f p p
#fix_modify AtC control momentum flux
fix_modify AtC mesh create_elementset all -INF INF -INF INF -INF INF
# bcs/ics conditions
fix_modify AtC fix displacement x all 0.0
fix_modify AtC fix displacement y all 0.0
fix_modify AtC fix displacement z all 0.0
fix_modify AtC fix velocity x all 0.0
fix_modify AtC fix velocity y all 0.0
fix_modify AtC fix velocity z all 0.0
fix_modify AtC mesh create_nodeset lbc -10 -10 INF INF INF INF
fix_modify AtC mesh create_nodeset rbc 10 10 INF INF INF INF
fix_modify AtC fix electric_potential lbc linear 0 0 0 $E 0 0 $V
fix_modify AtC fix electric_potential rbc linear 0 0 0 $E 0 0 $V
# run
thermo_style custom step cpu etotal ke
thermo $s
fix_modify AtC output electron_densityFE $s text
fix_modify AtC output index step
log electron_density.log
# run default material
print "default material - table linear"
run $s
# run CNT1 material
print "CNT1 material - analytical linear"
fix_modify AtC material all CNT1
run $s
# run CNT2 material
print "CNT2 material - analytical exponetial"
fix_modify AtC material all CNT2
run $s
# run CNT material
print "CNT material - table DOS"
fix_modify AtC material all CNT
#variable E equal 10*$E
#fix_modify AtC fix electric_potential lbc linear 0 0 0 $E 0 0 $V
#fix_modify AtC fix electric_potential rbc linear 0 0 0 $E 0 0 $V
run $s
c4science · Help