Skip to content

Commit

Permalink
Added HoppingMatrix, CrystalFieldSplitting to properties
Browse files Browse the repository at this point in the history
Added WignerSeitz to variables

Removed common.py module
  • Loading branch information
JosePizarro3 committed Apr 30, 2024
1 parent a7412b2 commit 4e3cfab
Show file tree
Hide file tree
Showing 4 changed files with 124 additions and 74 deletions.
74 changes: 0 additions & 74 deletions src/nomad_simulations/common.py

This file was deleted.

1 change: 1 addition & 0 deletions src/nomad_simulations/properties/__init__.py
Original file line number Diff line number Diff line change
Expand Up @@ -24,3 +24,4 @@
ElectronicDensityOfStates,
XASSpectra,
)
from .hopping_matrix import HoppingMatrix, CrystalFieldSplitting
99 changes: 99 additions & 0 deletions src/nomad_simulations/properties/hopping_matrix.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,99 @@
#
# Copyright The NOMAD Authors.
#
# This file is part of NOMAD. See https://nomad-lab.eu for further info.
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
#

# import typing
# from structlog.stdlib import BoundLogger
import numpy as np

from nomad.metainfo import Quantity, Section, Context

from ..physical_property import PhysicalProperty


class HoppingMatrix(PhysicalProperty):
"""
Transition probability between two atomic orbitals in a tight-binding model.
"""

iri = 'http://fairmat-nfdi.eu/taxonomy/HoppingMatrix'

n_orbitals = Quantity(
type=np.int32,
description="""
Number of orbitals in the tight-binding model.
""",
)

degeneracy_factors = Quantity(
type=np.int32,
shape=['*'],
description="""
Degeneracy of each real space point.
""",
)

value = Quantity(
type=np.complex128,
unit='joule',
description="""
Value of the hopping matrix in joules. The elements are complex numbers defined for each Wigner-Seitz point and
each pair of orbitals; thus, `rank = [n_orbitals, n_orbitals]`.
""",
)

def __init__(
self, m_def: Section = None, m_context: Context = None, **kwargs
) -> None:
super().__init__(m_def, m_context, **kwargs)
self.rank = [self.n_orbitals, self.n_orbitals]
self.name = self.m_def.name

def normalize(self, archive, logger) -> None:
super().normalize(archive, logger)


class CrystalFieldSplitting(PhysicalProperty):
"""
Energy difference between the degenerated orbitals of an ion in a crystal field environment.
"""

iri = 'http://fairmat-nfdi.eu/taxonomy/CrystalFieldSplitting'

n_orbitals = Quantity(
type=np.int32,
description="""
Number of orbitals in the tight-binding model.
""",
)

value = Quantity(
type=np.float64,
unit='joule',
description="""
Value of the crystal field splittings in joules.
""",
)

def __init__(
self, m_def: Section = None, m_context: Context = None, **kwargs
) -> None:
super().__init__(m_def, m_context, **kwargs)
self.rank = [self.n_orbitals]

def normalize(self, archive, logger) -> None:
super().normalize(archive, logger)
24 changes: 24 additions & 0 deletions src/nomad_simulations/variables.py
Original file line number Diff line number Diff line change
Expand Up @@ -129,3 +129,27 @@ def __init__(

def normalize(self, archive, logger) -> None:
super().normalize(archive, logger)


class WignerSeitz(Variables):
"""
Wigner-Seitz points in which the real space is discretized. This variable is used to define `HoppingMatrix(PhysicalProperty)` and
other inter-cell properties. See, e.g., https://en.wikipedia.org/wiki/Wigner–Seitz_cell.
"""

grid_points = Quantity(
type=np.float64,
shape=['n_grid_points', 3],
description="""
Wigner-Seitz points with respect to the origin cell, (0, 0, 0). These are 3D arrays stored in fractional coordinates.
""",
)

def __init__(
self, m_def: Section = None, m_context: Context = None, **kwargs
) -> None:
super().__init__(m_def, m_context, **kwargs)
self.name = self.m_def.name

def normalize(self, archive, logger) -> None:
super().normalize(archive, logger)

0 comments on commit 4e3cfab

Please sign in to comment.