vacspin 0.2.0__tar.gz → 0.3.0__tar.gz

{vacspin-0.2.0 → vacspin-0.3.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: vacspin
-Version: 0.2.0
+Version: 0.3.0
 Summary: Group-IV colour-centre spin-photon interfaces: strained spin Hamiltonians, cyclicity, transition-resolved Purcell budgets and exact single-shot readout statistics
 Author: Tanvir Mahmud Mahim
 License: Apache-2.0
@@ -128,8 +128,19 @@ No physical number is made up, and none is accepted without a source
   of Hepp et al., PRL 112, 036405 (2014). Strain and orbital
   quenching are sample-specific and deliberately not shipped; the
   reference string says so.
-
-For GeV-, PbV-, or your own sample, populate `SpinParameters` and
+- `gev_bhaskar2017()`: the measured GeV- splittings (152 and
+  981 GHz, zero-phonon line 602 nm) of Bhaskar et al., PRL 118,
+  223603 (2017), with the same deliberate omissions as SiV- --
+  strain and quenching are your sample's numbers, and `vacspin.lab`
+  exists to fit them.
+
+The lead-vacancy (PbV-) centre is documented but not shipped, on
+purpose: its ground splitting is measured (about 3900 GHz; Wang et
+al., ACS Photonics 8, 2947 (2021), confirmed in 2026), but no
+verified measured excited-state splitting was available, and this
+package does not ship half a parameter set.
+
+For PbV-, or your own sample of any centre, populate `SpinParameters` and
 `EmissionBudget` from your measurements or the literature; the
 provenance travels with every prediction.
 
@@ -174,7 +185,7 @@ round trip is exact.
 
 ## How it is checked
 
-56 tests (Python 3.9-3.13, run in CI on every push), every claim
+59 tests (Python 3.9-3.14 in CI on every push), every claim
 anchored to an exact result, a published measurement, or two
 independent code paths -- never a stored number. Highlights: exact
 Kramers doublets and the closed-form splitting against full

{vacspin-0.2.0 → vacspin-0.3.0}/README.md

@@ -105,8 +105,19 @@ No physical number is made up, and none is accepted without a source
   of Hepp et al., PRL 112, 036405 (2014). Strain and orbital
   quenching are sample-specific and deliberately not shipped; the
   reference string says so.
-
-For GeV-, PbV-, or your own sample, populate `SpinParameters` and
+- `gev_bhaskar2017()`: the measured GeV- splittings (152 and
+  981 GHz, zero-phonon line 602 nm) of Bhaskar et al., PRL 118,
+  223603 (2017), with the same deliberate omissions as SiV- --
+  strain and quenching are your sample's numbers, and `vacspin.lab`
+  exists to fit them.
+
+The lead-vacancy (PbV-) centre is documented but not shipped, on
+purpose: its ground splitting is measured (about 3900 GHz; Wang et
+al., ACS Photonics 8, 2947 (2021), confirmed in 2026), but no
+verified measured excited-state splitting was available, and this
+package does not ship half a parameter set.
+
+For PbV-, or your own sample of any centre, populate `SpinParameters` and
 `EmissionBudget` from your measurements or the literature; the
 provenance travels with every prediction.
 
@@ -151,7 +162,7 @@ round trip is exact.
 
 ## How it is checked
 
-56 tests (Python 3.9-3.13, run in CI on every push), every claim
+59 tests (Python 3.9-3.14 in CI on every push), every claim
 anchored to an exact result, a published measurement, or two
 independent code paths -- never a stored number. Highlights: exact
 Kramers doublets and the closed-form splitting against full

{vacspin-0.2.0 → vacspin-0.3.0}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "vacspin"
-version = "0.2.0"
+version = "0.3.0"
 description = "Group-IV colour-centre spin-photon interfaces: strained spin Hamiltonians, cyclicity, transition-resolved Purcell budgets and exact single-shot readout statistics"
 readme = "README.md"
 requires-python = ">=3.9"

{vacspin-0.2.0 → vacspin-0.3.0}/src/vacspin/__init__.py

@@ -14,8 +14,8 @@ diamond nanocavity on thin-film lithium niobate" (submitted to Optics
 Express, 2026), and the effective-Hamiltonian and readout formalism of
 Rosenthal et al., PRX 13, 031022 (2023) and arXiv:2403.13110.
 """
-from .params import (EmissionBudget, SpinParameters, siv_hepp2014,
-                     snv_emission, snv_rosenthal2023)
+from .params import (EmissionBudget, SpinParameters, gev_bhaskar2017,
+                     siv_hepp2014, snv_emission, snv_rosenthal2023)
 from .hamiltonian import (GAMMA_GHZ_PER_T, h_manifold, lab_to_spin,
                           qubit_frequency, qubit_frequency_perpendicular,
                           solve, zero_field_splitting)
@@ -31,10 +31,10 @@ from .lab import (SpinFit, design_fields, fit_spin_parameters,
                   load_observations_csv, parameter_information,
                   save_observations_csv)
 
-__version__ = "0.2.0"
+__version__ = "0.3.0"
 __all__ = [
     "SpinParameters", "EmissionBudget", "snv_rosenthal2023",
-    "snv_emission", "siv_hepp2014",
+    "snv_emission", "siv_hepp2014", "gev_bhaskar2017",
     "GAMMA_GHZ_PER_T", "h_manifold", "lab_to_spin", "solve",
     "zero_field_splitting", "qubit_frequency",
     "qubit_frequency_perpendicular",

{vacspin-0.2.0 → vacspin-0.3.0}/src/vacspin/params.py

@@ -24,7 +24,25 @@ without a source: both dataclasses refuse to exist without a
   string, and you supply your measured values for field-dependent
   work.
 
-For any other centre (GeV-, PbV-, or your own SnV sample), populate the
+* `gev_bhaskar2017`: the germanium-vacancy (GeV-) centre, measured
+  unstrained orbital splittings of M. K. Bhaskar et al., Phys. Rev.
+  Lett. 118, 223603 (2017) (arXiv:1612.03036): 152 GHz ground,
+  981 GHz excited, 602 nm zero-phonon line. As for SiV-, strain and
+  orbital quenching are sample-specific and deliberately NOT shipped
+  (Siyushev et al., PRB 96, 081201(R) (2017) measured 170 GHz on a
+  strained emitter -- about 20 GHz of transverse-strain contribution
+  -- which is exactly why); supply your measured values, or fit them
+  with `vacspin.lab`.
+
+The lead-vacancy (PbV-) centre is documented but NOT shipped, on
+purpose: its ground-state splitting is measured (about 3900 GHz;
+P. Wang et al., ACS Photonics 8, 2947 (2021), arXiv:2106.03413,
+confirmed by Chen et al., arXiv:2605.27841 (2026)), but no verified
+measured excited-state splitting was available to this package's
+maintainers, and this package does not ship half a parameter set.
+Populate `SpinParameters` yourself when your sample's numbers exist.
+
+For any other centre, or your own sample of any centre, populate the
 dataclasses from your measurements or the literature; the mandatory
 `reference` field keeps the provenance attached to every prediction.
 """
