Posts made by fbodo

undefined

esmuigors said in Unexpectedly cannot get the geometry after optimization: KEYWORD EXTPRT NOT ALLOWED:

About the manual extraction – do I understand correctly that I should take the primitive and not crystallographic cell parameters?! The input in the .d12 file is in crystallographic cell, isn't it?

Hi esmuigors,
If you want to proceed with the manual extraction of the geometry you should copy the Crystallographic cell lattice parameters required by the space group you are working with (e.g. Monoclinic will require you a, b, c and β, while a cubic only the a lattice vector), the only exception to this are P lattices where the primitive and crystallographic coincide. Accordingly, also the fractionary coordinate you will copy should be the one of the Crystallographic cell, but you should copy only the atoms present in the asymmetric unit of your system identified by a T in the FINAL OPTIMIZED GEOMETRY print in the .out file, I will include here below an example to better clarify.

 FINAL OPTIMIZED GEOMETRY - DIMENSIONALITY OF THE SYSTEM      3
 (NON PERIODIC DIRECTION: LATTICE PARAMETER FORMALLY SET TO 500)
 *******************************************************************************
 LATTICE PARAMETERS (ANGSTROMS AND DEGREES) - BOHR = 0.5291772083 ANGSTROM
 PRIMITIVE CELL - CENTRING CODE 6/0 VOLUME=    62.890767 - DENSITY 13.512 g/cm^3
         A              B              C           ALPHA      BETA       GAMMA 
     6.02167869     6.02167869     6.02167869   147.475653 147.475653  46.660012
 *******************************************************************************
 ATOMS IN THE ASYMMETRIC UNIT    2 - ATOMS IN THE UNIT CELL:    4
     ATOM                 X/A                 Y/B                 Z/C    
 *******************************************************************************
      1 T 273 TA   -2.420755613827E-03 -2.420755613827E-03  4.983867094700E-20
      2 F 273 TA   -2.524207556138E-01  2.475792443862E-01 -5.000000000000E-01
      3 T 233 AS    4.184207556138E-01  4.184207556138E-01 -2.220446049250E-16
      4 F 233 AS    1.684207556138E-01 -3.315792443862E-01 -5.000000000000E-01

 TRANSFORMATION MATRIX PRIMITIVE-CRYSTALLOGRAPHIC CELL
  0.0000  1.0000  1.0000  1.0000  0.0000  1.0000  1.0000  1.0000  0.0000

 *******************************************************************************
 CRYSTALLOGRAPHIC CELL (VOLUME=        125.78153367)
         A              B              C           ALPHA      BETA       GAMMA 
     3.37253732     3.37253732    11.05868170    90.000000  90.000000  90.000000

 COORDINATES IN THE CRYSTALLOGRAPHIC CELL
     ATOM                 X/A                 Y/B                 Z/C    
 *******************************************************************************
      1 T 273 TA   -1.956455588613E-19  2.296201495291E-19 -2.420755613827E-03
      2 F 273 TA   -5.000000000000E-01  2.084229456032E-17 -2.524207556138E-01
      3 T 233 AS    5.000000000000E-01  5.000000000000E-01 -8.157924438617E-02
      4 F 233 AS   -5.117515658395E-17 -5.000000000000E-01 -3.315792443862E-01

In this example you'll have to extract from the second geometry print (i.e. Crystallographic cell) the a and c lattice parameters, given the Tetragonal lattice, alongside the fractionary coordinate of atoms 1 and 3.

esmuigors said in Unexpectedly cannot get the geometry after optimization: KEYWORD EXTPRT NOT ALLOWED:

Are there any workarounds to persuade the script to actually copy the fort.34 file?

In this regards, I don't have your script, but it should be sufficient to add to your script something along these lines in the section where all the files are copied back to your working folder:

if [ -e fort.34 ] then
    cp fort.34 $HERE/$INPUTFILE.f34
fi

in this case $HERE=$PWD, while $INPUTFILE is the variable corresponding to your input name.

I hope this helps.

undefined

Hi job314,
I have tried to run the input file and geometry you provided us as you can see in the enclosed output file it runned without any hickups for us. Please let us know if you were to encounter this issue once more

undefined

