Search found 75 matches

by roginovicci
Thu Nov 29, 2018 6:52 am
Forum: Input file
Topic: Reduce memory demand in rf
Replies: 3
Views: 781

Re: Reduce memory demand in rf

Dear Eric, The nonlinear optical calculations by "optic" utility (SHG) require rather hi k-mesh and ecut to converge (which is also explained in tutorial m). I've found the convergence by ecut is 30Ha and k-mesh is 12x12x8 for my system. I know that the optical properties calculation (in d...
by roginovicci
Wed Nov 21, 2018 11:03 am
Forum: Input file
Topic: Reduce memory demand in rf
Replies: 3
Views: 781

Reduce memory demand in rf

Hi, To perform optical calculations one needs to perform ddk response function calculations ( rfelfd > 0 ) using a big set of empty bands ( nband > 300), which leads to huge amount of memory demand. For crystal with tetragonal symmetry and 10 atoms in unit cell the required memory (RAM) is more than...
by roginovicci
Tue Nov 06, 2018 6:04 pm
Forum: Configuration
Topic: Compilation fails with intel 19 compilator
Replies: 3
Views: 1141

Re: Compilation fails with intel 19 compilator

Dear gmatteo,
The answer is positive. without bigdft the compilations do not falls into compiler crush.
by roginovicci
Fri Nov 02, 2018 2:57 pm
Forum: Configuration
Topic: Compilation fails with intel 19 compilator
Replies: 3
Views: 1141

Re: Compilation fails with intel 19 compilator

Thank you gmatteo, I will try and give you know. Also I would like to let you know, that with vesion 2018 update 1 of intel compilator there is no such problem.
by roginovicci
Sun Oct 28, 2018 1:49 pm
Forum: Configuration
Topic: Compilation fails with intel 19 compilator
Replies: 3
Views: 1141

Compilation fails with intel 19 compilator

Hello, I've found that abinit source code could not be compiled with the recent version of Intel compilator. The following error emerges: /opt/intel//impi/2019.0.117/intel64/bin/mpiifort -I. -I. -I../../includes -I../../libABINIT/src -I/distrib/abinit/abinit-8.10.1/fallbacks/exports/include -free -f...
by roginovicci
Fri Jun 23, 2017 1:37 pm
Forum: GW
Topic: How to calculate full dielectric tensor within GWA
Replies: 1
Views: 2075

Re: How to calculate full dielectric tensor within GWA

just to make it clear. The components of dielectric tensor could be separately calculated by setting gw_qlwl variable ( for instance 0.001 0 0 for \varepsilon_xx) and inspecting spectra in _EM1_LF _EM1_NLF files
by roginovicci
Thu May 12, 2016 4:10 pm
Forum: GW
Topic: possible bug with icutcoul=2
Replies: 9
Views: 4849

Re: possible bug with icutcoul=2

All my hugs for your suggestion Yannick! But while changing ecuteps/ecutsigx parameters do influence on contour deformation method (gwcalctyp = 2) it fails with Hartree-Fock calculation (gwcalctyp = 5). I've found doubled FFT grid (fftgw =31) do cover all g-points, and that is probably the solution ...
by roginovicci
Thu May 12, 2016 12:16 pm
Forum: GW
Topic: possible bug with icutcoul=2
Replies: 9
Views: 4849

Re: possible bug with icutcoul=2

Dear Yannick! Thank you for the info, I'll check it out in a while. But, we still have a possibility to use a spheric Coulomb cut-off icutcoul=0 and rcut=a/2 (a is a lattice constant) in rhombohedral lattices. Why it is not possible to use icutcoul=0 with hex lattice? Sorry for (probably) stupid que...
by roginovicci
Thu May 12, 2016 8:52 am
Forum: GW
Topic: possible bug with icutcoul=2
Replies: 9
Views: 4849

Re: possible bug with icutcoul=2

I'v got the same Error in Sigma calculation when treating hexagonal graphene surface: --- !ERROR message: | Found q-points with non-zero component along non-periodic direction This is not allowed, see Notes in cutoff_surface.F90 ACTION : Modify the q-point sampling src_file: m_vcoul.F90 The q-point ...
by roginovicci
Fri Apr 15, 2016 1:53 pm
Forum: Response calculations
Topic: [BUG] Effective charge calculation with PAW fails.
Replies: 0
Views: 1854

[BUG] Effective charge calculation with PAW fails.

Dear developers, Calculation of effective charges always fails with PAW approach. Example of Effective charge tensors calculated with DFPT approach with norm conserving pseudo potentials: atom displacement 1 1 -1.362653E+00 -7.005936E-17 2.734595E-15 1 2 -7.005936E-17 -1.362653E+00 2.734595E-15 1 3 ...
by roginovicci
Fri Apr 08, 2016 12:34 pm
Forum: GW
Topic: How to calculate full dielectric tensor within GWA
Replies: 1
Views: 2075

How to calculate full dielectric tensor within GWA

I know it is possible to get the macroscopic dielectric constant at the calculation of the screening step. That is supposed to be an average value, right? But how to extract each component (xx,yy,zz,xy,xz,yz) of dielectric tensor within GW approximation? Is it possible?
by roginovicci
Thu Apr 07, 2016 5:34 pm
Forum: GW
Topic: self-consistent GW + PAW crash
Replies: 5
Views: 4284

Re: self-consistent GW + PAW crash

I was using datasets from GPAW project. The single shot G0W0 approximation shows similar results as we we obtained with VASP package, thus it seems these datasets is suitable, while there is no (as far as I know) special NC pseudopotential database except several potentials for selected atoms from O...
by roginovicci
Mon Mar 21, 2016 10:42 am
Forum: GW
Topic: self-consistent GW + PAW crash
Replies: 5
Views: 4284

