Molecular Dynamics Calculations for the Temperature Response of Poly(alkylated tri(ethylene oxide)isocyanate) Aqueous Solution

: Aqueous solutions of conventional temperature-responsive amphiphilic polymers undergo a coil–globule conformational transition around the lower critical solution temperature (LCST) that causes the polymer surfaces to become hydrophobic and the polymers to aggregate together. Iso-cyanate polymers with alkylated oligo(ethylene oxide) side chains are expected to have rigid main chains and, thus, do not undergo the coil–globule structural transition, but they have recently been reported to exhibit temperature-responsive properties. In this study, molecular dynamics was used to calculate the agglomeration tendencies of two chains of poly(alkylated tri(ethylene oxide)isocyanate) (PRTEOIC, where R = methyl (Me) or ethyl (Et)) in aqueous solution to elucidate the LCST phenomenon in the absence of coil–globule conformational transition. Our MD simulations showed that aggregation also occurs in rod polymers. Furthermore, we found that both (PMeTEOIC) 2 and (PEtTEOIC) 2 showed parallel agglomeration of the two molecular chains with increasing temperature, but only (PMeTEOIC) 2 showed a metastable T-shaped agglomeration in the middle temperature range. The crossing-point temperature ( T CRP ) at which the density of the ﬁrst hydrophobic hydration shell around the sidechain alkyl group equals the bulk water density is a useful indicator for predicting the LCST of rod polymers with dense side chains terminated by alkyl groups.


Introduction
Aqueous solutions of amphiphilic polymers possessing a flexible main chain exhibit a phase transition at the lower critical solution temperature (LCST) [1].Such polymers become soluble below the LCST and insoluble above it.A representative example is poly(N-isopropylacrylamide) (PNIPAM), which has been extensively studied both experimentally [2,3] and theoretically [4].Below the LCST, the side-chain amide group is hydrogen bonded to water molecules, which results in a "coil" state in which the polymer chain expands.Above the LCST, the side-chain amide group forms hydrogen bonds within the polymer, and the polymer chain adopts a collapsed "globule" state.Therefore, a coil-globule structural transition is considered a necessary driving force for LCST-type phase transitions to occur.Experimentally, the coil-globule transition at the LCST has been observed by fluorescence measurements [5][6][7], dynamic light scattering (DLS) measurements [8][9][10], small-angle X-ray scattering [11], and small-angle neutron scattering [12].Molecular dynamics (MD) can be used to explain the correlation between the phase transition and coil-globule transition.Deshmukh and Mancini et al. reported that 30 mer PNIPAM in water undergoes the coil-globule structural transition in 3 ns at 310 K, which is higher than the LCST of 305 K [13].Tavagnacco and Chiessi et al. performed MD calculations for 30 mer PNIPAM in solution and showed that a decrease in the first hydration shell around the isopropyl groups was accompanied by an increase in hydrogen bonding between the polymer and water molecules, indicating that the coil-globule structural transition leads to a significant rearrangement of the hydrogen-bonding pattern [14].
Our group recently reported a novel temperature-responsive water-soluble isocyanate polymer called poly(alkylated tri(ethylene oxide)isocyanate) (PRTEOIC, where R = methyl (Me) or ethyl (Et)), with a degree of polymerization of approximately 20 and polydispersity index of 1.1, that is expected to undergo a LCST phase transition without a coil-globule structural transition [15].Increasing the bulkiness of the alkyl groups at the end of the side chains should lower the LCST of isocyanate polymers.DLS, dipole moment measurements, and quantum chemical calculations have indicated that polyisocyanates with bulky alkyl side chains generally have rigid helical rod structures [16][17][18][19][20].In particular, poly(hexyl isocyanate) (PHIC) becomes a liquid crystal material because of its stiff and helical structure with chirality [21].DLS has shown that PHIC with a degree of polymerization of ~1600 has the relatively long persistence length of 20-40 nm in solvents such as hexane and chloroform [22][23][24][25].The change in the persistence length with solvent polarity can be attributed to changes in the torsional oscillation around the main chain with polar solvents, as shown by Cook et al., who took linewidth measurements of nuclear magnetic resonance (NMR) spectra [26].
Interestingly, PROEOIC, which has a small polymerization degree of 20, also undergoes a LCST-type phase transition, even though it should have a rigid structure, unlike PNIPAM.Turbidity measurements have shown that the cloud point temperature (T CLP ) is strongly dependent on the hydrophobicity of the end group of the side chain.For example, PEtTEOIC showed off T CLP = 308 K, while PMeTEOIC T CLP = 353 K (Supplementary Figure S1).
In the present study, we performed all-atom classical MD calculations on the temperature response for the agglomeration tendencies of two chains of PRTEOIC in aqueous solution.It has been reported that the solvent exclusion volume effect could be an indicator of LCST phase transition [27,28].Herein, we investigated the mechanism for the LCST-type phase transition of rod-type polymers with a small degree of polymerization and found that the hydrophobic hydration shell strength at the end of side chain could be an indicator of the phase-transition temperature in the absence of the coil-globule structural transition.

