.. 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_15_RHL2.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_15_RHL2.py: RHL2 **** Rhombohedral cell is defined by two parameters :math:`a` and :math:`\alpha`. RHL lattice has variation RHL2 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 = 110^{\circ}`. .. GENERATED FROM PYTHON SOURCE LINES 37-65 .. code-block:: Python import wulfric cell = wulfric.cell.sc_get_example("RHL2") 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 == "RHL2" print(variation) .. rst-class:: sphx-glr-script-out .. code-block:: none RHL2 .. 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-P-Z-Q-GAMMA-F-P1-Q1-L-Z .. 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 F 0.0000 0.5000 -0.5000 0.8717 0.6104 -0.5982 L 0.0000 0.5000 -0.0000 0.8717 0.6104 0.6475 P 0.3726 0.3726 -0.6274 1.2991 -0.0000 -0.5982 P1 -0.3726 0.6274 -0.3726 0.4443 1.2208 -0.5982 Q 0.2451 0.2451 0.2451 0.8548 -0.0000 1.2457 Q1 -0.2451 0.7549 -0.2451 0.8887 1.2208 0.0493 Z 0.5000 0.5000 -0.5000 1.7434 -0.0000 0.0493 .. 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.956 seconds) .. _sphx_glr_download_user-guide_conventions_bravais-lattices_2_sc_plot_15_RHL2.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_15_RHL2.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_15_RHL2.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_15_RHL2.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_