def get_task(task_name):
cell = vasp.parse_poscar(Si_str.split('\n'))
# Vasp settings
incar = vasp.incar()
incar.set_structure_optimization()
incar.set_nsw(20)
incar.set_encut(300)
incar_ph_rx = vasp.incar()
incar_ph_rx.set_structure_optimization()
incar_ph_rx.set_nsw(20)
incar_ph_rx.set_encut(300)
incar_ph_rx.set_isif(4) # volume constant
incar_ph_dsp = vasp.incar()
incar_ph_dsp.set_electronic_structure()
incar_ph_dsp.set_encut(300)
# Queue
job = ge.job(script="mpirun vasp5212mpi",
shell="/bin/zsh",
jobname=task_name,
pe="mpi* 8",
stdout="std.log",
stderr="err.log")
# Task
task = vasp.quasiharmonic_phonon(max_iteration=10,
min_iteration=1,
sampling_mesh=[20, 20, 20],
supercell_matrix=np.diag([2, 2, 2]),
cell=cell,
pseudo_potential_map=ps_map,
incar=[incar,
incar_ph_rx,
incar_ph_dsp],
k_mesh=[[8, 8, 8],
[8, 8, 8],
[4, 4, 4]],
job=[job,
job.copy(task_name + "-ph_relax"),
job.copy(task_name + "-ph")],
traverse=False)
return task
def get_task(task_name):
cell = vasp.parse_poscar(Si_str.split('\n'))
# Vasp settings
incar = vasp.incar()
incar.set_structure_optimization()
incar.set_nsw(20)
incar.set_encut(300)
incar_ph_rx = vasp.incar()
incar_ph_rx.set_structure_optimization()
incar_ph_rx.set_nsw(20)
incar_ph_rx.set_encut(300)
incar_ph_rx.set_isif(4) # volume constant
incar_ph_dsp = vasp.incar()
incar_ph_dsp.set_electronic_structure()
incar_ph_dsp.set_encut(300)
# Queue
job = ge.job(script="mpirun vasp5212mpi",
shell="/bin/zsh",
jobname=task_name,
pe="mpi* 8",
stdout="std.log",
stderr="err.log")
# Task
task = vasp.quasiharmonic_phonon(max_iteration=10,
min_iteration=1,
sampling_mesh=[20, 20, 20],
supercell_matrix=np.diag([2, 2, 2]),
cell=cell,
pseudo_potential_map=ps_map,
incar=[incar, incar_ph_rx, incar_ph_dsp],
k_mesh=[[8, 8, 8], [8, 8, 8], [4, 4, 4]],
job=[
job,
job.copy(task_name + "-ph_relax"),
job.copy(task_name + "-ph")
],
traverse=False)
return task
incar.set_structure_optimization()
incar.set_nsw(20)
incar_ph_rx = vasp.incar()
incar_ph_rx.set_structure_optimization()
incar_ph_rx.set_nsw(20)
incar_ph_rx.set_isif(4) # volume constant
# incar_ph_rx.set_isif(2) # lattice vectors constant
incar_ph_dsp = vasp.incar()
incar_ph_dsp.set_electronic_structure()
# Queue
job = ge.job(script="mpirun vasp5212mpi",
shell="/bin/zsh",
jobname=task_name,
pe="mpi* 4",
stdout="std.log",
stderr="err.log")
# Task
task = vasp.mode_gruneisen(max_iteration=10,
min_iteration=1,
delta_strain=0.002,
supercell_matrix=np.diag([2, 2, 2]),
cell=cell,
pseudo_potential_map=ps_map,
incar=[incar,
incar_ph_rx,
incar_ph_dsp],
k_mesh=[[4, 4, 6],
[4, 4, 6],
cell = cogue.cell(lattice=lattice,
points=points,
symbols=symbols)
# Vasp settings
ps_map = {'Sn': 'Sn_PBE',
'O': 'O_PBE'}
incar = vasp.incar()
incar.set_structure_optimization()
incar.set_encut(400)
incar.set_prec("Normal")
# Queue
job = ge.job(script="vasp5212serial",
shell="/bin/zsh",
jobname=task_name,
stdout="std.log",
stderr="err.log")
# Task
task = vasp.bulk_modulus(max_iteration=2,
cell=cell,
pseudo_potential_map=ps_map,
k_mesh=[4, 4, 6],
incar=incar,
job=job)
# Automatic calculation
calc = cogue.autocalc()
calc.append(task_name, task) # More tasks can be appended.
calc.set_queue(ge.queue())
cell = cogue.cell(lattice=lattice, points=points, symbols=symbols)
# Vasp settings
ps_map = {'Si': 'Si_PBE', 'O': 'O_PBE'}
incar = vasp.incar()
incar.set_structure_optimization()
incar.set_nsw(20)
incar.set_pstress(800)
incar_phonon = vasp.incar()
incar_phonon.set_electronic_structure()
# Queue
job = ge.job(script="mpirun vasp5212mpi",
shell="/bin/zsh",
jobname=task_name,
pe="mpi* 4",
stdout="std.log",
stderr="err.log")
# Task
task = vasp.phonon(max_iteration=10,
min_iteration=1,
supercell_matrix=np.diag([2, 2, 2]),
cell=cell,
pressure_target=80,
pseudo_potential_map=ps_map,
incar=[incar, incar_phonon],
k_mesh=[[4, 4, 6], [2, 2, 3]],
k_shift=[[0, 0, 0], [0, 0, 0.5]],
job=job)
