.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "user-guide/conventions/bravais-lattices/2_sc/plot_14_RHL1.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_user-guide_conventions_bravais-lattices_2_sc_plot_14_RHL1.py: RHL1 **** Rhombohedral cell is defined by two parameters :math:`a` and :math:`\alpha`. RHL lattice has variation RHL1 when :math:`\alpha < 90^{\circ}`. Cell constructor ================ To get an example of the cell use :py:func:`wulfric.cell.SC_RHL`. :py:func:`wulfric.cell.sc_get_example` returns an example where :math:`a = \pi` and :math:`\alpha = 70^{\circ}`. .. GENERATED FROM PYTHON SOURCE LINES 37-65 .. code-block:: Python import wulfric cell = wulfric.cell.sc_get_example("RHL1") atoms = dict(positions=[[0, 0, 0]], spglib_types=[1]) # To avoid multiple calls to spglib one can do it once and then pass spglib_data # to the functions where it is needed spglib_data = wulfric.get_spglib_data(cell=cell, atoms=atoms) kp = wulfric.Kpoints.from_crystal(cell=cell, atoms=atoms, convention="SC") conv_cell, conv_atoms = wulfric.crystal.get_conventional( cell=cell, atoms=atoms, convention="SC", spglib_data=spglib_data ) prim_cell, prim_atoms = wulfric.crystal.get_primitive( cell=cell, atoms=atoms, convention="SC", spglib_data=spglib_data ) variation = wulfric.crystal.sc_get_variation( cell=cell, atoms=atoms, spglib_data=spglib_data ) assert variation == "RHL1" print(variation) .. rst-class:: sphx-glr-script-out .. code-block:: none RHL1 .. GENERATED FROM PYTHON SOURCE LINES 66-68 K-path ====== .. GENERATED FROM PYTHON SOURCE LINES 68-71 .. code-block:: Python print(kp.path_string) .. rst-class:: sphx-glr-script-out .. code-block:: none GAMMA-L-B1|B-Z-GAMMA-X|Q-F-P1-Z|L-P .. GENERATED FROM PYTHON SOURCE LINES 72-74 High-symmetry points ==================== .. GENERATED FROM PYTHON SOURCE LINES 74-77 .. code-block:: Python print(kp.hs_table(decimals=4)) .. rst-class:: sphx-glr-script-out .. code-block:: none Name rel_b1 rel_b2 rel_b3 k_x k_y k_z GAMMA 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 B 0.7031 0.5000 0.2969 1.4687 -0.3541 -0.0214 B1 0.5000 0.2969 -0.2969 0.9728 -0.3541 -1.1005 F 0.5000 0.5000 -0.0000 1.2208 0.0000 -0.5609 L 0.5000 -0.0000 -0.0000 0.6104 -0.8717 -0.2805 L1 -0.0000 -0.0000 -0.5000 -0.0000 -0.0000 -1.1005 P 0.7031 0.3985 0.3985 1.3447 -0.5311 0.2591 P1 0.6015 0.6015 0.2969 1.4687 0.0000 -0.0214 P2 0.3985 0.3985 -0.2969 0.9728 0.0000 -1.1005 Q 0.6015 0.3985 -0.0000 1.2208 -0.3541 -0.5609 X 0.3985 -0.0000 -0.3985 0.4864 -0.6947 -1.1005 Z 0.5000 0.5000 0.5000 1.2208 0.0000 0.5396 .. GENERATED FROM PYTHON SOURCE LINES 78-80 Brillouin zone and default k-path ================================= .. GENERATED FROM PYTHON SOURCE LINES 80-95 .. code-block:: Python pe = wulfric.PlotlyEngine(_sphinx_gallery_fix=True) pe.plot_brillouin_zone( cell=prim_cell, color="red", legend_label="Brillouin zone of the primitive cell" ) pe.plot_brillouin_zone( cell=cell, color="chocolate", legend_label="Brillouin zone of the original cell" ) pe.plot_kpath(kp=kp) pe.plot_kpoints(kp=kp, only_from_kpath=True) pe.show(axes_visible=False) .. raw:: html

.. GENERATED FROM PYTHON SOURCE LINES 96-101 Cells of real space =================== .. hint Click on the legend to hide some of the cells .. GENERATED FROM PYTHON SOURCE LINES 101-113 .. code-block:: Python pe = wulfric.PlotlyEngine(_sphinx_gallery_fix=True) pe.plot_cell(cell=cell, legend_label="Original cell", color="Chocolate") pe.plot_cell(cell=prim_cell, legend_label="Primitive cell", color="Black") pe.plot_cell(cell=conv_cell, legend_label="Conventional cell", color="Blue") pe.plot_wigner_seitz_cell( cell=prim_cell, legend_label="Wigner-Seitz cell", color="green" ) pe.show(axes_visible=False) .. raw:: html

.. GENERATED FROM PYTHON SOURCE LINES 114-128 Edge cases ========== In rhombohedral lattice :math:`a = b = c` and :math:`\alpha = \beta = \gamma`, thus three edge cases exist: If :math:`\alpha = 60^{\circ}`, then the lattice is :ref:`sphx_glr_user-guide_conventions_bravais-lattices_2_sc_plot_02_FCC.py`. If :math:`\alpha \approx 109.47122063^{\circ}` (:math:`\cos(\alpha) = -1/3`), then the lattice is :ref:`sphx_glr_user-guide_conventions_bravais-lattices_2_sc_plot_03_BCC.py`. If :math:`\alpha = 90^{\circ}`, then the lattice is :ref:`sphx_glr_user-guide_conventions_bravais-lattices_2_sc_plot_01_CUB.py`. .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 1.605 seconds) .. _sphx_glr_download_user-guide_conventions_bravais-lattices_2_sc_plot_14_RHL1.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_14_RHL1.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_14_RHL1.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_14_RHL1.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_