diff --git a/PW/Doc/INPUT_PW.def b/PW/Doc/INPUT_PW.def index 2e29a1cbc..08266383c 100644 --- a/PW/Doc/INPUT_PW.def +++ b/PW/Doc/INPUT_PW.def @@ -2090,11 +2090,21 @@ input_description -distribution {Quantum Espresso} -package PWscf -program pw.x var tqr -type LOGICAL { default { .FALSE. } info { - If .true., use the real-space algorithm for augmentation - charges in ultrasoft pseudopotentials. - Must faster execution of ultrasoft-related calculations, - but numerically less accurate than the default algorithm. - Use with care and after testing! + If .true., use a real-space algorithm for augmentation + charges of ultrasoft pseudopotentials and PAWsets. + Faster but numerically less accurate than the default + G-space algorithm. Use with care and after testing! + } + } + + var real_space -type LOGICAL { + default { .FALSE. } + info { + If .true., exploit real-space localization to compute + matrix elements for nonlocal projectors. Faster and in + principle better scaling than the default G-space algorithm, + but numerically less accurate, may lead to some loss of + translational invariance. Use with care and after testing! } } } diff --git a/PW/Doc/INPUT_PW.html b/PW/Doc/INPUT_PW.html index c599a9417..7904177e4 100644 --- a/PW/Doc/INPUT_PW.html +++ b/PW/Doc/INPUT_PW.html @@ -55,59 +55,59 @@
-

INTRODUCTION

-

&CONTROL

+

INTRODUCTION

+

&CONTROL

-calculation | title | verbosity | restart_mode | wf_collect | nstep | iprint | tstress | tprnfor | dt | outdir | wfcdir | prefix | lkpoint_dir | max_seconds | etot_conv_thr | forc_conv_thr | disk_io | pseudo_dir | tefield | dipfield | lelfield | nberrycyc | lorbm | lberry | gdir | nppstr | lfcpopt | gate +calculation | title | verbosity | restart_mode | wf_collect | nstep | iprint | tstress | tprnfor | dt | outdir | wfcdir | prefix | lkpoint_dir | max_seconds | etot_conv_thr | forc_conv_thr | disk_io | pseudo_dir | tefield | dipfield | lelfield | nberrycyc | lorbm | lberry | gdir | nppstr | lfcpopt | gate
-

&SYSTEM

+

&SYSTEM

-ibrav | celldm | A | B | C | cosAB | cosAC | cosBC | nat | ntyp | nbnd | tot_charge | starting_charge | tot_magnetization | starting_magnetization | ecutwfc | ecutrho | ecutfock | nr1 | nr2 | nr3 | nr1s | nr2s | nr3s | nosym | nosym_evc | noinv | no_t_rev | force_symmorphic | use_all_frac | occupations | one_atom_occupations | starting_spin_angle | degauss | smearing | nspin | noncolin | ecfixed | qcutz | q2sigma | input_dft | exx_fraction | screening_parameter | exxdiv_treatment | x_gamma_extrapolation | ecutvcut | nqx1 | nqx2 | nqx3 | lda_plus_u | lda_plus_u_kind | Hubbard_U | Hubbard_J0 | Hubbard_alpha | Hubbard_beta | Hubbard_J | starting_ns_eigenvalue | U_projection_type | edir | emaxpos | eopreg | eamp | angle1 | angle2 | lforcet | constrained_magnetization | fixed_magnetization | lambda | report | lspinorb | assume_isolated | esm_bc | esm_w | esm_efield | esm_nfit | fcp_mu | vdw_corr | london | london_s6 | london_c6 | london_rvdw | london_rcut | dftd3_version | dftd3_threebody | ts_vdw_econv_thr | ts_vdw_isolated | xdm | xdm_a1 | xdm_a2 | space_group | uniqueb | origin_choice | rhombohedral | zgate | relaxz | block | block_1 | block_2 | block_height +ibrav | celldm | A | B | C | cosAB | cosAC | cosBC | nat | ntyp | nbnd | tot_charge | starting_charge | tot_magnetization | starting_magnetization | ecutwfc | ecutrho | ecutfock | nr1 | nr2 | nr3 | nr1s | nr2s | nr3s | nosym | nosym_evc | noinv | no_t_rev | force_symmorphic | use_all_frac | occupations | one_atom_occupations | starting_spin_angle | degauss | smearing | nspin | noncolin | ecfixed | qcutz | q2sigma | input_dft | exx_fraction | screening_parameter | exxdiv_treatment | x_gamma_extrapolation | ecutvcut | nqx1 | nqx2 | nqx3 | lda_plus_u | lda_plus_u_kind | Hubbard_U | Hubbard_J0 | Hubbard_alpha | Hubbard_beta | Hubbard_J | starting_ns_eigenvalue | U_projection_type | edir | emaxpos | eopreg | eamp | angle1 | angle2 | lforcet | constrained_magnetization | fixed_magnetization | lambda | report | lspinorb | assume_isolated | esm_bc | esm_w | esm_efield | esm_nfit | fcp_mu | vdw_corr | london | london_s6 | london_c6 | london_rvdw | london_rcut | dftd3_version | dftd3_threebody | ts_vdw_econv_thr | ts_vdw_isolated | xdm | xdm_a1 | xdm_a2 | space_group | uniqueb | origin_choice | rhombohedral | zgate | relaxz | block | block_1 | block_2 | block_height
-