Hi job314,
I gave a quick look at the files you have uploaded, and I don't seem to find the RESTART option in the INPUT file. Could you please check if this is, indeed, the input associated with your issue?

undefined

CRYSTALClear is an open source project that provides an easy Python interface with CRYSTAL. The package allows you to quickly extract information from the CRYSTAL output files and to easily generate customizable plots of computed quantities. In particular, the package supports plotting functionalities for:

Band Structures and DOSS	Mechanical Properties	IR and Raman Spectra

QTAIM	Transport Properties	Density Maps

For a more comprehensive review of the package capabilities we refer you to the Documentation Website.

How to Install

Prerequisites:

Python: Ensure you have Python installed on your system.
conda (Optional but Recommended): conda is a powerful package and environment manager that simplifies the process of managing Python dependencies. If you don't have conda, install it from the official Anaconda distribution: https://www.anaconda.com/

Installation Steps:

Create a Dedicated Environment (Recommended):
Creating a separate conda environment for CRYSTALClear helps isolate its dependencies from other projects and avoid potential conflicts.

conda create -n env_name python=3.9  # Replace 3.9 with your desired Python version

Activate the Environment (Reccommended):

conda activate env_name

Install CRYSTALClear using pip:

pip install CRYSTALClear

This command will install the CRYSTALClear library and its dependencies.

Verification:

Check Installation:
After installation, you can verify the installation by importing the library in a Python script or interactive session:

import CRYSTALClear

If the import is successful without any errors, the installation was successful.

Additional Notes:

Dependencies: The CRYSTALClear library might have dependencies on other libraries. pip will automatically install these dependencies during the installation process.
Updating: To update the CRYSTALClear library to the latest version:

pip install --upgrade CRYSTALClear

Troubleshooting:

If you encounter any issues during installation:
- Check for any error messages and try to resolve them based on the error descriptions.
- Ensure you have an active internet connection.
- Try updating pip to the latest version: pip install --upgrade pip
- If the issue persists, consider creating a new conda environment and reinstalling the library.

If you the issues persist please feel free to comment here below we will try our best to help you.

How to Use it

In order to use any of the feature available in the package the user need to generate
a CRYSTAL object through the crystal_io module. Three classes are available to the user
for different use cases:

Crystal_output: To extract informations and generate the object required to plot
from any .out file of the CRYSTAL package.
Properties_output: To generate the objects required to plot any of quantities
stored in the .DAT and .f25 files generated by the PROPERTIES module.
External_unit: To generate the objects required to plot any of the quantites stored
in .DAT and .f25 file generated by CRYSTAL.

Once the object is generated the user can generate the desired plot using the functions
in the plot module

Example

   # Example Band Structure plot
   from CRYSTALClear.crystal_io import Properties_output
   import CRYSTALClear.plot as CCplt
   import matplotlib.pyplot as plt

   data = Properties_output().read_electron_band()
   CCplt.plot_electron_band(data, **kwargs) # see the documentation for the available **kwargs 

   # The plotting function returns a matplotlib object that can be visualized as follows
   plt.show()

The plotting function will return one (or a list of) matplotlib object(s) that the user can easily modify as any other plot produced with the library.

For a more comprehensive set of example we refer you to the following jupyter notebook repository.

How to Contribute

Any contribution to the project is more than welcome and greatly appreciated, given the open source nature of the project. If you'd like to do so, please check our guidelines on GitHub

How to Cite

If you use this package in any of your work we kindly ask you to cite the following pubblication:

Camino, Bruno, Huanyu Zhou, Eleonora Ascrizzi, Alberto Boccuni, Filippo Bodo, Alessandro Cossard, Davide Mitoli, Anna Maria Ferrari, Alessandro Erba, and Nicholas M. Harrison, Comput. Phys. Commun. 292, 108853 (2023).

@article{camino2023crystalpytools,
  title={CRYSTALpytools: A Python infrastructure for the CRYSTAL code},
  author={Camino, Bruno and Zhou, Huanyu and Ascrizzi, Eleonora and Boccuni, Alberto and Bodo, Filippo and Cossard, Alessandro and Mitoli, Davide and Ferrari, Anna Maria and Erba, Alessandro and Harrison, Nicholas M},
  journal={Computer Physics Communications},
  volume={292},
  pages={108853},
  year={2023},
  publisher={Elsevier}
}