@@ -35,7 +53,7 @@ import dataclasses
 import numpy as np
 
 __all__ = ["SpinParameters", "EmissionBudget", "snv_rosenthal2023",
-           "snv_emission", "siv_hepp2014"]
+           "snv_emission", "siv_hepp2014", "gev_bhaskar2017"]
 
 
 def _check_ref(reference):
@@ -186,3 +204,32 @@ def siv_hepp2014(theta_rad=None, phi_rad=None) -> SpinParameters:
                   "splittings 50/260 GHz, ZPL 737 nm; strain and "
                   "orbital quenching deliberately not shipped (sample-"
                   "specific) -- supply measured values")
+
+
+def gev_bhaskar2017(theta_rad=None, phi_rad=None) -> SpinParameters:
+    """GeV- measured unstrained values of Bhaskar et al., PRL 118,
+    223603 (2017): orbital splittings 152 GHz (ground) and 981 GHz
+    (excited), zero-phonon line 602 nm. Strain and orbital quenching
+    are sample-specific and NOT shipped: this set carries ups = 0 and
+    f = delta = 0, stated here on purpose (Siyushev et al., PRB 96,
+    081201(R) (2017) measured 170 GHz on a strained emitter, about
+    20 GHz of it transverse strain) -- supply your measured values,
+    or fit them from your own frequencies with `vacspin.lab`. Default
+    orientation: the <111> axis of a (001)-oriented sample, exact
+    crystallography, overridable. Millisecond-scale coherence at
+    millikelvin temperatures is measured for this centre (Senkalla et
+    al., PRL 132, 026901 (2024): 24.1 ms CPMG below 300 mK)."""
+    if theta_rad is None:
+        theta_rad = float(np.arccos(1.0 / np.sqrt(3.0)))
+    if phi_rad is None:
+        phi_rad = float(np.pi / 4)
+    return SpinParameters(
+        lam_g=152.0, ups_g=0.0, lam_e=981.0, ups_e=0.0,
+        f_g=0.0, f_e=0.0, delta_g=0.0, delta_e=0.0,
+        theta_rad=theta_rad, phi_rad=phi_rad,
+        reference="Bhaskar et al., PRL 118, 223603 (2017): measured "
+                  "unstrained splittings 152/981 GHz, ZPL 602 nm; "
+                  "strain and orbital quenching deliberately not "
+                  "shipped (sample-specific; cf. Siyushev et al., "
+                  "PRB 96, 081201(R) (2017)) -- supply measured "
+                  "values")