Simulation Details
The Amber (ver.14 or 20) software package [29] was used to carry out all-atom MD simulations for aqueous solutions of one or two chains of PMeTEOIC and PEtTEOIC at six temperatures from 288 K to 370 K.The degree of polymerization was set to 20 to be close to the experimental molecular weight reported by Sakai and Kakuchi et al. [30].MD simulations were also performed for two molecules of alkyl tri(ethylene oxide) (MeTEO, EtTEO) in aqueous solution for comparison with their combined state when attached to the isocyanate polymer (Scheme 1).Scheme 1.Chemical structures of PRTEOIC and RTEO (R = CH3, C2H5).
These molecules were structurally optimized by the semi-empirical molecular orbital method (AM1), and their atomic charges were determined by the AM1-BCC method [31] by using the Antechamber tool of the AMBER package.The general AMBER force field (GAFF) was employed as the molecular force field [32,33].The two polymer chains were placed 12 Å apart as an initial arrangement by using the XLeap tool of the AMBER package.The TIP3P water model was used as a solvent molecule, and the basic simulation cell was solvated to form a cleaved octahedron under the periodic boundary condition.Approximately, 7000 water units were placed in the cell.The concentrations were about 60 g/L for the double polymer solutions (Supplementary Materials Table S1).Energy minimization of the whole system consisting of polymers and solvent water molecules was carried out in two stages, at 0 K.The first energy optimization was carried out in 1000 steps by the steepest descent method, in which two solvated polymer molecules were constrained by a force of 500 kcal/mol Å 2 .The structure was then allowed to move slightly in the direction of the greatest change on the energy surface.In the second energy optimization, the constraint was removed, and 1000 steps of the steepest descent method were carried out, followed by 1500 steps of the conjugate gradient method.The latter method efficiently found the minimum point by determining the vector through which the structure moves to satisfy the conjugacy property.For temperature elevation followed by equilibration, MD calculations were performed under a weak constraint of 10 kcal/mol Å 2 during the first 40 ps.Subsequently, the system temperature was raised from 0 K to the target temperature (i.e., 288 K, 323 K, 333 K, 343 K, 363 K, or 370 K).The system was then given an additional 100 ps for equilibration, which was a total equilibration time of 140 ps.The temperature of the system was equilibrated by using the Langevin temperature control method.A 20 ns simulation was performed for the PRTEOIC dimer.The calculations were carried out 10 times, and the average value at each temperature was taken.A 200 ns simulation was performed for the evaluation of the coil-globule structural transition of PRTE-OIC, and 40 ns was performed for the alkylated tri(ethylene oxide) (RTEO) dimer.The analysis used the last 30 ns.All calculations used a time step of 2 fs.
The hydrogen atoms were constrained by the SHAKE method.The NPT ensemble (constant pressure periodic boundary condition) was used, where the number of molecules, temperature, and pressure were kept constant.The temperature and pressure were controlled by using the Langevin temperature control method and Berendsen method, with a mean pressure of 1 atom and relaxation time of 2 ps.Long-range electrostatic interactions were calculated by using the particle mesh Ewald method [34], with a cutoff distance of 9 Å.Trajectories were plotted every 2.0 ps.The collision frequency was set to 1.0 ps −1 .These molecules were structurally optimized by the semi-empirical molecular orbital method (AM1), and their atomic charges were determined by the AM1-BCC method [31] by using the Antechamber tool of the AMBER package.The general AMBER force field (GAFF) was employed as the molecular force field [32,33].The two polymer chains were placed 12 Å apart as an initial arrangement by using the XLeap tool of the AMBER package.The TIP3P water model was used as a solvent molecule, and the basic simulation cell was solvated to form a cleaved octahedron under the periodic boundary condition.Approximately, 7000 water units were placed in the cell.The concentrations were about 60 g/L for the double polymer solutions (Supplementary Materials Table S1).Energy minimization of the whole system consisting of polymers and solvent water molecules was carried out in two stages, at 0 K.The first energy optimization was carried out in 1000 steps by the steepest descent method, in which two solvated polymer molecules were constrained by a force of 500 kcal/mol Å 2 .The structure was then allowed to move slightly in the direction of the greatest change on the energy surface.In the second energy optimization, the constraint was removed, and 1000 steps of the steepest descent method were carried out, followed by 1500 steps of the conjugate gradient method.The latter method efficiently found the minimum point by determining the vector through which the structure moves to satisfy the conjugacy property.For temperature elevation followed by equilibration, MD calculations were performed under a weak constraint of 10 kcal/mol Å 2 during the first 40 ps.Subsequently, the system temperature was raised from 0 K to the target temperature (i.e., 288 K, 323 K, 333 K, 343 K, 363 K, or 370 K).The system was then given an additional 100 ps for equilibration, which was a total equilibration time of 140 ps.The temperature of the system was equilibrated by using the Langevin temperature control method.A 20 ns simulation was performed for the PRTEOIC dimer.The calculations were carried out 10 times, and the average value at each temperature was taken.A 200 ns simulation was performed for the evaluation of the coil-globule structural transition of PRTEOIC, and 40 ns was performed for the alkylated tri(ethylene oxide) (RTEO) dimer.The analysis used the last 30 ns.All calculations used a time step of 2 fs.
The hydrogen atoms were constrained by the SHAKE method.The NPT ensemble (constant pressure periodic boundary condition) was used, where the number of molecules, temperature, and pressure were kept constant.The temperature and pressure were controlled by using the Langevin temperature control method and Berendsen method, with a mean pressure of 1 atom and relaxation time of 2 ps.Long-range electrostatic interactions were calculated by using the particle mesh Ewald method [34], with a cutoff distance of 9 Å.Trajectories were plotted every 2.0 ps.The collision frequency was set to 1.0 ps −1 .

