Sure. I will try it. Thankyou so much.

Thank you so much for your response. I tried it, but it shows an error

GCALCO - MAX INDICES DIRECT LATTICE VECTOR 29 29 3
NO.OF VECTORS CREATED 6999 STARS 3500 RMAX 122.19169 BOHR
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTT gordsh1 TELAPSE 0.05 TCPU 0.02

MATRIX SIZE: P(G) 377103, F(G) 76624, P(G) IRR 193585, F(G) IRR 78287
MAX G-VECTOR INDEX FOR 1- AND 2-ELECTRON INTEGRALS 321

CAPPA:IS1 8;IS2 8;IS3 8; K PTS MONK NET 260; SYMMOPS:K SPACE 2;G SPACE 1

PLANE AVERAGED ELECTROSTATIC POTENTIAL
N. OF POINTS 201 ZMIN -20.000 (AU) - ZMAX 20.000 (AU)

**** LDIM, ICASO= 4 2
ERROR **** POTCAL **** CASE NOT AVAILABLE
Abort(1) on node 2 (rank 2 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 2

Does this method only work for SLAB systems created using CRYSTAL17’s SLABCUT utility or SLAB keyword used in input geometry block? My present work involves a 2D material slab, but I generated the slab structure using VESTA and modeled it in CRYSTAL keyword in input block with a 15 Å vacuum along the z-direction

Dear CRYSTAL developers and users,

I would like to perform a band edge alignment calculation by referencing the valence/conduction band edges to the vacuum level in the context of photocatalysis. For this purpose, I need the planar-averaged electrostatic potential along the slab normal direction.

I have already generated the electrostatic potential in Gaussian cube format using the POTCUBE keyword (file: POT_CUBE.DAT).

My questions are:

Does CRYSTAL provide any direct tool or utility to compute the planar average of the electrostatic potential from POT_CUBE.DAT (or other related output files)?

If not, is it recommended to post-process the cube file externally (e.g., with scripts, VESTA, or custom Python routines)?

Any guidance, example workflows, or scripts would be greatly appreciated.

Best regards,
Chhatra Subba

Thankyou so much. This will help me

Dear Developer,

I have referred to the user manual and tutorial materials, but I am still having difficulty understanding how to correctly define the following sections when generating 2D electron density and electrostatic potential plots using CRYPLOT (Specifically following section)
'''''''''''''
ATOMS
1 1 1 -1
1 0 0 0
1 0 0 2
MARGINS
2. 2. 2. 2.
END
''''''''

OR

'''''''''''
COORDINA
-2.498 0.000 1.696
-2.498 0.000 -1.696
-1.249 -2.164 -1.696
RECTANGU
MARGINS
3. 3. 3. 3.
'''''''''

Could you kindly explain how to determine or generate these inputs properly?

Thank you very much for your assistance.

Best regards,
Chhatra Subba

Thankyou so much. It worked perfectly. This is a big help

Thank you for your response.
The procedure worked perfectly for the 3D structure. However, I'm encountering difficulties when applying it to a 2D system. Specifically, after using the SLABCUT option and optimizing the structure, the same plotting method doesn't seem to work as expected.

Could you kindly guide me on how to correctly generate the 3D charge density for a 2D slab model?

I got this error

DENSITY MATRIX - CONTRIBUTION FROM BANDS
22 23
WARNING **** PDIBAN **** K POINTS WEIGHTS MODIFIED - COMPTON PROFILES NOT ALLOWED
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTT PBAN TELAPSE 0.02 TCPU 0.02
ERROR **** GRID3D **** OPTION NOT ALLOWED
TTTTTTTTTTTTTTTTTTTTTTTTTTTTTT ERR TELAPSE 2.02 TCPU 0.02

Regards
Chhatra Subba

I am a new user of CRYSTAL17. I successfully completed a properties calculation and obtained the file fort.25, which I can visualize using CRYPLOT. However, I would like to plot the charge density in VESTA, and currently I'm not sure how to convert the output to a compatible format.

Could you please guide me on how to:
Convert fort.25 to a Gaussian CUBE file or any other format supported by VESTA?
Specify the correct options using runprop17 (or any other utility) to export the 3D total charge density?

Any help would be greatly appreciated.

Best regards,
Chhatra Subba