Molecular Approach To Understand the Tacticity Effects on the Hydrophilicity of Poly(N-isopropylacrylamide): Solubility of Dimer Model Compounds in Water

RISMiCal: A software package to perform fast RISM/3D-RISM calculations

Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.

Effect of Tacticity Sequence of the Poly(N-isopropylacrylamide) Oligomer on Phase Transition Behavior in Aqueous Solution

The tacticity of poly(N-isopropylacrylamide) (PNIPAM) has a strong impact on the lower critical solution temperature (LCST) in aqueous solution. The sequence of meso diads (m) and racemo diads (r) further contributes to such an effect. In this work, the phase transition behaviors of poly(N-isopropylacrylamide) pentamers with four kinds of sequences, i.e., rrmm, rmmr, mrrm, and rmrm, in water were studied applying replica exchange molecular dynamics with a modified OPLS/AA force field. The difference in local component concentration in the system was used as an order parameter to quantitatively describe the phase separation extent. It was found that the phase separation degree of rrmm and rmmr is higher than that of mrrm and rmrm at the same temperature. The LCSTs of rrmm and rmmr are lower than those of mrrm and rmrm. The radial distribution function and hydrogen bond analysis revealed that the average values of hydrogen bonds between pentamers for rrmm and rmmr are greater than those of mrrm and rmrm, whereas the average values of hydrogen bonds between pentamers and water for rrmm and rmmr are less than those of mrrm and rmrm. It was demonstrated that the isotactic triad (mm) plays an important role in the thermosensitive behaviors of the PNIPAM pentamer. The increase of isotactic triad (mm) content in the PNIPAM chain promotes the formation of intermolecular hydrogen bonds between amide and amide and leads to a higher aggregation of the pentamer with the sequence of rrmm or rmmr. Finally, the effect of the isotactic triad was qualitatively explained with the mean-field theory.

Phase-Separated Nanodroplets Formed below the Cloud Point for the Aqueous Solution of Stereo-Controlled Poly(N-isopropylacrylamide)

The relationship between dehydration of polymer chains and the nanodroplet formation through the macroscopic liquid-liquid phase separation (LLPS) has been investigated for the aqueous solution of stereocontrolled poly(N-isopropylacrylamide) (PNiPAm) and poly(N-diethylacrylamide) (PNdEAm). The fluorescent probe method reveals that the temperature range of dehydration for PNiPAm chains is much narrower than that for PNdEAm. The sharp dehydration of polymer chains may give rise to the characteristic thermoresponsive behavior of PNiPAm in water. For meso-rich PNiPAms, the dehydration point (T_dh), which is defined as the temperature where the single chains start assembling in the solution, locates far from the cloud point (T_c). That is, the dehydration of the chain occurs antecedently before the system undergoes a macroscopic LLPS. For PNdEAm, however, the dissociation between T_dh and T_c is not found. For the aqueous solution of PNiPAm with 52% of the meso content, the fluorescence correlation spectroscopy has revealed that nano-order droplets (ca. 45 nm of the hydration radius) are stabilized in the intermediate state between T_dh and T_c. The sharp dehydration of PNiPAm chains may enable an acute condensation of polymers in droplets, causing a viscoelastic hindrance in the coalescence of droplets.

Molecular Chirality and Cloud Activation Potentials of Dimeric α-Pinene Oxidation Products

The surface activity of ten atmospherically relevant α-pinene-derived dimers having varying terminal functional groups and backbone stereochemistry is reported. We find ∼10% differences in surface activity between diastereomers of the same dimer, demonstrating that surface activity depends upon backbone stereochemistry. Octanol-water (K_OW) and octanol-ammonium sulfate partitioning coefficient (K_OAS) measurements of our standards align well with the surface activity measurements, with the more surface-active dimers exhibiting increased hydrophobicity. Our findings establish a link between molecular chirality and cloud activation potential of secondary organic aerosol particles. Given the diurnal variations in enantiomeric excess of biogenic emissions, possible contributions of such a link to biosphere:atmosphere feedbacks as well as aerosol particle viscosity and phase separation are discussed.

