Modelización de los Procesos de Cristalización en Poliolefinas

“Por lo tanto, la cristalización de polímeros se supone, y en las teorías se describe a menudo, como un proceso de múltiples pasos con muchos aspectos físico-químicos y estructurales influyendo en él. Debido a la propia estructura de la cadena, es fácil entender que un proceso que es termodinámicamente forzado a aumentar su ordenamiento local, se vea obstaculizado geométricamente y, por tanto, no puede conducirse a un estado de equilibrio final. Como resultado, se forman habitualmente estructuras de no equilibrio con diferentes características dependiendo de la temperatura, presión, cizallamiento y otros parámetros físico-químicos del sistema”. Estas palabras, pronunciadas recientemente por el profesor Bernhard Wunderlich, uno de los mas relevantes fisico-quimicos que han abordado en las ultimas décadas el estudio del estado físico de las macromoléculas, adelantan lo que de alguna manera se explicita en esta memoria y constituyen el “leitmotiv” de este trabajo de tesis. El mecanismo de la cristalización de polímeros esta aun bajo debate en la comunidad de la física de polímeros y la mayoría de los abordajes experimentales se explican a través de la teoría LH. Esta teoría clásica debida a Lauritzen y Hoffman (LH), y que es una generalización de la teoría de cristalización de una molécula pequeña desde la fase de vapor, describe satisfactoriamente muchas observaciones experimentales aunque esta lejos de explicar el complejo fenómeno de la cristalización de polímeros. De hecho, la formulación original de esta teoría en el National Bureau of Standards, a comienzos de la década de los 70, sufrió varias reformulaciones importantes a lo largo de la década de los 80, buscando su adaptación a los hallazgos experimentales. Así nació el régimen III de cristalización que posibilita la creacion de nichos moleculares en la superficie y que dio pie al paradigma ofrecido por Sadler y col., para justificar los experimentos que se obtenian por “scattering” de neutrones y otras técnicas como la técnica de “droplets” o enfriamiento rapido. Por encima de todo, el gran éxito de la teoría radica en que explica la dependencia inversa entre el tamaño del plegado molecular y el subenfriamiento, definido este ultimo como el intervalo de temperatura que media entre la temperatura de equilibrio y la temperatura de cristalización. El problema concreto que aborda esta tesis es el estudio de los procesos de ordenamiento de poliolefinas con distinto grado de ramificacion mediante simulaciones numéricas. Los copolimeros estudiados en esta tesis se consideran materiales modelo de gran homogeneidad molecular desde el punto de vista de la distribución de tamaños y de ramificaciones en la cadena polimérica. Se eligieron estas poliolefinas debido al gran interes experimental en conocer el cambio en las propiedades fisicas de los materiales dependiendo del tipo y cantidad de comonomero utilizado. Además, son modelos sobre los que existen una ingente cantidad de información experimental, que es algo que preocupa siempre al crear una realidad virtual como es la simulación. La experiencia en el grupo Biophym es que los resultados de simulación deben de tener siempre un correlato mas o menos próximo experimental y ese argumento se maneja a lo largo de esta memoria. Empíricamente, se conoce muy bien que las propiedades físicas de las poliolefinas, en suma dependen del tipo y de la cantidad de ramificaciones que presenta el material polimérico. Sin embargo, tal como se ha explicado no existen modelos teóricos adecuados que expliquen los mecanismos subyacentes de los efectos de las ramas. La memoria de este trabajo es amplia por la complejidad del tema. Se inicia con una extensa introducción sobre los conceptos básicos de una macromolecula que son relevantes para entender el contenido del resto de la memoria. Se definen los conceptos de macromolecula flexible, distribuciones y momentos, y su comportamiento en disolución y fundido con los correspondientes parametros caracteristicos. Se pone especial énfasis en el concepto de “entanglement” o enmaranamiento por considerarse clave a la hora de tratar macromoléculas con una longitud superior a la longitud critica de enmaranamiento. Finaliza esta introducción con una reseña sobre el estado del arte en la simulación de los procesos de cristalización. En un segundo capitulo del trabajo se expone detalladamente la metodología usada en cada grupo de casos. En el primer capitulo de resultados, se discuten los estudios de simulación en disolución diluida para sistemas lineales y ramificados de cadena única. Este caso mas simple depende claramente del potencial de torsión elegido tal como se discute a lo largo del texto. La formación de los núcleos “babys” propuestos por Muthukumar parece que son consecuencia del potencial de torsión, ya que este facilita los estados de torsión mas estables. Así que se propone el análisis de otros potenciales que son igualmente utilizados y los resultados obtenidos sobre la cristalización, discutidos en consecuencia. Seguidamente, en un segundo capitulo de resultados se estudian moleculas de alcanos de cadena larga lineales y ramificados en un fundido por simulaciones atomisticas como un modelo de polietileno. Los resultados atomisticos pese a ser de gran detalle no logran captar en su totalidad los efectos experimentales que se observan en los fundidos subenfriados en su etapa previa al estado ordenado. Por esta razon se discuten en los capítulos 3 y 4 de resultados sistemas de cadenas cortas y largas utilizando dos modelos de grano grueso (CG-PVA y CG-PE). El modelo CG-PE se desarrollo durante la tesis. El uso de modelos de grano grueso garantiza una mayor eficiencia computacional con respecto a los modelos atomísticos y son suficientes para mostrar los fenómenos a la escala relevante para la cristalización. En todos estos estudios mencionados se sigue la evolución de los procesos de ordenamiento y de fusión en simulaciones de relajación isoterma y no isoterma. Como resultado de los modelos de simulación, se han evaluado distintas propiedades fisicas como la longitud de segmento ordenado, la cristalinidad, temperaturas de fusion/cristalizacion, etc., lo que permite una comparación con los resultados experimentales. Se demuestra claramente que los sistemas ramificados retrasan y dificultan el orden de la cadena polimérica y por tanto, las regiones cristalinas ordenadas decrecen al crecer las ramas. Como una conclusión general parece mostrarse una tendencia a la formación de estructuras localmente ordenadas que crecen como bloques para completar el espacio de cristalización que puede alcanzarse a una temperatura y a una escala de tiempo determinada. Finalmente hay que señalar que los efectos observados, estan en concordancia con otros resultados tanto teoricos/simulacion como experimentales discutidos a lo largo de esta memoria. Su resumen se muestra en un capitulo de conclusiones y líneas futuras de investigación que se abren como consecuencia de esta memoria. Hay que mencionar que el ritmo de investigación se ha acentuado notablemente en el ultimo ano de trabajo, en parte debido a las ventajas notables obtenidas por el uso de la metodología de grano grueso que pese a ser muy importante para esta memoria no repercute fácilmente en trabajos publicables. Todo ello justifica que gran parte de los resultados esten en fase de publicación. Abstract “Polymer crystallization is therefore assumed, and in theories often described, to be a multi step process with many influencing aspects. Because of the chain structure, it is easy to understand that a process which is thermodynamically forced to increase local ordering but is geometrically hindered cannot proceed into a final equilibrium state. As a result, nonequilibrium structures with different characteristics are usually formed, which depend on temperature, pressure, shearing and other parameters”. These words, recently written by Professor Bernhard Wunderlich, one of the most prominent researchers in polymer physics, put somehow in value the "leitmotiv "of this thesis. The crystallization mechanism of polymers is still under debate in the physics community and most of the experimental findings are still explained by invoking the LH theory. This classical theory, which was initially formulated by Lauritzen and Hoffman (LH), is indeed a generalization of the crystallization theory for small molecules from the vapor phase. Even though it describes satisfactorily many experimental observations, it is far from explaining the complex phenomenon of polymer crystallization. This theory was firstly devised in the early 70s at the National Bureau of Standards. It was successively reformulated along the 80s to fit the experimental findings. Thus, the crystallization regime III was introduced into the theory in order to explain the results found in neutron scattering, droplet or quenching experiments. This concept defines the roughness of the crystallization surface leading to the paradigm proposed by Sadler et al. The great success of this theory is the ability to explain the inverse dependence of the molecular folding size on the supercooling, the latter defined as the temperature interval between the equilibrium temperature and the crystallization temperature. The main scope of this thesis is the study of ordering processes in polyolefins with different degree of branching by using computer simulations. The copolymers studied along this work are considered materials of high molecular homogeneity, from the point of view of both size and branching distributions of the polymer chain. These polyolefins were selected due to the great interest to understand their structure– property relationships. It is important to note that there is a vast amount of experimental data concerning these materials, which is essential to create a virtual reality as is the simulation. The Biophym research group has a wide experience in the correlation between simulation data and experimental results, being this idea highly alive along this work. Empirically, it is well-known that the physical properties of the polyolefins depend on the type and amount of branches presented in the polymeric material. However, there are not suitable models to explain the underlying mechanisms associated to branching. This report is extensive due to the complexity of the topic under study. It begins with a general introduction to the basics concepts of macromolecular physics. This chapter is relevant to understand the content of the present document. Some concepts are defined along this section, among others the flexibility of macromolecules, size distributions and moments, and the behavior in solution and melt along with their corresponding characteristic parameters. Special emphasis is placed on the concept of "entanglement" which is a key item when dealing with macromolecules having a molecular size greater than the critical entanglement length. The introduction finishes with a review of the state of art on the simulation of crystallization processes. The second chapter of the thesis describes, in detail, the computational methodology used in each study. In the first results section, we discuss the simulation studies in dilute solution for linear and short chain branched single chain models. The simplest case is clearly dependent on the selected torsion potential as it is discussed throughout the text. For example, the formation of baby nuclei proposed by Mutukhumar seems to result from the effects of the torsion potential. Thus, we propose the analysis of other torsion potentials that are also used by other research groups. The results obtained on crystallization processes are accordingly discussed. Then, in a second results section, we study linear and branched long-chain alkane molecules in a melt by atomistic simulations as a polyethylene-like model. In spite of the great detail given by atomistic simulations, they are not able to fully capture the experimental facts observed in supercooled melts, in particular the pre-ordered states. For this reason, we discuss short and long chains systems using two coarse-grained models (CG-PVA and CG-PE) in section 3 and 4 of chapter 2. The CG-PE model was developed during the thesis. The use of coarse-grained models ensures greater computational efficiency with respect to atomistic models and is enough to show the relevant scale phenomena for crystallization. In all the analysis we follow the evolution of the ordering and melting processes by both isothermal and non isothermal simulations. During this thesis we have obtained different physical properties such as stem length, crystallinity, melting/crystallization temperatures, and so on. We show that branches in the chains cause a delay in the crystallization and hinder the ordering of the polymer chain. Therefore, crystalline regions decrease in size as branching increases. As a general conclusion, it seems that there is a tendency in the macromolecular systems to form ordered structures, which can grown locally as blocks, occupying the crystallization space at a given temperature and time scale. Finally it should be noted that the observed effects are consistent with both, other theoretical/simulation and experimental results. The summary is provided in the conclusions chapter along with future research lines that open as result of this report. It should be mentioned that the research work has speeded up markedly in the last year, in part because of the remarkable benefits obtained by the use of coarse-grained methodology that despite being very important for this thesis work, is not easily publishable by itself. All this justify that most of the results are still in the publication phase.