&ELECTRONS

+

&ELECTRONS

-electron_maxstep | scf_must_converge | conv_thr | adaptive_thr | conv_thr_init | conv_thr_multi | mixing_mode | mixing_beta | mixing_ndim | mixing_fixed_ns | diagonalization | ortho_para | diago_thr_init | diago_cg_maxiter | diago_david_ndim | diago_full_acc | efield | efield_cart | efield_phase | startingpot | startingwfc | tqr +electron_maxstep | scf_must_converge | conv_thr | adaptive_thr | conv_thr_init | conv_thr_multi | mixing_mode | mixing_beta | mixing_ndim | mixing_fixed_ns | diagonalization | ortho_para | diago_thr_init | diago_cg_maxiter | diago_david_ndim | diago_full_acc | efield | efield_cart | efield_phase | startingpot | startingwfc | tqr | real_space
-

&IONS

+

&IONS

-ion_dynamics | ion_positions | pot_extrapolation | wfc_extrapolation | remove_rigid_rot | ion_temperature | tempw | tolp | delta_t | nraise | refold_pos | upscale | bfgs_ndim | trust_radius_max | trust_radius_min | trust_radius_ini | w_1 | w_2 +ion_dynamics | ion_positions | pot_extrapolation | wfc_extrapolation | remove_rigid_rot | ion_temperature | tempw | tolp | delta_t | nraise | refold_pos | upscale | bfgs_ndim | trust_radius_max | trust_radius_min | trust_radius_ini | w_1 | w_2
-

&CELL

+

&CELL

-cell_dynamics | press | wmass | cell_factor | press_conv_thr | cell_dofree +cell_dynamics | press | wmass | cell_factor | press_conv_thr | cell_dofree
-

ATOMIC_SPECIES

+

ATOMIC_SPECIES

-X | Mass_X | PseudoPot_X +X | Mass_X | PseudoPot_X
-

ATOMIC_POSITIONS

+

ATOMIC_POSITIONS

-X | x | y | z | if_pos(1) | if_pos(2) | if_pos(3) +X | x | y | z | if_pos(1) | if_pos(2) | if_pos(3)
-

K_POINTS

+

K_POINTS

-nks | xk_x | xk_y | xk_z | wk | nk1 | nk2 | nk3 | sk1 | sk2 | sk3 +nks | xk_x | xk_y | xk_z | wk | nk1 | nk2 | nk3 | sk1 | sk2 | sk3
-

CELL_PARAMETERS

+

CELL_PARAMETERS

-v1 | v2 | v3 +v1 | v2 | v3
-

CONSTRAINTS

+

CONSTRAINTS

-nconstr | constr_tol | constr_type | constr(1) | constr(2) | constr(3) | constr(4) | constr_target +nconstr | constr_tol | constr_type | constr(1) | constr(2) | constr(3) | constr(4) | constr_target
-

OCCUPATIONS

+

OCCUPATIONS

-f_inp1 | f_inp2 +f_inp1 | f_inp2
-

ATOMIC_FORCES

+

ATOMIC_FORCES

-X | fx | fy | fz +X | fx | fy | fz
-

INTRODUCTION

+

INTRODUCTION

 Input data format: { } = optional, [ ] = it depends, | = or
 
@@ -190,11 +190,11 @@ if (not automatic)
    label_n Fx(n) Fy(n) Fz(n) ]
    
- +

Namelist: &CONTROL

- +
@@ -249,7 +249,7 @@ A string describing the task to be performed. Options are:
calculation CHARACTER
- +
@@ -264,7 +264,7 @@ reprinted on output.
title CHARACTER
- +
@@ -295,7 +295,7 @@ Currently two verbosity levels are implemented:
verbosity CHARACTER
- +
@@ -329,7 +329,7 @@ and startingpot.
restart_mode CHARACTER
- +
@@ -356,7 +356,7 @@ Note that this flag has no effect on reading, only on writing.
wf_collect LOGICAL
- +
@@ -374,7 +374,7 @@ performed in this run
nstep INTEGER
- +
@@ -389,7 +389,7 @@ band energies are written every iprint iterations
iprint INTEGER
- +
@@ -405,7 +405,7 @@ calculate stress. It is set to .TRUE. automatically if
tstress LOGICAL
- +
@@ -416,7 +416,7 @@ calculate forces. It is set to .TRUE. automatically if
tprnfor LOGICAL
- +
@@ -433,7 +433,7 @@ time step for molecular dynamics, in Rydberg atomic units
dt REAL
- +
@@ -451,7 +451,7 @@ see also wfcdir
outdir CHARACTER
- +
@@ -473,7 +473,7 @@ may need to copy files to outdir. Works only for pw.x.
wfcdir CHARACTER
- +
@@ -489,7 +489,7 @@ prefix.wfc, prefix.rho, etc.
prefix CHARACTER
- +
@@ -507,7 +507,7 @@ doesn't work together with wf_collect
lkpoint_dir LOGICAL
- +
@@ -525,7 +525,7 @@ fit into your batch queues.
max_seconds REAL
- +
@@ -545,7 +545,7 @@ See also forc_conv_thr - both criteria must be sati
etot_conv_thr REAL
- +
@@ -563,7 +563,7 @@ See also etot_conv_thr - both criteria must be sati
forc_conv_thr REAL
- +
@@ -615,7 +615,7 @@ but you cannot restart in disk_io=='none'
disk_io CHARACTER
- +
@@ -632,7 +632,7 @@ directory containing pseudopotential files
pseudo_dir CHARACTER
- +
@@ -650,7 +650,7 @@ the added potential.
tefield LOGICAL
- +
@@ -670,7 +670,7 @@ with the discontinuity FALLING IN THE EMPTY SPACE.
dipfield LOGICAL
- +
@@ -687,7 +687,7 @@ This is different from tefield == .true. !
lelfield LOGICAL
- +
@@ -705,7 +705,7 @@ external iteration on the charge density
nberrycyc INTEGER
- +
@@ -727,7 +727,7 @@ Works ONLY with norm-conserving pseudopotentials.
lorbm LOGICAL
- +
@@ -743,7 +743,7 @@ See the header of PW/src/bp_c_phase.f90 for documentation.
lberry LOGICAL
- +
@@ -757,7 +757,7 @@ For calculations with finite electric fields
gdir INTEGER
- +
@@ -770,7 +770,7 @@ The same for calculation with finite electric fields
nppstr INTEGER
- +
@@ -798,7 +798,7 @@ NB:
lfcpopt LOGICAL
- +
@@ -832,11 +832,11 @@ Note: this option was called "monopole" in v6.0 and 6.1 of pw.x
gate LOGICAL
- +

