diff --git a/four_levels_model.py b/four_levels_model.py
index 4829e1d..f05c54e 100644
--- a/four_levels_model.py
+++ b/four_levels_model.py
@@ -236,13 +236,13 @@ class HeatingGas:
         self.pah_cross_dc = self.pah_cross_dc[self.energy_range]
         
         ''' yield from anion to neutral '''
-        part = np.where(self.energy_negative_charged >= 0)[0]
+        part = np.where(self.energy_negative_charged <=13.6)[0]
 
         self.detachment_yield = np.full(len(part), 1)
         #the anion being unstable, any energy is enough to detach the electron
 
         ''' yield from neutral to the first photoionization '''
-        first_part = np.where(self.energy_neutral<self.ip_neutral)[0]
+        first_part = np.where(self.energy_neutral<=self.ip_neutral)[0]
 
         second_part = np.where( (self.energy_neutral>=self.ip_neutral) & (self.energy_neutral<(self.ip_neutral+9.2) ) )[0]
         third_part = np.where((self.energy_neutral>self.ip_neutral+9.2))[0]
@@ -367,8 +367,8 @@ class HeatingGas:
         #erg s-1 /(molecule of size n_c)
         
         ''' heating rate of the molecule itself '''
-        heating_pah_a = np.trapz(photo_absorption_a,\
-                                        self.energy) #erg s-1
+        heating_pah_a = np.trapz(photo_absorption_a[energy_range_power_absorbed],\
+                                        self.energy[energy_range_power_absorbed]) #erg s-1
         heating_pah_n = np.trapz(photo_absorption_n[energy_range_power_absorbed],\
                                         self.energy[energy_range_power_absorbed]) #erg s-1
         heating_pah_c = np.trapz(photo_absorption_c[energy_range_power_absorbed],\
--
libgit2 0.21.2