The path towards type V deep eutectic solvents: inductive effect and steric hindrance in the system tert-butanol + perfluoro tert-butanol

The solid-liquid phase behaviour of two tertiary alcohols, perfluoro-tert -butanol and tert - butanol, was here studied using experimental (ITC, DSC and density measurements) and theoretical (MD simulations) approaches. The phase diagram of the binary mixture reveals highly negative deviations to ideality at low concentrations, as well as the formation of co-crystals and is characterized by two eutectic points, and a congruent melting point and a peritectic reaction corresponding to TBH:TBF stoichiometries of 2:1 and 1:1 respectively. Measured excess molar enthalpies and volumes were calculated, showing negative and positive deviations from ideality, respectively. The effect of acidity, stereochemical hindrance and phobic effects and how they affect intermolecular interactions in these binary mixtures is discussed, with the aim of designing and fine-tuning type V deep eutectic solvents. The


Introduction
2][3] Because DESs are mixtures, their physicochemical properties can be easily tailored by carefully selecting their precursors and varying mole ratios.Together with the sustainable character of common DES constituents and the absence of synthesis steps in their preparation, this led to extensive of research on the use of DESs as green solvents. 3,4However, most DESs proposed in the past contain at least one ionic precursor (types I -IV), limiting their application scope.Recently, though, a new type of DESs based solely on non-ionic compounds has been proposed (type V DESs), opening the road to create a whole range of non-ionic systems that better mimic the properties of conventional organic solvents. 5,6e liquid phase of DESs arises from the establishment of eutectic-type solid-liquid equilibrium (SLE) between the solid precursors, leading to a melting temperature depression of the mixture and, consequently, a liquid phase, at a given operational temperature.Because this depression is heightened by negative deviations from thermodynamic ideality, 2,6 the design of DESs entails selecting precursors that display strong cross intermolecular interactions, usually hydrogen bond donors (HBDs) and hydrogen bond acceptors (HBAs). 2 Please note the use of the term "cross interactions": negative deviations are favoured by strong HBD-HBA interactions and weak HBD-HBD or HBA-HBA association.For type V DESs where electrostatic interactions are not present, this is commonly achieved by selecting as precursors HBDs with weak HBA capability (also known as asymmetric HBDs) or HBAs without any HBD sites (also known as lone HBAs). 6For example, the prototypical type V DESs is the thymol-menthol system.Despite the structural similarity between thymol and menthol, the aromatic ring in thymol acts as an electron withdrawing agent to its hydroxyl group, leading to a stronger HBD but weaker HBA site. 5This results in a weaker thymol-thymol hydrogen bond (or, in other words, less thymol-thymol association) and a stronger thymol-menthol hydrogen bond, leading to a DES that is liquid at room temperature in a wide composition range.
Tertiary alcohols are a particularly interesting class of compounds in this area.By changing the substituents around the tertiary carbon, it is possible to drastically influence the molecule´s electronic distribution.In fact, akin to the aromaticity in the case of thymol explained above, fluorination can be used to increase the acidity of the alcohol hydroxyl group, with tert-butanol (TBH) and perfluoro tert-butanol (TBF) being an excellent example of this behaviour.Previous studies [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23] have shown the structure of TBH to be formed by coexisting monomers, dimers, and cyclic tetramers and hexamers.The population of these cyclic oligomers represents more than 65% of the total species in the liquid phase.However, the structure of TBF is markedly different. 24,25Its liquid phase is essentially constituted by monomers (70 -80%), with a minor fraction of dimers and almost no trimers.This weak tendency for self-association is due to the electron-withdrawing effect of the -CF 3 groups 26 that confers to TBF a strong acidic character that is reflected on its pKa=5.4, 27which is significantly lower than TBH, which has a pKa close to 19. 28 Additionally, given the bulky size of -CF 3 groups, steric hindrance also affects its hydrogen bonding capabilities preventing aggregation beyond the dimer.All these properties of TBF suggest that it could be an excellent candidate to act as an asymmetric HBD in the formulation of type V DES.
The behaviour of TBF as an asymmetric HBD is further supported by previous studies of the liquid phase of mixtures of TBF-TBH (note the resemblance between TBF-TBH and thymolmenthol). 29,30Using a variety of techniques including vibrational spectroscopy (Raman and Infrared) and X-ray diffraction, supported by DFT and MD simulations, it was possible to show that TBF acts as a structure breaker of TBH alcohol networks and that the 1:1 TBF-TBH hetero-dimer is the dominant species in the mole fraction concentration range of 0.3-07.For mixtures richer in TBF the characteristic monomers of the pure TBF liquid phase become dominant, while for mixtures richer in TBH oligomers start to appear with some relevance of the 2:1 trimer TBH-TBF-TBH.Once again, these characteristics of the TBF-TBH system are identical to those observed for the thymol-menthol system. 5taining a rigorous understanding of the behaviour of the TBF-TBH system connects not only to the formulation of DESs, but also to how the simultaneous presence of hydrogenated and perfluorinated chains affect the organization of a fluid.It is well known that hydrogenated and perfluorinated chains are mutually phobic, despite interacting through similar dispersion forces. 31Their mixtures display positive deviations from Raoult's law, liquid-liquid immiscibility, large positive excess enthalpies and volumes, and nanosegregation, [31][32][33] all clear signs of weak unlike intermolecular forces.Anomalies in transport, 34 surface, 35,36 and conformational 37 properties have also been reported.Liquids in which the two types of chains coexist are thus prone to organization 38 and self-assembling. 391][42] A network of hydrogen bonds are formed, zigzagging throughout the liquid, while the hydrogenated and perfluorinated segments occupy the remaining space while trying to avoid each other.The TBH-TBF mixture studied in this work is a further example of this type of mixtures.In this case, the presence of tertiary groups on both alcohols eliminates the influence of flexibility vs. stiffness and gauche conformers on the organization of the liquid, which was shown to be important in mixtures involving linear alcohols.
In order to probe the organization of the liquid mixture and its impact on the liquid phase nonideality and the solid-liquid equilibrium of the system, the excess enthalpies and excess volumes of TBH-TBF mixtures were measured across the entire composition range.These results were interpreted in terms of intermolecular interactions with the aid of molecular dynamics simulations, while gaining molecular level information on the organization and energetics of the tertiary alcohol binary system.This combination of excess properties, molecular dynamics, and hydrogen bonding analysis has been extensively and successfully used in the field of DESs to probe their liquid structures. 43

