[QE-users] 2D correction for using pw.x with Environ

Rabet, Sahar s.rabet at tu-braunschweig.de
Tue Jun 20 17:21:56 CEST 2023


Dear QE community,

I have a question regarding using of continuum solvation model using ENVIRON with pw.x.
I am just running a relaxation simulation for an approximately large hydrocarbon (15 atoms) on a metallic surface. I did this before in a vacuum. I am now wondering by using the ENVIRON, does the dipole correction still work (tefield     = .true.,  dipfield    = .true.)? or just the implemented PBC correction in ENVIRON (PCC) is enough (pbc_correction = 'parabolic') and no need for any dipole correction? shall I consider any other correction for 2D system by using ENVIRON?

Many thanks,

Sahar Rabet

Technische Universität Braunschweig
Institut für Thermodynamik

Deutschland / Germany

p.s. following is the input file that I am using.

&ENVIRON
  !
  verbose = 0
  env_surface_tension = 11.5
  env_electrostatic = .true.
  env_pressure = 0.0
  environ_thr = 0.1
  env_static_permittivity = 5.7
  env_external_charges = 2
  !
/
&BOUNDARY
  !
  rhomax = 0.005
  rhomin = 0.0001
  solvent_mode = 'full'
  !
/
&ELECTROSTATIC
  !
  pbc_dim = 2
  solver = 'cg'
  pbc_correction = 'parabolic'
  mix = 0.5
  pbc_axis = 3
  tol = 1D-11
  auxiliary = 'none'
  !
/
EXTERNAL_CHARGES (bohr)
  2 0 0 17.1178   0.5 2 3
  2 0 0 -8.3808   0.5 2 3
-----------------------------------
&CONTROL
 calculation = 'relax'
 prefix = 'tmp',
 pseudo_dir ='',
 outdir = "./"
 restart_mode = "from_scratch"
 tprnfor=.true
 disk_io='low',
 tefield     = .true.,
 dipfield    = .true.
  /
 &SYSTEM
 ibrav = 0,
 nat=90
 ntyp=4
 vdw_corr = 'DFT-D3'
 occupations = "smearing",
 smearing    = "gaussian",
 degauss     = 0.01,
 emaxpos = 0.80
 eopreg = 0.01
 edir = 3
 ecutwfc = 50.0
 ecutrho = 400.0
 tot_charge = -4.0
  /
 &ELECTRONS
 mixing_beta=0.1,
 conv_thr=1e-06,
 mixing_mode = 'local-TF'
 electron_maxstep = 500
  /
 &IONS
  ion_dynamics="bfgs",
 /
ATOMIC_SPECIES
Cu 63.546 Cu.pbe-dn-kjpaw_psl.1.0.0.UPF
H  1.008 H.pbe-kjpaw_psl.1.0.0.UPF
O  15.9994 O.pbe-n-kjpaw_psl.1.0.0.UPF
C  12.0107 C.pbe-n-kjpaw_psl.1.0.0.UPF

CELL_PARAMETERS {angstrom}
12.69256670000000 0.00000000000000 0.00000000000000
0.00000000000000 12.69256670000000 0.00000000000000
0.00000000000000 0.00000000000000 33.59000000000000

K_POINTS {automatic}
 3 3 1 0 0 0

ATOMIC_POSITIONS {angstrom}



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20230620/c2d83a84/attachment.html>


More information about the users mailing list