TB2J 0.9.9.13__py3-none-any.whl → 0.9.9.15__py3-none-any.whl

TB2J/MAEGreen.py

@@ -94,12 +94,11 @@ class MAEGreen(ExchangeNCL):
                 self.thetas.append(i * np.pi / 180)
                 self.phis.append(j * np.pi / 180)
 
-    def set_angels_ztox(self, n=16):  
+    def set_angels_ztox(self, n=16):
         """Set angles for a scan from z to x"""
         self.thetas = np.linspace(0, np.pi, n)
         self.phis = np.zeros(n)
 
-
     def set_angles_random(self, n=16):
         # n random pairs of theta, phi
         self.thetas = np.random.random(n) * np.pi

TB2J/__init__.py

@@ -1 +1 @@
-__version__ = "0.9.9.12"
+__version__ = "0.9.9.14"

TB2J/exchange.py

@@ -48,7 +48,7 @@ class Exchange(ExchangeParams):
         self.emin = self.G.find_energy_ingap(rbound=self.efermi - 15.0) - self.efermi
         # self.emin = self.G.find_energy_ingap(rbound=self.efermi - 15.0) - self.efermi
         # self.emin = -42.0
-        #print(f"A gap is found at {self.emin}, set emin to it.")
+        # print(f"A gap is found at {self.emin}, set emin to it.")
 
     def set_tbmodels(self, tbmodels):
         pass

TB2J/interfaces/abacus/stru_api.py

@@ -7,21 +7,28 @@
 Modified on Wed Aug 01 11:44:51 2022
 @author: Ji Yu-yang
 """
-import warnings
-warnings.simplefilter(action='ignore', category=FutureWarning)
-
 
 import os
 import re
 import shutil
+import warnings
+from copy import deepcopy
+from math import sqrt
 from pathlib import Path
 
 import numpy as np
 from ase import Atoms
+from ase.calculators.calculator import kpts2ndarray, kpts2sizeandoffsets
 from ase.calculators.singlepoint import SinglePointDFTCalculator, arrays_to_kpoints
+from ase.constraints import FixAtoms, FixCartesian
+from ase.dft.kpoints import bandpath
+from ase.stress import full_3x3_to_voigt_6_stress
 from ase.units import Bohr, GPa, Hartree, Rydberg, _me, mol
 from ase.utils import lazymethod, lazyproperty, reader, writer
 
+warnings.simplefilter(action="ignore", category=FutureWarning)
+
+
 _re_float = r"[-+]?\d+\.*\d*(?:[Ee][-+]\d+)?"
 AU_to_MASS = mol * _me * 1e3
 UNIT_V = np.sqrt(Hartree / AU_to_MASS)
@@ -42,7 +49,6 @@ def write_input(fd, parameters=None):
     parameters: dict
         The dictionary of all paramters for the calculation.
     """
