-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Add XY Model with coupling distribution.
- Loading branch information
Showing
3 changed files
with
138 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 1,87 @@ | ||
| #include "xy_inhom.hpp" | ||
| #include <cmath> | ||
|
|
||
| namespace XYInhom | ||
| { | ||
| XYInhom::XYInhom(int L, const Real &J0, const CouplingDistribution &coupling_distribution, const Real &hx, | ||
| const Real &hy, const Real &hz) | ||
| : SpinHalf1D(L, false, false, false), m_J0(J0), m_coupling_distribution(coupling_distribution), m_hx(hx), m_hy(hy), | ||
| m_hz(hz) | ||
| { | ||
| m_name = "XYInhom"; | ||
| if (m_L % 2) | ||
| { | ||
| throw std::domain_error("L must be even"); | ||
| } | ||
|
|
||
| std::function<Real(int)> dist; | ||
| switch (m_coupling_distribution) | ||
| { | ||
| case Uniform: | ||
| dist = [](int r) { | ||
| (void)r; | ||
| return 1.0; | ||
| }; | ||
| break; | ||
| case Exponential: | ||
| dist = [](int r) { return std::exp(-static_cast<Real>(r - 1)); }; | ||
| break; | ||
| case Gaussian: | ||
| dist = [](int r) { return std::exp(-static_cast<Real>((r - 1) * (r - 1))); }; | ||
| break; | ||
| case Logarithmic: | ||
| dist = [](int r) { return 1.0 / std::log2(static_cast<Real>(r 1)); }; | ||
| break; | ||
| default: | ||
| throw std::domain_error("invalid coupling distribution"); | ||
| break; | ||
| } | ||
|
|
||
| for (int i = 0; i < m_L; i) | ||
| { | ||
| for (int j = 0; j < i; j) | ||
| { | ||
| const auto J = m_J0 * dist(std::abs(i - j)); | ||
| m_two_body_terms.emplace_back(4.0 * J, "S ", i, "S-", j); | ||
| m_two_body_terms.emplace_back(4.0 * J, "S-", i, "S ", j); | ||
| } | ||
| } | ||
|
|
||
| for (int i = 0; i < m_L; i) | ||
| { | ||
| m_one_body_terms.emplace_back(-2.0 * m_hx, "Sx", i); | ||
| m_one_body_terms.emplace_back(-2.0 * m_hy, "Sy", i); | ||
| m_one_body_terms.emplace_back(-2.0 * m_hz, "Sz", i); | ||
| } | ||
| } | ||
|
|
||
| XYInhom::XYInhom(const json &js) | ||
| : XYInhom(js.at("L").get<int>(), jsonGetDefault<Real>(js, "J0", 1.0), | ||
| jsonGetDefault<CouplingDistribution>(js, "coupling_distribution", CouplingDistribution::Uniform), | ||
| jsonGetDefault<Real>(js, "hx", 0.0), jsonGetDefault<Real>(js, "hy", 0.0), | ||
| jsonGetDefault<Real>(js, "hz", 0.0)) | ||
| { | ||
| } | ||
|
|
||
| const Real &XYInhom::getJ0() const | ||
| { | ||
| return m_J0; | ||
| } | ||
| const CouplingDistribution &XYInhom::getCouplingDistribution() const | ||
| { | ||
| return m_coupling_distribution; | ||
| } | ||
| const Real &XYInhom::getHx() const | ||
| { | ||
| return m_hx; | ||
| } | ||
| const Real &XYInhom::getHy() const | ||
| { | ||
| return m_hy; | ||
| } | ||
| const Real &XYInhom::getHz() const | ||
| { | ||
| return m_hz; | ||
| } | ||
|
|
||
| } // namespace XYInhom |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 1,50 @@ | ||
| #ifndef MPSKIT_MODELS_SPIN_HALF_1D_XY_INHOM | ||
| #define MPSKIT_MODELS_SPIN_HALF_1D_XY_INHOM | ||
|
|
||
| #include "../../../json.hpp" | ||
| #include "../../../types.hpp" | ||
| #include "spin_half_1d.hpp" | ||
|
|
||
| namespace XYInhom | ||
| { | ||
| enum CouplingDistribution : int | ||
| { | ||
| Uniform, | ||
| Exponential, | ||
| Gaussian, | ||
| Logarithmic, | ||
| Invalid | ||
| }; | ||
|
|
||
| NLOHMANN_JSON_SERIALIZE_ENUM(CouplingDistribution, { | ||
| {XYInhom::Invalid, nullptr}, | ||
| {XYInhom::Uniform, "uniform"}, | ||
| {XYInhom::Exponential, "exponential"}, | ||
| {XYInhom::Gaussian, "gaussian"}, | ||
| {XYInhom::Logarithmic, "logarithmic"}, | ||
| }) | ||
|
|
||
| class XYInhom : public SpinHalf1D | ||
| { | ||
| private: | ||
| Real m_J0; | ||
| CouplingDistribution m_coupling_distribution; | ||
| Real m_hx; | ||
| Real m_hy; | ||
| Real m_hz; | ||
|
|
||
| public: | ||
| explicit XYInhom(int L, const Real &J0, const CouplingDistribution &coupling_distribution, const Real &hx = 0.0, | ||
| const Real &hy = 0.0, const Real &hz = 0.0); | ||
| explicit XYInhom(const json &js); | ||
|
|
||
| const Real &getJ0() const; | ||
| const CouplingDistribution &getCouplingDistribution() const; | ||
| const Real &getHx() const; | ||
| const Real &getHy() const; | ||
| const Real &getHz() const; | ||
| }; | ||
|
|
||
| } // namespace XYInhom | ||
|
|
||
| #endif /* MPSKIT_MODELS_SPIN_HALF_1D_XY_INHOM */ |