Atomistic simulations of cosolvent/water mixtures

In this thesis, atomistic simulations are performed to investigate hydrophobic solvation and hydrophobic interactions in cosolvent/water binary mixtures. Many cosolvent/water binary mixtures exhibit non-ideal behavior caused by aggregation at the molecular scale level although they are stable and homogenous at the macroscopic scale. Force-field based atomistic simulations provide routes to relate atomistic-scale structure and interactions to thermodynamic solution properties. The predicted solution properties are however sensitive to the parameters used to describe the molecular interactions. In this thesis, a force field for tertiary butanol (TBA) and water mixtures is parameterized by making use of the Kirkwood-Buff theory of solution. The new force field is capable of describing the alcohol-alcohol, water-water and alcohol-water clustering in the solution as well as the solution components’ chemical potential derivatives in agreement with experimental data. With the new force field, the preferential solvation and the solvation thermodynamics of a hydrophobic solute in TBA/water mixtures have been studied. First, methane solvation at various TBA/water concentrations is discussed in terms of solvation free energy-, enthalpy- and entropy- changes, which have been compared to experimental data. We observed that the methane solvation free energy varies smoothly with the alcohol/water composition while the solvation enthalpies and entropies vary nonmonotonically. The latter occurs due to structural solvent reorganization contributions which are not present in the free energy change due to exact enthalpy-entropy compensation. It is therefore concluded that the enthalpy and entropy of solvation provide more detailed information on the reorganization of solvent molecules around the inserted solute. Hydrophobic interactions in binary urea/water mixtures are next discussed. This system is particularly relevant in biology (protein folding/unfolding), however, changes in the hydrophobic interaction induced by urea molecules are not well understood. In this thesis, this interaction has been studied by calculating the free energy (potential of mean force), enthalpy and entropy changes as a function of the solute-solute distance in water and in aqueous urea (6.9 M) solution. In chapter 5, the potential of mean force in both solution systems is analyzed in terms of its enthalpic and entropic contributions. In particular, contributions of solvent reorganization in the enthalpy and entropy changes are studied separately to better understand what are the changes in interactions in the system that contribute to the free energy of association of the nonpolar solutes. We observe that in aqueous urea the association between nonpolar solutes remains thermodynamically favorable (i.e., as it is the case in pure water). This observation contrasts a long-standing belief that clusters of nonpolar molecules dissolve completely in the presence of urea molecules. The consequences of our observations for the stability of proteins in concentrated urea solutions are discussed in the chapter 6 of the thesis.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions?

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions? The collapse of polyions in a poor solvent is a complex system and is an active research subject in the theoretical polyelectrolyte community. The complexity is due to the subtle interplay between hydrophobic effects, electrostatic interactions, entropy elasticity, intrinsic excluded volume as well as specific counter-ion and co-ion properties. Long range Coulomb forces can obscure single molecule properties. The here presented approach is to use just a small amount of screening salt in combination with a very high sample dilution in order to screen intermolecular interaction whereas keeping intramolecular interaction as much as possible (polyelectrolyte concentration cp ≤ 12 mg/L, salt concentration; Cs = 10^-5 mol/L). This is so far not described in literature. During collapse, the polyion is subject to a drastic change in size along with strong reduction of free counterions in solution. Therefore light scattering was utilized to obtain the size of the polyion whereas a conductivity setup was developed to monitor the proceeding of counterion collection by the polyion. Partially quaternized PVP’s below and above the Manning limit were investigated and compared to the collapse of their uncharged precursor. The collapses were induced by an isorefractive solvent/non-solvent mixture consisting of 1-propanol and 2-pentanone, with nearly constant dielectric constant. The solvent quality for the uncharged polyion could be quantified which, for the first time, allowed the experimental investigation of the effect of electrostatic interaction prior and during polyion collapse. Given that the Manning parameter M for QPVP4.3 is as low as lB / c = 0.6 (lB the Bjerrum length and c the mean contour distance between two charges), no counterion binding should occur. However the Walden product reduces with first addition of non solvent and accelerates when the structural collapse sets in. Since the dielectric constant of the solvent remains virtually constant during the chain collapse, the counterion binding is entirely caused by the reduction in the polyion chain dimension. The collapse is shifted to lower wns with higher degrees of quaternization as the samples QPVP20 and QPVP35 show (M = 2.8 respectively 4.9). The combination of light scattering and conductivity measurement revealed for the first time that polyion chains already collect their counter ions well above the theta-dimension when the dimensions start to shrink. Due to only small amounts of screening salt, strong electrostatic interactions bias dynamic as well as static light scattering measurements. An extended Zimm formula was derived to account for this interaction and to obtain the real chain dimensions. The effective degree of dissociation g could be obtained semi quantitatively using this extrapolated static in combination with conductivity measurements. One can conclude the expansion factor a and the effective degree of ionization of the polyion to be mutually dependent. In the good solvent regime g of QPVP4.3, QPVP20 and QPVP35 exhibited a decreasing value in the order 1 > g4.3 > g20 > g35. The low values of g for QPVP20 and QPVP35 are assumed to be responsible for the prior collapse of the higher quaternized samples. Collapse theory predicts dipole-dipole attraction to increase accordingly and even predicts a collapse in the good solvent regime. This could be exactly observed for the QPVP35 sample. The experimental results were compared to a newly developed theory of uniform spherical collapse induced by concomitant counterion binding developed by M. Muthukumar and A. Kundagrami. The theory agrees qualitatively with the location of the phase boundary as well as the trend of an increasing expansion with an increase of the degree of quaternization. However experimental determined g for the samples QPVP4.3, QPVP20 and QPVP35 decreases linearly with the degree of quaternization whereas this theory predicts an almost constant value.

