diff --git a/.gitignore b/.gitignore
index fa9452f6..213f8a44 100644
--- a/.gitignore
+++ b/.gitignore
@@ -123,3 +123,4 @@ tests/TEMP_test_addModule_materials.py
 /tutorials/04_scenario/pvd_job_two-material-demo
 /tutorials/04_scenario/pvd_job_stack
 /tutorials/04_scenario/pvd_job_perov-module-stack
+/_si_imgs
diff --git a/pvdeg/data/DegradationDatabase.json b/pvdeg/data/DegradationDatabase.json
index 3e1eaf00..015f5162 100644
--- a/pvdeg/data/DegradationDatabase.json
+++ b/pvdeg/data/DegradationDatabase.json
@@ -1186,17 +1186,37 @@
         "SourceTitle": "Accelerated aging of all-inorganic, interface-stabilized perovskite solar cells",
         "Authors": "Xiaoming Zhao, Tianran Liu, Quinn C. Burlingame, Tianjun Liu, Rudolph Holley III, Guangming Cheng, Nan Yao, Feng Gao, Yueh-Lin Loo",
         "Reference": "Science 377, 307-310 (2022)",
-        "Material": "Halide Perovskite (CsPbI3, uncapped device)",
+        "Material": "Halide Perovskite (CsPbI3, uncapped device, ISOS-L-3)",
         "EquationType": "Perovskite_biexponential",
