A Theoretical Study on the Phosgenation of 2,4-Toluenediamine (2,4-TDA)

Industrially relevant phosgenation mechanisms of 2,4-toluenediamine (2,4-TDA) were investigated using G3MP2B3 model chemistry. Six reaction pathways had been explored, which resulted in the formation of toluene diisocyanate (2,4-TDI) including different scenarios of the ‘phosgenations first’ and ‘consecutive phosgenations’ mechanisms in both gas and condensed phases. Two possible ‘phosgenations first’ mechanisms show superior to the others in terms of energy, regardless of which phases are considered. Due to the o-dichlorobenzene (ODCB) solvation, the reaction barriers are dramatically reduced compared to the gas-phase reaction mechanism and the solvent effect can be described by linear relationship. Standard enthalpy of formation value was also recommended for 2,4-TDA (59.3 kJ/mol) and 2,4-TDI (−94.1 kJ/mol), as well as for the gas-phase intermediates (IM).


Introduction
Toluene diisocyanate (TDI) is an organic compound, which is a colorless to pale yellowish liquid with a sharp pungent odor used predominantly in the production of polyurethanes (primarily flexible foams) and consumer products, such as coatings, elastomers, adhesives, and sealants [1]. It is also used as curing agent in solid propellant formulations [2]. It is insoluble in water and miscible in most common organic solvents [3]. The most important isomers of TDI are toluene-2,4-diisocyanate (2,4-TDI) and toluene-2,6-diisocyanate (2,6-TDI), commercially produced in different weight ratios. The most common ratio is an 80/20 mixture of 2,4-TDI and 2,6-TDI [3]. The global demand for TDI was 2.37 million tons in 2020, and it is forecasted to increase to about 2.77 million tons by end of 2022 [4].
As shown in Figure 1, the 2,4-TDI is primarily produced by the reaction of toluenediamine (4-methylbenzene-1,3-diamine, abbreviated as 2,4-TDA) with phosgene (COCl 2 ) in a continuous reactor, and the resultant crude 2,4-TDI is purified by distillation. The distillation residue formed is hydrolyzed at <230 • C under absolute pressures of less than 30 bar [5]. For these process, the reactive phosgene compound is produced by passing purified carbon monoxide and chlorine gas through a bed of porous activated carbon [6,7]. The industrial phosgenation of TDA is carried out in a chemically inert orga vent, maintaining a suitable temperature from 40 to 150 °C [8]. The gas-phase ph tion of TDA also can be carried out in the presence of inert gas or vapors of an inert and subsequently heated to a temperature of 300-400 °C [9]. However, the rea amine and phosgene in the liquid phase is very rapid at all industrial temperatu pressures, only a good mixing is desired to suppress the secondary reactions [10]. solvents such as benzene, toluene, xylene, chlorinated aromatic hydrocarbons, an leum hydrocarbons, which are inert with hydrogen chloride, phosgene, and iso can be employed [11]. While using ortho-dichlorobenzene (ODCB) as a solvent, m erating conditions (31.4-100 °C) [12] can be applied, which are mainly due to t boiling point of ODCB (Tb = 180.5 °C). The unreacted phosgene and hydrochloric recycled and utilized in the technology of TDI production [13]. Performing the ph tion reaction in the gas phase is more economical and environmentally benign tha liquid phase due to low residence time, and reduced usage of solvent and ener However, there are side reactions in the gas phase because of the evaporation of leading to the formation of deposits with unresolved structures in the reactors [15] a thorough understanding of the fundamental reactions is necessary. In general, th genation of diamines [15] takes place by two reaction mechanisms: (i) the transfo of both amino groups into carbamoyl chloride followed by the two consecutive H ination (termed as "phosgenations first"), and (ii) one of the amine groups turns cyanate while the other amine group remains intact and the same occurs with th amine group (termed as "stepwise phosgenations").
On one hand, there is an increasing request to understand industrial process molecular level in order to use this information for the further process devel Many theoretical investigations are going on in the field of urethane industry, in phosgene production in presence of [NEt3Me]Cl/[NEt3Me][Cl3] catalyst system [1 mal decomposition of phosgene [17], the formation of other diamine intermediate the phosgenation of MDA [15], the liquid structure of different isocyanates [18,1 thane formation with an excess of isocyanate or alcohol [20], and in presence of c [21], just to name a few. On the other hand, the chemical industries largely use th genation of 2,4-TDA. Its reaction mechanism and the corresponding energy pr less explored. To the best of our knowledge, the thermochemical properties of 2,4their intermediates have not yet been reported. Therefore, the current work aims t tigate the reaction mechanism of the phosgenation of 2,4-TDA using quantum c calculations. The thermodynamic profiles of the phosgenation reaction of 2,4-TD solvent ODCB and their gas-phase reaction were compared. Moreover, the effect perature on the thermodynamic properties of gas-phase and liquid-phase reaction The industrial phosgenation of TDA is carried out in a chemically inert organic solvent, maintaining a suitable temperature from 40 to 150 • C [8]. The gas-phase phosgenation of TDA also can be carried out in the presence of inert gas or vapors of an inert solvent and subsequently heated to a temperature of 300-400 • C [9]. However, the reaction of amine and phosgene in the liquid phase is very rapid at all industrial temperatures and pressures, only a good mixing is desired to suppress the secondary reactions [10]. Organic solvents such as benzene, toluene, xylene, chlorinated aromatic hydrocarbons, and petroleum hydrocarbons, which are inert with hydrogen chloride, phosgene, and isocyanate can be employed [11]. While using ortho-dichlorobenzene (ODCB) as a solvent, mild operating conditions (31.4-100 • C) [12] can be applied, which are mainly due to the high boiling point of ODCB (T b = 180.5 • C). The unreacted phosgene and hydrochloric acid are recycled and utilized in the technology of TDI production [13]. Performing the phosgenation reaction in the gas phase is more economical and environmentally benign than in the liquid phase due to low residence time, and reduced usage of solvent and energy [14]. However, there are side reactions in the gas phase because of the evaporation of amines leading to the formation of deposits with unresolved structures in the reactors [15]. Hence, a thorough understanding of the fundamental reactions is necessary. In general, the phosgenation of diamines [15] takes place by two reaction mechanisms: (i) the transformation of both amino groups into carbamoyl chloride followed by the two consecutive HCl elimination (termed as "phosgenations first"), and (ii) one of the amine groups turns into isocyanate while the other amine group remains intact and the same occurs with the other amine group (termed as "stepwise phosgenations").
On one hand, there is an increasing request to understand industrial processes at the molecular level in order to use this information for the further process development. Many theoretical investigations are going on in the field of urethane industry, including phosgene production in presence of [NEt 3 Me]Cl/[NEt 3 Me][Cl 3 ] catalyst system [16], thermal decomposition of phosgene [17], the formation of other diamine intermediate (MDA), the phosgenation of MDA [15], the liquid structure of different isocyanates [18,19], urethane formation with an excess of isocyanate or alcohol [20], and in presence of catalysts [21], just to name a few. On the other hand, the chemical industries largely use the phosgenation of 2,4-TDA. Its reaction mechanism and the corresponding energy profile are less explored. To the best of our knowledge, the thermochemical properties of 2,4-TDI and their intermediates have not yet been reported. Therefore, the current work aims to investigate the reaction mechanism of the phosgenation of 2,4-TDA using quantum chemical calculations. The thermodynamic profiles of the phosgenation reaction of 2,4-TDA in the solvent ODCB and their gas-phase reaction were compared. Moreover, the effect of temperature on the thermodynamic properties of gas-phase and liquid-phase reactions for the formation of 2,4-TDI was studied. [22] is considered a reasonable compromise between the accuracy and the computational costs of the system of interest; therefore, G3MP2B3 was applied to compute thermochemical properties of the species involved in the title reaction, as implemented in the Gaussian09 program package [23]. As part of this protocol, geometry optimizations were carried out using the B3LYP/6-31G(d) level of theory [24] by applying the "tight" convergence criterion. Normal mode analysis was performed on the optimized structures at the same level of theory to characterize the identities of the computed species on the potential energy surface (PES). Furthermore, the transition-state (TS) structures also were confirmed to locate first-order saddle points on the PES by visual inspection of the intramolecular motions corresponding to the imaginary wavenumber using GaussView 6 [25]. Intrinsic reaction coordinate (IRC) calculations [26] were carried out to map the minimal energy pathways (MEP). As a first proxy, the zero-point corrected relative B3LYP/6-31G(d) energy (∆E 0,B3LYP/6-31G(d) ) had been calculated. To refine the electronic energy, additional single point calculations of the critical points of the PES were carried out using QCISD(T)/6-31G(d) (including MP2/6-31G(d) level of theory) and MP2/GTMP2 levels of theories based on B3LYP/6-31G(d) geometries according to the G3MP2B3 composite method. Harmonic wavenumbers obtained at B3LYP/6-31G(d) level of theory were scaled by a factor of 0.96 [22] for refining the accuracy of the thermodynamic properties, such as zero-point corrected relative energy (∆E 0 , G3MP2B3 ), relative enthalpy (∆H • G3MP2B3 ), the standard enthalpy of formation (∆ f,298.15K H • (g)), relative Gibbs free energy (∆G • G3MP2B3 ), and entropy (S • ) under standard condition (T = 298.15K and P = 1 atm). The ∆ f,298.15K H • (g) were obtained by group additivity (GA) role using online NIST tool [27], as well as by atomization scheme (AS) using the standard enthalpy of formation for atoms from CCCBDB database [28].

G3MP2B3 model chemistry
The integral equation formalism of the polarizable continuum model (PCM) with the radii and non-electrostatic terms of Truhlar and co-workers (SMD) [29] was used to mimic implicitly the surrounding ortho-dichloro-benzene (ODCB) environment of the phosgenation. The SMD model had been applied for the re-optimization of all the gasphase structures meaning the SMD treatment of ODCB was fully integrated into the G3MP2B3 protocol.

Gas-Phase Reaction Mechanism
Due to the asymmetric arrangement of the amine group in 2,4-TDA, six different phosgenation reaction pathways can be considered and they can be classified into two types, such as "stepwise phosgenations" and "phosgenations first", as shown in Figure 2. The "phosgenations first" mechanism means that both amines transform to carbamic chloride groups in a consecutive manner and then HCl-elimination takes place to give the isocyanate groups. In the "stepwise phosgenations" mechanism, one of the amine groups turns into isocyanate while the other amine group remains intact, and the same occurs with the other amine group. Each type has different reaction pathways, "stepwise phosgenations" have two reaction pathways (noted as Pathway 1 and Pathway 2 in Figures 3 and 4), and "phosgenations first" have four reaction pathways (noted as Pathway 3-6 in Figures 3 and 4).
In the first elementary step, one of the amine groups turns into carbamic chloride via a four-membered transition state, in which synchronous phosgene addition and HCl elimination take place. If the reaction center is close to the methyl group, that is TSa1, the reaction barrier is roughly 5 kJ/mol higher (∆E 0 = 44.2 kJ/mol) than in the case of amine in position 3 (TSb1, ∆E 0 = 39.4 kJ/mol), as shown in Figure 3, for the gas-phase reaction (as discussed later, the potential energy surface for the ODCB phase reaction shows the same trends shown in Figure 4), and in Table 1. This is also manifested in the transition-state structure in such a way that the critical bond distance for the C-Cl bond being broken is significantly larger in the case of Tsa1 (2.602 Å) than for TSb1 (2.512 Å) and the released chlorine approaches the hydrogen amine closer (2.007 Å vs. 2.069 Å), as shown in Figures 5 and 6. Both transition states correlate with exothermic products, namely (5-amino-2-methylphenyl)carbamic chloride (IMa1, ∆ r H 0 = −65.5 kJ/mol) and (3-amino-4-methylphenyl)carbamic chloride (Imb1, ∆ r H 0 = −67.8 kJ/mol).      In the first elementary step, one of the amine groups turns into carbamic chloride via a four-membered transition state, in which synchronous phosgene addition and HCl elimination take place. If the reaction center is close to the methyl group, that is TSa1, the reaction barrier is roughly 5 kJ/mol higher (ΔE0 = 44.2 kJ/mol) than in the case of amine in position 3 (TSb1, ΔE0 = 39.4 kJ/mol), as shown in Figure 3, for the gas-phase reaction (as discussed later, the potential energy surface for the ODCB phase reaction shows the same trends shown in Figure 4), and in Table 1. This is also manifested in the transition-state structure in such a way that the critical bond distance for the C-Cl bond being broken is significantly larger in the case of Tsa1 (2.602 Å ) than for TSb1 (2.512 Å ) and the released chlorine approaches the hydrogen amine closer (2.007 Å vs. 2.069 Å ), as shown in Figures 5 and 6. Both transition states correlate with exothermic products, namely (5-amino-2methylphenyl)carbamic chloride (IMa1, ΔrH 0 = −65.5 kJ/mol) and (3-amino-4methylphenyl)carbamic chloride (Imb1, ΔrH 0 = −67.8 kJ/mol). Table 1. G3MP2B3 thermochemical properties calculated in the gas phase and ortho-dichlorobenzene (ODCB), including zero-point corrected relative energies (ΔE0,G3MP2B3), relative enthalpies (ΔH 0 G3MP2B3), and relative Gibbs free energies (ΔG 0 G3MP2B3).        In the case of "stepwise phosgenations", the carbamic chloride intermediate IMa1 or Imb1 undergoes HCl-elimination resulting in the formation of the first isocyanato group, namely 3-isocyanato-4-methylaniline (IMa2) and 5-isocyanato-methylaniline (Imb2), respectively. Their relative energies (IMa2 (∆E 0 = −51.3 kJ/mol), Imb2 (∆E 0 = −51.8 kJ/mol)) are higher than the intermediate of the previous step (IMa1 (∆E 0 = −69.2 kJ/mol), Imb1 (∆E 0 = −71.1 kJ/mol)). Therefore, the corresponding transition states (Tsa2 and TSb2) lie high (∆E 0 = 43.6 kJ/mol for Tsa2 and ∆E 0 = 44.6 kJ/mol for TSb2), which can be understood by looking at the drastically elongated C-Cl bond in the transition-state structures in Figure 5 (2.753 Å for Tsa2) and Figure 6 (2.750 Å for TSb2). This HCl elimination step is the most energy-consuming one in the "stepwise phosgenations" mechanism. Considering the further step either from IMa2 (→ Tsa3 → IMa3 → 2,4-TDI) or Imb2 (→ TSb3 → Imb3 → 2,4-TDI), it turned out that the transformation of the second amine group into isocyanate showed similar energy profile as the first one. In the case of "stepwise phosgenations", the carbamic chloride intermediate IMa1 or Imb1 undergoes HCl-elimination resulting in the formation of the first isocyanato group, namely 3-isocyanato-4-methylaniline (IMa2) and 5-isocyanato-methylaniline (Imb2), respectively. Their relative energies (IMa2 (ΔE0 = −51.3 kJ/mol), Imb2 (ΔE0 = −51.8 kJ/mol)) are higher than the intermediate of the previous step (IMa1 (ΔE0 = −69.2 kJ/mol), Imb1 (ΔE0 = −71.1 kJ/mol)). Therefore, the corresponding transition states (Tsa2 and TSb2) lie high (ΔE0 = 43.6 kJ/mol for Tsa2 and ΔE0 = 44.6 kJ/mol for TSb2), which can be understood by looking at the drastically elongated C-Cl bond in the transition-state structures in Figure  5 (2.753 Å for Tsa2) and Figure 6 (2.750 Å for TSb2). This HCl elimination step is the most energy-consuming one in the "stepwise phosgenations" mechanism. Considering the further step either from IMa2 (→ Tsa3 → IMa3 → 2,4-TDI) or Imb2 (→ TSb3 → Imb3 → 2,4-TDI), it turned out that the transformation of the second amine group into isocyanate showed similar energy profile as the first one. For instance, IMa2 reacts with excess phosgene and forms carbamic chloride compound (IMa3, 3-isocyanato-4-methylphenyl)carbamic chloride) and the exoergicity of this elementary step is −65.6 kJ/mol, which is only slightly different from that of IMa1 (∆E 0 = −69.2 kJ/mol). However, the Imb3 pathway is less exoergic by 7.4 kJ/mol (∆E 0 = −63.7 kJ/mol vs. ∆E 0 = −71.1 kJ/mol for Imb1). Therefore, one can state that the reactivity of these amine groups only slightly interferes, which is also confirmed by the high similarity between corresponding transition states and intermediate structures (therefore these structures are not discussed here thoroughly). Regardless of the small deviations, the energy pattern for an isocyanate group formation with the "stepwise phosgenations" mechanism can be described as follows: the first step is exothermic with a transition state lying below the transition state of the endothermic second step. Due to the high exothermicity of the synchronous phosgene addition and HCl-elimination step the overall reaction is strongly exothermic with the formation of the TDI product and four HCl molecules (∆H 0 = −81.4 kJ/mol). Regardless of IMa1 or Imb1 formed in the first reaction step of the "phosgenations first", mechanism can take place via Tsa5 or TSb5 transition state, respectively. The resulted intermediate is the same, that is 4-methyl-1,3-phenylene dicarbamic chloride (IM4) and HCl. The relative energies of Tsa5 and TSb5 are close to each other (−11.2 kJ/mol for Tsa5 and −12.4 kJ/mol for TSb5), and both are below the entrance level even in the case of gas-phase reactions. They are also similar in structure to each other, as well as Tsa1, Tsa3, TSb1, and TSb3. The largest variation in the parameters can be found in the case of the C-Cl distance (a bond being broken), which is the smallest for TSb3 (2.488 Å) and the largest for Tsa1 (2.602 Å).
IM4 with the two carbamic chloride groups can turn into either IMa3 or Imb3, regardless of forming through either Tsa2 or TSb2, making four possible pathways for the "phosgenations first" mechanism. These HCl-eliminations (Tsa6 and TSb6) are structurally similar to the transition states of final HCl-eliminations (Tsa4 and TSb4); therefore, they are all discussed at once. The C-Cl distance is elongated (2.723-2.753 Å) significantly, which means that the C-Cl is already partially being broken and this liberated chlorine atom synchronously approaches the amine hydrogen in the distance of 1.751-1.774 Å in a planar four-centered transition state. These HCl-elimination reactions are all endoergic by 17.8 kJ/mol and 19.2 kJ/mol.
Comparing the two types of reaction mechanisms from an enthalpic perspective, "stepwise phosgenations" have an exothermic and then an endothermic step in consecutive manner, whereas the "phosgenations first" type has two exothermic steps at the beginning and two endothermic steps at the end.

Thermochemical Properties of the Intermediates
The thermodynamic properties also can be calculated for all the intermediates, reactant and product in the gas phase at thermodynamic standard conditions and all of them are tabulated in Table 2. Table 2. Gas-phase thermochemical properties for intermediates of the phosgenation of the 2,4-TDA. The standard enthalpy of formation is calculated from G3MP2B3 enthalpies through group additivity (GA) role obtained by the means of online NIST tool [27] and through atomization scheme (AS). Standard molar enthalpy of formation (∆ f,298.15K H 0 (g)), standard molar entropy (S 0 (g)), and molar heat capacity (C v (g)) were obtained at 1 atm pressure at 298.15 K. Entropy and heat capacity values are obtained from scaled B3LYP/6-31G(d) frequencies. Standard molar enthalpy of formation (∆ f,298.15K H 0 (g)) was directly derived from our G3MP2B3 results for the atomization scheme (AS), calculated by using the accurate atomization enthalpy values (∆ atom H 0 ( 3 C) = 716.68 ± 0.45 kJ/mol, ∆ atom H 0 ( 2 H) = 218.00 ± 0.01 kJ/mol, ∆ atom H 0 ( 4 N) = 472.68 ± 0.40 kJ/mol, ∆ atom H 0 ( 3 O) = 249.23 ± 0.00 kJ/mol, ∆ atom H 0 ( 2 Cl)= 121.31 ± 0.01 kJ/mol) from CCBDB database and the reaction enthalpy of the following hypothetical reaction with the stochiometric numbers x, y, z, q and w:

Species
x 3 C + y 2 H + z 4 N + q 3 O + w 2 Cl = 1 C x H y N z O q Cl w However, the NIST group additivity tool (GA) was used to generate values for TDA and TDI. In the GA model, the group increment for methyl group was −42.7 kJ/mol, while that of the phenyl ring is 0 kJ/mol. The enthalpy contribution from the amine group bonded to aromatic carbon (C a-ring ) is 20 kJ/mol. Finally, enthalpy increments of C a-ring -C, C a-ring -H and C a-ring -N were 23 kJ/mol, 14 kJ/mol and −2 kJ/mol, respectively. All these additivity values were obtained from the NIST tool [27]. In case of TDI, the aromatic isocyanate group increment (−61.2 kJ/mol) was adopted from our previous work [15,30]. In both cases, the AS and GA values are within 3 kJ/mol, which also suggest that the energetics of the title reaction obtained with G3MP2B3 are consistent with the literature. Surprisingly, there is no reported standard enthalpy of formation for any of the compounds listed in Table 2 to the best of our knowledge. As also can be seen in Table 2, only TDA has endothermic formation and formation of IMa2 and IMb2 have slightly exothermic while the remaining intermediates (IMa1, IMb1, IMa3, IMb1, IM4) and TDI are all highly exothermic product of the phosgenation. However, the overall reaction of the 2,4-TDI formation is exothermic, 2,4-TDI is reactive enough to form highly thermally and mechanically stable frameworks, such as isocyanate dimers and isocyanurates [31], but such condensed phase reactions is out of the scope of the current study. For completeness, the entropy content (S • (g)) and heat capacity values (C v (g)) are also tabulated in Table 2 for all species studied.

Effect of Solvent and Temperature
As shown earlier, the first activation barrier of gas-phase phosgenation (both TSa1 and TSb1) lies relatively high, and this is a crucial step in this industrially important system. To reduce the activation energy, ODCB is used in the chemical industry as an inert, but polar solvent (the permanent electric dipole moment of ODCB is 2.27 D). Furthermore, to get closer to the industrial relevant condition, two different temperature effects were studied (273.15 K, and 423.15 K) and their thermochemical properties are shown in Table 1. Due to the polar environment provided by ODCB, the relative energies of the transition states and intermediates are all drastically reduced. Relative energies of the first transition states in both directions become submerged (∆E 0 = −0.2 kJ/mol for TSa1 and ∆E 0 = −3.7 kJ/mol for TSb1) and are kinetically easily accessible.
Comparing Figures 3 and 4, the potential energy surfaces in gas phase and ODCB remained similar, therefore the relative energies of gas phase against those in ODCB are plotted in Figure 7, and they show a very strong linear relationship with the R 2 value of 0.9634, and 0.9995 for the transition-state (TS) and the intermediate (IM) structures, respectively. Analyzing the fitted expression for TS relative energies, the gas-phase relative energies are all multiplied by a factor of 1.3. Then, this effect is overcompensated for the positive relative energy state by a general strong shift (−53.7 kJ/mol) toward smaller energies as the intercept indicates. Since all intermediates are below the energy level of the reactants, the 1.18 slope value resulted in a significant decrease in their energies.
The relative enthalpies are not significantly affected by changing the temperature from 298.15 K to 423.15 K. The largest absolute deviation was 2.6 kJ/mol, regardless whether gas-phase or ODCB phase values were considered.

Comparison of the Relative Energies from G3MP2B3 and B3LYP/6-31G(d)
Phosgenation of amines can be studied for larger analogous systems, although the G3MP2B3 computation is not suitable for a larger system. To provide a trade-off alternative according to accuracy and computation time, we studied the correlation between the G3MP2B3 and B3LYP/6-31G(d) relative energies. As shown in Figure 8 for gas-phase energies and in Figure 9 for the ODCB phase, G3MP2B3 and B3LYP/6-31G(d) energies are in a strong linear relationship for transition states (TS), and intermediates (IM) with the R 2 value of 0.9177, and 0.9953, respectively, for gas phase and 0.9786, and 0.9945, respectively, for ODCB. For both surrounding media, slope values for TS and IM agree with the uncertainty of the fit, but they differ according to gas-phase (roughly 0.88) or ODCB (roughly 0.95) results. The magnitude of the intercepts provides a relatively small shift in the relative-energy values, from 0.02 kJ/mol to 6.2 kJ/mol. Their deviations for IM and TS are close to or somewhat above their uncertainties.  The relative enthalpies are not significantly affected by changing the temperature from 298.15 K to 423.15 K. The largest absolute deviation was 2.6 kJ/mol, regardless whether gas-phase or ODCB phase values were considered.

Comparison of the Relative Energies from G3MP2B3 and B3LYP/6-31G(d)
Phosgenation of amines can be studied for larger analogous systems, although the G3MP2B3 computation is not suitable for a larger system. To provide a trade-off alternative according to accuracy and computation time, we studied the correlation between the G3MP2B3 and B3LYP/6-31G(d) relative energies. As shown in Figure 8 for gas-phase energies and in Figure 9 for the ODCB phase, G3MP2B3 and B3LYP/6-31G(d) energies are in a strong linear relationship for transition states (TS), and intermediates (IM) with the R 2 value of 0.9177, and 0.9953, respectively, for gas phase and 0.9786, and 0.9945, respectively, for ODCB. For both surrounding media, slope values for TS and IM agree with the uncertainty of the fit, but they differ according to gas-phase (roughly 0.88) or ODCB (roughly 0.95) results. The magnitude of the intercepts provides a relatively small shift in the relative-energy values, from 0.02 kJ/mol to 6.2 kJ/mol. Their deviations for IM and TS are close to or somewhat above their uncertainties.

Thermochemistry of the Studied Reaction at Different Conditions
The phosgenation of 2,4-TDA in the gas-phase and ODCB environments was studied at two different temperatures, 298.15 K and 423.15 K, to explore the temperature dependence on the thermodynamics of the studied reaction mechanism. Relative enthalpies and Gibbs free energies obtained at the mentioned temperature and 1 atm pressure are given in Table 1.
No significant change in the relative enthalpies (∆H G3MP2B3 (T)) (0.0 and 2.6 kJ/mol for gas phase and −0.2 and 2.6 kJ/mol for ODCB) was observed in both reaction environments while increasing the temperature from 298.15 K to 423.15 K. However, the relative Gibbs free energies (∆G 0 ) were substantially increased from −33.5 to 22.9 kJ/mol for gas phase and −32.3 to 24.2 kJ/mol for ODCB.

Conclusions
In this work, the energy profile of the six different possible reaction mechanisms for the 2,4-TDA phosgenation were discussed based on the results obtained by G3MP2B3 composite method and B3LYP/6-31G(d) level of theory, which can be classified as follows: (a) "Phosgenations first": first amine groups turn into carbamic chloride groups and then the isocyanate groups form; (b) "Stepwise phosgenations": amine group presented in ortho, or para, position converts as a carbamic chloride and then it loses HCl and forms isocyanate. The same step repeats for the other position and forms another isocyanate group.
Our key findings of these reaction mechanisms are the followings: 1.
The addition of the first COCl 2 is the rate limiting step for the reaction, regardless of which mechanism is taking place; 2.
The activation barriers for the gas-phase reaction are relatively high and the ODCB has reduced the activation barriers for all the reaction pathways; 3.
For transition states and intermediates, the gas phase and ODCB phase relative energies show strong linear correlation, and the presence of ODCB medium causes drastic shrinking in energy level of the transition state as one compares them with the gas-phase values. For the intermediates, the relative energies lowered by a factor of 1.18; 4.
Comparing the energy profile of the "phosgenations first" and "stepwise phosgenations", the former has two exothermic consecutive elementary steps and then two consecutive endothermic ones, while the latter has consecutive exothermic and endothermic steps. Therefore, the "phosgenation first" seems to be the dominant channel due to both thermodynamic and kinetic points of view; 5.
Relative energies obtained by G3MP2B3 and B3LYP/6-31G(d) has linear relation in both the gas phase and ODCB for transition states and intermediates. Therefore, B3LYP/6-31G(d) computation can be an excellent compromise between accuracy and computation time for phosgenation reactions. Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the policy of the University of Miskolc.