Namelist: &SYSTEM

- +
@@ -968,7 +968,7 @@ ibrav structure celldm(2)-celldm(6)
ibrav INTEGER

Either:

- +
@@ -992,7 +992,7 @@ cell vectors are read from card CELL_PARAMETERS
celldm(i), i=1,6 REAL
+A, B, C, cosAB, cosAC, cosBC @@ -1019,7 +1019,7 @@ cell vectors are read from card CELL_PARAMETERS.
-A, B, C, cosAB, cosAC, cosBC REAL
- +
@@ -1035,7 +1035,7 @@ space_group is set, in which case, INEQUIVALENT atoms)
nat INTEGER
- +
@@ -1050,7 +1050,7 @@ number of types of atoms in the unit cell
ntyp INTEGER
- +
@@ -1070,7 +1070,7 @@ k-point, not the number of bands per k-point, is doubled
nbnd INTEGER
- +
@@ -1091,7 +1091,7 @@ inserted to remove divergences if the cell is not neutral.
tot_charge REAL
- +
@@ -1107,7 +1107,7 @@ to create starting potential with startingpot = 'atom
starting_charge(i), i=1,ntyp REAL
- +
@@ -1126,7 +1126,7 @@ the self-consistent cycle.
tot_magnetization REAL
- +
@@ -1157,7 +1157,7 @@ kept zero (the internal variable domag is .FALSE.).
starting_magnetization(i), i=1,ntyp REAL
- +
@@ -1172,7 +1172,7 @@ kinetic energy cutoff (Ry) for wavefunctions
ecutwfc REAL
- +
@@ -1198,7 +1198,7 @@ to be accurately converged.
ecutrho REAL
- +
@@ -1222,7 +1222,7 @@ to instabilities.
ecutfock REAL
+nr1, nr2, nr3
-nr1, nr2, nr3 INTEGER
@@ -1238,7 +1238,7 @@ be used.
       
+nr1s, nr2s, nr3s
-nr1s, nr2s, nr3s INTEGER
@@ -1250,7 +1250,7 @@ be used.
          
- +
@@ -1282,7 +1282,7 @@ and what you get. May be useful in the following cases:
nosym LOGICAL
- +
@@ -1308,7 +1308,7 @@ in all cases to have the full symmetry of the Bravais lattice
nosym_evc LOGICAL
- +
@@ -1324,7 +1324,7 @@ if (.TRUE.) disable the usage of k => -k symmetry
noinv LOGICAL
- +
@@ -1340,7 +1340,7 @@ that consist in a rotation + time reversal.
no_t_rev LOGICAL
- +
@@ -1356,7 +1356,7 @@ symmetry operations having an associated fractionary translation
force_symmorphic LOGICAL
- +
@@ -1377,7 +1377,7 @@ in real space.
use_all_frac LOGICAL
- +
@@ -1439,7 +1439,7 @@ with the value of tot_charge.
occupations CHARACTER
- +
@@ -1481,7 +1481,7 @@ can be forced to be spin-angle functions by setting
one_atom_occupations LOGICAL
- +
@@ -1508,7 +1508,7 @@ When lspinorb is .FALSE. this flag is not used.
starting_spin_angle LOGICAL
- +
@@ -1524,7 +1524,7 @@ integration in metals.
degauss REAL
- +
@@ -1567,7 +1567,7 @@ smearing with Fermi-Dirac function
smearing CHARACTER
- +
@@ -1590,7 +1590,7 @@ nspin = 4 : spin-polarized calculation, noncollinear
nspin INTEGER
- +
@@ -1605,7 +1605,7 @@ if .true. the program will perform a noncollinear calculation.
noncolin LOGICAL
- +
@@ -1621,7 +1621,7 @@ if .true. the program will perform a noncollinear calculation.
ecfixed REAL
- +
@@ -1637,7 +1637,7 @@ if .true. the program will perform a noncollinear calculation.
qcutz REAL
- +
@@ -1660,7 +1660,7 @@ See: M. Bernasconi et al, J. Phys. Chem. Solids 56, 501 (1995),
q2sigma REAL
- +
@@ -1678,7 +1678,7 @@ Use with care and if you know what you are doing!
input_dft CHARACTER
- +
@@ -1695,7 +1695,7 @@ the exx_fraction default value is 0.20.
exx_fraction REAL
- +
@@ -1713,7 +1713,7 @@ J. Chem. Phys. 124, 219906 (2006),
screening_parameter REAL
- +
@@ -1751,7 +1751,7 @@ for treating the Coulomb potential divergencies at small q vectors.
exxdiv_treatment CHARACTER
- +
@@ -1768,7 +1768,7 @@ Set this to .false. for GAU-PBE.
x_gamma_extrapolation LOGICAL
- +
@@ -1791,7 +1791,7 @@ divergencies at small q vectors.
ecutvcut REAL
+nqx1, nqx2, nqx3
-nqx1, nqx2, nqx3 INTEGER
@@ -1804,7 +1804,7 @@ In QE =< 5.0.2 it defaulted to nqx1=nqx2=nqx3=1.
          
