Relating side chain organization of PNIPAm with its conformation in aqueous methanol

Combining nuclear magnetic resonance (NMR), dynamic light scattering (DLS), and m s long all-atom simulations with two million particles, we establish a delicate correlation between increased side chain organization of PNIPAm and its collapse in aqueous methanol mixtures. We find that the preferential binding of methanol with PNIPAm side chains, bridging distal monomers along the polymer backbone, results in increased organization. Furthermore, methanol–PNIPAm preferential binding is dominated by hydrogen bonding. Our findings reveal that the collapse of PNIPAm is dominated by enthalpic interactions and that the standard poor solvent (entropic) effects play no major role.


Introduction
Polymer conformations in solvent mixtures often exhibit puzzling and paradoxical behavior.One such phenomenon is co-nonsolvency that occurs when two competing (miscible) good solvents for a polymer are mixed together, as a result the same polymer collapses within intermediate solvent-cosolvent mixing ratios in bulk [1][2][3][4][5][6][7][8][9] and near surfaces. 10,11One popular system that shows co-non-solvency is the conformational behavior of poly(N-isopropylacrylamide) (PNIPAm) in aqueous alcohol mixtures.Even when the phenomenon is usually associated with PNIPAm, [1][2][3][4] the name co-non-solvency was first coined when polystyrene chains were dissolved in a mixture of cyclohexane and DMF solution. 124][15][16][17][18] Furthermore, these systems present both LCST and UCST temperature effects.This suggests that the effect, unlike the common chemical believe, is not ''only'' restricted to polymers exhibiting LCST behavior and thus is independent of specific chemical details.
The microscopic origin of the puzzling coil-globule-coil transition is a matter of intense debate.Here extensive experimental, 1-5,9 theoretical 6,8,[19][20][21] and computer simulations 4,7,8,22 studies have been performed.3][4][5]9 In some cases, also in conjunction with analytical theory 1,5 and parameter estimation for interaction strength. 9On the theory side, ever since the first theoretical work employing the Flory-Huggins theory, 1 several analytical works have been proposed to explain co-non-solvency based on the cooperativity effect, 6 particle based theory, 8 Flory-Huggins lattice model at high polymer concentrations, 20 and the off-lattice statistical model, 21 to name a few.Furthermore, former simulations are mostly limited to a few studies. 4,7,8In this context, using a semigrand canonical molecular dynamics approach, two of us have previously shown that PNIPAm has a significantly higher affinity (or preferential binding) towards methanol than towards water by a factor of B4k B T per monomer. 7This is further exemplified by potential of mean force (PMF) calculations, where a clear preferential interaction of methanol with PNIPAm was observed. 23This indicates that the polymer collapse in miscible good solvents, such as water and methanol for PNIPAm, is dictated by the relative (enthalpic) interaction strengths between methanol-PNIPAm and water-PNIPAm.Recently, it has been shown that the effect of co-non-solvency can also be observed in tertiary butyl alcohol in methanol-water mixtures, that is driven by enthalpic interactions. 24Furthermore, the enthalpically driven collapse of PNIPAm is against the common understanding of standard poor solvent collapse of LCST polymers, where solvent entropy gain plays a crucial role in polymer collapse.Furthermore, because a wide variety of polymers show the co-non-solvency phenomenon when dissolved in appropriate mixtures of solvents, a unified concept of co-non-solvency was proposed within a generic approach. 8,19Therefore, simple Lennard-Jones (LJ) interactions between monomer-(co)solvent are sufficient to explain the co-non-solvency effect at constant temperature 8,19 and more effects, 25,26 while ignoring all chemical details that often only contribute to a mere numerical prefactor.However, LJ representation of a (co)solvent bead may represent several methanol and/or water molecules. 8his journal is © The Royal Society of Chemistry 2016 What makes the co-non-solvency of PNIPAm an interesting and puzzling effect is that concepts known from the conventional polymer science often are insufficient to describe the phenomenon.For example: (1) a polymer collapses when the solvent quality remains good or even gets increasingly better by the addition of better cosolvent, 7,8 (2) standard poor solvent collapse or entropic effects are irrelevant, (3) the phenomenon is independent of any underlying LCST or UCST temperature behavior of polymers, [13][14][15][16][17][18][19] and ( 4) is driven by large (local) concentration fluctuations of different solvent components near the polymer, 19 making the mean-field description unsuitable.In this work, we revisit the phenomenon of co-non-solvency of PNIPAm in aqueous methanol mixtures by combining nuclear magnetic resonance (NMR), dynamic light scattering (DLS) and ms long, two million particles all-atom simulations.Note that we have used all-atom simulation instead of a semi-grand canonical setup. 7This is because in the case of a good solvent chain, the chain extension covers almost a full simulation domain consisting of B30 nm box boundary.To better correlate the different approaches, we match the polymer length N l to be similar, in the unit of persistence length l p , which was chosen as N l B 100l p .Our results provide experimental support for the claims presented in the four points mentioned above.
The remainder of the paper is organized as follows: in Section 2 we briefly state the methodology for simulations, material synthesis and experimental measurements.Section 3 presents results and discussion.Finally we draw our conclusions in Section 4.