Re: self-consistent GW + PAW crash

I'm sorry I've posted wrong input in previous post. Actually problem arises when full set of k-points required. In my case Abinit crashes when array kptgw contains point 2.50000000E-01 2.50000000E-01 2.50000000E-01 In the screening step calculation that point do exist in q-mesh: ==== Q-mesh for the ...
by roginovicci
Thu Mar 17, 2016 5:14 pm
Forum: GW
Topic: self-consistent GW + PAW crash
Replies: 5
Views: 4284

self-consistent GW + PAW crash

Abinit crashes when performing self-consistent GW calculations within PAW formalism. If gwcalctyp > 10 it fails if <10 goes clear. With NCPP all goes fine Input file: ecut 30 # PAW pawecutdg 60 usexcnhat=-1 pawovlp 20 # DATASET 1 : Sigma calculation # optdriver 4 # Sigma run. irdkss 1 irdscr 1 gwcal...
by roginovicci
Wed Mar 09, 2016 1:08 pm
Forum: GW
Topic: GW nkptgw variable bigger than one possible? [SOLVED]
Replies: 2
Views: 2478

GW nkptgw variable bigger than one possible? [SOLVED]

Hello, I'm trying to perform GW calculations for alpha-TeO2 crystal. My input file is as follows: ecut 25 ndtset 3 # Definition of parameters for the calculation of the KSS file kssform1 3 nbandkss1 298 # Number of bands in KSS file (-1 means the maximum possible) nband1 304 # Number of (occ and emp...
by roginovicci
Fri Nov 06, 2015 11:20 am
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

OK, I've performed calculations with higher ecut=50Ha. Now all acoustic modes is negative and of an order of tenth cm-1 (-45cm-1 max) Phonon frequencies in cm-1 : - -4.589082E+01 -1.724093E+01 -1.724093E+01 5.334798E+01 1.029513E+02 After applying ASR the negative frequencies wanishes and the differ...
by roginovicci
Thu Nov 05, 2015 3:24 pm
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

That's sound strange as my result (for JTH datasets) gives energy convergence of 0.1mHa at ecut=40Ha (k-grid is ngkpt ) ecut1 2.50000000E+01 Hartree ecut2 3.00000000E+01 Hartree ecut3 3.50000000E+01 Hartree ecut4 4.00000000E+01 Hartree ecut5 4.50000000E+01 Hartree etotal1 -1.4590089049E+02 etotal2 -...
by roginovicci
Mon Nov 02, 2015 3:59 pm
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

The paw dataset I've used in calculations above was generated with gpaw project the following way: For oxigen: gpaw-setup -l -f LDA -r 1.1,1.05,1.05 -c "[He]" -w -p --compensation-charge-radius=1. -e "0.78;0.97;0.0" O for silicium: gpaw-setup -l -f LDA -r 1.8,1.9,1.8 -c "[Ne...
by roginovicci
Fri Oct 30, 2015 7:51 pm
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

Some additional info. I'm using JTH PAW datasets, ABINIT complies on high values in Psp strength Dij matrix pawio_print_ij: WARNING - The matrix seems to have high value(s) ! ( 3 components have a value greater than 50.0). It can cause instabilities during SCF convergence. Action: you should check y...
by roginovicci
Fri Oct 30, 2015 4:21 pm
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

Dear Dr. Boris Dorado! Thank you very much for the detailed answer! But concerning your suggestion I've found all my input data is correct. Namely: 1. I've checked in output file each 'Fourteen' section and found the values is pretty close: displacement of atom 1 along direction 1 2DEtotal= 0.414978...
by roginovicci
Mon Oct 26, 2015 4:51 pm
Forum: Response calculations
Topic: PAW+phonons
Replies: 44
Views: 24134

Re: PAW+phonons

Dear Dr. Boris Dorado, Thanks a lot for working on DFPT+PAW problem! I may assume the problem partly (or fully) resolved. Talking about subject structure using LDA PAW datasets taken from JTH project and using acoustic sum rule the negative frequencies vanishes in 7.10.5 version. But the value of ac...
by roginovicci
Mon Sep 15, 2014 1:52 pm
Forum: Response calculations
Topic: symmetry of phonon modes
Replies: 13
Views: 7287

Re: symmetry of phonon modes

I understand from this tutorial that the list of character numbers are given in the same order as the list of symmetry operations given in the log file. I wonder why the symmetry operations are only given in the log file, in my opinion this is important enough information to figure in the out file....
by roginovicci
Tue Jul 15, 2014 11:00 pm
Forum: Platform specific questions
Topic: Problem with cuda [SOLVED]
Replies: 18
Views: 14152

Re: Problem with cuda [SOLVED]

Is it possible to talk about fast configuration when code compiled with gnu C? I've found PGI or Intel compilators produces boost in performance up to 70%. That was a couple years thought. I'm not sure topic starter use linux platform, but I can assume based on my experience that best linux distros ...
by roginovicci
Sat Jul 12, 2014 9:45 pm
Forum: Configuration
Topic: Abinit +intel compilators wierd issue while opt geom [SOLVED]
Replies: 8
Views: 20752

Re: Abinit +intel compilators wierd issue while opt geom [SOLVED]

But if you were not sourcing the intel compiler, what compiler did you use ??? So I guess there was a conflict between compilation time and runtime dynamic libraries... Using ldd command I have found the binary was linked to intel libraries only. no libgfortran in the list. what fallback fails to c...