- +
@@ -1834,7 +1834,7 @@ see lda_plus_u_kind, Hubbard
lda_plus_u LOGICAL
- +
@@ -1855,7 +1855,7 @@ Specifies the type of DFT+U calculation:
lda_plus_u_kind INTEGER
- +
@@ -1870,7 +1870,7 @@ Hubbard_U(i): U parameter (eV) for species i, DFT+U calculation
Hubbard_U(i), i=1,ntyp REAL
- +
@@ -1886,7 +1886,7 @@ see P
Hubbard_J0(i), i=1,ntype REAL
- +
@@ -1904,7 +1904,7 @@ Cococcioni and de Gironcoli, [Back to Top] -
Hubbard_alpha(i), i=1,ntyp REAL
+
@@ -1923,7 +1923,7 @@ Cococcioni and de Gironcoli, [Back to Top] -
Hubbard_beta(i), i=1,ntyp REAL
+
@@ -1946,7 +1946,7 @@ calculated from J using atomic ratios.
Hubbard_J(i,ityp), (i,ityp) = (1,1) ... (3,ntyp)
- +
@@ -1970,7 +1970,7 @@ when the default choice takes another path.
starting_ns_eigenvalue(m,ispin,ityp), (m,ispin,ityp) = (1,1,1) ... (2*lmax+1,nspin or npol,ntyp)
- +
@@ -2035,7 +2035,7 @@ NB: forces and stress currently implemented only for the
U_projection_type CHARACTER
- +
@@ -2048,7 +2048,7 @@ potential is constant in planes defined by FFT grid points;
edir INTEGER
- +
@@ -2065,7 +2065,7 @@ Used only if tefield is .TRUE.
emaxpos REAL
- +
@@ -2081,7 +2081,7 @@ Zone in the unit cell where the saw-like potential decreases.
eopreg REAL
- +
@@ -2103,7 +2103,7 @@ unphysical forces will result.
eamp REAL
- +
@@ -2115,7 +2115,7 @@ only; index i runs over the atom types.
angle1(i), i=1,ntyp REAL
- +
@@ -2127,7 +2127,7 @@ For noncollinear calculations only.
angle2(i), i=1,ntyp REAL
- +
@@ -2140,7 +2140,7 @@ When starting a non collinear calculation using an existing density
lforcet LOGICAL
- +
@@ -2225,7 +2225,7 @@ configurations or use the nosym flag (only as a last remedy)
constrained_magnetization CHARACTER
- +
@@ -2245,7 +2245,7 @@ fixed when constrained_magnetization=='
fixed_magnetization(i), i=1,3 REAL
- +
@@ -2266,7 +2266,7 @@ N.B.: if the scf calculation does not converge, try to reduce lambda
lambda REAL
- +
@@ -2282,7 +2282,7 @@ writes all the atomic magnetic moments.
report INTEGER
- +
@@ -2293,7 +2293,7 @@ spin-orbit.
lspinorb LOGICAL
- +
@@ -2411,7 +2411,7 @@ NB:
assume_isolated CHARACTER
- +
@@ -2457,7 +2457,7 @@ See also esm_efield.
esm_bc CHARACTER
- +
@@ -2479,7 +2479,7 @@ z = +/- [L_z/2 + esm_w] ).
esm_w REAL
- +
@@ -2500,7 +2500,7 @@ between semi-infinite ESM electrodes.
esm_efield REAL
- +
@@ -2520,7 +2520,7 @@ for the polynomial fit along the cell edge.
esm_nfit INTEGER
- +
@@ -2540,7 +2540,7 @@ with appropriate total charge of the system by giving 'tot_charge'.
fcp_mu REAL
- +
@@ -2600,7 +2600,7 @@ A. Otero de la Roza et al., J. Chem. Phys. 136, 174109 (2012),
vdw_corr CHARACTER
- +
@@ -2618,7 +2618,7 @@ OBSOLESCENT, same as vdw_corr='DFT-D'
london LOGICAL
- +
@@ -2633,7 +2633,7 @@ global scaling parameter for DFT-D. Default is good for PBE.
london_s6 REAL
- +
@@ -2651,7 +2651,7 @@ atomic C6 coefficient of each atom type
london_c6(i), i=1,ntyp REAL
- +
@@ -2669,7 +2669,7 @@ atomic vdw radii of each atom type
london_rvdw(i), i=1,ntyp REAL
- +
@@ -2684,7 +2684,7 @@ cutoff radius (a.u.) for dispersion interactions
london_rcut REAL
- +
@@ -2701,7 +2701,7 @@ parametrized.
dftd3_version integer
- +
@@ -2718,7 +2718,7 @@ If dftd3_version=2, three-body contribution is always disabled.
dftd3_threebody LOGICAL
- +
@@ -2734,7 +2734,7 @@ is a safe choice, likely too safe, but you do not gain much in increasing it
ts_vdw_econv_thr REAL
- +
@@ -2750,7 +2750,7 @@ for an isolated (non-periodic) system.
ts_vdw_isolated LOGICAL
- +
@@ -2768,7 +2768,7 @@ OBSOLESCENT, same as vdw_corr='xdm'
xdm LOGICAL
- +
@@ -2788,7 +2788,7 @@ in this list, the coefficients are given in:
xdm_a1 REAL
- +
@@ -2808,7 +2808,7 @@ in this list, the coefficients are given in:
xdm_a2 REAL
- +
@@ -2829,7 +2829,7 @@ are of type crystal_sg. See also uniqueb,
space_group INTEGER
- +
@@ -2848,7 +2848,7 @@ b axis. If .FALSE. it is parallel to the c axis.
uniqueb LOGICAL
- +
@@ -2866,7 +2866,7 @@ second origin.
origin_choice INTEGER
- +
@@ -2889,7 +2889,7 @@ to the rhombohedral axes and ibrav=5 is used in both cases.
rhombohedral LOGICAL