Microanalysis of Single Poly(N-isopropylacrylamide) Droplet Produced by an Optical Tweezer in Water: Isotacticity Dependence of Growth and Chemical Structure of the Droplet

Thermoresponsive phase separation mechanisms of aqueous poly(N-isopropylacrylamide) (PNIPAM) solutions were investigated using an optical tweezer combined with a Raman microspectroscope. A near-infrared laser beam (λ = 1064 nm) was focused into the solution to produce and trap a single polymer microdroplet under an optical microscope. The laser beam played two important roles: The first role is to locally heat the solution to induce phase separation in which numerous polymer microdroplets are generated around the focus, while the second one is to collect these microdroplets. Eventually, a single polymer droplet was stably produced and trapped at the focus. Our method enabled us to perform two types of microanalysis for the droplet. Analysis I is real-time monitoring the growth of the polymer droplets by which we can determine the growth rate of droplets. Analysis II is Raman microspectroscopy to reveal chemical components of the droplets. By means of these two analyses, we revealed important phase separation mechanisms in terms of stereoregularity (isotacticity) dependence. From analysis I, we show that droplet growth is governed by the Ostwald ripening mechanism and the growth is accelerated by increasing the isotacticity. From analysis II, we show that the gelation is promoted in the droplet (physical gel formation) with increasing isotacticity. Our technique should be a versatile tool to explore liquid-liquid phase separation mechanisms for various binary solution systems.

Stereoselectivity in Atmospheric Autoxidation

We show that the diastereomers of hydroxy peroxy radicals formed from OH and O₂ addition to C2 and C3, respectively, of crotonaldehyde (CH₃CHCHCHO) undergo gas-phase unimolecular aldehydic hydrogen shift (H-shift) chemistry with rate coefficients that differ by an order of magnitude. The stereospecificity observed here for crotonaldehyde is general and will lead to a significant diastereomeric-specific chemistry in the atmosphere. This enhancement of specific stereoisomers by stereoselective gas-phase reactions could have widespread implications given the ubiquity of chirality in nature. The H-shift rate coefficients calculated using multiconformer transition state theory (MC-TST) agree with those determined experimentally using stereoisomer-specific gas-chromatography chemical ionization mass spectroscopy (GC-CIMS) measurements.

Solution behaviour of poly(N-isopropylacrylamide) stereoisomers in water: a molecular dynamics simulation study

The water affinity of poly(N-isopropylacrylamide), PNIPAM, is tuned by tacticity, since the hydrophobicity rises with the increase of the degree of isotacticity. On the basis of this experimental evidence, atomistic molecular dynamics simulations of pairs of PNIPAM stereoisomers in 1.6% w/w polymer aqueous solution, a condition intermediate between the dilute and semidilute regimes, were carried out to comparatively investigate the solution behaviour and hydration of atactic and isotactic-rich PNIPAMs, both below and above the lower critical solution temperature, LCST. 30-mers with contents of meso dyads, m, of 45% and 59%, built assuming a Bernoullian dyad distribution, are used as models since their stereochemical composition corresponds to that of experimentally characterized PNIPAM stereoisomers. The simulation results at 283 K, below the LCST, show a slight influence of tacticity on the chain size, but a higher propensity for inter-chain association of the meso-dyad-rich system, in agreement with the experimental results. Junctions between chains are formed because of hydrophobic interactions and are stabilized by a layer of hydrogen bonded water molecules, whose mobility is reduced as compared to that observed for the same meso-dyad-rich stereoisomer at infinite dilution. At 323 K, above the LCST, simulations detect both the coil-globule transition and the aggregation of chains. Under these conditions, the influence of tacticity on the characteristics of PNIPAM aggregate is negligible.

Development and Application of a Coarse-Grained Model for PNIPAM by Iterative Boltzmann Inversion and Its Combination with Lattice Boltzmann Hydrodynamics

