preparations for implementing outer mass transfer

Author: Silvia Toonen

interactions.py

@@ -869,6 +869,18 @@ def triple_stable_mass_transfer(bs, donor, accretor, self):
         bs.bin_type = bin_type['stable_mass_transfer']                
 
     #implementation is missing
+
+def triple_common_envelope_phase(bs, donor, accretor, self):
+    # mass transfer of both inner and outer orbit is not yet considered here
+    
+    # orbital evolution is being taken into account in secular_code        
+    if REPORT_FUNCTION_NAMES:
+        print('Triple common envelope')
+
+    bs.bin_type = bin_type['common_envelope']                
+    self.save_snapshot()       
+    
+    #implementation is missing
 
 #when the tertiary star transfers mass to the inner binary
 def outer_mass_transfer(bs, donor, accretor, self):
@@ -884,13 +896,11 @@ def outer_mass_transfer(bs, donor, accretor, self):
         # possible the outer binary needs part_dt_mt as well. 
         #adjusting triple is done in secular evolution code
     else:        
-        if REPORT_FUNCTION_NAMES:
-            print('outer_mass_transfer: unstable mass transfer in outer binary')
-        bs.bin_type = bin_type['common_envelope_energy_balance']                
-        self.save_snapshot()        
-        
-        #implementation is missing
-        #snapshot
+        triple_common_envelope_phase(bs, donor, accretor, self)
+
+
+    #stopping condition 0:False, 1:True, -1: calculate through outer mass transfer - effect on inner & outer orbit is taken care off here. 
+    return -1 
 
 #-------------------------
 
@@ -1044,7 +1054,7 @@ def perform_stellar_interaction(bs, self):
 
             if bs.child1.is_donor:
                 if bs.child2.child1.is_donor or bs.child2.child2.is_donor:
-                    outer_mass_transfer(bs, bs.child1, bs.child2, self)
+                    stopping_condition = outer_mass_transfer(bs, bs.child1, bs.child2, self)
             else:
                 detached(bs, self)
 

triple_class.py

@@ -863,10 +863,6 @@ def check_RLOF(self):
             if star.radius >= Rl3 - (1.0 * small_numerical_error|units.RSun):
                 star.is_donor = True
 
-            if star.is_donor and (bin.child1.is_donor or bin.child2.is_donor):
-                print(Rl1, bin.child1.radius, Rl2, bin.child2.radius)
-                print(Rl3, star.radius)
-                sys.exit('RLOF in inner and outer binary')
 
         else:
             sys.exit('check_RLOF: structure stellar system unknown')
@@ -1842,29 +1838,33 @@ def resolve_stellar_interaction(self, stellar_system = None):
             return
         elif self.is_binary(stellar_system):
             if REPORT_TRIPLE_EVOLUTION:
-                print('\n perform stellar interaction: binary - double star')
+                print('\n perform stellar interaction: binary')
 #            stellar_system = perform_stellar_interaction(stellar_system, self)
             stopping_condition = perform_stellar_interaction(stellar_system, self)
             return stopping_condition #stellar interaction
         else:
-            if not stellar_system.child1.is_star: #child1 is a multiple
-                stopping_condition = self.resolve_stellar_interaction(stellar_system.child1)  
-                if not stopping_condition: #stellar interaction
-                    return False                                     
-            elif not stellar_system.child2.is_star: #child2 is a multiple
-                stopping_condition = self.resolve_stellar_interaction(stellar_system.child2)        
-                if not stopping_condition: #stellar interaction
-                    return False                                     
-            else:
-                sys.exit('resolve_stellar_interaction: structure stellar system unknown, both children are binaries')
-            
             if REPORT_TRIPLE_EVOLUTION:
-                print('\n perform stellar interaction: binary')
-#            stellar_system = perform_stellar_interaction(stellar_system, self)            
+                print('\n perform stellar interaction')
+      
             stopping_condition = perform_stellar_interaction(stellar_system, self)            
             if not stopping_condition: #stellar interaction
                 return False                                     
+
+            if stopping_conditions>0:
+                #in case of outer mass transfer, skip inner evolution, is taken care of in outer mass transfer
+                if not stellar_system.child1.is_star: #child1 is a multiple
+                    stopping_condition = self.resolve_stellar_interaction(stellar_system.child1)  
+                    if not stopping_condition: #stellar interaction
+                        return False                                     
+                elif not stellar_system.child2.is_star: #child2 is a multiple
+                    stopping_condition = self.resolve_stellar_interaction(stellar_system.child2)        
+                    if not stopping_condition: #stellar interaction
+                        return False                                     
+                else:
+                    sys.exit('resolve_stellar_interaction: structure stellar system unknown, both children are binaries')
+
             return True                  
+            
 
     def determine_mass_transfer_timescale(self, stellar_system = None):