Lipid biomarkers in the Enol Lake (Asturias, Northern Spain): coupled natural and human induced environmental history

The lipid content of three cores from Lake Enol (Picos de Europa National Park, Asturias, Northern Spain) was studied. The n-alkane profiles indicated a major input from terrigenous plants [predominance of high molecular weight (HMW) alkanes] since ca. 1695 AD to the water body, although the uppermost cm revealed a predominance of organic matter (OM) derived from algae, as the most abundant alkane was C17. Three units revealing different environmental conditions were defined. Unit A (ca. 1695–1860 AD) in the lowermost parts of ENO13-10 (< 12 cm) and ENO13-15 (< 28 cm) was identified and was characterized by higher OM input and evidence of minimal degradation (high CPI values, predominance of HMW n-alkanoic acids and good correspondence between the predominant n-alkane and n-alkanoic acid chains). These findings could be linked to the Little Ice Age, when cold and humid conditions may have favored an increase in total organic carbon (TOC) and n-alkane and n-alkanoic acid content (greater terrigenous OM in-wash), and may have also reduced bacterial activity. In Unit B (ca. 1860–1980 AD) the lack of correspondence between the n-alkane and n-alkanoic acid profiles of ENO13-10 (12–4 cm) and ENO13-15 (28–8 cm) suggested a certain preferential microbial synthesis of long chain saturated fatty acids from primary OM and/or bacterial activity, coinciding with a decrease in OM input, which could be linked to the global warming that started in the second half of the 19th century. In ENO13-7 the low OM input (low TOC) was accompanied by some bacterial degradation (predominance ofHMWn-alkanoic acids but with a bimodal distribution) in the lowermost 16–5 cm. Evidence of considerable phytoplankton productivity and microbial activity was especially significant in Unit C (ca. 1980–2013 AD) identified in the uppermost part of all three cores (5 cm in ENO13-7, 4 cm in ENO13-10 and 8 cm in ENO13-15), coinciding with higher concentrations of n-alkanes and n-alkanoic acids, which were considered to be linked to warmer and drier conditions, as well as to greater anthropogenic influence in modern times. Plant sterols, such as b-sitosterol, campesterol and stigmasterol, were significantly present in the cores. In addition, fecal stanols, such as 24-ethylcoprostanol from herbivores, were present, thereby indicating a continuous and significant pollution input derived from these animals since the 17th century, being more important in the last 20 years.