{vacspin-0.2.0 → vacspin-0.3.0}/src/vacspin.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: vacspin
-Version: 0.2.0
+Version: 0.3.0
 Summary: Group-IV colour-centre spin-photon interfaces: strained spin Hamiltonians, cyclicity, transition-resolved Purcell budgets and exact single-shot readout statistics
 Author: Tanvir Mahmud Mahim
 License: Apache-2.0
@@ -128,8 +128,19 @@ No physical number is made up, and none is accepted without a source
   of Hepp et al., PRL 112, 036405 (2014). Strain and orbital
   quenching are sample-specific and deliberately not shipped; the
   reference string says so.
-
-For GeV-, PbV-, or your own sample, populate `SpinParameters` and
+- `gev_bhaskar2017()`: the measured GeV- splittings (152 and
+  981 GHz, zero-phonon line 602 nm) of Bhaskar et al., PRL 118,
+  223603 (2017), with the same deliberate omissions as SiV- --
+  strain and quenching are your sample's numbers, and `vacspin.lab`
+  exists to fit them.
+
+The lead-vacancy (PbV-) centre is documented but not shipped, on
+purpose: its ground splitting is measured (about 3900 GHz; Wang et
+al., ACS Photonics 8, 2947 (2021), confirmed in 2026), but no
+verified measured excited-state splitting was available, and this
+package does not ship half a parameter set.
+
+For PbV-, or your own sample of any centre, populate `SpinParameters` and
 `EmissionBudget` from your measurements or the literature; the
 provenance travels with every prediction.
 
@@ -174,7 +185,7 @@ round trip is exact.
 
 ## How it is checked
 
-56 tests (Python 3.9-3.13, run in CI on every push), every claim
+59 tests (Python 3.9-3.14 in CI on every push), every claim
 anchored to an exact result, a published measurement, or two
 independent code paths -- never a stored number. Highlights: exact
 Kramers doublets and the closed-form splitting against full

{vacspin-0.2.0 → vacspin-0.3.0}/src/vacspin.egg-info/SOURCES.txt

@@ -15,6 +15,7 @@ src/vacspin.egg-info/dependency_links.txt
 src/vacspin.egg-info/requires.txt
 src/vacspin.egg-info/top_level.txt
 tests/test_cavity.py
+tests/test_gev.py
 tests/test_hamiltonian.py
 tests/test_integration.py
 tests/test_lab.py

vacspin-0.3.0/tests/test_gev.py

@@ -0,0 +1,42 @@
+"""GeV- anchors: the shipped values are locked to their source
+(Bhaskar et al., PRL 118, 223603 (2017): 152 GHz ground, 981 GHz
+excited, measured); at zero strain the closed-form zero-field
+splitting IS the spin-orbit constant, exactly; the full solver
+reproduces both splittings; the dipole sum rule holds for this centre
+as an identity; an aligned field keeps the flip channel exactly dark
+(the same protection SnV and SiV show); and the provenance string
+carries its source and its deliberate omissions."""
+import numpy as np
+
+from vacspin import (cyclicity, gev_bhaskar2017, solve,
+                     transition_table, zero_field_splitting)
+
+
+def test_values_locked_to_source():
+    p = gev_bhaskar2017()
+    assert (p.lam_g, p.lam_e) == (152.0, 981.0)
+    assert (p.ups_g, p.ups_e) == (0.0, 0.0)
+    assert (p.f_g, p.f_e, p.delta_g, p.delta_e) == (0.0,) * 4
+    assert "Bhaskar" in p.reference and "223603" in p.reference
+    assert "not" in p.reference        # the omissions are stated too
+
+
+def test_zero_field_splittings_exact():
+    p = gev_bhaskar2017()
+    assert zero_field_splitting(p.lam_g, p.ups_g) == 152.0
+    assert zero_field_splitting(p.lam_e, p.ups_e) == 981.0
+    eg, _, ee, _ = solve(p, [0.0, 0.0, 0.0])
+    assert abs((eg[2] - eg[0]) - 152.0) < 1e-9
+    assert abs((ee[2] - ee[0]) - 981.0) < 1e-9
+
+
+def test_sum_rule_and_aligned_protection():
+    p = gev_bhaskar2017()
+    for b in ([0, 0, 0], [0.03, -0.05, 0.09]):
+        t = transition_table(p, b)
+        assert np.allclose(t["strengths"].sum(axis=1), 3.0, rtol=0,
+                           atol=1e-10)
+    mu = np.array([np.sin(p.theta_rad) * np.cos(p.phi_rad),
+                   np.sin(p.theta_rad) * np.sin(p.phi_rad),
+                   np.cos(p.theta_rad)])
+    assert cyclicity(p, 0.1 * mu) > 1e12