pvlib · kandersolar · Feb 8, 2024 · Dec 22, 2023 · Dec 22, 2023 · Dec 22, 2023
diff --git a/docs/sphinx/source/whatsnew/v0.10.4.rst b/docs/sphinx/source/whatsnew/v0.10.4.rst
@@ -7,7 +7,7 @@ v0.10.4 (Anticipated March, 2024)
 
 Enhancements
 ~~~~~~~~~~~~
-* Added the Huld PV model (:pull:`1940`)
+* Added the Huld PV model used by PVGIS (:pull:`1940`)
 
 
 Bug fixes
@@ -30,3 +30,4 @@ Requirements
 Contributors
 ~~~~~~~~~~~~
 * Cliff Hansen (:ghuser:`cwhanse`)
+* :ghuser:`matsuobasho`
diff --git a/pvlib/pvarray.py b/pvlib/pvarray.py
@@ -296,6 +296,10 @@ def huld(effective_irradiance, temp_mod, pdc0, k=None, cell_type=None):
         P_{dc} = G' P_{dc0} (1 + k'_1 \log(G') + k'_2 \log^2 (G') + k'_3 T' +
                  k'_4 T' \log(G') + k'_5 T' \log^2 (G') + k'_6 T'^2)
 
+
+    Users should be aware that at very low irradiance, i.e.,
+    :math:`G' < 20 W/m^2`, :math:`P_{dc}` may be negative.
-    :math:`G' < 20 W/m^2`, :math:`P_{dc}` may be negative.
+    :math:`G' < 20 W/m^2`, :math:`P_{dc}` may be negative.
+    This is caused by the :math:`\log(G')` term, which tends toward 
+    negative infinity as irradiance approaches zero.
-    :math:`G' < 20 W/m^2`, :math:`P_{dc}` may be negative.
+    :math:`G' < 20 W/m^2`, :math:`P_{dc}` may be negative.
+    This is caused by the :math:`\log(G')` term, which tends toward 
+    negative infinity as irradiance approaches zero.
+
     PVGIS documentation shows a table of default parameters :math:`k'` for
     different cell types. The parameters :math:`k'` differ from the parameters
     :math:`k` expected by :py:func:`huld` by the factor ``pdc0``, that is,