variables used only if gate = .TRUE.

- +
@@ -2909,7 +2909,7 @@ T. Brumme, M. Calandra, F. Mauri; [Back to Top] -
zgate REAL
+
@@ -2929,7 +2929,7 @@ correction (dipfield=.TRUE.).
relaxz LOGICAL
- +
@@ -2951,7 +2951,7 @@ decrease of the potential barrier.
block LOGICAL
- +
@@ -2968,7 +2968,7 @@ unit cell size along z,block_1 in ]0,1[
block_1 REAL
- +
@@ -2985,7 +2985,7 @@ unit cell size along z,block_2 in ]0,1[
block_2 REAL
- +
@@ -3004,11 +3004,11 @@ Height of the potential barrier in Rydberg.
block_height REAL
- +

Namelist: &ELECTRONS

- +
@@ -3023,7 +3023,7 @@ maximum number of iterations in a scf step
electron_maxstep INTEGER
- +
@@ -3039,7 +3039,7 @@ when electron_maxstep is reached. Use with care.
scf_must_converge LOGICAL
- +
@@ -3060,7 +3060,7 @@ iterative diagonalizazion: see diago_thr_init
conv_thr REAL
- +
@@ -3076,7 +3076,7 @@ the inner scf loops when using EXX.
adaptive_thr LOGICAL
- +
@@ -3092,7 +3092,7 @@ used for the first scf cycle.
conv_thr_init REAL
- +
@@ -3109,7 +3109,7 @@ max( conv_thr, conv_thr_multi<
conv_thr_multi REAL
- +
@@ -3144,7 +3144,7 @@ as above, with local-density-dependent TF screening
mixing_mode CHARACTER
- +
@@ -3159,7 +3159,7 @@ mixing factor for self-consistency
mixing_beta REAL
- +
@@ -3175,7 +3175,7 @@ If you are tight with memory, you may reduce it to 4 or so.
mixing_ndim INTEGER
- +
@@ -3191,7 +3191,7 @@ atomic density appearing in the Hubbard term ).
mixing_fixed_ns INTEGER
- +
@@ -3236,7 +3236,7 @@ on linear algebra.
diagonalization CHARACTER
- +
@@ -3253,7 +3253,7 @@ on linear algebra.
ortho_para INTEGER
- +
@@ -3272,7 +3272,7 @@ For non-scf calculations: default is (conv_thr/N elec)/1
diago_thr_init REAL
- +
@@ -3282,7 +3282,7 @@ For conjugate gradient diagonalization: max number of iterations
diago_cg_maxiter INTEGER
- +
@@ -3304,7 +3304,7 @@ is significant compared to the time spent in h_psi
diago_david_ndim INTEGER
- +
@@ -3322,7 +3322,7 @@ total energy, forces, and other ground-state properties).
diago_full_acc LOGICAL
- +
@@ -3339,7 +3339,7 @@ and if k-points (K_POINTS card) are not automatic.
efield REAL
- +
@@ -3356,7 +3356,7 @@ k-points (K_POINTS card) are automatic.
efield_cart(i), i=1,3 REAL
- +
@@ -3392,7 +3392,7 @@ none of the above points
efield_phase CHARACTER
- +
@@ -3419,7 +3419,7 @@ and the only sensible possibility.
startingpot CHARACTER
- +
@@ -3469,7 +3469,7 @@ directory specified by variables prefix and [Back to Top] -
startingwfc CHARACTER
+
@@ -3480,24 +3480,42 @@ directory specified by variables prefix and
-If .true., use the real-space algorithm for augmentation
-charges in ultrasoft pseudopotentials.
-Must faster execution of ultrasoft-related calculations,
-but numerically less accurate than the default algorithm.
-Use with care and after testing!
+If .true., use a real-space algorithm for augmentation
+charges of ultrasoft pseudopotentials and PAWsets.
+Faster but numerically less accurate than the default
+G-space algorithm. Use with care and after testing!
+         
+
tqr LOGICAL
+ + + + + + + + + + +
real_spaceLOGICAL
Default: .FALSE. +
+If .true., exploit real-space localization to compute
+matrix elements for nonlocal projectors. Faster and in
+principle better scaling than the default G-space algorithm,
+but numerically less accurate, may lead to some loss of
+translational invariance. Use with care and after testing!
          
