Skip to content

CRYSTAL

Discuss features, updates, and general use of the CRYSTAL module

40 Topics 187 Posts

Subcategories


  • SCF, Convergence, Thresholds, Density Functionals, Spin

    7 Topics
    36 Posts

    Hi,

    The MP2 option is no longer supported in recent versions of the CRYSTAl program. If you are interested in a periodic MP2 calculation, my suggestion is to contact Lorenzo Maschio ([email protected]) and Denis Usvyat ([email protected]) directly, who may provide guidance in the use of the CRYSCOR program.

  • Input Format, Pseudopotentials

    1 Topics
    9 Posts

    thank you

  • Internal Coordinates, Constraints, Convergence

    5 Topics
    17 Posts

    Thank you very much. Now it works without any issue. I was not aware of the necessity to change this parameter as I did not encounter any issue with other calculations I have done due to it being 0 0 0.

  • Hessian, Phonons, Quasi-Harmonic Approximation, Anharmonic Force Constants

    3 Topics
    12 Posts

    Hello Dr. Erba,

    Thank you very much for your thorough explanation! This clarifies all my questions. I will make sure to re-calculate using a larger supercell.

    Best,
    Danny

  • Harmonic and Anharmonic Vibrational Spectra, Born Tensor, Raman Activities, Phonon Density-of-States

    10 Topics
    39 Posts

    Sure, I've uploaded the files input files in this link because it exceeded the file size limit to upload them in this comment:
    https://otagouni-my.sharepoint.com/:f:/g/personal/rempe782_student_otago_ac_nz/EmtG3VsEJYJArghuK_IF1x0BvJNBiYsE1mrbQELFWOE1-A?e=NhXhQs

    Used the same inputs (.d12, fort.9, fort.13, fort.20, FREQINFO.DAT) for to do the frequency calculation restarts for CRYSTAL17 and 23

    Thanks heaps

  • Strain, Elastic Tensor, Seismic Waves Velocities

    2 Topics
    5 Posts

    Hi,

    Direct piezoelectric constants of a 3D lattice in CRYSTAL are defined and computed as:
    $$
    e_{ci}^{3D} = \left( \frac{\partial P_c}{\partial \eta_i}\right) = \frac{1}{V}\left( \frac{\partial^2 E}{\partial E_c\partial \eta_i}\right)
    $$
    that is as first derivatives of Cartesian components of the polarization (c=x,y,z) with respect to strain components, or, equivalently as second derivatives of the energy density (V is the volume of the 3D lattice cell) with respect to Cartesian components of an electric field \(E_c\) and strain components, where the strain \( \eta \) is dimensionless and thus the direct piezoelectric constants have units of \( \textup{charge/length}^2 \).

    For 1D and 2D periodic lattices, as the volume (V) is not uniquely defined (or not defined at all in some cases), one may divide by the length \(l \) and area \( A\) of the lattice cell instead:
    $$
    e_{ci}^{1D} = \frac{1}{l}\left( \frac{\partial^2 E}{\partial E_c\partial \eta_i}\right) \quad \textup{and} \quad e_{ci}^{2D} = \frac{1}{A}\left( \frac{\partial^2 E}{\partial E_c\partial \eta_i}\right)
    $$
    that would thus be expressed in units of \( \textup{charge} \) or \( \textup{charge/length} \) for 1D and 2D lattices, respectively.

    However, in CRYSTAL for 1D and 2D lattices we do not divide by \(l \) or \( A\) , and just define and compute the piezoelectric constants as:
    $$
    e_{ci}^\textup{1D and 2D} = \left( \frac{\partial^2 E}{\partial E_c\partial \eta_i}\right)
    $$
    with units of \( \textup{charge}\cdot\textup{length} \).

    Yes, these constants are physically meaningful for 1D and 2D systems. For a 2D monolayer system, for instance, depending on what you need to compare with, you can do one of two things:

    keep them as they are printed in the CRYSTAL output (units of \( \textup{charge}\cdot\textup{length} \))

    divide the values you get in the CRYSTAL output by the area of the 2D cell (and thus express them in units of \( \textup{charge/length} \))

    I would not divide by a volume because I would not know the physical meaning of the volume of a 2D monolayer system.

    Hope this helps,

  • 2 Topics
    24 Posts

    For the purpose of finding the minimum energy structure to then do Raman calculations, it is.

    EOS gives you much more than that of course: the p(V) or, equivalently, V(p) relation (i.e. structure as a function of pressure), the bulk modulus K(p), and allows to compute the enthalpy H(p).

  • Two-Component Spinors, Non-Collinear Magnetisation, Spin and Particle Currents

    4 Topics
    18 Posts

    Thanks Jacques! These are very informative.
    The information for specific basis sets is helpful, I will try them.

  • Input Format, Symmetry, Manipulation, Slabs, Nanotubes, Fullerenes, Helices

    3 Topics
    16 Posts

    Good! The clean way is:

    CRYSTAL
    1 0 0
    P 21/A

    I guess the extra "1"s you put are safe as they correspond to the identity operator.

  • Electric Field, Polarizability, Dielectric Tensor, Hyper-Polarizabilities

    2 Topics
    8 Posts

    Hi,
    Thank you very much for your reply!

    Best regards,
    Masoud

  • Questions that do not fit in other categories

    1 Topics
    3 Posts

    Thank you very much, that worked (but not restart)