Experimental and Computational Details
The alcohols used in this study, pure (99.5%, anhydrous) tert-butanol and pure (99%) perfluoro-tert-butanol, were obtained from Thermo Scientifc and Alfa Aesar, respectively.To ensure the dryness of the alcohol, activated molecular sieves of pore size 0.3 nm with a humidity indicator, obtained from Metrohm AG, were used.

Excess Enthalpies
A high sensitivity heat conduction isothermal titration microcalorimeter from Thermometric AB/TA (now Waters TM ), together with a prototype water bath from Lund University, Sweden, was used to measure the partial excess molar enthalpies of mixtures of TBF in TBH, at atmospheric pressure and at 298.15 ± 0.01 K. 46 A 7½ digit Agilent nanovoltmeter (model 34420A) is connected to the calorimetric channel and interfaced with the computer, where data acquisition and the syringe pump control are performed through a modified version of the program LABTERMO. 46,47The system was electrically calibrated by release of controlled amounts of heat through a 'homemade' insertion heater with accurately known resistance, as well as chemically by use of the test reaction of dilution of a 10% propan-1-ol solution in water. 48Hastelloy® made titration cell of 1 mL volume equipped with a gold propeller (Lund University) was used for the ITC experiments, together with a stainless-steel solid cell as reference cell.Previous to the titration run the titration cell was warmed up to 310 K to prevent the crystallization of the tert-butanol.The sample cell was charged with solutions of different molar fraction of the two alcohols and the perfluoro-tert-butanol was titrated in aliquots of 10-40 L through a modified gas-tight Hamilton syringe.
The system was allowed to equilibrate between injections, re-establishing the baseline.After each injection, a peak corresponding to the heat generated or absorbed by the mixing process is recorded.The signal obtained was corrected using Tian's equation, and the corresponding peaks integrated numerically. 49The area of each peak is proportional to the heat involved in the mixing experiment and is thus used to calculate the partial excess molar enthalpy of TBF, .The experimental data was fitted using the derivative of the Redlich−Kister equation and   2 from the obtained parameters the excess molar enthalpy of mixing ( ) was calculated.

𝐴 𝑖 Δ mix 𝐻
To reproduce closely the composition dependence of the partial excess molar enthalpy while avoid introducing unphysical behaviour into the fitting equation, five (n=4) fitting parameters were used in equations 1-4 (1)