Biomarkers and inorganic proxies in the paleoenvironmental reconstruction of mires: The importance of landscape in Las Conchas (Asturias, Northern Spain)

We determined the lipid distributions (n-alkanes, n-alkan-2-ones, n-alkanoic acids), total organic carbon (TOC), total nitrogen (TN), Ca/Mg and ash content in Las Conchas mire, a 3.2 m deep bryophyte-dominated mire in Northern Spain covering 8000 cal yr BP. Bog conditions developed in the bottom 20 cm of the profile, and good preservation of organic matter (OM) was inferred from n-alkanoic acid distribution, with the exception of the uppermost 20 cm (last ca. 200 yr). Microbial synthesis of long chain saturated fatty acids from primary OM likely produced a dominance of short chain n-alkanoic acids with a bimodal distribution, as well as the lack of correspondence between the n-alkane and n-alkanoic acid profiles in the upper 20 cm. This was accompanied by an increase in ash content, a decrease in TOC and variation in n-alkane ratios, thereby suggesting significant changes in the mire, namely drainage and transformation to a meadow, in the last ca. 200 yr. The distribution of n-alkan-2-ones indicated an increase in bacterial source from the bottom of the record to 94 cm, whereas their distribution in the upper part could be attributed mainly to plant input and/or the microbial oxidation of n-alkanes. The different n-alkane proxies showed variations, which we interpreted in terms of changes in vegetation (Sphagnum vs. non-Sphagnum dominated phases) during the last 8000 cal yr BP. C23 was the most abundant homolog throughout most of the record, thereby suggesting dominant humid conditions alternating with short drier phases. However, such humid conditions were not linked to paleoclimatic variation but rather to geomorphological characteristics: Las Conchas mire, at the base of the Cuera Range, receives continuous runoff—even during drier periods—which is not necessarily accompanied by additional mineral input to peat, producing the development of Sphagnum moss typical of waterlogged ecotopes and damp habitats. Thus, although geochemical proxies indicated an ombrotrophic regime in the mire, geomorphological characteristics may make a considerable contribution to environmental conditions.