-        "Comments": "Parameters for the uncapped CsPbI3 device from Zhao et al. 2022. Ea values read directly from the Arrhenius plot (Fig. 3C, uncapped fast/slow dashed lines). k0 values are calibrated to reproduce T90,Agg ~ 1440 h under ISOS-L2 (85 C, 1000 W/m^2) as reported by Orooji et al. 2026 (doi:10.1039/d6el00021e) using the same Zhao parameters. A1, A2, B are estimated from the biexponential shape in Fig. 3A (uncapped). Exact fitted coefficients are in the Zhao supplementary (doi:10.1126/science.abn5679); users should re-fit k0, A1, A2, B to their own device data for quantitative predictions.",
-        "Ea_fast": {"value": 0.248, "units": "eV"},
-        "Ea_slow": {"value": 0.243, "units": "eV"},
-        "k0_fast": {"value": 0.56, "units": "h^-1 (pre-exponential at reference irradiance I_ref)"},
-        "k0_slow": {"value": 0.48, "units": "h^-1 (pre-exponential at reference irradiance I_ref)"},
-        "A1": {"value": 0.25, "units": "dimensionless (fast-process amplitude, A1+A2+B=1)"},
-        "A2": {"value": 0.70, "units": "dimensionless (slow-process amplitude)"},
-        "B": {"value": 0.05, "units": "dimensionless (residual stable fraction)"},
-        "gamma": {"value": 1.0, "units": "dimensionless (light-intensity exponent)"},
-        "I_ref": {"value": 1200.0, "units": "W/m^2 (reference irradiance used in Zhao experiments)"}
+            "Comments": "Uncapped CsPbI3 device from Zhao et al. 2022 under ISOS-L-3-style light-soaking conditions. Ea_fast and Ea_slow are taken directly from Fig. 3C; k0_fast, k0_slow, A1, A2, and B were obtained by AI-assisted digitization of Fig. S8 followed by a global biexponential fit with Ea fixed to the paper values. The paper does not tabulate numerical k0 or amplitude coefficients. Amplitude terms are only moderately constrained because the two activation energies are similar.",
+        "Ea_fast": {"value": 0.248, "units": "eV", "source": "Zhao 2022 Fig. 3C"},
+        "Ea_slow": {"value": 0.243, "units": "eV", "source": "Zhao 2022 Fig. 3C"},
+        "k0_fast": {"value": 4.32, "units": "h^-1 (pre-exponential at reference irradiance I_ref)", "source": "AI-digitized Zhao 2022 Fig. S8 + global biexp fit (Ea fixed to paper)"},
+        "k0_slow": {"value": 0.243, "units": "h^-1 (pre-exponential at reference irradiance I_ref)", "source": "AI-digitized Zhao 2022 Fig. S8 + global biexp fit (Ea fixed to paper)"},
+        "A1": {"value": 0.222, "units": "dimensionless (fast-process amplitude, A1+A2+B=1)", "source": "AI-digitized Zhao 2022 Fig. S8 global biexp fit"},
+        "A2": {"value": 0.716, "units": "dimensionless (slow-process amplitude)", "source": "AI-digitized Zhao 2022 Fig. S8 global biexp fit"},
+        "B": {"value": 0.062, "units": "dimensionless (residual stable fraction)", "source": "AI-digitized Zhao 2022 Fig. S8 global biexp fit"},
+        "gamma": {"value": 1.0, "units": "dimensionless (light-intensity exponent)", "source": "Orooji et al. 2026 (doi:10.1039/d6el00021e); Zhao paper does not specify"},
+        "I_ref": {"value": 1200.0, "units": "W/m^2 (reference irradiance used in Zhao experiments)", "source": "Zhao 2022 Methods"}
+    },
+    "D047": {
+        "DataEntryPerson": "Rajiv (Dax) Daxini",
+        "DateEntered": "5/8/2026",
+        "doi": "10.1126/science.abn5679",
+        "SourceTitle": "Accelerated aging of all-inorganic, interface-stabilized perovskite solar cells",
+        "Authors": "Xiaoming Zhao, Tianran Liu, Quinn C. Burlingame, Tianjun Liu, Rudolph Holley III, Guangming Cheng, Nan Yao, Feng Gao, Yueh-Lin Loo",
+        "Reference": "Science 377, 307-310 (2022)",
+        "Material": "Halide Perovskite (CsPbI3 with 2D Cs2PbI2Cl2 capping layer, ISOS-L-3)",
+        "EquationType": "Perovskite_biexponential",
+            "Comments": "Capped CsPbI3 device with a 2D Cs2PbI2Cl2 interface stabilizer from Zhao et al. 2022. Ea_fast and Ea_slow are taken directly from Fig. 3C; k0_fast, k0_slow, A1, A2, and B were obtained by AI-assisted digitization of Fig. S8 followed by a global biexponential fit with Ea fixed to the paper values. The paper does not tabulate numerical k0 or amplitude coefficients. The capped data is effectively single-pathway over the measured window, so A2 is near zero and k0_slow is weakly constrained.",
+        "Ea_fast": {"value": 0.487, "units": "eV", "source": "Zhao 2022 Fig. 3C"},
+        "Ea_slow": {"value": 0.432, "units": "eV", "source": "Zhao 2022 Fig. 3C"},
+        "k0_fast": {"value": 1712.0, "units": "h^-1 (pre-exponential at reference irradiance I_ref)", "source": "AI-digitized Zhao 2022 Fig. S8 (capped) + global biexp fit (Ea fixed to paper)"},
+        "k0_slow": {"value": 1.0e-06, "units": "h^-1 (placeholder; A2~0 so k0_slow is unconstrained)", "source": "AI-digitized Zhao 2022 Fig. S8 (capped) global biexp fit"},
+        "A1": {"value": 0.225, "units": "dimensionless (fast-process amplitude, A1+A2+B=1)", "source": "AI-digitized Zhao 2022 Fig. S8 (capped) global biexp fit"},
+        "A2": {"value": 0.0, "units": "dimensionless (slow-process amplitude; statistically zero in fit)", "source": "AI-digitized Zhao 2022 Fig. S8 (capped) global biexp fit"},
+        "B": {"value": 0.775, "units": "dimensionless (residual stable fraction)", "source": "AI-digitized Zhao 2022 Fig. S8 (capped) global biexp fit"},
+        "gamma": {"value": 1.0, "units": "dimensionless (light-intensity exponent)", "source": "Orooji et al. 2026 (doi:10.1039/d6el00021e); Zhao paper does not specify"},
+        "I_ref": {"value": 1200.0, "units": "W/m^2 (reference irradiance used in Zhao experiments)", "source": "Zhao 2022 Methods"}
     }
 }