Excess Volumes
The density of TBF-TBH liquid mixtures was measured at 298.15 K, using an Anton Paar DMA 5000 vibrating tube densimeter.The densimeter was previously calibrated with atmospheric air and freshly boiled Millipore water at several temperatures, and the maximum deviation from literature values was found to be less than 0.00002 g•cm -3 .Between each measurement, the cell was thoroughly cleaned with isopropanol and dried with air and the density of air checked.Samples of the binary mixtures were prepared gravimetrically in 4 mL screw-cap glass vials, using a Mettler AE 240 analytical balance.The composition uncertainty was estimated as better than ± 0.0005 mole fraction units, leading to uncertainties in the individual values of the excess volumes under 0.05 cm 3 •mol -1 .

Differential Scanning Calorimetry
Thermal The TBH-TBF mixtures at different compositions were prepared gravimetrically in a glove box under dry nitrogen in batches of 1 g and kept in sealed glass bottles under nitrogen, to avoid contamination and sample evaporation.The batches of each mixture were stored for at least 24 hours (mixture homogenization) prior to the DSC crucibles preparation under dry N 2 (g).All measurements were performed under a constant flow of nitrogen (50 mL•min -1 ), using varying heating rates, 50 L sealed aluminium crucibles and samples were used.The mass of each experiment was selected according to the needs of the studied process as described previously. 54e composition uncertainty was estimated as better than ± 0.01 mole fraction.Each sample was cooled until 173 K and maintained at that temperature to ensure complete crystallization.
A final scan at 5 K•min −1 was performed to determine the temperatures of the solid-solid and solid-liquid phase transitions.

Molecular Dynamics Simulations
Atomistic molecular dynamics simulations of the studied TBF-TBH mixtures were performed using the GROMACS 5.0.7 simulation package. 55The tertiary alcohols were modelled with the OPLS-AA force field, 56 except for the partial charges which were obtained from DFT calculations.Full details on the used models were described in a previous publication. 25e unlike dispersive parameters were calculated using geometric combining rules, as prescribed by the OPLS-AA force field, with the exception of the interactions between hydrogen and fluorine atoms, where factors of 1.035 and 0.8 were applied to the crossed size ( ) and energy ( ) parameters, respectively.These deviations to the combining rules   -   - account for the peculiar volumetric and enthalpic effects that occur when perfluorinated and hydrogenated compounds are mixed, and were determined in a previous work for the n-butanol + 1H,1H-perfluorobutanol system. 40bic simulation boxes with periodic boundary conditions, containing 300 molecules were generated with the GROMACS "insert-molecules" program and pre-equilibrated in NpT using the Berendsen thermostat and barostat until the density reached an approximately constant value.The systems were then run in the NpT ensemble, using the leapfrog integrator with a 2 fs time step, for at least 35 ns, of which the first 5 ns were discarded as equilibration time and the remainder used to compute the averaged properties.Temperature and pressure were held constant during the production run, using the Nosé-Hoover thermostat 57,58 and Parrinello-Rahman barostat 59 with coupling times of 0.5 and 10 ps, respectively.The LINCS algorithm was applied to constrain all bonds containing hydrogen atoms to the respective equilibrium length. 60A cut-off distance of 1.4 nm was used to calculate both dispersive and electrostatic interactions.Beyond the cut-off distance, standard analytical corrections to the energy and pressure dispersive terms were applied and the particle-mesh Ewald method was used to account for long-range electrostatic interactions. 61,62ysical properties obtained from molecular simulation models are well known to depend significantly on the electrostatic characteristics, especially in the case of highly polar and hydrogen bonding compounds such as the alcohols studied in this work.On the other hand, charge distributions obtained by ab initio calculations vary significantly with the details of the calculation, namely the level of theory, basis sets and algorithm used to attribute the point charges from the full electronic density distribution.In this work, a sensitivity analysis was performed on the partial charges of the hydroxyl groups, to evaluate its effect on the shape and magnitude of the calculated excess enthalpies and volumes.Several combinations of the charge distributions of both compounds were tested, and those that yielded the most relevant results are detailed in Table 1.
In this table, the molecular models with the original point charges previously determined 25