The polymer poly(N-isopropylacrylamide) (PNIPAM) is studied using a novel combination of multiscale modeling methodologies. We develop an iterative Boltzmann inversion potential of concentrated PNIPAM solutions and combine it with lattice Boltzmann as a Navier-Stokes equation solver for the solvent. We study in detail the influence of the methodology on statics and dynamics of the system. The combination is successful and significantly simpler and faster than other mapping techniques for polymer solution while keeping the correct hydrodynamics. The model can semiquantitatively describe the correct phase behavior and polymer dynamics.

Tacticity-Dependent Interchain Interactions of Poly(N-Isopropylacrylamide) in Water: Toward the Molecular Dynamics Simulation of a Thermoresponsive Microgel

The discovery that the lower critical solution temperature (LCST) of poly(N-Isopropylacrylamide) (PNIPAM) in water is affected by the tacticity opens the perspective to tune the volume phase transition temperature of PNIPAM microgels by changing the content of meso dyads in the polymer network. The increased hydrophobicity of isotactic-rich PNIPAM originates from self-assembly processes in aqueous solutions also below the LCST. The present work aims to detect the characteristics of the pair interaction between polymer chains, occurring in a concentration regime close to the chain overlap concentration, by comparing atactic and isotactic-rich PNIPAM solutions. Using atomistic molecular dynamics simulations, we successfully modelled the increased association ability of the meso-dyad-rich polymer in water below the LCST, and gain information on the features of the interchain junctions as a function of tacticity. Simulations carried out above the LCST display the PNIPAM transition to the insoluble state and do not detect a relevant influence of stereochemistry on the structure of the polymer ensemble. The results obtained at 323 K provide an estimate of the swelling ratio of non-stereocontrolled PNIPAM microgels which is in agreement with experimental findings for microgels prepared with low cross-linker/monomer feed ratios. This study represents the first step toward the atomistic modelling of PNIPAM microgels with a controlled tacticity.

Assessing the Efficacy of Poly(N-isopropylacrylamide) for Drug Delivery Applications Using Molecular Dynamics Simulations

All-atom molecular dynamic simulations (AA-MD) are performed for aqueous solutions of hydrophobic drug molecules (phenytoin) with model polymer excipients, namely, (1) N-isopropylacrylamide, (pNIPAAm), (2) pNIPAAm-co-acrylamide (Am), and (3) pNIPAAm-co-dimethylacrylamide (DMA). After validating the force field parameters using the well-known lower critical solution behavior of pNIPAAm, we simulate the polymer-drug complex in water and its behavior at temperatures below (295 K) and above the LCST (310 K). Using radial distribution functions, we find that there is an optimum comonomer molar fraction of around 20-30% DMA at which interaction with phenytoin drug molecules is strongest, consistent with recent experimental findings. The results provide evidence that molecular simulations are able to provide guidance in the optimization of novel polymer excipients for drug release.

A review of stereochemical implications in the generation of secondary organic aerosol from isoprene oxidation

The atmospheric reactions leading to the generation of secondary organic aerosol (SOA) from the oxidation of isoprene are generally assumed to produce only racemic mixtures, but aspects of the chemical reactions suggest this may not be the case. In this review, the stereochemical outcomes of published isoprene-degradation mechanisms contributing to high amounts of SOA are evaluated. Despite evidence suggesting isoprene first-generation oxidation products do not contribute to SOA directly, this review suggests the stereochemistry of first-generation products may be important because their stereochemical configurations may be retained through to the second-generation products which form SOA. Specifically, due to the stereochemistry of epoxide ring-opening mechanisms, the outcome of the reactions involving epoxydiols of isoprene (IEPOX), methacrylic acid epoxide (MAE) and hydroxymethylmethyl-α-lactone (HMML) are, in principle, stereospecific which indicates the stereochemistry is predefined from first-generation precursors. The products from these three epoxide intermediates oligomerise to form macromolecules which are proposed to form chiral structures within the aerosol and are considered to be the largest contributors to SOA. If conditions in the atmosphere such as pH, aerosol water content, relative humidity, pre-existing aerosol, aerosol coatings and aerosol cation/anion content (and other) variables acting on the reactions leading to SOA affect the tacticity (arrangement of chiral centres) in the SOA then they may influence its physical properties, for example its hygroscopicity. Chamber studies of SOA formation from isoprene encompass particular sets of controlled conditions of these variables. It may therefore be important to consider stereochemistry when upscaling from chamber study data to predictions of SOA yields across the range of ambient atmospheric conditions. Experiments analysing the stereochemistry of the reactions under varying conditions of the above variables would help elucidate whether there is stereoselectivity in SOA formation from isoprene and if the rates of SOA formation are affected.