Analysis
The radius of gyration, R g , represents the extension degree of the molecular chain from the center of gravity.R g was employed to investigate the agglomeration degree of two polymer chains, as well as the chain extension of a single polymer chain.R g is expressed as follows: where N is the total number of atoms, r i (t) is the atomic coordinate vector, and r c (t) is the center of gravity.Considering the time at which the double polymer chains reaches the thermal equilibrium state, the distribution function of R g (t) was obtained by excluding the first 10 ns of the simulation.
To analyze the structure of the hydrated shell around the polymer, we obtained a radial distribution function (RDF) for the terminal carbon atoms of the alkylated tri(ethylene oxide) side chain to the oxygen atoms of water molecules.During the hydration of nonpolar molecules, hydrogen bonds do not occur between the solute and solvent.However, the water molecules surrounding the nonpolar solute stabilize into a network structure, which is called a hydrophobic hydrated shell.The RDF can be calculated as follows: where n(r) is the number of particles between spherical shells at distances r and dr from a certain particle, and ρ average is the average density of the system.When the inter-particle distance, r, is very small, the particle cannot exist because of the repulsion between particles, so g(r) = 0.When r → ∞ , the relation of g Ce−Ow (r) = 1 holds because water molecules exist at the average density.After the RDF was obtained, the peak intensity of g Ce−Ow (r) for the first hydrated shell was extracted and was plotted against the temperature.Then, the crossing-point temperature (T CRP ), which is defined as the temperature at which the average density of bulk water (g Ce−Ow (r) = 1) and peak intensity of the RDF intersect, was determined for comparison with the T CLP .