Modeling and simulations of some anisotropic soft-matter systems: from biaxial to chiral materials

We have modeled various soft-matter systems with molecular dynamics (MD) simulations. The first topic concerns liquid crystal (LC) biaxial nematic (Nb) phases, that can be possibly used in fast displays. We have investigated the phase organization of biaxial Gay-Berne (GB) mesogens, considering the effects of the orientation, strength and position of a molecular dipole. We have observed that for systems with a central dipole, nematic biaxial phases disappear when increasing dipole strength, while for systems characterized by an offset dipole, the Nb phase is stabilized at very low temperatures. In a second project, in view of their increasing importance as nanomaterials in LC phases, we are developing a DNA coarse-grained (CG) model, in which sugar and phosphate groups are represented with Lennard-Jones spheres, while bases with GB ellipsoids. We have obtained shape, position and orientation parameters for each bead, to best reproduce the atomistic structure of a B-DNA helix. Starting from atomistic simulations results, we have completed a first parametrization of the force field terms, accounting for bonded (bonds, angles and dihedrals) and non-bonded interactions (H-bond and stacking). We are currently validating the model, by investigating stability and melting temperature of various sequences. Finally, in a third project, we aim to explain the mechanism of enantiomeric discrimination due to the presence of a chiral helix of poly(gamma-benzyl L-glutamate) (PBLG), in solution of dimethylformamide (DMF), interacting with chiral or pro-chiral molecules (in our case heptyl butyrate, HEP), after tuning properly an atomistic force field (AMBER). We have observed that DMF and HEP molecules solvate uniformly the PBLG helix, but the pro-chiral solute is on average found closer to the helix with respect to the DMF. The solvent presents a faster isotropic diffusion, twice as HEP, also indicating a stronger interaction of the solute with the helix.

Thermodynamic and mechanical properties of polymer-solvent systems

This work presents the results of theoretical and experimental characterization of thermodynamic, mechanical and transport properties in polymer solvent systems. The polymer solvent pairs considered ranged to those in which the polymer is rubbery, to those in which the initially glassy polymeric matrix is plasticized by the action of the low molecular weight species. Advanced Equation of State models have been adopted for thermodynamic modeling,along with a rigorous procedure that enables to extend their applicability to the non equilibrium, glassy region. Mass sorption kinetics had been modeled with phenomenological models and with advanced kinetic models.

Complex colloids by the emulsion solvent evaporation process