-    from copy import deepcopy
 
     params = deepcopy(parameters)
     params["dft_functional"] = (
@@ -105,8 +111,6 @@ def write_kpt(fd=None, parameters=None, atoms=None):
         return
     elif kpts is not None:
         if isinstance(kpts, dict) and "path" not in kpts:
-            from ase.calculators.calculator import kpts2sizeandoffsets
-
             kgrid, shift = kpts2sizeandoffsets(atoms=atoms, **kpts)
             koffset = []
             for i, x in enumerate(shift):
@@ -121,8 +125,6 @@ def write_kpt(fd=None, parameters=None, atoms=None):
         koffset = [koffset] * 3
 
     if isinstance(kgrid, dict) or hasattr(kgrid, "kpts"):
-        from ase.calculators.calculator import kpts2ndarray
-
         kmode = "Direct"
         kgrid = kpts2ndarray(kgrid, atoms=atoms)
     elif isinstance(kgrid, str) and (kgrid == "gamma"):
@@ -272,8 +274,6 @@ def judge_exist_stru(stru=None):
 
 
 def read_ase_stru(stru=None, coordinates_type="Cartesian"):
-    from ase.constraints import FixAtoms, FixCartesian
-
     fix_cart = np.ones((len(stru), 3), dtype=int).tolist()
     for constr in stru.constraints:
         for i in constr.index:
@@ -567,8 +567,6 @@ def read_kpt(fd, cell=None):
         else:
             knumbers = klines[:, 3].astype(int)
             if cell is not None:
-                from ase.dft.kpoints import bandpath
-
                 return bandpath(kpts, cell, npoints=knumbers.sum())
             else:
                 return {
@@ -683,8 +681,6 @@ def read_abacus(fd, latname=None, verbose=False):
     If `verbose` is True, pseudo-potential, basis and other information along with the Atoms object will be output as a dict.
     """
 
-    from ase.constraints import FixCartesian
-
     contents = fd.read()
     title_str = r"(?:LATTICE_CONSTANT|NUMERICAL_DESCRIPTOR|NUMERICAL_ORBITAL|ABFS_ORBITAL|LATTICE_VECTORS|LATTICE_PARAMETERS|ATOMIC_POSITIONS)"
 
@@ -887,8 +883,6 @@ def read_abacus(fd, latname=None, verbose=False):
 
 
 def get_lattice_from_latname(lines, latname=None):
-    from math import sqrt
-
     if lines:
         lines = lines.group(1).split(" ")
 
@@ -1488,7 +1482,6 @@ class AbacusOutCalcChunk(AbacusOutChunk):
     @lazymethod
     def get_stress(self):
         """Get the stress from the output file according to index"""
-        from ase.stress import full_3x3_to_voigt_6_stress
 
         try:
             stress = (

TB2J/interfaces/siesta_interface.py

@@ -165,7 +165,7 @@ Warning: The DMI component parallel to the spin orientation, the Jani which has
             # thetas = [0, np.pi / 2, np.pi, 3 * np.pi / 2]
             # phis = [0, 0, 0, 0]
             # MAE.set_angles(thetas=thetas, phis=phis)
-            #MAE.set_xyz_angles()
+            # MAE.set_xyz_angles()
             MAE.run(output_path=f"{output_path}_anisotropy", with_eigen=False)
             # print(
             #    f"MAE calculation finished. The results are in {output_path} directory."

TB2J/io_exchange/io_exchange.py

@@ -340,20 +340,19 @@ Generation time: {now.strftime("%y/%m/%d %H:%M:%S")}
         """
         i = self.i_spin(i)
         j = self.i_spin(j)
-        J=D=Ja = None
+        J = D = Ja = None
         if Jiso:
-            J= self.get_Jiso(i, j, R)
+            J = self.get_Jiso(i, j, R)
         if DMI:
             D = self.get_DMI(i, j, R)
             if D is not None:
-                D*=1
+                D *= 1
         if Jani:
             Ja = self.get_Jani(i, j, R)
             if Ja is not None:
-                Ja*=1
+                Ja *= 1
         Jtensor = combine_J_tensor(Jiso=J, D=D, Jani=Ja)
 
-
         # if iso_only:
         #    J = self.get_Jiso(i, j, R)
         #    if J is not None:

TB2J/io_exchange/io_vampire.py

@@ -60,7 +60,7 @@ def write_vampire_unitcell_file(cls, fname):
                     val = np.real(Jtensor[i, j] * 2.0 / J)
                     if np.abs(val) < 1e-30:
                         val = 0.0
-                    myfile.write(f"{val:<012.5e} ")
+                    myfile.write(f"{val:< 12.5e} ")
             myfile.write("\n")
 
 

TB2J/io_merge.py

@@ -1,31 +1,35 @@
-import os
 import copy
+import os
 import warnings
+from itertools import combinations_with_replacement
+
 import numpy as np
-from itertools import combinations_with_replacement, product
+
 from TB2J.io_exchange import SpinIO
 
 u0 = np.zeros(3)
 uy = np.array([0.0, 1.0, 0.0])
 uz = np.array([0.0, 0.0, 1.0])
 
-def get_Jani_coefficients(a, R=np.eye(3)):
 
+def get_Jani_coefficients(a, R=np.eye(3)):
     if len(a) == 1:
         u = a
         v = a
     else:
         u = a[[0, 0, 1]]
         v = a[[0, 1, 1]]
-    
+
     ur = u @ R.T
     vr = v @ R.T
-    coefficients = np.hstack([ur*vr, np.roll(ur, -1, axis=-1)*vr + np.roll(vr, -1, axis=-1)*ur])
+    coefficients = np.hstack(
+        [ur * vr, np.roll(ur, -1, axis=-1) * vr + np.roll(vr, -1, axis=-1) * ur]
+    )
 
     return coefficients, u, v
 
-def get_projections(a, b, tol=1e-2):
 
+def get_projections(a, b, tol=1e-2):
     projections = np.empty((2, 3))
     if np.linalg.matrix_rank([a, b], tol=tol) == 1:
         if np.linalg.matrix_rank([a, uy], tol=tol) == 1:
@@ -41,13 +45,13 @@ def get_projections(a, b, tol=1e-2):
 
     return projections
 
+
 class SpinIO_merge(SpinIO):
     def __init__(self, *args, **kwargs):
         super(SpinIO_merge, self).__init__(*args, **kwargs)
         self.projv = None
 
     def _set_projection_vectors(self):
-
         norm = np.linalg.norm(self.spinat, axis=-1).reshape(-1, 1)
         spinat = self.spinat / norm
         idx = [self.ind_atoms[i] for i in self.index_spin if i >= 0]
@@ -60,27 +64,29 @@ class SpinIO_merge(SpinIO):
         self.projv = projv
 
     @classmethod
-    def load_pickle(cls, path='TB2J_results', fname='TB2J.pickle'):
+    def load_pickle(cls, path="TB2J_results", fname="TB2J.pickle"):
         obj = super(SpinIO_merge, cls).load_pickle(path=path, fname=fname)
         obj._set_projection_vectors()
 
         return obj
 
+
 def read_pickle(path):
-    p1 = os.path.join(path, 'TB2J_results', 'TB2J.pickle')
-    p2 = os.path.join(path, 'TB2J.pickle')
+    p1 = os.path.join(path, "TB2J_results", "TB2J.pickle")
+    p2 = os.path.join(path, "TB2J.pickle")
     if os.path.exists(p1) and os.path.exists(p2):
         print(f" WARNING!: Both file {p1} and {p2} exist. Use default {p1}.")
     if os.path.exists(p1):
-        ret = SpinIO_merge.load_pickle(os.path.join(path, 'TB2J_results'))
+        ret = SpinIO_merge.load_pickle(os.path.join(path, "TB2J_results"))
     elif os.path.exists(p2):
         ret = SpinIO_merge.load_pickle(path)
     else:
         raise FileNotFoundError(f"Cannot find either file {p1} or {p2}")
     return ret
 
-class Merger():
-    def __init__(self, *paths, main_path=None): 
+
+class Merger:
+    def __init__(self, *paths, main_path=None):
         self.dat = [read_pickle(path) for path in paths]
 
         if main_path is None:
@@ -90,32 +96,34 @@ class Merger():
             self.dat.append(copy.deepcopy(self.main_dat))
 
         self._set_projv()
-        
+
     def _set_projv(self):
         cell = self.main_dat.atoms.cell.array
-        rotated_cells = np.stack(
-            [obj.atoms.cell.array for obj in self.dat], axis=0
-        )
+        rotated_cells = np.stack([obj.atoms.cell.array for obj in self.dat], axis=0)
         R = np.linalg.solve(cell, rotated_cells)
         indices = range(len(self.dat))
 
-        proju = {}; projv = {}; coeff_matrix = {}; projectors = {};
+        proju = {}
+        projv = {}
+        coeff_matrix = {}
+        projectors = {}
         for key in self.main_dat.projv.keys():
             vectors = [obj.projv[key] for obj in self.dat]
-            coefficients, u, v = zip(*[get_Jani_coefficients(vectors[i], R=R[i]) for i in indices])
+            coefficients, u, v = zip(
+                *[get_Jani_coefficients(vectors[i], R=R[i]) for i in indices]
+            )
             projectors[key] = np.vstack([u[i] @ R[i].T for i in indices])
             coeff_matrix[key] = np.vstack(coefficients)
             proju[key] = np.stack(u)
             projv[key] = np.stack(v)
             if np.linalg.matrix_rank(coeff_matrix[key], tol=1e-2) < 6:
-                warnings.warn('''
+                warnings.warn("""
                     WARNING: The matrix of equations to reconstruct the exchange tensors is
                     close to being singular. This happens when the magnetic moments between
                     different configurations are cloes to being parallel. You need to consider 
                     more rotated spin configurations, otherwise the results might have a large 
-                    error.'''
-                )
-        
+                    error.""")
+
         self.proju = proju
         self.projv = projv
         self.coeff_matrix = coeff_matrix
@@ -123,71 +131,78 @@ class Merger():
 
     def merge_Jani(self):
         Jani_dict = {}
-        proju = self.proju; projv = self.projv; coeff_matrix = self.coeff_matrix;
+        proju = self.proju
+        projv = self.projv
+        coeff_matrix = self.coeff_matrix
         for key in self.main_dat.Jani_dict.keys():
             try:
                 R, i, j = key
                 u = proju[i, j]
                 v = projv[i, j]
                 Jani = np.stack([sio.Jani_dict[key] for sio in self.dat])
-                projections = np.einsum('nmi,nij,nmj->nm', u, Jani, v).flatten()
+                projections = np.einsum("nmi,nij,nmj->nm", u, Jani, v).flatten()
             except KeyError as err:
                 raise KeyError(
                     "Can not find key: %s, Please make sure the three calculations use the same k-mesh and same Rcut."
-                    % err)
+                    % err
+                )
             newJani = np.linalg.lstsq(coeff_matrix[i, j], projections, rcond=1e-2)[0]
-            Jani_dict[key] = np.array([
-                [newJani[0], newJani[3], newJani[5]],
-                [newJani[3], newJani[1], newJani[4]],
-                [newJani[5], newJani[4], newJani[2]]
-            ])
+            Jani_dict[key] = np.array(
+                [
+                    [newJani[0], newJani[3], newJani[5]],
+                    [newJani[3], newJani[1], newJani[4]],
+                    [newJani[5], newJani[4], newJani[2]],
+                ]
+            )
         self.main_dat.Jani_dict = Jani_dict
 
     def merge_Jiso(self):
-        Jdict={}
+        Jdict = {}
         for key in self.main_dat.exchange_Jdict.keys():
             try:
-               J = np.mean([obj.exchange_Jdict[key] for obj in self.dat])
+                J = np.mean([obj.exchange_Jdict[key] for obj in self.dat])
             except KeyError as err:
                 raise KeyError(
-                        "Can not find key: %s, Please make sure the three calculations use the same k-mesh and same Rcut."
-                        % err)
+                    "Can not find key: %s, Please make sure the three calculations use the same k-mesh and same Rcut."
+                    % err
+                )
             Jdict[key] = J
         self.main_dat.exchange_Jdict = Jdict
-               
 
     def merge_DMI(self):
         dmi_ddict = {}
         if all(obj.has_dmi for obj in self.dat):
-            projectors = self.projectors; proju = self.proju;
+            projectors = self.projectors
+            proju = self.proju
             for key in self.main_dat.dmi_ddict.keys():
                 try:
                     R, i, j = key
                     u = proju[i, j]
                     DMI = np.stack([sio.dmi_ddict[key] for sio in self.dat])
-                    projections = np.einsum('nmi,ni->nm', u, DMI).flatten()
+                    projections = np.einsum("nmi,ni->nm", u, DMI).flatten()
                 except KeyError as err:
                     raise KeyError(
                         "Can not find key: %s, Please make sure the three calculations use the same k-mesh and same Rcut."
-                        % err)
+                        % err
+                    )
                 newDMI = np.linalg.lstsq(projectors[i, j], projections, rcond=4e-1)[0]
                 dmi_ddict[key] = newDMI
             self.main_dat.dmi_ddict = dmi_ddict
 
     def standardize(self):
         # make sure that the Jani has the trace of zero
-        Jdict=self.main_dat.exchange_Jdict
-        Jani_dict=self.main_dat.Jani_dict
+        Jdict = self.main_dat.exchange_Jdict
+        Jani_dict = self.main_dat.Jani_dict
         for key in self.main_dat.Jani_dict.keys():
             Jani = self.main_dat.Jani_dict[key]
-            shift = np.trace(Jani)/3.0
-            Jani_dict[key]  -= shift * np.eye(3)
+            shift = np.trace(Jani) / 3.0
+            Jani_dict[key] -= shift * np.eye(3)
             Jdict[key] += shift
         self.main_dat.Jani_dict = Jani_dict
         self.main_dat.exchange_Jdict = Jdict
-        
 
-def merge(*paths, main_path=None, save=True, write_path='TB2J_results'):
+
+def merge(*paths, main_path=None, save=True, write_path="TB2J_results"):
     m = Merger(*paths, main_path=main_path)
     m.merge_Jiso()
     m.merge_DMI()

TB2J/magnon/magnon3.py

@@ -2,12 +2,16 @@ from dataclasses import asdict, dataclass
 from pathlib import Path
 from typing import List, Optional, Tuple, Union
 
+import matplotlib.pyplot as plt
 import numpy as np
 import tomli
 import tomli_w
+from ase.dft.dos import DOS
+from ase.units import J, eV
 from scipy.spatial.transform import Rotation
 
 from TB2J.io_exchange import SpinIO
+from TB2J.kpoints import monkhorst_pack
 from TB2J.magnon.magnon_band import MagnonBand
 from TB2J.magnon.magnon_math import get_rotation_arrays
 from TB2J.mathutils.auto_kpath import auto_kpath
@@ -162,10 +166,10 @@ class Magnon:
                 phase = 2 * np.pi * R @ qpt
                 Jq[iqpt] += np.exp(1j * phase) * JRprime[iR]
 
-        #Jq_copy = Jq.copy()
-        #Jq.swapaxes(-1, -2)  # swap xyz
-        #Jq.swapaxes(-3, -4)  # swap ij
-        #Jq = (Jq.conj() + Jq_copy) / 2.0
+        # Jq_copy = Jq.copy()
+        # Jq.swapaxes(-1, -2)  # swap xyz
+        # Jq.swapaxes(-3, -4)  # swap ij
+        # Jq = (Jq.conj() + Jq_copy) / 2.0
         return Jq
 
     def Hq(self, kpoints):
@@ -187,7 +191,6 @@ class Magnon:
         magmoms = self.magmom.copy()
         magmoms /= np.linalg.norm(magmoms, axis=-1)[:, None]
 
-        
         U, V = get_rotation_arrays(magmoms, u=self._uz)
 
         J0 = -self.Jq(np.zeros((1, 3)))[0]
@@ -195,7 +198,6 @@ class Magnon:
         # Jq = -Hermitize(self.Jq(kpoints, anisotropic=anisotropic))
 
         Jq = -self.Jq(kpoints)
-        print(f"J0 shape: {J0.shape}")
 
         C = np.diag(np.einsum("ix,ijxy,jy->i", V, 2 * J0, V))
         B = np.einsum("ix,kijxy,jy->kij", U, Jq, U)
@@ -203,7 +205,6 @@ class Magnon:
         A2 = np.einsum("ix,kijxy,jy->kij", U.conj(), Jq, U)
 
         H = np.block([[A1 - C, B], [B.swapaxes(-1, -2).conj(), A2 - C]])
-        print(f"H shape: {H.shape}")
         return H
 
     def _magnon_energies(self, kpoints, u=None):
@@ -801,12 +802,11 @@ def main():
     )
 
     parser.add_argument(
-        "--show", 
+        "--show",
         action="store_true",
         default=False,
         help="show figure on screen.",
-        )
-
+    )
 
     args = parser.parse_args()
 
@@ -834,7 +834,7 @@ def main():
             Q=args.Q if args.Q is not None else None,
             uz_file=args.uz_file,
             n=args.n if args.n is not None else None,
-            show=args.show
+            show=args.show,
         )
 
     plot_magnon_bands_from_TB2J(params)
@@ -842,3 +842,97 @@ def main():
 
 if __name__ == "__main__":
     main()
+
+
+class MagnonASEWrapper:
+    def __init__(self, magnon: Magnon):
+        self.magnon = magnon
+        self.atoms = None
+        self.kpts = None
+        self.weights = None
+        self.dos_args = {}
+
+    def set(self, atoms=None, kmesh=[9, 9, 9], gamma=True, **kwargs):
+        self.atoms = atoms
+        self.dos_args = {
+            "kmesh": kmesh,
+            "gamma": gamma,
+        }
+        self.kpts = monkhorst_pack(
+            self.dos_args["kmesh"], gamma_center=self.dos_args["gamma"]
+        )
+        self.weights = np.ones(len(self.kpts)) / len(self.kpts)
+
+    def get_k_points_and_weights(self):
+        return self.kpts, self.weights
+
+    def get_k_point_weights(self):
+        return self.weights
+
+    def get_number_of_spins(self):
+        return 1
+
+    def get_eigenvalues(self, kpt, spin=0):
+        """
+        return the eigenvalues at a given k-point. The energy unit is eV
+        args:
+            kpt: k-point index.
+            spin: spin index.
+        """
+        kpoint = self.kpts[kpt]
+        # Magnon energies are already in eV, convert to meV for consistency with plot
+        evals = self.magnon._magnon_energies(np.array([kpoint]))[0]
+        evals = evals * J / eV  # Convert to eV
+        return evals
+
+    def get_fermi_level(self):
+        return 0.0
+
+    def get_bz_k_points(self):
+        return self.kpts
+
+    def get_dos(self, width=0.1, window=None, npts=401):
+        dos = DOS(self, width=width, window=window, npts=npts)
+        energies = dos.get_energies()
+        tdos = dos.get_dos()
+        return energies, tdos
+
+    def plot_dos(
+        self,
+        smearing_width=0.0001,
+        window=None,
+        npts=401,
+        output="magnon_dos.pdf",
+        ax=None,
+        show=True,
+        dos_filename="magnon_dos.txt",
+    ):
+        """
+        plot total DOS.
+        :param width: width of Gaussian smearing
+        :param window: energy window
+        :param npts: number of points
+        :param output: output filename
+        :param ax: matplotlib axis
+        :return: ax
+        """
+        if ax is None:
+            _fig, ax = plt.subplots()
+        energies, tdos = self.get_dos(width=smearing_width, window=window, npts=npts)
+        energies = energies * 1000  # Convert to meV
+        tdos = tdos / 1000  # Convert to states/meV
+        if dos_filename is not None:
+            np.savetxt(
+                dos_filename,
+                np.array([energies, tdos]).T,
+                header="Energy(meV) DOS(state/meV)",
+            )
+        ax.plot(energies, tdos)
+        ax.set_xlabel("Energy (meV)")
+        ax.set_ylabel("DOS (states/meV)")
+        ax.set_title("Total DOS")
+        if output is not None:
+            plt.savefig(output)
+        if show:
+            plt.show()
+        return ax

TB2J/magnon/magnon_band.py

@@ -40,7 +40,7 @@ class MagnonBand:
         if self.xcoords is None:
             self.xcoords = np.arange(len(self.kpoints))
 
-    def plot(self, ax=None, filename=None, show=False, shift=0.0,**kwargs):
+    def plot(self, ax=None, filename=None, show=False, shift=0.0, **kwargs):
         """Plot the magnon band structure.
 
         Parameters
@@ -79,7 +79,7 @@ class MagnonBand:
                 for band in segment_bands:
                     ax.plot(
                         x,
-                        band[start_idx : start_idx + nbands]+shift,
+                        band[start_idx : start_idx + nbands] + shift,
                         linewidth=linewidth,
                         color=color,
                         linestyle=linestyle,
@@ -91,7 +91,7 @@ class MagnonBand:
             for band in self.energies.T:
                 ax.plot(
                     self.xcoords,
-                    band+shift,
+                    band + shift,
                     linewidth=linewidth,
                     color=color,
                     linestyle=linestyle,

TB2J/magnon/magnon_dos.py

@@ -0,0 +1,308 @@
+"""Module for magnon density of states calculations and plotting."""
+
+import json
+from dataclasses import dataclass
+from pathlib import Path
+from typing import Optional
+
+import matplotlib.pyplot as plt
+import numpy as np
+from ase.dft.dos import DOS
+
+from TB2J.kpoints import monkhorst_pack
+
+
+@dataclass
+class MagnonDOS:
+    """Data class for storing magnon DOS data"""
+
+    energies: np.ndarray  # DOS energy points in meV
+    dos: np.ndarray  # DOS values in states/meV
+    weights: Optional[np.ndarray] = None  # k-point weights
+    kpoints: Optional[np.ndarray] = None  # k-points used for DOS
+
+    def save(self, filename: str):
+        """Save DOS data to a JSON file.
+
+        Parameters
+        ----------
+        filename : str
+            Output filename (should end in .json)
+        """
+        # Convert numpy arrays to lists for JSON serialization
+        data = {
+            "energies": self.energies.tolist(),
+            "dos": self.dos.tolist(),
+        }
+        if self.weights is not None:
+            data["weights"] = self.weights.tolist()
+        if self.kpoints is not None:
+            data["kpoints"] = self.kpoints.tolist()
+
+        with open(filename, "w") as f:
+            json.dump(data, f)
+
+    @classmethod
+    def load(cls, filename: str) -> "MagnonDOS":
+        """Load DOS data from a JSON file.
+
+        Parameters
+        ----------
+        filename : str
+            Input JSON filename
+
+        Returns
+        -------
+        MagnonDOS
+            Loaded DOS object
+        """
+        with open(filename) as f:
+            data = json.load(f)
+
+        # Convert lists back to numpy arrays
+        data["energies"] = np.array(data["energies"])
+        data["dos"] = np.array(data["dos"])
+        if "weights" in data:
+            data["weights"] = np.array(data["weights"])
+        if "kpoints" in data:
+            data["kpoints"] = np.array(data["kpoints"])
+
+        return cls(**data)
+
+    def plot(self, ax=None, color="blue", show=True, filename=None, **plot_kwargs):
+        """Plot the magnon DOS.
+
+        Parameters
+        ----------
+        ax : matplotlib.axes.Axes, optional
+            Axis to plot on. If None, creates new figure
+        color : str, optional
+            Color for DOS line
+        show : bool, optional
+            Whether to show plot
+        filename : str, optional
+            If provided, saves plot to file
+        **plot_kwargs : dict
+            Additional keyword arguments passed to plot
+
+        Returns
+        -------
+        matplotlib.axes.Axes
+            The plotting axes
+        """
+        if ax is None:
+            _, ax = plt.subplots()
+
+        ax.plot(self.energies, self.dos, color=color, **plot_kwargs)
+        ax.set_xlabel("Energy (meV)")
+        ax.set_ylabel("DOS (states/meV)")
+        ax.set_title("Magnon DOS")
+
+        if filename:
+            plt.savefig(filename)
+        if show:
+            plt.show()
+
+        return ax
+
+
+class MagnonDOSCalculator:
+    """Calculator for magnon density of states"""
+
+    def __init__(self, magnon):
+        """Initialize DOS calculator
+
+        Parameters
+        ----------
+        magnon : Magnon
+            Magnon object containing exchange parameters
+        """
+        self.magnon = magnon
+        self.kpts = None
+        self.weights = None
+        self.dos_args = {}
+
+    def estimate_energy_range(self, padding_factor=1.2):
+        """Estimate the energy range of eigenvalues.
+
+        Computes eigenvalues at zone center and high-symmetry points at zone boundaries
+        to estimate the full range of magnon energies.
+
+        Parameters
+        ----------
+        padding_factor : float, optional
+            Factor to extend the energy window beyond min/max values.
+            Default is 1.2 (20% padding).
+
+        Returns
+        -------
+        tuple
+            (min_energy, max_energy) in eV
+        """
+        # Generate high-symmetry points
+        kpoints = np.array(
+            [
+                [0.0, 0.0, 0.0],  # Γ (zone center)
+                [0.5, 0.0, 0.0],  # X
+                [0.5, 0.5, 0.0],  # M
+                [0.5, 0.5, 0.5],  # R (zone corner)
+                [0.0, 0.5, 0.0],  # Y
+                [0.0, 0.0, 0.5],  # Z
+            ]
+        )
+
+        # Calculate eigenvalues at these points
+        evals = self.magnon._magnon_energies(kpoints)
+        min_energy = evals.min()
+        max_energy = evals.max()
+
+        # Add padding and convert to eV
+        window_size = max_energy - min_energy
+        min_energy = min_energy - (padding_factor - 1) * window_size
+        max_energy = max_energy + (padding_factor - 1) * window_size
+
+        return min_energy, max_energy
+
+    def set_kmesh(self, kmesh=[9, 9, 9], gamma=True):
+        """Set k-point mesh for DOS calculation.
+
+        Parameters
+        ----------
+        kmesh : list, optional
+            Number of k-points along each direction
+        gamma : bool, optional
+            Whether to include Gamma point
+        """
+        self.kpts = monkhorst_pack(kmesh, gamma_center=gamma)
+        self.weights = np.ones(len(self.kpts)) / len(self.kpts)
+
+    def get_fermi_level(self):
+        return 0.0  # Fermi energy is not used in magnon calculations
+
+    def get_eigenvalues(self, kpt, spin=0):
+        """Get eigenvalues at a k-point.
+
+        Parameters
+        ----------
+        kpt : int
+            K-point index
+        spin : int, optional
+            Spin index (unused)
+
+        Returns
+        -------
+        numpy.ndarray
+            Eigenvalues in eV
+        """
+        kpoint = self.kpts[kpt]
+        evals = self.magnon._magnon_energies(np.array([kpoint]))[0]
+        return evals
+
+    def get_dos(self, width=0.1, window=None, npts=1001):
+        """Calculate DOS using ASE's DOS module.
+
+        Parameters
+        ----------
+        width : float, optional
+            Gaussian smearing width in eV
+        window : tuple, optional
+            Energy window (min, max) in eV
+        npts : int, optional
+            Number of energy points
+
+        Returns
+        -------
+        MagnonDOS
+            Calculated DOS object
+        """
+        if self.kpts is None:
+            self.set_kmesh()
+
+        # Estimate energy window if not provided
+        if window is None:
+            window = self.estimate_energy_range()
+
+        dos_calc = DOS(self, width=width, window=window, npts=npts)
+        energies = dos_calc.get_energies()
+        dos_vals = dos_calc.get_dos()
+
+        # Convert to meV
+        energies = energies * 1000  # eV to meV
+        dos_vals = dos_vals / 1000  # states/eV to states/meV
+
+        return MagnonDOS(
+            energies=energies,
+            dos=dos_vals,
+            weights=self.weights,
+            kpoints=self.kpts,
+        )
+
+    def get_number_of_spins(self):
+        """Required by ASE DOS calculator."""
+        return 1
+
+    def get_k_point_weights(self):
+        """Required by ASE DOS calculator."""
+        return self.weights
+
+    def get_bz_k_points(self):
+        """Required by ASE DOS calculator."""
+        return self.kpts
+
+
+def plot_magnon_dos(
+    magnon,
+    kmesh=[9, 9, 9],
+    gamma=True,
+    width=0.0005,
+    window=None,
+    xlim=None,
+    npts=1001,
+    filename=None,
+    save_data=True,
+    show=True,
+):
+    """Convenience function to calculate and plot magnon DOS.
+
+    Parameters
+    ----------
+    magnon : Magnon
+        Magnon object containing exchange parameters
+    kmesh : list, optional
+        Number of k-points along each direction
+    gamma : bool, optional
+        Whether to include Gamma point
+    width : float, optional
+        Gaussian smearing width in eV
+    window : tuple, optional
+        Energy window (min, max) in eV
+    npts : int, optional
+        Number of energy points
+    filename : str, optional
+        Output filename for plot
+    save_data : bool, optional
+        Whether to save DOS data to JSON
+    show : bool, optional
+        Whether to show plot
+
+    Returns
+    -------
+    MagnonDOS
+        The calculated DOS object
+    """
+    calculator = MagnonDOSCalculator(magnon)
+    calculator.set_kmesh(kmesh=kmesh, gamma=gamma)
+    dos = calculator.get_dos(width=width, window=window, npts=npts)
+
+    # Plot DOS
+    dos.plot(filename=filename, show=show)
+
+    # Save data if requested
+    if save_data:
+        data_file = (
+            Path(filename).with_suffix(".json") if filename else Path("magnon_dos.json")
+        )
+        dos.save(data_file)
+        print(f"DOS data saved to {data_file}")
+
+    return dos

TB2J/magnon/magnon_math.py

@@ -32,7 +32,7 @@ def get_rotation_arrays(magmoms, u=uz):
     v = magmoms
     n = np.cross(u, v)
     n /= np.linalg.norm(n, axis=-1).reshape(dim, 1)
-    #z = u #np.repeat(u, dim, axis=0)
+    # z = u #np.repeat(u, dim, axis=0)
     z = np.repeat(u, dim, axis=0)
     A = np.stack([z, np.cross(n, z), n], axis=1)
     B = np.stack([v, np.cross(n, v), n], axis=1)

TB2J/magnon/plot_magnon_dos_cli.py

@@ -0,0 +1,99 @@
+#!/usr/bin/env python3
+"""Command-line tool for plotting magnon DOS from TB2J results."""
+
+import argparse
+from pathlib import Path
+
+from TB2J.magnon.magnon3 import Magnon
+from TB2J.magnon.magnon_dos import plot_magnon_dos
+
+
+def main():
+    parser = argparse.ArgumentParser(
+        description="Calculate and plot magnon DOS from TB2J results"
+    )
+    parser.add_argument(
+        "--path",
+        default="TB2J_results",
+        help="Path to TB2J results directory (default: TB2J_results)",
+    )
+    parser.add_argument(
+        "--kmesh",
+        type=int,
+        nargs=3,
+        default=[20, 20, 20],
+        metavar=("nx", "ny", "nz"),
+        help="k-point mesh dimensions (default: 20, 20, 20)",
+    )
+    parser.add_argument(
+        "--no-gamma",
+        action="store_false",
+        dest="gamma",
+        help="Exclude Gamma point from k-mesh",
+    )
+    parser.add_argument(
+        "--width",
+        type=float,
+        default=0.001,
+        help="Gaussian smearing width in eV (default: 0.001)",
+    )
+    parser.add_argument(
+        "--window",
+        type=float,
+        nargs=2,
+        metavar=("emin", "emax"),
+        help="Energy window in meV (optional)",
+    )
+    parser.add_argument(
+        "--npts",
+        type=int,
+        default=401,
+        help="Number of energy points (default: 401)",
+    )
+    parser.add_argument(
+        "--output",
+        default="magnon_dos.png",
+        help="Output filename for plot (default: magnon_dos.png)",
+    )
+    parser.add_argument(
+        "-show",
+        action="store_true",
+        dest="show",
+        help="Show plot window",
+    )
+
+    args = parser.parse_args()
+
+    # Check if TB2J results exist
+    if not Path(args.path).exists():
+        raise FileNotFoundError(f"TB2J results not found at {args.path}")
+
+    # Load magnon calculator
+    print(f"Loading exchange parameters from {args.path}...")
+    magnon = Magnon.from_TB2J_results(path=args.path)
+
+    # Convert window from meV to eV if provided
+    window = None
+    if args.window is not None:
+        window = (args.window[0] / 1000, args.window[1] / 1000)
+
+    # Calculate and plot DOS
+    print("\nCalculating magnon DOS...")
+    _dos = plot_magnon_dos(
+        magnon,
+        kmesh=args.kmesh,
+        gamma=args.gamma,
+        width=args.width,
+        window=window,
+        npts=args.npts,
+        filename=args.output,
+        show=args.show,
+    )
+
+    print(f"\nPlot saved to {args.output}")
+    data_file = Path(args.output).with_suffix(".json")
+    print(f"DOS data saved to {data_file}")
+
+
+if __name__ == "__main__":
+    main()

{tb2j-0.9.9.13.data → tb2j-0.9.9.15.data}/scripts/TB2J_merge.py

@@ -1,8 +1,8 @@
 #!python
 import argparse
-import os
-import sys
+
 from TB2J.io_merge import merge
+from TB2J.version_info import print_license
 
 
 def main():
@@ -32,14 +32,18 @@ def main():
         "--main_path",
         help="The path containning the reference structure.",
         type=str,
-        default=None
+        default=None,
     )
 
     args = parser.parse_args()
     # merge(*(args.directories), args.type.strip().lower(), path=args.output_path)
     # merge(*(args.directories), method=args.type.strip().lower(), path=args.output_path)
-    #merge2(args.directories, args.type.strip().lower(), path=args.output_path)
+    # merge2(args.directories, args.type.strip().lower(), path=args.output_path)
+    print_license()
+    print("Merging the TB2J results from the following directories: ", args.directories)
     merge(*args.directories, main_path=args.main_path, write_path=args.output_path)
+    print("Merging completed. The results are saved in:", args.output_path)
 
 
-main()
+if __name__ == "__main__":
+    main()

tb2j-0.9.9.15.data/scripts/TB2J_plot_magnon_bands.py

@@ -0,0 +1,27 @@
+#!python
+
+
+if __name__ == "__main__":
+    # add a warning messege that this functionality is under development and should not be used in production.
+    # make it visually distinct, e.g. with a different color or formatting.
+    import warnings
+
+    from TB2J.magnon.magnon3 import main
+
+    warnings.warn(
+        """ 
+        # !!!!!!!!!!!!!!!!!! WARNING: =============================
+        # 
+        # This functionality is under development and should not be used in production.
+        # It is provided for testing and development purposes only.
+        # Please use with caution and report any issues to the developers.
+        #
+        # This warning will be removed in future releases.
+        # =====================================
+
+        """,
+        UserWarning,
+        stacklevel=2,
+    )
+    # Call the main function from the magnons module
+    main()

{tb2j-0.9.9.13.dist-info → tb2j-0.9.9.15.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: TB2J
-Version: 0.9.9.13
+Version: 0.9.9.15
 Summary: TB2J: First principle to Heisenberg exchange J using tight-binding Green function method
 Author: Xu He
 Author-email: mailhexu@gmail.com

{tb2j-0.9.9.13.dist-info → tb2j-0.9.9.15.dist-info}/RECORD

@@ -1,17 +1,17 @@
 TB2J/Jdownfolder.py,sha256=n5BeQCYP4mD9JsAPeE1F3ZKKR3SUxADfDbaG_rzi77k,9658
 TB2J/Jtensor.py,sha256=Wi06AAbfKFU6f2a2jkFr9zU2cwwfVroajTp-dwCtkTE,3160
 TB2J/MAE.py,sha256=fM8U-Dgp9HcQOEeC_kyZV1oVrygBvcux9BraUXVouvY,10994
-TB2J/MAEGreen.py,sha256=DHoU6nlM0Fk42T_SYd_s-VafB6wCZNzVzRXhuv-hCwM,15762
+TB2J/MAEGreen.py,sha256=zVLjQtFDFcRmC8PYm7MOgKnkOuFJEMKr-eoNBABUsMc,15759
 TB2J/Oiju.py,sha256=cNGv8N5uH_swGq7cnAt2OyiDfqtjLlLrwseGu0E4iaM,3383
 TB2J/Oiju_epc.py,sha256=oytM3NYW7nWmklrGgNlqwIpI_JYv_hb7ZnR4o9nYNog,6809
-TB2J/__init__.py,sha256=TcMW0FTsZFR9McJTh8TBwZl2lylHc2uwsxVV9PsgSwY,25
+TB2J/__init__.py,sha256=VIRI2sifGP0sSwxkmWEp9IvuVn2fBhvrNpmqbwLrG8s,25
 TB2J/anisotropy.py,sha256=0zmvXkmDmakbBOwGYLa3IIkv5cE99SHLAQJsGoZz7JQ,25463
 TB2J/basis.py,sha256=DFo6_QUwjBwisP6zGxvoO0lpGTMDPAOkiL9giNCjOjA,1558
 TB2J/citation.py,sha256=gcQeyJZaT1Qrtsl8Y3s4neOH3-vvgmIcCvXeV2o3vj0,2891
 TB2J/contour.py,sha256=zLHQZZ3hhgLkLFPATCraLOJyLJKLC0fba_L_5sRz23o,3246
 TB2J/density_matrix.py,sha256=D5k8Oe21OCiLVORNYbo4TZOFG0slrQSbj91kJ3TMFjs,1514
 TB2J/epc.py,sha256=zLbtqZJhDr8DnnGN6YENcXwrMb3Qxu6KB08mLy9Pw20,3474
-TB2J/exchange.py,sha256=rJ7XeHOo6j9A0alauM2XQ5DvTO741dq_pDQQxyyn8mE,26981
+TB2J/exchange.py,sha256=-BEXSQqx374cRbJWjvjyJAZkchzNK1uPUDvTYox3Axs,26982
 TB2J/exchangeCL2.py,sha256=P7bklMXVYX_tn9DbjEPqeTir1SeZyfPBIP1fhWUzLmY,11419
 TB2J/exchange_params.py,sha256=VW9nGVio6M_Ub9-36L_LExhjgdD1E_joYpI8AxmM360,8029
 TB2J/exchange_pert.py,sha256=jmFMtQbYa_uczM4VAeS6TijkIHRFIqEzZJswzE9Wfuo,8523
@@ -19,7 +19,7 @@ TB2J/exchange_qspace.py,sha256=ZL68qBGFUaQ9BsSPsJaaoWOr9RssPiqX34R_9I3nk_8,8436
 TB2J/gpaw_wrapper.py,sha256=aJ--9Dtyq7jOP1Hkh-Sh1nWcfXm6zKcljOCO0DNCAr0,6890
 TB2J/green.py,sha256=ySXjoV3Cj_E2C41dRfc3TkS7M4vmM6qcHXkt-_-jxsY,17071
 TB2J/greentest.py,sha256=2ISSfhor9ecSEOi_E6b4Cv26wEIQlwlzca0ru8z44_E,1603
-TB2J/io_merge.py,sha256=xj5sIo5o62HYwSnijGycem_rAc8kSzY-HI2IJDRhzyE,7405
+TB2J/io_merge.py,sha256=mgoc8ZUP_a_g-9KvWWISp25qo6aBJ8TQsDsBrMSB5nk,7509
 TB2J/kpoints.py,sha256=9L7tBarFBHoIhpuc9zuwA6HdnlgH834SQrPek4yRoWk,3191
 TB2J/myTB.py,sha256=ok_B4my29bOIghMSZfx0Es6G8FaXaIiLP4gPxTdSj00,17659
 TB2J/mycfr.py,sha256=ZF1PEE2khlKd_4gPyMkoNXepX3XqwWAL2kDbRJNVX-Y,3908
@@ -43,31 +43,33 @@ TB2J/interfaces/__init__.py,sha256=4tiLoXQ73Nlyi9L4j8jJXOYzXluVNPxQZkwfkQZEGHg,3
 TB2J/interfaces/gpaw_interface.py,sha256=GCDlJ-hRWfChvWwsgBDYSmVqO4sH9HAuGZTV9GqgN6c,1504
 TB2J/interfaces/lawaf_interface.py,sha256=PieLnmppdafOYsgeHznqOou1g9L1sam5jOm3KaObdqo,4408
 TB2J/interfaces/manager.py,sha256=PQMLEfMCT5GnDWSl2nI4JOgRPm_fysyR-6Y6l97xWcw,860
-TB2J/interfaces/siesta_interface.py,sha256=Y4rRn_MRZUOI4y4btgCKTtVy8iXrSpaBuK7SMCJ1C2o,7594
+TB2J/interfaces/siesta_interface.py,sha256=Te9-Au0ZGAt-LfEnG8exUFD4lMj6uZu8RmHjopb0XLU,7595
 TB2J/interfaces/wannier90_interface.py,sha256=qzRgXUBb7t1Aiegrl_RV51BB8csdtVM0EP0Z4pjmTcs,4467
 TB2J/interfaces/abacus/__init__.py,sha256=leas71oCvM_HxrF4gnO5A_VKcJmDAgsI1BUctLU3OBw,177
 TB2J/interfaces/abacus/abacus_api.py,sha256=lNV4LNkLcKw7Zux4MQYM9wnh3eFTlcSqbf4Pb7pqhrk,7243
 TB2J/interfaces/abacus/abacus_wrapper.py,sha256=bQFnvvy-0hruTavH_VspMk1wD32t2UFybEkCgwkwle0,11941
 TB2J/interfaces/abacus/gen_exchange_abacus.py,sha256=v-AUHGkJWeMNf4D5A4wOpCjM6DygQsFB6SWu-BGdTxM,3262
 TB2J/interfaces/abacus/orbital_api.py,sha256=QEAyy4y_uM5GnwRyLvmShjd_z5A2fwqk7L8yj0Y2Mgc,1577
-TB2J/interfaces/abacus/stru_api.py,sha256=UsIs55Y_1AM31N2oCjwmgHIMgRIcusKA6g3K2tjr3uM,67945
+TB2J/interfaces/abacus/stru_api.py,sha256=Hb9MCo_4d_DVEEhntvBNAFw5KQ46TkAADn4z7CT6qOk,67801
 TB2J/interfaces/abacus/test_density_matrix.py,sha256=bMWWJYtDS57SpPZ-eZXZ9Hr_UK4mv8ZHM7SzItG3IVA,774
 TB2J/interfaces/abacus/test_read_HRSR.py,sha256=W1oO_yigT50Yb5_u-KB_IfTpM7kArGkBuMSMs0H4CTs,1235
 TB2J/interfaces/abacus/test_read_stru.py,sha256=hoKPHVco8vwzC7Gao4bOPCdAPhh29x-9DTJJqRr7AYM,788
 TB2J/io_exchange/__init__.py,sha256=LqEnG67qDVKt4hCUywDEQvUIJ7jsQjmtueOW_J16NOE,54
-TB2J/io_exchange/io_exchange.py,sha256=KnR2UyvMSoBye8dD6ETqZW6nld-CcYchcCrAmHbVsBA,22735
+TB2J/io_exchange/io_exchange.py,sha256=TQPsxpV0vyFq-HFRRt9GQbpnmRw-5KwyM26eb9fKaoY,22743
 TB2J/io_exchange/io_multibinit.py,sha256=8PDmWxzGuv-GwJosj2ZTmiyNY_duFVWJ4ekCuSqGdd8,6739
 TB2J/io_exchange/io_tomsasd.py,sha256=NqkAC1Fl-CUnFA21eBzSy_S5F_oeQFJysw4UukQbN8o,4173
 TB2J/io_exchange/io_txt.py,sha256=BMr1eSILlKpgtjvDx7uw2VMAkEKSvGEPNxpaT_zev0I,10547
 TB2J/io_exchange/io_uppasd.py,sha256=bI4iPEgnK4TvCZNvb6x2xYXgjW7pEehCqmcizy2pqFU,3301
-TB2J/io_exchange/io_vampire.py,sha256=u4NZhqoC_JBcVq19WZfBkyM5p8u-OMuRoUrtFIcPp5E,5762
+TB2J/io_exchange/io_vampire.py,sha256=vOStLmtCiWLp9GPhZpsAmrtaRHg9KSmtOM2Fky6yCQA,5762
 TB2J/magnon/__init__.py,sha256=Q69duroGIIotgW_71ZdHYarSiJLGAu9NPYgEacUk6eI,117
 TB2J/magnon/io_exchange2.py,sha256=EcC3x6H13qq61WBsr__xKzHDtSvy_dMz1tEdUaqSG2I,23265
-TB2J/magnon/magnon3.py,sha256=75YmSTT652q7lLEmJmn4-TbIkwJnjjelMfu0eyEzjA0,27401
-TB2J/magnon/magnon_band.py,sha256=8HJWl1zcYSxJQflFUj4qaITnTWE-rhF5zQ1YgQN78oU,6098
+TB2J/magnon/magnon3.py,sha256=PFZFiyfYTIBscVu8SzvGxtrXhT-QM66EwIhWOqgSYQE,30219
+TB2J/magnon/magnon_band.py,sha256=ZLHK1qf2Clu9jlcvBwBCBEi1MrPfy2BWrow8JJ2v5pk,6103
+TB2J/magnon/magnon_dos.py,sha256=uztzsxCmFZpDHl-JziaUb_NNwBeZw9-L2ozd3XNTvxI,8506
 TB2J/magnon/magnon_io.py,sha256=H4bmzCcnh1D3Sb6UBIIKWa6jIrA20dg9lX4wfLXHEjo,1241
-TB2J/magnon/magnon_math.py,sha256=b4Po_aZm8aIhr5u2nbjjtDDMTKV_ecI2d6Fc_AhiuNc,1358
+TB2J/magnon/magnon_math.py,sha256=0mVQY5tXe6fHsfijivXf-1-Tk-Ku02Gv7-TdAd8qUpc,1359
 TB2J/magnon/plot.py,sha256=oNavax15jRfW4Sxl6FgwLXSrMntyVz2cgVDbmgkWehw,3193
+TB2J/magnon/plot_magnon_dos_cli.py,sha256=XVozodrIO_w4FBhp-ZcpIZbyGWLdegUFLdvY5zL3XOI,2592
 TB2J/magnon/structure.py,sha256=rKefzXgQlEjVvV-A7E2IogVDWwf5egZr3Wxxu6HuJU4,7685
 TB2J/mathutils/__init__.py,sha256=E70Mx8mLXh3DJGmdN3TLWmGYMcsbvKxmgxUbAI6Mzmo,49
 TB2J/mathutils/auto_kpath.py,sha256=bDXDyQzRByVpPSOpz7NU7F5hYqS6IX-47lwWoMU89lM,4978
@@ -89,21 +91,21 @@ TB2J/spinham/supercell.py,sha256=y17uUC6r3gQb278FhxIW4CABihfLTvKFj6flyXrCPR8,122
 TB2J/wannier/__init__.py,sha256=7ojCbM84PYv1X1Tbo4NHI-d3gWmQsZB_xiYqbfxVV1E,80
 TB2J/wannier/w90_parser.py,sha256=dbd63LuKyv2DVUzqRINGsbDzEsOxsQyE8_Ear_LQIRg,4620
 TB2J/wannier/w90_tb_parser.py,sha256=qt8pnuprmPp9iIAYwPkPbmEzk6ZPgMq2xognoQp7vwc,4610
-tb2j-0.9.9.13.data/scripts/TB2J_downfold.py,sha256=i4BVqnpDdgrX_amookVWeLGefGBn-qeAutWiwuY9SfQ,2099
-tb2j-0.9.9.13.data/scripts/TB2J_eigen.py,sha256=Qs9v2hnMm2Tpfoa4h53muUKty2dZjwx8948MBoQooNg,1128
-tb2j-0.9.9.13.data/scripts/TB2J_magnon.py,sha256=q7UwAmorRcFNk4tfE7gl_ny05l6p7pbD9Wm_LkIpKEw,3101
-tb2j-0.9.9.13.data/scripts/TB2J_magnon2.py,sha256=tMa7Fg_Wd2UytnrH3C_AsgGM7BciUW0iy6NiPlWvar8,1920
-tb2j-0.9.9.13.data/scripts/TB2J_magnon_dos.py,sha256=TMXQvD2dIbO5FZ4tUMmxJgCgH2O2hDAPUNfEKO4z-x4,110
-tb2j-0.9.9.13.data/scripts/TB2J_merge.py,sha256=y834SF4rIRn1L1ptkhczvavQpC-8Px6DTmDOOSaq_DE,1854
-tb2j-0.9.9.13.data/scripts/TB2J_plot_magnon_bands.py,sha256=qUbisPFLvW_j1IdTpc_jU0WiT-7CT-WErbkhNvAPM0Y,91
-tb2j-0.9.9.13.data/scripts/TB2J_rotate.py,sha256=zgiDFuYZNmzKK0rwDmTaYD2OpRlmKA_VGeBx83w2Xwc,873
-tb2j-0.9.9.13.data/scripts/TB2J_rotateDM.py,sha256=kCvF7gotuqAX1VnJ06cwfVm7RrhrdtiV5v7d9P2Pn_E,567
-tb2j-0.9.9.13.data/scripts/abacus2J.py,sha256=ozYI7qZyja1WEs9oCYVpeYAfVj5PGhehhmdZFcvrd3g,1795
-tb2j-0.9.9.13.data/scripts/siesta2J.py,sha256=QJ6c0DbqxaqYEesxiL5R9nK9-flNLrr7hajKfCwirYc,2318
-tb2j-0.9.9.13.data/scripts/wann2J.py,sha256=OA31VHEXbQMD-JozoLUHDF6vB9Sr62d804OApSKtSnU,3240
-tb2j-0.9.9.13.dist-info/licenses/LICENSE,sha256=CbZI-jyRTjiqIcWa244cRSHJdjjtUNqGR4HeJkgEwJw,1332
-tb2j-0.9.9.13.dist-info/METADATA,sha256=XJtwVxFDDV27HFqNVsgHrjlRbY3g-dKnsc0FiLjJ3kg,1768
-tb2j-0.9.9.13.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-tb2j-0.9.9.13.dist-info/entry_points.txt,sha256=Hdz1WC9waUzyFVmowKnbbZ6j-J4adHh_Ko6JpxGYAtE,131
-tb2j-0.9.9.13.dist-info/top_level.txt,sha256=whYa5ByLYhl5XnTPBHSWr-IGD6VWmr5Ql2bye2qwV_s,5
-tb2j-0.9.9.13.dist-info/RECORD,,
+tb2j-0.9.9.15.data/scripts/TB2J_downfold.py,sha256=i4BVqnpDdgrX_amookVWeLGefGBn-qeAutWiwuY9SfQ,2099
+tb2j-0.9.9.15.data/scripts/TB2J_eigen.py,sha256=Qs9v2hnMm2Tpfoa4h53muUKty2dZjwx8948MBoQooNg,1128
+tb2j-0.9.9.15.data/scripts/TB2J_magnon.py,sha256=q7UwAmorRcFNk4tfE7gl_ny05l6p7pbD9Wm_LkIpKEw,3101
+tb2j-0.9.9.15.data/scripts/TB2J_magnon2.py,sha256=tMa7Fg_Wd2UytnrH3C_AsgGM7BciUW0iy6NiPlWvar8,1920
+tb2j-0.9.9.15.data/scripts/TB2J_magnon_dos.py,sha256=TMXQvD2dIbO5FZ4tUMmxJgCgH2O2hDAPUNfEKO4z-x4,110
+tb2j-0.9.9.15.data/scripts/TB2J_merge.py,sha256=WQVrS2FvrjxkKukiPLY-naQxEK-d2T-bAb5n7e6ml2I,2096
+tb2j-0.9.9.15.data/scripts/TB2J_plot_magnon_bands.py,sha256=wLlueNI-wU79cmPFu4LmPDwVBMerKSb5sUlTDSPQK6I,874
+tb2j-0.9.9.15.data/scripts/TB2J_rotate.py,sha256=zgiDFuYZNmzKK0rwDmTaYD2OpRlmKA_VGeBx83w2Xwc,873
+tb2j-0.9.9.15.data/scripts/TB2J_rotateDM.py,sha256=kCvF7gotuqAX1VnJ06cwfVm7RrhrdtiV5v7d9P2Pn_E,567
+tb2j-0.9.9.15.data/scripts/abacus2J.py,sha256=ozYI7qZyja1WEs9oCYVpeYAfVj5PGhehhmdZFcvrd3g,1795
+tb2j-0.9.9.15.data/scripts/siesta2J.py,sha256=QJ6c0DbqxaqYEesxiL5R9nK9-flNLrr7hajKfCwirYc,2318
+tb2j-0.9.9.15.data/scripts/wann2J.py,sha256=OA31VHEXbQMD-JozoLUHDF6vB9Sr62d804OApSKtSnU,3240
+tb2j-0.9.9.15.dist-info/licenses/LICENSE,sha256=CbZI-jyRTjiqIcWa244cRSHJdjjtUNqGR4HeJkgEwJw,1332
+tb2j-0.9.9.15.dist-info/METADATA,sha256=4MA4Id_pqZPASA3ig_9xiICywGLB06rlUT6bPiccdW4,1768
+tb2j-0.9.9.15.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+tb2j-0.9.9.15.dist-info/entry_points.txt,sha256=S81PATMTaKCIznjd3pei3Kb3q3xOy2dU-x_cAFSCS7k,194
+tb2j-0.9.9.15.dist-info/top_level.txt,sha256=whYa5ByLYhl5XnTPBHSWr-IGD6VWmr5Ql2bye2qwV_s,5
+tb2j-0.9.9.15.dist-info/RECORD,,

{tb2j-0.9.9.13.dist-info → tb2j-0.9.9.15.dist-info}/entry_points.txt

@@ -1,3 +1,4 @@
 [console_scripts]
+TB2J_plot_magnon_dos.py = TB2J.magnon.plot_magnon_dos_cli:main
 TB2J_symmetrize.py = TB2J.symmetrize_J:symmetrize_J_cli
 lawaf2J.py = TB2J.interfaces.lawaf_interface:lawaf2J_cli

tb2j-0.9.9.13.data/scripts/TB2J_plot_magnon_bands.py

@@ -1,7 +0,0 @@
-#!python
-
-
-if __name__ == "__main__":
-    from TB2J.magnon.magnon3 import main
-
-    main()