Rigidity of the Main Chain
Isocyanates with bulky alkyl side chains are expected to be rigid.Figure 1a shows the change in R g over time at three temperatures selected according to the experimental turbidity measurements (Supplementary Materials Figure S1) [30]-288 K (at which both PMeTEOIC and PEtTEOIC dissolve), 343 K (at which PMeTEOIC dissolves and PEtTEOIC aggregates), and 370 K (at which both PMeTEOIC and PEtTEOIC aggregate).No significant changes in R g over time were observed for either PMeTEOIC or PEtTEOIC at any temperature (Figure 1b).These R g results suggest that PRTEOIC's are rigid.
To further examine the rigidity of PRTEOIC during the MD simulations, the distributions of the dihedral angles formed by C-N-C-O in the main chain at 363 K were investigated based on the obtained trajectories (Supplementary Materials Figures S2-S4).For the dihedral angle distributions, most (85%) had only one peak.Therefore, little isomerization occurred between the cisoids and transoids.This means that the structure of the main chain was mostly maintained even at high temperatures.The proportions of cisoids and transoids remained at 32% and 68%, respectively.The C-N-C-O bond near the H end of the main chain underwent almost no cisoid-transoid isomerization, while the C-N-C-O bond near the OC 4 H 9 end of the main chain was more likely to undergo cisoid-transoid isomerization.Therefore, the hydrophilic H end of the main chain was less likely to isomerize, while the hydrophobic OC 4 H 9 end was more likely to isomerize.
The MD simulations with the GAFF indicated that the degree of chain extension of PRTEOIC is independent of the temperature and that there is no coil-globule transition above and below the experimentally obtained LCST.
for polymers with methacrylate main chains and OEG side chains [35,36].To confirm the rigidity of the PRTEOIC main chain, the rotational barrier was estimated by applying the density functional theory (DFT) B3LYP/6-31G(d) for MeTEOIC3 (Supplementary Materials Figure S6), using Gaussian 16 [37].The rotational barrier obtained by subtracting the energy of the stable structure from the energy of the transition structure was at least 0.664 eV, which is about 26 times higher than the kinetic energy at room temperature (i.e., 0.0257 eV at 298 K).An analysis of the side-chain parts indicated that the averaged side-chain end-to-end distances also remained constant regardless of the temperature (Supplementary Materials Figure S5).This suggests that the folding of the PRTEOIC side chain is not related to its LCST phase transition.This result is consistent with the calculations of Dalgakiran et al. for polymers with methacrylate main chains and OEG side chains [35,36].To confirm the rigidity of the PRTEOIC main chain, the rotational barrier was estimated by applying the density functional theory (DFT) B3LYP/6-31G(d) for MeTEOIC 3 (Supplementary Materials Figure S6), using Gaussian 16 [37].The rotational barrier obtained by subtracting the energy of the stable structure from the energy of the transition structure was at least 0.664 eV, which is about 26 times higher than the kinetic energy at room temperature (i.e., 0.0257 eV at 298 K).
Ute and Green et al. performed temperature-dependent NMR and experimentally obtained a free-energy activation barrier of 0.85 eV for poly(2-butylhexyl isocyanate) (PB-HIC) [20].Young and Cook performed MD calculations by using the Merck molecular force field and obtained a rotational barrier of 0.76 eV at maximum for 21 mer poly(methyl isocyanate) [18].Our calculation results are comparable with these values in the literature.These quantum calculations also revealed that the main chain of PRTEOIC is rigid and that the LCST phase transition is not coupled with the structural transition of the isocyanate main chain.