In this thesis, different complex colloids were prepared by the process of solvent evaporation from emulsion droplets (SEED). The term “complex” is used to include both an addressable functionality as well as the heterogeneous nature of the colloids.Firstly, as the SEED process was used throughout the thesis, its mechanism especially in regard to coalescence was investigated,. A wide variety of different techniques was employed to study the coalescence of nanodroplets during the evaporation of the solvent. Techniques such as DLS or FCS turned out not to be suitable methods to determine droplet coalescence because of their dependence on dilution. Thus, other methods were developed. TEM measurements were conducted on mixed polymeric emulsions with the results pointing to an absence of coalescence. However, these results were not quantifiable. FRET measurements on mixed polymeric emulsions also indicated an absence of coalescence. Again the results were not quantifiable. The amount of coalescence taking place was then quantified by the application of DC-FCCS. This method also allowed for measuring coalescence in other processes such as the miniemulsion polymerization or the polycondensation reaction on the interface of the droplets. By simulations it was shown that coalescence is not responsible for the usually observed broad size distribution of the produced particles. Therefore, the process itself, especially the emulsification step, needs to be improved to generate monodisperse colloids.rnThe Janus morphology is probably the best known among the different complex morphologies of nanoparticles. With the help of functional polymers, it was possible to marry click-chemistry to Janus particles. A large library of functional polymers was prepared by copolymerization and subsequent post-functionalization or by ATRP. The polymers were then used to generate Janus particles by the SEED process. Both dually functionalized Janus particles and particles with one functionalized face could be obtained. The latter were used for the quantification of functional groups on the surface of the Janus particles. For this, clickable fluorescent dyes were synthesized. The degree of functionality of the polymers was found to be closely mirrored in the degree of functionality of the surface. Thus, the marriage of click-chemistry to Janus particles was successful.Another complex morphology besides Janus particles are nanocapsules. Stimulus-responsive nanocapsules that show triggered release are a highly demanding and interesting system, as nanocapsules have promising applications in drug delivery and in self-healing materials. To achieve heterogeneity in the polymer shell, the stimulus-responsive block copolymer PVFc-b-PMMA was employed for the preparation of the capsules. The phase separation of the two blocks in the shell of the capsules led to a patchy morphology. These patches could then be oxidized resulting in morphology changes. In addition, swelling occurred because of the hydrophobic to hydrophilic transition of the patches induced by the oxidation. Due to the swelling, an encapsulated payload could diffuse out of the capsules, hence release was achieved.The concept of using block copolymers responsive to one stimulus for the preparation of stimulus-responsive capsules was extended to block copolymers responsive to more than one stimulus. Here, a block copolymer responsive to oxidation and a pH change as well as a block copolymer responsive to a pH change and temperature were studied in detail. The release from the nanocapsules could be regulated by tuning the different stimuli. In addition, by encapsulating stimuli-responsive payloads it was possible to selectively release a payload upon one stimulus but not upon the other one.In conclusion, the approaches taken in the course of this thesis demonstrate the broad applicability and usefulness of the SEED process to generate complex colloids. In addition, the experimental techniques established such as DC-FCCS will provide further insight into other research areas as well.

The contribution of the Finite Volume Point Dilution Method to the determination of solute mass flux in an alluvial aquifer.

Groundwater represents one of the most important resources of the world and it is essential to prevent its pollution and to consider remediation intervention in case of contamination. According to the scientific community the characterization and the management of the contaminated sites have to be performed in terms of contaminant fluxes and considering their spatial and temporal evolution. One of the most suitable approach to determine the spatial distribution of pollutant and to quantify contaminant fluxes in groundwater is using control panels. The determination of contaminant mass flux, requires measurement of contaminant concentration in the moving phase (water) and velocity/flux of the groundwater. In this Master Thesis a new solute flux mass measurement approach, based on an integrated control panel type methodology combined with the Finite Volume Point Dilution Method (FVPDM), for the monitoring of transient groundwater fluxes, is proposed. Moreover a new adsorption passive sampler, which allow to capture the variation of solute concentration with time, is designed. The present work contributes to the development of this approach on three key points. First, the ability of the FVPDM to monitor transient groundwater fluxes was verified during a step drawdown test at the experimental site of Hermalle Sous Argentau (Belgium). The results showed that this method can be used, with optimal results, to follow transient groundwater fluxes. Moreover, it resulted that performing FVPDM, in several piezometers, during a pumping test allows to determine the different flow rates and flow regimes that can occurs in the various parts of an aquifer. The second field test aiming to determine the representativity of a control panel for measuring mass flus in groundwater underlined that wrong evaluations of Darcy fluxes and discharge surfaces can determine an incorrect estimation of mass fluxes and that this technique has to be used with precaution. Thus, a detailed geological and hydrogeological characterization must be conducted, before applying this technique. Finally, the third outcome of this work concerned laboratory experiments. The test conducted on several type of adsorption material (Oasis HLB cartridge, TDS-ORGANOSORB 10 and TDS-ORGANOSORB 10-AA), in order to determine the optimum medium to dimension the passive sampler, highlighted the necessity to find a material with a reversible adsorption tendency to completely satisfy the request of the new passive sampling technique.