Excess Properties
The study of the binary system TBF+TBH should be informative as to how a significant difference in acidity between two compounds is reflected in the non-ideality of their liquid mixture.As such, to probe the non-ideality of this system, excess enthalpies were measured by isothermal titration calorimetry (ITC).The experimental partial excess molar enthalpies H̅ 2 E , (TBH being defined as the component (1) and TBF component (2)) were fitted using the derivative of the Redlich−Kister equation, with five (n=4) fitting parameters.The very large negative deviations seen in Figure 1, around -7 kJ•mol -1 , confirm that, upon mixing, there is a considerable release of energy conferred by the preferred cross-association interaction.The perfluorinated compound TBF, which in its pure liquid phase is mostly found as isolated monomers, is able to form a very strong hydrogen bond in the presence of TBH. 30 Due to the lower ability of TBF to associate in the pure form, which can be attributed to the steric hindrance and inductive effects created by the -CF 3 substituents, the magnitude of the excess molar enthalpy is comparable to the energy of a single hydrogen bond.As such, the very negative H E values in this mixture derive essentially from the "excess" of hydrogen bonds formed, in comparison with the average of the pure compounds.The enthalpy of solution at infinite dilution of TBF in TBH and TBH in TBF was derived from the from the extrapolation of TBF partial excess molar enthalpies H̅ 2 E and H̅ 1 E to infinite dilution were found to be -(20.7 ± 0.1) kJ•mol -1 and -(19 ± 2) kJ•mol -1 respectively, which is related to the strong hydrogen bond interaction between TBF and TBH.The MD simulations, using the original charge distributions, are able to reproduce the sign, shape and order of magnitude of the H E curve, predicting a value that is around 50% of the experimental for the equimolar composition.This agreement can be considered remarkable, as the MD excess enthalpies are very sensitive to small variations in the partial charges of the hydroxyl groups.A sensitivity analysis was performed, and it was found that very small corrections of 0.015 to the charges of the two OH groups (Model 2), to the effect of promoting the preferred cross-association, increase the (absolute) calculated excess enthalpy enough to quantitatively reproduce the experimental values.It should be noted that these corrections are fully within the uncertainty of the ab initio scheme used to attribute the point charges.
The origin of the very negative values of H E in this system was further interpreted by analysing the MD trajectories to calculate the average distribution of hydrogen bonds in the studied mixtures, which has been used in the past to understand the liquid structure of DESs.The results are shown in figure 3, where the distribution of the four possible types of H-bonds and their sum is represented as a function of composition.For this calculation, any two molecules were counted as hydrogen bonded if the intermolecular O-H distance was under 0.27 nm, and 100% corresponds to a total limit number of two hydrogen bonds per molecule.The results shown were calculated from the simulations with the original charges (Model 1) for both compounds, and the alternative models tested present very similar results.As can be seen, in the mixtures the total number of hydrogen bonds formed is always larger than the average of the pure compounds and, moreover there is always a clear preference for formation of the most energetically favourable cross hydrogen bond between the oxygen of TBH and the hydrogen of TBF.From the results reported thus far, the enthalpic behaviour of the TBF-TBH system is clearly dominated by the changes in hydrogen bonding behaviour from the pure compounds to the mixtures.On the contrary, the well-known phobicity between hydrogenated and fluorinated chains, which normally results in positive excess enthalpies for mixtures of alkanes + perfluoroalkanes and hydrogenated alcohols + fluorinated alcohols, 40,63,64 plays a secondary role in the case of these tertiary alcohols and, instead, large negative deviations are observed due to the increased association in the mixture.The volumetric behaviour of the TBH+TBF system was also studied at 298.15 K, and the respective excess volumes are presented as Table SI.6 and in figure 4, along with the results obtained using MD simulations.The experimental data have been fitted to a four parameters (n=3) Redlich-Kister equation (analogous to Eq. 1), and the obtained coefficients are reported in table 2. As can be seen, the excess volume curve is very asymmetrical, with positive values on the TBH-rich side which peak at a value of 0.5 cm 3 •mol -1 at the 1(TBF) : 2(TBH) proportion and then decrease, becoming very close to zero for mixtures at the 1(TBF):1(TBH) proportion.
This unusual V E shape should be linked both with the variety of associating behaviours observed along the composition range 8,26,27 and with the volume expansion that is typically observed when perfluorinated and hydrogenated chains are mixed. 65,66At low TBF compositions, the volumetric behaviour follows the characteristic expansion of hydrogenated + perfluorinated systems, even though clearly attenuated by the enthalpically favourable mixing process.On the opposite side of the diagram, for TBF-rich mixtures, this characteristic expansion is countered, possibly because the hydrogen-bonded aggregates (especially the very stable heterodimer) are able to accommodate more efficiently.The MD simulation results are also shown in Figure 4.As can be seen, the simulations with Models 1 and 2 predict positive excess volume curves of the correct order of magnitude but fail to reproduce the peculiar shape of the experimental curve, predicting V E values that are essentially symmetric with composition (as is observed in mixtures of linear alcohols).Thus, a sensitivity analysis was again performed on the charge distribution of the OH group of the pure alcohols, to evaluate its influence in the shape of the simulated V E curve.It was found that, if the difference between O and H partial charges of TBH is slightly increased (keeping the original model for TBF), the simulated V E curve reproduces the experimental behaviour, albeit with smaller absolute values.These changes correspond to increasing the self-association of TBH and the intensity of the O(TBH)-H(TBF) cross-association.This combination of Model 3 of TBH and Model 1 of TBF slightly overpredicts the experimental excess enthalpy values, as seen in Figure 2. In conclusion, these results show that the peculiar properties of the TBH+TBF mixture result from a very subtle balance between the hydrogen-bonding characteristics of the pure compounds, combined with the unfavorable dispersive interactions between their hydrogenated and perfluorinated moieties.The expanded uncertainties of the experimental results, including the calibration uncertainty, were assigned based on independent experiments as t  s  , where t is obtained from student's t-distribution (0.95 level of n confidence), s is the standard deviation and n is the number of independent experiments."a": sum of the enthalpic contributions (SS + isotropization) ; "b ": sum of the entropic contributions (SS + isotropization) The asymmetry in melting temperatures of two seemingly similar compounds is related with differentiation of the hydrogen bonding capabilities in the pure form.Due to their low cohesive interaction in liquid phase and high conformational rigidity, fluorinated compounds usually have lower boiling points and higher melting points than the hydrogenated counterparts.

