Pt表面の液滴接触角(先行事例)
分子動力学法による白金表面上の水液滴構造分子動力学法による白金表面上の水液滴構造(東京大学)
https://www.photon.t.u-tokyo.ac.jp/~maruyama/papers/02/DS39_Kimura.pdf
lammpsで界面作成
液体:水(tip3p力場)
固体:Pt(INTERFACE力場)
固体の作成
-
INTERFACE FF からPtの構造データ(.car, .mdf)をダウンロード
-
lammpsのメインフォルダ内にあるtoolsにmsi2lmpというツールがあり、これによって.carファイル、.mdfファイルからlammpsファイルを作成できる。
液体をPt上に乗せる
units real
atom_style full
bond_style harmonic
angle_style harmonic
dihedral_style zero
improper_style zero
pair_style lj/cut/coul/long 12 12
pair_modify mix arithmetic
special_bonds lj/coul 0.0 0.0 0.5
kspace_style pppm 1e-4
neigh_modify every 1 delay 0 check yes
variable N_A index 6.022140857e+23
variable temp index 300
variable nstep index 100000
#____________SiO2読み込み___________#
read_data pt111_surface_3.data &
extra/atom/types 3 &
extra/bond/types 2 &
extra/angle/types 1 &
extra/dihedral/types 0 &
extra/improper/types 0 &
extra/bond/per/atom 2 &
extra/angle/per/atom 1 &
extra/dihedral/per/atom 0 &
extra/improper/per/atom 0 &
extra/special/per/atom 100
variable lx index $(lx)
variable ly index $(ly)
variable nwall index $(count(all))
group bottom_wall id 1:${nwall}
#________Pt 固定相と熱浴の設定____________#
group Pt type 1
displace_atoms all move 0 0 3
region fixed_Pt plane 0 0 $(bound(Pt,zmin)+0.5) 0 0 -1
group fixed_Pt region fixed_Pt
region fixed_Pt delete
print $(count(fixed_Pt))
region hbath_Pt plane 0 0 $(bound(Pt,zmin)+5) 0 0 -1
group deleteme region hbath_Pt
group hbath_Pt subtract deleteme fixed_Pt
region hbath_Pt delete
group deleteme delete
print $(count(hbath_Pt))
#__________液体分子を格子状に配置___________#
省略
#__________界面の力場パラメータを設定___________#
#Pair Coeffs # lj/cut/coul/long/gpu
#1 7.8005 2.53457
#2 0.152073 3.15066
#3 0 0
#4 0 0
variable para index 0.5
variable surface_eps index $(v_para*sqrt(7.8005*0.152073))
variable surface_sig index $((2.53457+3.15066)/2)
pair_coeff 1 2 ${surface_eps} ${surface_sig}
write_data data.Pt_liquidCA
#_________________________________________________________#
#___________________________実行_____________________________#
minimize 1.0e-4 1.0e-6 100 1000
fix hbath_Pt hbath_Pt langevin ${temp} ${temp} 100 1 tally yes
fix vliquid liquid temp/csvr ${temp} ${temp} 100 12345
group integ subtract all fixed_Pt
fix integ integ nve
compute comliquid liquid temp
compute comSiO2 integ temp
dump methDump all atom 1000 dump.lammpstrj
reset_timestep 0
compute Dens_liquid2D liquid chunk/atom bin/2d z lower 0.5 y lower 0.5
compute Dens_liquid liquid chunk/atom bin/1d z lower 0.005 units reduced
fix 2 liquid ave/chunk 1 ${nstep} ${nstep} Dens_liquid density/mass ave running file liquid_density.dat
fix 3 liquid ave/chunk 1 ${nstep} ${nstep} Dens_liquid2D density/mass ave running file liquid_density2D.dat
timestep 1
#variable GamNsurf equal lz*(pzz-0.5*(pxx+pyy))*101325/1e10*1e3
#fix ave/G all ave/time 100 5 1000 v_GamNsurf file surface_tension.dat
thermo_style custom step time temp press vol density c_comliquid c_comSiO2
thermo 1000
run ${nstep}
write_data data.last
液滴の密度可視化
きれいな液滴形になっていないが、これはPt表面が平らなせいか液滴が左右に動いてしまったため。界面に密度が高くなっている場所が確認できる。これは界面からの相互作用によって水の構造化が起きているのだと考えられる。