Influence of Tacticity on Hydrophobicity of Poly(N-isopropylacrylamide): A Single Chain Molecular Dynamics Simulation Study

Several pieces of experimental evidence show that the water affinity of poly(N-isopropylacrylamide), PNIPAM, decreases with the increase of the degree of isotacticity. To give a rationale to such effect we carried out atomistic molecular dynamics simulations of PNIPAM 30-mers with content of meso dyads, m, of 45% and 59%, assuming a Bernoullian dyad distribution. The single chain behavior of these stereoisomers in water was modeled at 283 and 323 K, i.e., below and above the lower critical solution temperature, LCST. Simulation results show that the dyad composition has influence on size and conformation of the oligomer below the LCST, the isotactic-rich stereoisomer preferring conformations with lower radius of gyration. With reference to the cooperative hydration model of PNIPAM, we analyzed the organization and the dynamics of water in the surroundings of the polymer. Below the LCST the number of hydrogen bonds per water molecule in the first hydration shell and the water surface concentration around PNIPAM are affected by the meso/racemo dyad ratio, showing the higher hydrophobicity of the isotactic-rich system. Above the LCST the subtle effects leading to the modulation of conformation and hydration by means of tacticity are overtaken, and the chain collapse is observed for both systems up to a similar globular state. The comparison of simulation findings of the m 45% stereoisomer with available experimental results of atactic PNIPAM highlights additional molecular details of this macromolecule in aqueous environment. The characteristic time for motion of water molecules in the PNIPAM first hydration shell at 283 K is about 34 ps, in agreement with the experimental value.

phase instability and molecular kinetics provoked by repeated crossing of the demixing transition of PNIPAM solutions

The demixing process of aqueous poly(N-isopropylacrylamide) (PNIPAM) solutions can occur either via a nucleation and growth process or via spinodal decomposition. The ensuing self-assembly, leading to heterogeneous morphologies within the PNIPAM solution, is codetermined by kinetic processes caused by molecular transport. By subjecting PNIPAM solutions to cyclic changes in temperature leading to repeated crossing of the demixing transition, we are able to assess the importance of kinetics as well as of overheating and supercooling of the phase transition within the metastable range delimited by the binodal and spinodal lines. First indications about the location of these stability limits for the low- and high-temperature phases, separated by about 1.6 K, could be gained by detailed kinetic studies of the refractive index. These investigations are made possible due to the novel technique of temperature-modulated optical refractometry.

Infrared spectroscopic study of stereo-controlled poly(N-isopropylacrylamide) with an extended chain conformation induced by adsorption on a gold surface

Poly(N-isopropylacrylamide) (PNiPAM) compounds with various diad tacticities were prepared, and the molecular interaction properties in a thin film deposited on a gold surface were analyzed using infrared spectroscopy. The intramolecular and intermolecular interactions were found to depend on the tacticity, and only atactic (diad ratio 46 %) PNiPAM exhibits poor molecular interaction even in the bulk sample. On the other hand, the same series of compounds dissolved in an acetone solution were spread on a gold surface to form a thin film. In the dissolution process, the polymer molecules are relaxed via solvation, and they are bound to the gold surface by a molecular interaction to form a submonolayer thin film. In the thin film, the molecular interaction with the gold surface via the N-H group was monitored in the infrared spectra only for a nearly isotactic (m = 90) PNiPAM by an apparent shift of the N-H stretching vibration band. This shift was confirmed by changing the degree of hydrophilicity of the gold surface: a larger shift is found on a gold surface with stronger hydrophilicity. As a result, the conformation of a nearly isotactic molecule is found to be extended by the interaction with the gold surface, which works to immobilize the molecule.