- +

Namelist: &IONS

input this namelist only if calculation == 'relax', 'md', 'vc-relax', or 'vc-md'

- +
@@ -3582,7 +3600,7 @@ Newton's equation
ion_dynamics CHARACTER
- +
@@ -3613,7 +3631,7 @@ from standard input, even if restarting.
ion_positions CHARACTER
- +
@@ -3659,7 +3677,7 @@ only for molecular dynamics calculations
pot_extrapolation CHARACTER
- +
@@ -3697,7 +3715,7 @@ only for molecular dynamics calculations
wfc_extrapolation CHARACTER
- +
@@ -3728,7 +3746,7 @@ but the relaxation itself should be faster.

variables used for molecular dynamics

-
remove_rigid_rot LOGICAL
+
@@ -3801,7 +3819,7 @@ and leave uncontrolled further on
ion_temperature CHARACTER
- +
@@ -3817,7 +3835,7 @@ target temperature for most thermostats.
tempw REAL
- +
@@ -3833,7 +3851,7 @@ the run-averaged and target temperature differ more than tolp.
tolp REAL
- +
@@ -3863,7 +3881,7 @@ equipartition between all quadratic degrees of freedom).
delta_t REAL
- +
@@ -3896,7 +3914,7 @@ if ion_temperature == 'andersen' :
nraise INTEGER
- +
@@ -3918,7 +3936,7 @@ the supercell.

keywords used only in BFGS calculations

-
refold_pos LOGICAL
+
@@ -3936,7 +3954,7 @@ but conv_thr will not be reduced to less that [Back to Top] -
upscale REAL
+
@@ -3956,7 +3974,7 @@ used.
bfgs_ndim INTEGER
- +
@@ -3972,7 +3990,7 @@ Maximum ionic displacement in the structural relaxation.
trust_radius_max REAL
- +
@@ -3989,7 +4007,7 @@ BFGS is reset when trust_radius < [Back to Top] -
trust_radius_min REAL
+
@@ -4005,7 +4023,7 @@ Initial ionic displacement in the structural relaxation.
trust_radius_ini REAL
- +
@@ -4021,7 +4039,7 @@ Initial ionic displacement in the structural relaxation.
w_1 REAL
- +
@@ -4040,14 +4058,14 @@ Parameters used in line search based on the Wolfe conditions.
w_2 REAL
- +

Namelist: &CELL

input this namelist only if calculation == 'vc-relax' or 'vc-md'

- +
@@ -4112,7 +4130,7 @@ BFGS quasi-newton algorithm (default)
cell_dynamics CHARACTER
- +
@@ -4127,7 +4145,7 @@ Target pressure [KBar] in a variable-cell md or relaxation run.
press REAL
- +
@@ -4145,7 +4163,7 @@ Fictitious cell mass [amu] for variable-cell simulations
wmass REAL
- +
@@ -4162,7 +4180,7 @@ cell during a simulation.
cell_factor REAL
- +
@@ -4179,7 +4197,7 @@ relaxation ('vc-relax' : note that the other convergence
press_conv_thr REAL
- +
@@ -4278,7 +4296,7 @@ BEWARE: if axis are not orthogonal, some of these options do not
cell_dofree CHARACTER
- +

Card: ATOMIC_SPECIES

@@ -4286,22 +4304,22 @@ BEWARE: if axis are not orthogonal, some of these options do not

Syntax:

ATOMIC_SPECIES
-
+
- - - + + + - - - + + + - - - + + +
 X(1)  Mass_X(1)  PseudoPot_X(1)  X(1)  Mass_X(1)  PseudoPot_X(1) 
 X(2)  Mass_X(2)  PseudoPot_X(2)  X(2)  Mass_X(2)  PseudoPot_X(2) 
 . . .
 X(ntyp)  Mass_X(ntyp)  PseudoPot_X(ntyp)  X(ntyp)  Mass_X(ntyp)  PseudoPot_X(ntyp) 
@@ -4312,7 +4330,7 @@ BEWARE: if axis are not orthogonal, some of these options do not
- +
XX CHARACTER
@@ -4326,7 +4344,7 @@ max total length cannot exceed 3 characters)
 
 
-
+
Mass_XMass_X REAL
@@ -4341,7 +4359,7 @@ these values unless other values are provided)
 
 
-
+
PseudoPot_XPseudoPot_X CHARACTER
@@ -4361,7 +4379,7 @@ none of the above  old PWscf norm-conserving format
 
- +

Card: ATOMIC_POSITIONS { alat | bohr | angstrom | crystal | crystal_sg }