Nuclear magnetic resonance measurements
The H-NMR experiments 27,28 were measured with a 5 mm triple resonance TXI 1 H/ 13 C/ 15 N probe equipped with a z-gradient on an 850 MHz Bruker AVANCE III system.For proton spectra, 128 transients were used with a 9.5 ms long 901 pulse and a 17 600 Hz spectral width together with a recycling delay of 5 s.The temperature was regulated at 298.3 K and calibrated with a standard 1 H methanol NMR sample using the Topspin 3.1 software (Bruker).Temperature was controlled with a VTU (variable temperature unit) and an accuracy of AE0.1 K.
Diffusion Ordered NMR Spectroscopy (DOSY-NMR) experiments were performed with a gradient strength of 5.350 G mm À1 on a Bruker Avance-III 850 NMR Spectrometer.The gradient strength of probes was calibrated by using a sample of 2 H 2 O/ 1 H 2 O at a defined temperature and compared with the theoretical diffusion coefficient of 2 H 2 O/ 1 H 2 O (values taken from Bruker diffusion manual) at 298.3 K.
The diffusion delay time (D), Bruker term D 2 O was optimized for the TXI probe for 60 ms, while the gradient pulse length was kept at 1.6 ms.The optimization was realized by comparing the remaining intensity of the signals at 2% and 95% gradient strengths.The intensity loss of the echo was in the range of 90%.Using longer diffusion time causes a loss of signal intensity (from the echo) due to a short spin lattice relation time T 1 , which was determined with the inversion recovery method 27 before diffusion measurements.
The diffusion measurements were done with a 2D DOSY sequence 29 by incrementing in 16 linear steps from 2% to 100% with the TXI probe.The calculation of the diffusion value was automatically done with the mono exponential function: 30 ln where I(G) and I(0) are the intensities of the signals with and without gradient, g the gyromagnetic ratio of the nucleus ( 1 H in these measurements), G is the gradient strength, d the duration of the pulse field gradient (PFG), D the diffusion value in m 2 s À1 and D the ''diffusion time'' between the beginning of the two gradient pulses.The relaxation delay between the scans was 3 s.The 2D sequence for diffusion measurement used double stimulated echo with three spoil gradients for convection compensation and with an eddy current delay of 5 ms for reduction 31,32 (acronym Bruker pulse program: dstebpgp3s).The spin-spin relaxation times T 2 were obtained via the CPMG method (Carr-Purcell-Meiboom-Gill) 33,34 using eight different increments.In the experiment, the time between the inversion (1801 pulse) and the read pulse (901 pulse) is incremented in eight steps and the quantitative analysis of the eight integrals (from the eight increments) exponentially fitted.
2.2 Size exclusion chromatography and light scattering 2.2.1 Size exclusion chromatography.For relative molecular weight determinations, a PSS SECcurity Agilent 1260 Infinity Setup (Polymer Standards Service GmbH (PSS)) was used, including a column set from PSS (2Â GRAM 1000, 1Â GRAM 100, particle size 10 mm) maintained at 60 1C, a UV (270 nm) and an RI detector.The eluent was DMF (containing 1 gL À1 LiBr) with a flow rate of 1 mL min À1 .The RI detector signal was used and the molecular weights are relative to linear polystyrene (PS) standards provided by Polymer Standards Service (PSS).Relative M w values agreed within 5% (high M w range) to 10% (low M w range) with absolute M w values obtained by static light scattering (SLS).
2.2.2Light scattering measurements.All light scattering experiments were performed on a commercially available instrument from ALV GmbH consisting of an electronically controlled goniometer and an ALV-5004 multiple tau full-digital correlator (320 channels).A HeNe laser with a wavelength of 632.8 nm and an output power of 25 mW (JDS Uniphase, Type 1145P) was utilized as a light source.All solutions were filtered through Millex-LCR 0.45 mm filters (Merck Millipore), directly into quartz light scattering cuvettes (inner diameter 18 mm), which were cleaned before in a Thurmond apparatus with distilled acetone.All light scattering measurements were carried out, similar to the NMR measurements, at a temperature of 25 1C.For the dynamic light scattering (DLS) experiments, PNIPAm samples were dissolved at a concentration of 1 gL À1 in pure water and in pure methanol.The z-average diffusion coefficients were determined after angular dependent measurements and averaging the apparent diffusion coefficients.The refractive index increment of each sample in methanol was determined with a Michelson interferometer.