High-resolution NMR of encapsulated proteins dissolved in low-viscosity fluids

The majority of known proteins are too large to be comprehensively examined by solution NMR methods, primarily because they tumble too slowly in solution. Here we introduce an approach to making the NMR relaxation properties of large proteins amenable to modern solution NMR techniques. The encapsulation of a protein in a reverse micelle dissolved in a low-viscosity fluid allows it to tumble as fast as a much smaller protein. The approach is demonstrated and validated with the protein ubiquitin encapsulated in reverse micelles prepared in a variety of alkane solvents.

Enhancing the catalytic repertoire of nucleic acids: a systematic study of linker length and rigidity

The incorporation of potentially catalytic groups in DNA is of interest for the in vitro selection of novel deoxyribozymes. A series of 10 C5-modified analogues of 2′-deoxyuridine triphosphate have been synthesised that possess side chains of differing flexibility and bearing a primary amino or imidazole functionality. For each series of nucleotide analogues differing degrees of flexibility of the C5 side chain was achieved through the use of alkynyl, alkenyl and alkyl moieties. The imidazole function was conjugated to these C5-amino-modified nucleotides using either imidazole 4-acetic acid or imidazole 4-acrylic acid (urocanic acid). The substrate properties of the nucleotides (fully replacing dTTP) with Taq polymerase during PCR have been investigated in order to evaluate their potential applications for in vitro selection experiments. 5-(3-Aminopropynyl)dUTP and 5-(E-3-aminopropenyl)dUTP and their imidazole 4-acetic acid- and urocanic acid-modified conjugates were found to be substrates. In contrast, C5-amino-modified dUTPs with alkane or Z-alkene linkers and their corresponding conjugates were not substrates. The incorporation of these analogues during PCR has been confirmed by inhibition of restriction enzyme digestion using XbaI and by mass spectrometry of the PCR products.

Cloning and characterization of ERG25, the Saccharomyces cerevisiae gene encoding C-4 sterol methyl oxidase.

We have cloned the Saccharomyces cerevisiae C-4 sterol methyl oxidase ERG25 gene. The sterol methyl oxidase performs the first of three enzymic steps required to remove the two C-4 methyl groups leading to cholesterol (animal), ergosterol (fungal), and stigmasterol (plant) biosynthesis. An ergosterol auxotroph, erg25, which fails to demethylate and concomitantly accumulates 4,4-dimethylzy-mosterol, was isolated after mutagenesis. A complementing clone consisting of a 1.35-kb Dra I fragment encoded a 309-amino acid polypeptide (calculated molecular mass, 36.48 kDa). The amino acid sequence shows a C-terminal endoplasmic reticulum retrieval signal KKXX and three histidine-rich clusters found in eukaryotic membrane desaturases and in a bacterial alkane hydroxylase and xylene monooxygenase. The sterol profile of an ERG25 disruptant was consistent with the erg25 allele obtained by mutagenesis.