@@ -4377,39 +4395,39 @@ previous scf calculation will be used instead !!!
ATOMIC_POSITIONS { alat | bohr | angstrom | crystal | crystal_sg }
-
+
- - - - + + + + - - - + + + - - - - + + + + - - - + + + - - - - + + + + - - - + + +
 X(1)  x(1)  y(1)  z(1)  X(1)  x(1)  y(1)  z(1)  {  if_pos(1)(1)  if_pos(2)(1)  if_pos(3)(1)  if_pos(1)(1)  if_pos(2)(1)  if_pos(3)(1)  }
 X(2)  x(2)  y(2)  z(2)  X(2)  x(2)  y(2)  z(2)  {  if_pos(1)(2)  if_pos(2)(2)  if_pos(3)(2)  if_pos(1)(2)  if_pos(2)(2)  if_pos(3)(2)  }
 . . .
 X(nat)  x(nat)  y(nat)  z(nat)  X(nat)  x(nat)  y(nat)  z(nat)  {  if_pos(1)(nat)  if_pos(2)(nat)  if_pos(3)(nat)  if_pos(1)(nat)  if_pos(2)(nat)  if_pos(3)(nat)  }
@@ -4420,7 +4438,7 @@ previous scf calculation will be used instead !!!

Description of items:

- +
Card's options: @@ -4487,7 +4505,7 @@ rhombohedral. - +
XX CHARACTER
 label of the atom as specified in ATOMIC_SPECIES
@@ -4497,7 +4515,7 @@ rhombohedral.
 
@@ -4548,7 +4566,7 @@ NOTE: each atomic coordinate can also be specified as a simple algebraic express
 
-x, y, z +x, y, z REAL
@@ -4572,7 +4590,7 @@ atoms.
-if_pos(1), if_pos(2), if_pos(3) +if_pos(1), if_pos(2), if_pos(3) INTEGER
- +

Card: K_POINTS { tpiba | automatic | crystal | gamma | tpiba_b | crystal_b | tpiba_c | crystal_c }

@@ -4581,25 +4599,25 @@ atoms.

Syntax:

K_POINTS tpiba | crystal | tpiba_b | crystal_b | tpiba_c | crystal_c
-nks  
+nks  
- - - - + + + + - - - - + + + + - - - - + + + +
 xk_x(1)  xk_y(1)  xk_z(1)  wk(1)  xk_x(1)  xk_y(1)  xk_z(1)  wk(1) 
 xk_x(2)  xk_y(2)  xk_z(2)  wk(2)  xk_x(2)  xk_y(2)  xk_z(2)  wk(2) 
 . . .
 xk_x(nks)  xk_y(nks)  xk_z(nks)  wk(nks)  xk_x(nks)  xk_y(nks)  xk_z(nks)  wk(nks) 
@@ -4609,7 +4627,7 @@ atoms.

Syntax:

K_POINTS automatic
-nk1  nk2  nk3  sk1  sk2  sk3  
+nk1  nk2  nk3  sk1  sk2  sk3  
@@ -4623,7 +4641,7 @@ atoms.

Description of items:

- +
Card's options: @@ -4716,7 +4734,7 @@ As tpiba_c, but k-points are in crystal coordinates.
- +
@@ -4728,7 +4746,7 @@ As tpiba_c, but k-points are in crystal coordinates.
nks INTEGER
@@ -4753,7 +4771,7 @@ for a band-structure plot), weights can be set to any value
-xk_x, xk_y, xk_z, wk +xk_x, xk_y, xk_z, wk REAL
+nk1, nk2, nk3
-nk1, nk2, nk3 INTEGER
@@ -4765,7 +4783,7 @@ These parameters specify the k-point grid
 
+sk1, sk2, sk3
-sk1, sk2, sk3 INTEGER
@@ -4779,7 +4797,7 @@ half a grid step in the corresponding direction ).
 
- +

Card: CELL_PARAMETERS { alat | bohr | angstrom }

@@ -4791,26 +4809,26 @@ Optional card, needed only if ibrav == 0 is specified, igno
CELL_PARAMETERS { alat | bohr | angstrom }
-
+
- - - + + + - - - + + + - - - + + +
 v1(1)  v1(2)  v1(3)  v1(1)  v1(2)  v1(3) 
 v2(1)  v2(2)  v2(3)  v2(1)  v2(2)  v2(3) 
 v3(1)  v3(2)  v3(3)  v3(1)  v3(2)  v3(3) 
@@ -4820,7 +4838,7 @@ Optional card, needed only if ibrav == 0 is specified, igno

Description of items:

- +
Card's options: @@ -4847,7 +4865,7 @@ If neither unit nor lattice parameter are specified, @@ -4863,7 +4881,7 @@ Crystal lattice vectors (in cartesian axis):
-v1, v2, v3 +v1, v2, v3 REAL
- +

Card: CONSTRAINTS

@@ -4878,11 +4896,11 @@ When this card is present the SHAKE algorithm is automatically used.

Syntax:

CONSTRAINTS
-nconstr   { constr_tol   }
+nconstr   { constr_tol   }
- - - + + + @@ -4890,13 +4908,13 @@ When this card is present the SHAKE algorithm is automatically used.   - + - - - + + + @@ -4904,14 +4922,14 @@ When this card is present the SHAKE algorithm is automatically used.   - + - - - + + + @@ -4919,7 +4937,7 @@ When this card is present the SHAKE algorithm is automatically used.   - +
 constr_type(1)  constr(1)(1)  constr(2)(1)  constr_type(1)  constr(1)(1)  constr(2)(1)  [  constr(3)(1)   ] {  constr_target(1)  constr_target(1)  }
 constr_type(2)  constr(1)(2)  constr(2)(2)  constr_type(2)  constr(1)(2)  constr(2)(2)  [  constr(3)(2)   ] {  constr_target(2)  constr_target(2)  }
 . . .
 constr_type(nconstr)  constr(1)(nconstr)  constr(2)(nconstr)  constr_type(nconstr)  constr(1)(nconstr)  constr(2)(nconstr)  [  constr(3)(nconstr)   ] {  constr_target(nconstr)  constr_target(nconstr)  }
@@ -4929,7 +4947,7 @@ When this card is present the SHAKE algorithm is automatically used.

Description of items:

- +
@@ -4938,7 +4956,7 @@ When this card is present the SHAKE algorithm is automatically used.
nconstr INTEGER
- +
@@ -4949,7 +4967,7 @@ When this card is present the SHAKE algorithm is automatically used.
constr_tol REAL
- +
constr_typeconstr_type CHARACTER
@@ -5013,7 +5031,7 @@ G. Roma, J.P. Crocombette: J. Nucl. Mater. 403, 32 (2010), @@ -5049,7 +5067,7 @@ These variables have different meanings for different constraint types:
-constr(1), constr(2), constr(3), constr(4) +constr(1), constr(2), constr(3), constr(4)
- +
constr_targetconstr_target REAL
@@ -5062,7 +5080,7 @@ This variable is optional.
 
- +

Card: OCCUPATIONS

@@ -5072,21 +5090,21 @@ This variable is optional.

Syntax:

OCCUPATIONS
-
+
- - + + - + - - + + - +
 f_inp1(1)  f_inp1(2)  f_inp1(1)  f_inp1(2)   . . . f_inp1(nbnd)  f_inp1(nbnd) 
[   f_inp2(1)  f_inp2(2)  f_inp2(1)  f_inp2(2)   . . . f_inp2(nbnd)  f_inp2(nbnd)    ]
@@ -5098,7 +5116,7 @@ This variable is optional.
- +
f_inp1f_inp1 REAL
@@ -5109,7 +5127,7 @@ For spin-polarized calculations, these are majority spin states.
 
 
-
+
f_inp2f_inp2 REAL
@@ -5122,7 +5140,7 @@ To be specified only for spin-polarized calculations.
 
- +

Card: ATOMIC_FORCES

@@ -5136,25 +5154,25 @@ BEWARE: if the sum of external forces is not zero, the center of mass of

Syntax:

ATOMIC_FORCES
-
+
- - - - + + + + - - - - + + + + - - - - + + + +
 X(1)  fx(1)  fy(1)  fz(1)  X(1)  fx(1)  fy(1)  fz(1) 
 X(2)  fx(2)  fy(2)  fz(2)  X(2)  fx(2)  fy(2)  fz(2) 
 . . .
 X(nat)  fx(nat)  fy(nat)  fz(nat)  X(nat)  fx(nat)  fy(nat)  fz(nat) 
@@ -5165,7 +5183,7 @@ BEWARE: if the sum of external forces is not zero, the center of mass of
- +
XX CHARACTER
 label of the atom as specified in ATOMIC_SPECIES
@@ -5175,7 +5193,7 @@ BEWARE: if the sum of external forces is not zero, the center of mass of
 
@@ -5191,7 +5209,7 @@ external force on atom X (cartesian components, Ry/a.u. units)
 
-fx, fy, fz +fx, fy, fz REAL
- This file has been created by helpdoc utility on Mon Jul 23 13:55:31 CEST 2018. + This file has been created by helpdoc utility on Tue Aug 21 17:08:54 CEST 2018. diff --git a/PW/Doc/INPUT_PW.txt b/PW/Doc/INPUT_PW.txt index 04c5ebbe2..837f82716 100644 --- a/PW/Doc/INPUT_PW.txt +++ b/PW/Doc/INPUT_PW.txt @@ -2193,11 +2193,22 @@ NAMELIST: &ELECTRONS Type: LOGICAL Default: .FALSE. - Description: If .true., use the real-space algorithm for augmentation - charges in ultrasoft pseudopotentials. - Must faster execution of ultrasoft-related calculations, - but numerically less accurate than the default algorithm. - Use with care and after testing! + Description: If .true., use a real-space algorithm for augmentation + charges of ultrasoft pseudopotentials and PAWsets. + Faster but numerically less accurate than the default + G-space algorithm. Use with care and after testing! + +-------------------------------------------------------------------- + + +-------------------------------------------------------------------- + Variable: real_space + + Type: LOGICAL + Default: .FALSE. + Description: If .true., exploit real-space localization to compute + matrix elements for nonlocal projectors. Faster and in + principle better scaling than the default G-space algorithm, + but numerically less accurate, may lead to some loss of + translational invariance. Use with care and after testing! +-------------------------------------------------------------------- ===END OF NAMELIST====================================================== @@ -3289,4 +3300,4 @@ CARD: ATOMIC_FORCES ===END OF CARD========================================================== -This file has been created by helpdoc utility on Mon Jul 23 13:55:31 CEST 2018 +This file has been created by helpdoc utility on Tue Aug 21 17:08:53 CEST 2018