PNIPAm samples and synthesis
Three different PNIPAm samples were used for our experiments.These include a commercial PNIPAm from Sigma-Aldrich [cat # 535311-10G].Using light scattering and size exclusion chromatography, we estimate absolute molecular weight M w and poly dispersity (PDI) to be M w (SLS) = 293 kDa, M w (GPC) = 309 kDa, and PDI = 2.69.N l value is estimated using the equation N l = (M w /PDI)/M 0 , where M 0 is the NIPAm formula weight.This corresponds to N l = 962 (sample PNIPAM-963).Furthermore, we have prepared two PNIPAm samples by RAFT polymerization giving chain lengths bracketing that corresponding to N l B 100l p to better correlate all-atom simulations and experimental observations (details will be described in the next section).

All-atom simulations
We employ all atom molecular dynamics simulations using GROMACS package. 35We use the Gromos96 force field 36 for methanol, the SPC/E water model 37 and the force field parameters for PNIPAm are taken from ref. 4. The temperature is set to 298 K using velocity rescaling with a coupling constant 0.5 ps. 38The electrostatics are treated using Particle Mesh Ewald. 39The interaction cutoff is chosen as 1.0 nm.The time step for the simulations is set to 2 fs.Initial equilibration of every configuration is performed at ambient pressure for 20 ns, where the pressure coupling is done using a Berendsen barostat 40 with a coupling time of 0.5 ps.The final configuration of constant pressure simulations, with equilibrated density, is used for canonical simulations.
We choose a PNIPAm chain of length N l = 256 that corresponds to B100l p , with l p being the persistence length of a PNIPAm chain.Note that an atactic PNIPAm chain has l p B 2-3 monomers.We choose five different methanol mole fractions x m in MD simulations.In Table 1 we present our system parameters.The production runs are performed for at least 1 ms long MD trajectory each.During the production run observables such as gyration radius R g and structure factor S(q) are calculated.In Fig. 1 we present the time evolution of R g for three different x m .It can be seen that the structure is rather stable over long simulation time scales.

PNIPAm conformation in aqueous methanol revisited
We first revisit PNIPAm conformation in water and methanol mixtures.In Fig. 2 we show R g as a function of x m .We have also included the data from the earlier generic simulations corresponding to N l = 100l p . 8It can be appreciated that, by matching l p between the all-atom and previous generic simulations, we obtain a very good (almost quantitative) agreement.Furthermore, Fig. 2 also suggests that 1s B 1 nm.This can be rationalized as follows:  R h is then translated into R g using the expressions R g = 1.5R h for coil conformations and R g = (3/5) 1/2 R h for globule. 43,44Note that we present experimental data of R g for the well defined collapsed structures and expanded chain conformations.Near the transition region 0.25 o x m o 1.0, we do not present data because of the nontrivial relation between R h and R g .
Having shown the results for the collapse-swelling-collapse transition, we now want to briefly explain the origin of this reentrant transition.In this context, it was previously shown that the initial collapse at lower x m values are due to preferential binding of methanol molecules with the PNIPAm chain. 7Therefore, when a small amount of methanol molecules is added, these molecules try to bind to more than one monomer acting as sticky contacts between distal monomers along the polymer backbone.This tendency leads to the formation of segmental loops along the backbone initiating the process of polymer collapse.When the x m exceeds a certain concentration, so that the system can overcome the solvent translational entropy, the polymer re-opens after complete decoration of the polymer with methanol molecules.In this case, conformational entropy contributes to a logarithmic correction 19 and, therefore, only leads to a weak effect on describing the overall phenomenon, which is otherwise dominated by enthalpy.The methanol molecules forming enthalpically driven sticky contacts between distal monomer units were termed as bridging cosolvent and their fractions defined as f B . 8 The analytical expression for f B 19 can be converted into gyration radius using an argument using the formulation 19 based on, 43,44 where R 0 g is the gyration radius for x m = 0.Here V is the magnitude of the negative excluded volume À|V|, which can be estimated from the simulations and analytical theory. 8Note that f B gives the direct measure of V using the relation V = 100f B (x c ).Furthermore, an analytical expression can be derived for f B using a Langmuir like adsorption isotherm taking into account the competitive displacement of both solvent and cosolvent. 8In Fig. 2 we also include R g estimated using eqn (2).It can be seen that the experimental, all-atom simulation and generic simulation data can be well described by the analytical theory.Suggesting a reasonable correlation between the phenomena of co-non-solvency and the bridging scenario proposed earlier. 8ig. 3 shows a scaling plot of R g obtained from all-atom simulations and from NMR and DLS measurements of the three polymer samples included in this work.For comparison we have also included data taken from the published literature. 3,7It can be appreciated that the data from different N l obtained from different methods falls within the universal scaling law, further showing a nice quantitative agreement between different methods.For the collapsed structure at x m = 0.1 and 0.25, we present the data calculated using nuclear magnetic resonance (NMR).For comparison, we also include data from our previous generic simulations 8 with N l B 100l p and an analytical expression presented in eqn (2).Arrows are the indicative of the corresponding y-axis of the corresponding data set.

Good solvent collapse
One of the most intriguing features of this phenomenon is that PNIPAm collapses in good solvent.This makes the polymer conformation completely decoupled from the thermodynamic solvent quality (see chemical potential data in ref. 7).Therefore, we expect the chains to maintain a self avoiding walk statistics between the bridging points, i.e. an exponent of n = 3/5.A quantity that best describes the polymer conformation is the static structure factor S(q), which usually requires very long simulation trajectories and also long chain lengths.In this context, thus far all-atom simulations are mostly limited to oligomer systems (consisting of N l = 20-40 to 8-16l p ).In this work, we calculate S(q) for a collapsed chain at x m = 0.1 over last 0.5 ms long all-atom MD run for N l = 256 B 100l p .In Fig. 4 we show S(q), which shows a cross over from an approximate q À5/3 scaling between 4 nm À1 r q r 10 nm À1 (or corresponding length scale of 1.6 nm Z l Z 0.62 nm) to a q À4 scaling for q r 4 nm À1 .This suggests that -while a PNIPAm chain remains globally collapsed, it consists of rather large good solvent blobs with typical sizes of l B 1.3 nm.

Folded PNIPAm structure and side chain organization
We have used solution NMR to monitor side chain dynamics of PNIPAm as a function of x m (predeuterated MeOD in D 2 O).Typical 1 H NMR spectra are presented in Fig. 5.All the signals decrease in area and broaden with increasing x m between 7.5% r x m r 37.6% and then again increases when x m Z 37.6%.The decrease in area and broadening of the signals indicate slowing dynamics, which can be reflected in spin-spin relaxation times T 2 , as also reported by others as a means to map out polymer collapse. 48The T 2 values for the H a signal plotted in Fig. 6, as a function of x m , give the expected clear signature of the coilglobule-coil scenario as a result of the side chains becoming less mobile when the polymer is in the collapsed state.
The data from the NMR experiments also suggests that the side group rigidity of a PNIPAm chain is because of the fact that a PNIPAm collapses in a way that the inner core of the folded structure is occupied by the side groups.MD simulations also support this claim for a collapsed structure.In this context, as described earlier, when there is preferential binding of PNIPAm with methanol, it is expected to observe methanol encapsulation within a collapsed structure of PNIPAm.Therefore, it is important to monitor if the side chain rich inner core also has methanol molecules sitting in between.For this purpose, we have performed a separate set of NMR experiments to identify preferential binding of the methanol molecules with PNIPAm chains.
Fig. 4 Static structure factor S(q) of a collapsed PNIPAm structure at x m = 0.1.For the calculations of S(q) we only consider the alkane backbone of a PNIPAm chain.

Preferential adsorption of methanol on PNIPAm and solvent intake
To experimentally observe the preferential interaction of the PNIPAm with methanol, we have used a concept where preferential PNIPAm-methanol binding 7 can lead to local aggregation of methanol near the polymer structure that leads to depletion away from the polymer, see a schematic in Fig. 7(a).For this purpose we have prepared a NMR tube separated into two compartments by a membrane filter, see Fig. 7(b).Both compartments contain 15% aqueous methanol, and 1.188 g PNI-PAM was added to the upper compartment.This essentially becomes an osmosis experiment, with quantitative 1 H-NMR spectra collected to monitor changes in solvent composition in the lower compartment, from which solvent uptake by the PNIPAM can be estimated.For all different time dependent proton spectra, the same phase correction, baseline and integration parameters (integral width for the proton and methyl group of methanol) were applied.
In Fig. 7(c) we present the amount of depleted methanol in the lower panel of the NMR tube in Fig. 7(b) that is engulfed by the PNIPAm sample in the upper panel of the NMR tube.The measurements were conducted over sixteen days.It can be seen that the methanol increases by B3% in the polymer system within the first 3-4 days.Beyond 4 days, solvent intake data shows a plateau suggesting no evaporation of methanol from the airtight experimental setup.Note that a rather large polymer concentration is needed in the upper panel of Fig. 7(b) to observe any significant solvent intake.The solvent intake, as observed in our NMR experiments, further supports our earlier claim that preferential adsorption of the methanol with PNIPAm drives the polymer collapse. 7This scenario can be further validated by looking into the potential of mean force (PMF) between PNIPAm-methanol and PNIPAm-water calculated using the umbrella sampling. 23,45,46n Fig. 8 we show PMF between different PNIPAm-(co)solvent pairs, which show a clear signature of preferential enthalpic interactions between PNIPAm-methanol.Furthermore, in a given polymer, the energy density within the solvation volume is dictated not only the interaction energy, but also by the sizes of the solvent (methanol and water in this case).Therefore, it should still be mentioned that the enthalpic interactions (or bridging) are usually not given by a single methanol molecule, rather a few collectively lead to sticky contacts. 8he window of PNIPAm collapse, or the LCST behavior, is strongly dependent on the temperature of the systems. 1,2Here, it is important to mention that -just because a polymer exhibits an LCST behavior, it does not suggest that the polymer collapse at a constant T should also be driven enhanced solvent entropy by the addition of methanol in water.Our arguments are based on the claim that the interaction asymmetry between PNIPAm-methanol and PNIPAm-water dictate PNIPAm collapse in aqueous methanol mixture. 8The smaller the asymmetry, the narrower the window of polymer collapse. 19To elucidate that a PNIPAm collapses because  of interaction asymmetry, we have calculated PMF between NIPAm-water and NIPAm-methanol for T = 278 K.The data is shown in Fig. 9.As expected, at a reduced T the asymmetry in interactions also reduces, making the background binary fluid homogeneous for the polymer.This can not facilitate a polymer to collapse in a binary mixture.Interestingly previous experiments 1,2 have shown that the PNIPAm remains in coil state for T = 278 K, thus showing a nice correlation capturing the temperature effects between the earlier experiments 1,2 and our simulations.

Mechanism of polymer collapse
Lastly we want to comment on the possible mechanism of cosolvent bridging leading to PNIPAm collapse in aqueous methanol mixtures.In this context, it is already presented in the previous section and consistent with previous work, 7 that there is a preferential binding of PNIPAm with methanol.To further illustrate the bridging scenario driven by preferential adsorption, we calculate the number of hydrogen bonds n H-bond /N l between a methanol molecules and a NIPAm monomer using all-atom simulations.In Fig. 10(a) we show n H-bond /N l as a function of x m .In the range 0.10 r x m r 0.25, NIPAm shows a strong tendency of hydrogen bonding with the methanol molecules that is evident from the maximum deviation from the expected linear behavior. 4Without attempting to describe any specific geometry or arrangement of H-bond donor and acceptors, we can assume the -OH end of the methanol points towards the PNIPAm amide linkage, such that the CH 3 of methanol now forms part of the local solvent accessible surface previously defined by the amide group.The sticky contacts could then be formed between these CH 3 groups and the isopropyl group of a distal NIPAm unit (see Fig. 10(b)), or between multiple CH 3 groups of bound methanol molecules (see Fig. 10(c)).Note that in Fig. 10(b) and (c), for simplicity of representation, we only highlight hydrogen bonds between hydrogen of methanol and the oxygen of NIPAm amide groups that we expect to be most dominant.However, it should be mentioned that there might be several more scenarios.For example, there is a possibilities of bonding between methanol oxygen and hydrogen of the amide group and methanol hydrogen with nitrogen of the amide groups.Furthermore, it should also be mentioned that a single CH 3 interaction with isopropyl group is Bk B T, which may not sound as a large enough interaction strength.However, these methanol mediated sticky contacts are each likely facilitated by a   few methanol molecules making the sticky contact attraction strength of the order of several Bk B T and not only by a single methanol molecule, as simplified in the schematics shown in Fig. 10(b) and (c).Here we want to emphasize that the solvation properties are intimately linked to the energy density within the solvation volume.Hence, not only dictated by the individual interaction strengths, but also related to number of (co)solvent particles within the solvation volume.Therefore, if one can reduces the energy density within the solvation volume, such that the solvent-cosolvent interaction contrast also reduces, one should expect to see narrower window of collapse.Indeed, it had been experimentally observed that for larger alcohols (such as ethanol or propanol) window of collapse reduces by B30-40% in comparison to methanol. 47

Conclusions
We have revisited the co-non-solvency of PNIPAm in aqueous methanol mixtures.For this purpose, we have combined nuclear magnetic resonance and dynamics light scattering experiments with the all-atom molecular dynamics simulations, complementing our earlier studies of generic simulations 8 and analytical theory. 19hese findings strongly support that the initial collapse at lower methanol concentration is due to the methanol-PNIPAm enthalpic (bridging) effects and the reopening at larger methanol concentrations is entropic.Furthermore, preferential PNIPAm-methanol binding leads to increased organization of the PNIPAm side chains, which is intimately linked to the global conformational behavior of PNIPAm in aqueous methanol mixtures.
While we study a specific system of PNIPAm in aqueous methanol mixture, our proposed mechanism of PNIPAm collapse at lower methanol concentrations provide a natural explanation to other phenomena, such as the initial collapse of PNIPAm in aqueous urea mixtures. 48,49Where the collapse was proposed is initiated by hydrogen bonded bridging of urea molecules between two NIPAm monomers that are far along the polymer backbone.
It should also be mentioned that when dealing with polymer physics and/or the thermodynamics of polymer solutions, two key considerations, out of several, are absolutely needed to make any reasonable comparison to experimental data: (1) a polymer should be studied and not an oligomer and (2) time scale of simulations compared to the polymer relaxation time.In this context, our simulations of ms long trajectory of a N l = 100l p chain consisting of two million particles is the largest all-atom simulations performed on this PNIPAm based systems.A reasonably good agreement between all-atom simulations and experiments complementing earlier generic simulations, suggests that the co-non-solvency phenomenon is indeed driven by enthalpy. 8

Fig. 1
Fig. 1 Time evolution of polymer gyration radii R g for three different methanol mole fractions x m .The results are shown for a chain length N l = 256 and at temperature T = 298 K.

Fig. 2
Fig. 2 Gyration radius R g as a function of methanol mole fraction x m .Results are shown for R g obtained from all-atom simulations for a chain of length N l = 256 B 100l p and experimental measurements for N l B 207 B 83l p (sample PNIPAm-207).For the pure water (x m = 0.0) and pure methanol (x m = 1.0) we use the data obtained from dynamic light scattering (DLS).For the collapsed structure at x m = 0.1 and 0.25, we present the data calculated using nuclear magnetic resonance (NMR).For comparison, we also include data from our previous generic simulations8 with N l B 100l p and an analytical expression presented in eqn(2).Arrows are the indicative of the corresponding y-axis of the corresponding data set.

Fig. 3
Fig. 3 Gyration radius R g as a function of chain length N l .Results are shown for coil (at x m = 1.0) and globule (at x m = 0.1) conformations.Data for R g are obtained from different experiments, simulations and also from published work from the literature, as described in the legend.Symbols are respective data and the lines are power law fits shown in the legend.For coil conformation R g p N l 3/5 and for globular conformations R g p N l 1/3 .

Fig. 5
Fig. 5 Left panel shows schematic representation of a monomer of PNIPAm.In NMR experiments we identify the rigidity of H a and H b indicated in this schematic.Right panel presents nuclear magnetic resonance (NMR) spectra highlighting H a and H b hydrogens as indicated in the left panel.Results are shown for different methanol mole fractions x m , starting from pure water x m = 0.0 (or 0%) to pure methanol x m = 1.0 (or 100%).The signal around 3.75 ppm corresponds to H a and H b peak appears around 1.10 ppm.Note that for clear representation of the data, we have aligned the peak positions of H a .We have used Sigma Aldrich sample of PNIPAm-962.

Fig. 6
Fig.6In-plane relaxation time T 2 of H a hydrogen as a function of methanol mole fractions x m .The data is obtained by integrating the intensity peak around 3.75 ppm in Fig.5.

Fig. 7
Fig. 7 Part (a) presents a schematic representation of the concept where local aggregation of cosolvents near a polymer lead to depletion away from the macromolecular structure.Part (b) shows a schematic showing the 5 mm tube with external reference capillary, a sealing function and a paper filter in the center that separates PNIPAM from bulk aqueous methanol mixture at the bottom.Part (c) presents the amount of excess methanol molecules encapsulated by the PNIPAm collapsed sample as measured from the depletion of methanol in the lower part of the NMR tube.Here 3% of methanol is with respect to the initial 15% of the methanol content.

Fig. 8
Fig. 8 Potential of mean force v(r) showing NIPAm-methanol and NIPAm-water interaction strengths for two different pressures.Simulations are performed at a temperature of 298 K and the data is taken from ref. 23.

Fig. 10
Fig. 10 Part (a) shows number of hydrogen bonds n H-bond /N l between a methanol molecule and a monomer of PNIPAm as a function of methanol mole fraction x m .The dashed red line is the linear extrapolation.Data is shown for all-atom simulations of chain length N l = 256 and for the temperature T = 298 K. Parts (b) and (c) present schematic of two possible scenarios of bridging methanol molecules, which we expect to be most relevant.

Table 1
System sizes for the simulations performed in this study.Here, N is the total number of solvent molecules, N w is the number of water molecules, N m is the number of methanol molecules, x m is the methanol mole fraction, L box the equilibrated box length, and chain gyration radius R g