Agglomeration of (PRTEOIC) 2
To investigate the agglomeration behavior of the bimolecular polymer chains, MD simulations were performed to determine the change in R g over time from the center of gravity of (PRTEOIC) 2 , which was used to calculate the R g distribution, as shown in Figures 2-4.Snapshots of the bimolecular chains corresponding to each peak are also shown in these figures.Figure 2 shows the R g distributions at 288 K, at which the two polyisocyanates are both soluble in water.The R g distribution was broad, with several peaks for both (PMeTEOIC) 2 and (PEtTEOIC) 2 .The peaks below 13 Å can be attributed to the bimolecular chains in bound states, while the peaks above 18 Å can be attributed to the bimolecular chains in dissociated states.Both polymers were in bond-dissociation equilibrium.Furthermore, (PMeTEOIC) 2 had an R g peak at about 20 Å that was absent in (PEtTEOIC) 2 , thus indicating that the polymers separated more in (PMeTEOIC) 2 than in (PEtTEOIC) 2 .
Macromol 2023, 3, FOR PEER REVIEW Ute and Green et al. performed temperature-dependent NMR and experime obtained a free-energy activation barrier of 0.85 eV for poly(2-butylhexyl isocya (PBHIC) [20].Young and Cook performed MD calculations by using the Merck mole force field and obtained a rotational barrier of 0.76 eV at maximum for 21 mer poly(m isocyanate) [18].Our calculation results are comparable with these values in the liter These quantum calculations also revealed that the main chain of PRTEOIC is rigi that the LCST phase transition is not coupled with the structural transition of the is nate main chain.

Agglomeration of (PRTEOIC)2
To investigate the agglomeration behavior of the bimolecular polymer chains simulations were performed to determine the change in Rg over time from the cen gravity of (PRTEOIC)2, which was used to calculate the Rg distribution, as shown i ures 2-4.Snapshots of the bimolecular chains corresponding to each peak are also s in these figures.Figure 2 shows the Rg distributions at 288 K, at which the two poly anates are both soluble in water.The Rg distribution was broad, with several pea both (PMeTEOIC)2 and (PEtTEOIC)2.The peaks below 13 Å can be attributed to the lecular chains in bound states, while the peaks above 18 Å can be attributed to the lecular chains in dissociated states.Both polymers were in bond-dissociation equilib Furthermore, (PMeTEOIC)2 had an Rg peak at about 20 Å that was absent in (PEtTEO thus indicating that the polymers separated more in (PMeTEOIC)2 than in (PEtTEO    Macromol 2023, 3, FOR PEER REVIEW Ute and Green et al. performed temperature-dependent NMR and experime obtained a free-energy activation barrier of 0.85 eV for poly(2-butylhexyl isocya (PBHIC) [20].Young and Cook performed MD calculations by using the Merck mole force field and obtained a rotational barrier of 0.76 eV at maximum for 21 mer poly(m isocyanate) [18].Our calculation results are comparable with these values in the liter These quantum calculations also revealed that the main chain of PRTEOIC is rigid that the LCST phase transition is not coupled with the structural transition of the is nate main chain.

Agglomeration of (PRTEOIC)2
To investigate the agglomeration behavior of the bimolecular polymer chains simulations were performed to determine the change in Rg over time from the cen gravity of (PRTEOIC)2, which was used to calculate the Rg distribution, as shown in ures 2-4.Snapshots of the bimolecular chains corresponding to each peak are also s in these figures.Figure 2 shows the Rg distributions at 288 K, at which the two polyi anates are both soluble in water.The Rg distribution was broad, with several peak both (PMeTEOIC)2 and (PEtTEOIC)2.The peaks below 13 Å can be attributed to the b lecular chains in bound states, while the peaks above 18 Å can be attributed to the b lecular chains in dissociated states.Both polymers were in bond-dissociation equilib Furthermore, (PMeTEOIC)2 had an Rg peak at about 20 Å that was absent in (PEtTEO thus indicating that the polymers separated more in (PMeTEOIC)2 than in (PEtTEO    distribution for (PMeTEOIC)2 had a tail around 17 Å.The agglomeration structures (PMeTEOIC)2 and (PEtTEOIC)2 were parallel around 13 Å, while that of (PMeTEOI was T-shaped around 17 Å. Figure 4 shows the Rg distributions at 370 K, at which both (PMeTEOIC)2 and (P TEOIC)2 are insoluble.The tail due to the T-shaped agglomeration of (PMeTEOIC)2 t was observed at 343 K disappeared, while the distribution of (PEtTEOIC)2 remained s ilar to that at 343 K.The snapshots show that both agglomerations had parallel shapes In summary, the Rg distributions were broad with multiple peaks; this suggests t the two molecules moved freely below the LCST and that the peaks sharpened near LCST.Above the LCST, the peaks remained sharp, which suggests that the agglomerat of (PRTEOIC)2 stabilized.The temperature responses of these Rg distributions are c sistent with an LCST-type temperature response in which a polymer aggregates and solution becomes cloudy at high temperatures.
Figure 5 shows the T-shaped aggregates observed only in (PMeTEOIC)2 at 343 K this temperature, PMeTEOIC was dispersed in the aqueous solution and formed sm aggregates that scattered light, which reduced the transmittance.The T-shaped dim may be the basic structure of this small metastable aggregate.This raises the question why the metastable T-shaped dimer was only observed in (PMeTEOIC)2.An interest study on the structure of the benzene dimer may help answer this question.Miyazaki a Fujii et al. performed detailed spectroscopic measurements of supersonic jets and show that benzene dimers have a T-shaped structure in the ground electronic state and a pa lel sandwich structure in the excimer state [38].The main difference between the grou state and excimer state is the charge distribution on each atom.These results suggest t a change in charge distribution affects the balance between the Coulomb force and L don dispersion force, which changes the dimer structure.In this study, we used the M simulation data to analyze the electrostatic energy and van der Waals energy, which related to intermolecular interactions between polymer chains.The distributions of electrostatic and van der Waals energies calculated for two polymer chains are shown Supplementary Materials Figures S7 and S8.
Figure 4 shows the R g distributions at 370 K, at which both (PMeTEOIC) 2 and (PEtTEOIC) 2 are insoluble.The tail due to the T-shaped agglomeration of (PMeTEOIC) 2 that was observed at 343 K disappeared, while the distribution of (PEtTEOIC) 2 remained similar to that at 343 K.The snapshots show that both agglomerations had parallel shapes.
In summary, the R g distributions were broad with multiple peaks; this suggests that the two molecules moved freely below the LCST and that the peaks sharpened near the LCST.Above the LCST, the peaks remained sharp, which suggests that the agglomeration of (PRTEOIC) 2 stabilized.The temperature responses of these R g distributions are consistent with an LCST-type temperature response in which a polymer aggregates and the solution becomes cloudy at high temperatures.
Figure 5 shows the T-shaped aggregates observed only in (PMeTEOIC) 2 at 343 K.At this temperature, PMeTEOIC was dispersed in the aqueous solution and formed small aggregates that scattered light, which reduced the transmittance.The T-shaped dimer may be the basic structure of this small metastable aggregate.This raises the question of why the metastable T-shaped dimer was only observed in (PMeTEOIC) 2 .An interesting study on the structure of the benzene dimer may help answer this question.Miyazaki and Fujii et al. performed detailed spectroscopic measurements of supersonic jets and showed that benzene dimers have a T-shaped structure in the ground electronic state and a parallel sandwich structure in the excimer state [38].The main difference between the ground state and excimer state is the charge distribution on each atom.These results suggest that a change in charge distribution affects the balance between the Coulomb force and London dispersion force, which changes the dimer structure.In this study, we used the MD simulation data to analyze the electrostatic energy and van der Waals energy, which are related to intermolecular interactions between polymer chains.The distributions of the electrostatic and van der Waals energies calculated for two polymer chains are shown in Supplementary Materials Figures S7 and S8.
Table 1 summarizes the ratio of the electrostatic energy to the van der Waals energy (|Eelec/EvdW|).Errors were estimated from the full width at half maximum for each energy distribution.The ratios of (PMeTEOIC)2 and (PEtTEOIC)2 were significantly different, which may explain why the metastable T-shaped agglomeration pair was only observed in (PMeTEOIC)2.

Radial Distribution Functions of Water Molecules around the PRTEOIC Dimer
To investigate the effect of temperature on the water molecules surrounding the polymer chains, we analyzed the RDF of the water oxygen atoms from the terminal carbon atoms of the side chain of PRTEOIC ( ()), which we used to estimate the peak intensity of the first hydrophobic hydration shell, as shown in Figure 6.This analysis was normalized to  () = 1 at 25 Å, where the presence of the polymer had a negligible effect on the density of water.Around the terminal carbon atoms of the polymer side chains, the first hydration shells formed at around 3.6 Å.At even longer distances, the multiple tri(ethylene oxide) side chains were densely attached to the main chain, which prevented water molecules from entering the periphery of the side chain and resulted in a density of less than unity.For the peaks of the first hydrophobic hydration shell of (PMeTEOIC)2 and (PEtTEOIC)2, both peak intensities decreased with the increasing temperature and fell below  () = 1 (corresponding to the bulk water density) at different temperatures.At all temperatures, the hydrophobic hydration shell of (PEtTEOIC)2 had a lower peak intensity than that of (PMeTEOIC)2.Figure 7 plots the peak intensities for the first hydrophobic hydration shells of (PRTEOIC)2 and (RTEO)2 against the temperature.TCRP (i.e., the temperature at which the peak intensity intersected the bulk water value of  () = 1) was 358 K for (PMeTEOIC)2 and 322 K for (PEtTEOIC)2.These temperatures are close to the cloud point temperatures of 353 K for PMeTEOIC and 308 K for PEtTEOIC.At 283 K, the electrostatic energy (E elec ) was about 1280 kcal/mol for (PMeTEOIC) 2 and about 1320 kcal/mol for (PEtTEOIC) 2 .The van der Waals energy (E vdW ) was about −370 kcal/mol for (PMeTEOIC) 2 and about −440 kcal/mol for (PEtTEOIC) 2 (Supplementary Materials Tables S2 and S3).The absolute values of E elec and E vdW were larger for (PEtTEOIC) 2 than for (PMeTEOIC) 2 , and this trend remained unchanged with the increasing temperature.Even at lower temperatures, (PEtTEOIC) 2 had greater intermolecular interactions than (PMeTEOIC) 2 , suggesting that it is more prone to agglomeration.
Table 1 summarizes the ratio of the electrostatic energy to the van der Waals energy (|E elec /E vdW |).Errors were estimated from the full width at half maximum for each energy distribution.The ratios of (PMeTEOIC) 2 and (PEtTEOIC) 2 were significantly different, which may explain why the metastable T-shaped agglomeration pair was only observed in (PMeTEOIC) 2 .

Radial Distribution Functions of Water Molecules around the PRTEOIC Dimer
To investigate the effect of temperature on the water molecules surrounding the polymer chains, we analyzed the RDF of the water oxygen atoms from the terminal carbon atoms of the side chain of PRTEOIC (g Ce−Ow (r)), which we used to estimate the peak intensity of the first hydrophobic hydration shell, as shown in Figure 6.This analysis was normalized to g Ce−Ow (r) = 1 at 25 Å, where the presence of the polymer had a negligible effect on the density of water.Around the terminal carbon atoms of the polymer side chains, the first hydration shells formed at around 3.6 Å.At even longer distances, the multiple tri(ethylene oxide) side chains were densely attached to the main chain, which prevented water molecules from entering the periphery of the side chain and resulted in a density of less than unity.For the peaks of the first hydrophobic hydration shell of (PMeTEOIC) 2 and (PEtTEOIC) 2 , both peak intensities decreased with the increasing temperature and fell below g Ce−Ow (r) = 1 (corresponding to the bulk water density) at different temperatures.At all temperatures, the hydrophobic hydration shell of (PEtTEOIC) 2 had a lower peak intensity than that of (PMeTEOIC) 2 .Figure 7 plots the peak intensities for the first hydrophobic hydration shells of (PRTEOIC) 2 and (RTEO) 2 against the temperature.T CRP (i.e., the temperature at which the peak intensity intersected the bulk water value of g Ce−Ow (r) = 1) was 358 K for (PMeTEOIC) 2 and 322 K for (PEtTEOIC) 2 .These  For (RTEO)2, which showed no agglomeration in the turbidity measurements, the peak intensity of the first hydrophobic hydration shell was much higher than  () = 1 at all temperatures because more water molecules were involved in the solvation around the solute.The Rg distributions of (RTEO)2 in aqueous solution at 288, 343, and 370 K are presented in Supplementary Materials Figure S1.For both (MeTEO)2 and (EtTEO)2, the broad shape of the Rg distribution was almost the same at all temperatures.There was a small peak near 5 Å, which means that the two molecules were temporarily close to each other.However, because the main peak was around 14 Å, the two molecules mostly remained apart.Therefore, no stable dimer structure formed in the (RTEO)2 bimolecular system at any temperature, thus indicating that RTEO is not temperature-responsive.
Aqueous solutions of poly(ethylene oxide) with a high molecular weight are known to exhibit an LCST around 373 K [39], whereas aqueous solutions of oligo(ethylene oxide) with a low molecular weight do not exhibit an LCST-type phase transition.In our MD simulations, the RTEO bimolecular chain did not form a stable dimer.This is because RTEO chains with low molecular weights form a strong hydrophobic hydration shell around their alkyl groups, and the density of this strong hydrophobic shell did not fall   For (RTEO)2, which showed no agglomeration in the turbidity measurements, the peak intensity of the first hydrophobic hydration shell was much higher than  () = 1 at all temperatures because more water molecules were involved in the solvation around the solute.The Rg distributions of (RTEO)2 in aqueous solution at 288, 343, and 370 K are presented in Supplementary Materials Figure S1.For both (MeTEO)2 and (EtTEO)2, the broad shape of the Rg distribution was almost the same at all temperatures.There was a small peak near 5 Å, which means that the two molecules were temporarily close to each other.However, because the main peak was around 14 Å, the two molecules mostly remained apart.Therefore, no stable dimer structure formed in the (RTEO)2 bimolecular system at any temperature, thus indicating that RTEO is not temperature-responsive.
Aqueous solutions of poly(ethylene oxide) with a high molecular weight are known to exhibit an LCST around 373 K [39], whereas aqueous solutions of oligo(ethylene oxide) with a low molecular weight do not exhibit an LCST-type phase transition.In our MD simulations, the RTEO bimolecular chain did not form a stable dimer.This is because RTEO chains with low molecular weights form a strong hydrophobic hydration shell around their alkyl groups, and the density of this strong hydrophobic shell did not fall For (RTEO) 2 , which showed no agglomeration in the turbidity measurements, the peak intensity of the first hydrophobic hydration shell was much higher than g Ce−Ow (r) = 1 at all temperatures because more water molecules were involved in the solvation around the solute.The R g distributions of (RTEO) 2 in aqueous solution at 288, 343, and 370 K are presented in Supplementary Materials Figure S1.For both (MeTEO) 2 and (EtTEO) 2 , the broad shape of the R g distribution was almost the same at all temperatures.There was a small peak near 5 Å, which means that the two molecules were temporarily close to each other.However, because the main peak was around 14 Å, the two molecules mostly remained apart.Therefore, no stable dimer structure formed in the (RTEO) 2 bimolecular system at any temperature, thus indicating that RTEO is not temperature-responsive.
Aqueous solutions of poly(ethylene oxide) with a high molecular weight are known to exhibit an LCST around 373 K [39], whereas aqueous solutions of oligo(ethylene oxide)

Figure 1 .
Figure 1.(a) Rg changes up to 200 ns of simulation time for the whole single polymer chain of PMeTEOIC (left) and PEtTEOIC (right) at 288 K, 343 K, and 370 K. (b) Rg values for the whole single polymer chain of PMeTEOIC (red) and PEtTEOIC (blue) vs. temperature.

Figure 1 .
Figure 1.(a) R g changes up to 200 ns of simulation time for the whole single polymer chain of PMeTEOIC (left) and PEtTEOIC (right) at 288 K, 343 K, and 370 K. (b) R g values for the whole single polymer chain of PMeTEOIC (red) and PEtTEOIC (blue) vs. temperature.

Figure 2 .
Figure 2. Rg distributions of (PMeTEOIC)2 (left) and (PEtTEOIC)2 (right) in water at 288 K dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 3 .
Figure 3. Rg distributions of (PMeTEOIC)2 (left) and (PEtTEOIC)2 (right) in water at 343 K dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 3
Figure3shows the Rg distributions at 343 K, at which (PMeTEOIC)2 is solubl (PEtTEOIC)2 is insoluble in water.Both distributions had a peak around 13 Å

Figure 2 .
Figure 2. R g distributions of (PMeTEOIC) 2 (left) and (PEtTEOIC) 2 (right) in water at 288 K.The dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 2 .
Figure 2. Rg distributions of (PMeTEOIC)2 (left) and (PEtTEOIC)2 (right) in water at 288 K dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 3 .
Figure 3. Rg distributions of (PMeTEOIC)2 (left) and (PEtTEOIC)2 (right) in water at 343 K dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 3
Figure3shows the Rg distributions at 343 K, at which (PMeTEOIC)2 is solubl (PEtTEOIC)2 is insoluble in water.Both distributions had a peak around 13 Å

Figure 3 .
Figure 3. R g distributions of (PMeTEOIC) 2 (left) and (PEtTEOIC) 2 (right) in water at 343 K.The dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 4 .
Figure 4. Rg distributions of (PMeTEOIC)2 (left) and (PEtTEOIC)2 (right) in water at 370 K. dashed line shows the borderline of the agglomerated state and the dissociated state.

Figure 4 .
Figure 4. Rg distributions of (PMeTEOIC) 2 (left) and (PEtTEOIC) 2 (right) in water at 370 K.The dashed line shows the borderline of the agglomerated state and the dissociated state.
1Value of Eelec averaged over six temperatures.2Standarddeviation of Eelec at six temperatures.

Table 1 .
Ratio of electrostatic energy to van der Waals energy (|E elec /E vdW |) at six temperatures.
1Value of E elec averaged over six temperatures.2Standarddeviation of E elec at six temperatures.