Acid, Silver, and Solvent-Free Gold-Catalyzed Hydrophenoxylation of Internal Alkynes

A range of arylgold compounds have been synthesized and investigated as single-component catalysts for the hydrophenoxylation of unactivated internal alkynes. Both carbene and phosphine-ligated compounds were screened as part of this work, and the most efficient catalysts contained either JohnPhos or IPr/SIPr. Phenols bearing either electron-withdrawing or electron-donating groups were efficiently added using these catalysts. No silver salts, acids, or solvents were needed for the catalysis, and either microwave or conventional heating afforded moderate to excellent yields of the vinyl ethers.

Excited state properties of formamide in water solution: An ab initio study

We present ab initio quantum calculation of the optical properties of formamide in vapor phase and in water solution. We employ time dependent density functional theory for the isolated molecule and many-body perturbation theory methods for the system in solution. An average over several molecular dynamics snapshots is performed to take into account the disorder of the liquid. We find that the excited stateproperties of the gas-phase formamide are strongly modified by the presence of the water solvent: the geometry of the molecule is distorted and the electronic and optical properties are severely modified. The important interaction among the formamide and the water molecules forces us to use fully quantum methods for the calculation of the excited stateproperties of this system. The excitonic wave function is localized both on the solute and on part of the solvent.

Scale Up of PLA Nanoparticle Synthesis by Solvent Displacement and Functionalization of the Particles with a Silica Shell for Targeted Drug Delivery Applications

Biodegradable polymer nanoparticles have the properties necessary to address many of the issues associated with current drug delivery techniques including targeted and controlled delivery. A novel drug delivery vehicle is proposed consisting of a poly(lactic acid) nanoparticle core, with a functionalized, mesoporous silica shell. In this study, the production of PLA nanoparticles is investigated using solvent displacement in both a batch and continuous manner, and the effects of various system parameters are examined. Using Pluronic F-127 as the stabilization agent throughout the study, PLA nanoparticles are produced through solvent displacement with diameters ranging from 200 to 250 nm using two different methods: dropwise addition and in an impinging jet mixer. The impinging jet mixer allows for easy scale-up of particle production. The concentration of surfactant and volume of quench solution is found to have minimal impact on particle diameter; however, the concentration of PLA is found to significantly impact the diameter mean and polydispersity. In addition, the stability of the PLA nanoparticles is observed to increase as residual THF is evaporated. Lastly, the isolated PLA nanoparticles are coated with a silica shell using the Stöber Process. It is found that functionalizing the silica with a phosphonic silane in the presence of excess Pluronic F-127 decreases coalescence of the particles during the coating process. Future work should be conducted to fine-tune the PLA nanoparticle synthesis process by understanding the effect of other system parameters and in synthesizing mesoporous silica shells.

Solvatochromic and Thermochromic Studies of Tryptanthrin

Solvatochromism and thermochromism describe how a solvent or environment affects the photophysical behavior of a photoluminescent solute. The most common use of solvatochromism is as a probe in which the polarity of a solvent in which a solvatochromic solute is dissolved can be spectroscopically measured. Solvatochromic and thermochromic studies of tryptanthrin in several different solvents are reported. Absorption and corrected emission spectra for tryptanthrin at ~10-6 M concentrations are reported in four aprotic and nine alcoholic solvents. The absorption spectra are relatively unaffected by changes in solvent polarity and by differences in the hydrogen bonding ability of the alcoholic solvents. The emission spectra are much more affected by changes in solvent polarity and hydrogen bonding ability. In aprotic solvents, emission energy decreases and emission intensity increases with increasing solvent polarity. In the alcoholic solvents, emission energy also decreases with increasing solvent polarity. However, emission intensity for the alcoholic solvents varies significantly from the aprotic solvents over similar polarity ranges. This suggests that in the alcoholic solvents, hydrogen bonding ability correlates better than polarity to emission energy and intensity trends. The absorption and emission data in the aprotic solvents were also used to estimate the ground and emitting excited state dipole moments for tryptanthrin. The value obtained for the ground state dipole moment (2.37 D) agrees with theoretical results (2.06 D) and a previously reported experimental value (2.0 D). Attempts to explain previously reported results and conclusions with respect to the solvatochromic behavior of the aromatic carbonyls fluorenone and benzo(b)fluorenone were explored in an attempt to understand the solvatochromic response of tryptanthrin. Such attempts include models dependent on non-radiative decay pathways like intersystem crossing, internal conversion, and hydrogen bonding interactions.