Physical Chemistry Chemical Physics Accepted Manuscript
However, TBF is an exception to the typical behaviour due to the lower ability of TBF to Another important aspect of this phase diagram is the formation of two co-crystals driven by the favourable interactions between two compounds and their favourable ability (size and/ shape) to pack in a co-crystal form.It is not surprising that the stoichiometries of solids and  are the most common aggregation states observed in the liquid mixtures; the hetero-dimer  TBH-TBF and the trimer TBH-TBF-TBH. 30The formation of stoichiometric compounds could be mitigated by decreasing the extent of fluorination and further decrease the melting temperature of eutectic mixtures.
The most relevant feature of the SLE diagram, is the considerable negative deviation from the ideality at compositions in-between the pure TBH and the eutectic points.This is particularly striking given that hydrocarbons and fluorocarbons mixtures are known to interact weakly as compared to their pure state, generally resulting in positive deviations to ideality and often leading to liquid-liquid phase separation.However, this effect is here heavily displaced by the strong hydrogen bonding cross-interactions, particularly the one between the proton of TBF and the oxygen of TBH.In fact, the cross-associative interactions are so strong that the two compounds react to form stoichiometric co-crystals.Thus, this clearly supports the idea that fluorinated compounds act as asymmetric hydrogen bond donors, with their liquid mixtures behaving as type V DESs, displaying severe negative deviations from ideality.

Conclusions
In this work, we covered the potential of fluorine atoms to increase the acidity of hydroxyl groups, in an attempt to understand and address the design of novel type V deep eutectic solvents.Besides the determination of the SLE diagram for the TBH+TBF system, energetic and volumetric characterizations of this system were also made, and interpretation was aided by the use of MD simulations.The results show significant negative deviations in the solidliquid transitions and excess molar enthalpies, as well as positive deviations in the excess molar volumes.The strong cross-association between the two tertiary alcohols is responsible for the non-ideality of these mixtures.Despite of the importance of the hetero-molecular interactions on the tendency of the two compounds of a binary mixture to display very low eutectic points, this aspect only reveals a part of the picture.The governing interactions of the pure constituents must not be overlooked, as well as their potential to neatly pack or form stable co-crystals.
Moreover, the results here reported showed that the degree of fluorination in tertiary alcohols, can be an interesting strategy for the tuning of the melting points, mixing properties, and control of co-crystallization.

53 Page
behaviour and temperatures of phase transition and fusion for TBF-TBH mixtures were evaluated in a power compensation differential scanning calorimeter, PERKIN ELMER model Pyris Diamond DSC, using hermetically sealed aluminium crucibles, samples of about 20 mg, and a constant flow of nitrogen (50 mL•min −1 ).The temperature and heat flux scales of the power compensation DSC were calibrated by measuring the temperature and the enthalpy of fusion of some reference materials namely, benzoic acid, 4-metoxibenzoic acid, triphenylene, naphthalene, anthracene, 1,3,5-triphenylbenzene, diphenylacetic acid, perylene, o-terphenyl, and 9,10-diphenylanthracene, at different scanning rates (2, 5 and 10 K•min -1 ). 50- Open Access Article.Published on 28 March 2023.Downloaded on 3/29/2023 9:45:02 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.View Article Online DOI: 10.1039/D3CP00701D are designated as "Model 1".In "Model 2", the difference between the partial charges of the O-H groups of each alcohol was increased relatively to one another, while keeping the O-H charge separation in each pure compound constant.This change increases the more probable crossassociation (the H-bond between the oxygen of TBF and the hydrogen of TBH).As will be shown below, this combination of charges reproduces the experimental H E curve but not the asymmetry of the V E curve.For "Model 3", the difference between the partial charges within the O-H group of TBH was increased, and this model was used together with the original (Model 1) of TBF.This combination increases both the self-association of TBH and the more probable cross-association in the mixture and is able to reproduce the highly asymmetric shape of the V E curve.

Figure 2 .
Figure 2. Excess molar enthalpies of the binary system TBF+TBH at T = 298.15K. Red crossesexperimental data derived from ITC; Solid curve -Second order Redlich-Kister equation fitted to the experimental points; Circles -Simulations with the original (Model 1) force field; Squares -Simulations with Model 2 for both TBF and TBH.Triangles -with the Model 3 of TBH and Model 1 of TBF.

Figure 3 .
Figure 3. Hydrogen-bond distribution in the (TBH+TBF) mixture as a function of composition, obtained from the MD trajectories for the original (Model 1) point charge models.Lines are visual guides.

Figure 4 .
Figure 4. Excess molar volumes of the binary system TBF+TBH at T = 298.15K. Crossesexperimental data; Dash curve -Redlich-Kister equation fitted to the experimental points; Circles -Simulations with the original (Model 1) force field; Squares -Simulations with Model 2 for TBH and TBF.Triangles -Simulations with the Model 3 of TBH and Model 1 of TBF.

Table 3 .
Thermal behaviour and temperatures of phase transition phase and fusion of tert-butanol and perfluoro-tert-butanol were evaluated with a DSC, PERKIN ELMER model Pyris Diamond DSC, using hermetically sealed aluminium crucibles.The obtained results (average of five independent experiments) are listed in table 3. The experimental results indicate an identical enthalpy and entropy of isotropization for both compounds (TBH and TBF).However, the perfluorinated tert-butanol sample shows a pre-melting solid-solid transition (SS) at 246.2 K with a strong enthalpic/entropic effect.Solid-solid transition temperature, T ss , standard (p o =10 5 Pa) molar enthalpy, (T fus ), and ∆ l s H o m entropy, (T fus ), of fusion for TBH and TBF.

OpenFigure 5 .
Figure 5. Experimental solid-liquid equilibrium phase diagram at atmospheric pressure for TBF+TBH.Green and grey circle symbols represent the data along the liquidus and solidus curves, respectively.Red circles represent the isotropization temperature of the pure samples, yellow square is the solidsolid transition of the TBF.The solid lines are guides for the eye following the experimental data, dotted curves represent the ideal behaviour (grey disregarding the solid-solid transition of TBF, yellow including it).The  and  correspond to the co-crystals of 1:2 and 1:1 stoichiometry, respectively, with their melting temperatures highlighted (blue and yellow circles, respectively).L represents the liquid phase.

Physical Chemistry Chemical Physics Accepted Manuscript Open Access
Article.Published on 28 March 2023.Downloaded on 3/29/2023 9:45:02 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Table 1 .
Atomic partial charges employed in the MD simulations.
Table 2 lists the obtained Redlich−Kister coefficients A i (n = 4), with additional information being provided in SI.Moreover, the Redlich−Kister equation was used to estimate the partial excess molar enthalpy of TBH, H̅ 1 E and the excess molar enthalpies of mixing, as described in SI.The excess molar enthalpies and partial excess molar enthalpies at 298.15 K are shown in Figure 1.