Les carbones désordonnés nanoporeux sont des supports efficaces pour le piégeage de polluants y compris à l’état de traces dans les eaux usées. Le phénomène d’adsorption à l’origine de la rétention des molécules est cependant complexe car dépendant d’une multitude de facteurs : structure, morphologie et charge de la surface carbonée d’une part,taille/forme et polarité de la molécule d’autre part, l’ensemble étant dépendant du pH et de la concentration. Pourune meilleure compréhension du phénomène, il est important de pouvoir étudier séparément certains paramètres.Dans la perspective d’étudier le phénomène d’adsorption en milieu aqueux sur des carbones nanoporeux à structure et morphologie modèle, des structures lamellaires nanoporeuses de type carbone turbostratique ont été générées numériquement en langage C++ avec le calcul de la fonction de distribution radiale ou de paires. L’adsorption gazeuse d’une molécule non polaire ou polaire puis de deux molécules polaires (H2O/CO2) et (H2O/C6H6O)a été simulée par la méthode Grand Canonique Monte Carlo sur ce support modèle (Isotherme d’adsorption,chaleur d’adsorption, densité des molécules adsorbées) en fonction de la température. Les temps de calcul ont été drastiquement diminués en développant des codes parallèles optimisés sous MPI C++. L’influence de la forme etde la distribution en taille des pores a été mise en évidence en simulant l’adsorption sur la structure d’un carbone activé déjà obtenue par reconstruction 3D de type RMC. Enfin, d’un point de vue expérimental, l’intercalation d’ions tetraalkylammonium par voie électrochimique dans des carbones lamellaires (HOPG et graphite) a été explorée en vue d’obtenir des carbones lamellaires nanoporeux (≈1 nm). La structure a été caractérisée par diffraction des rayons X. / Disordered nanoporous carbons are the good materials for capturing pollutants, including traces in wastewater. The phenomenon of adsorption at the origin of the retention of molecules is complex. However, depending on a multitude of factors : structure, morphology and loading of the carbonaceous surface, on the one hand, size/shapeand polarity of the molecule, on the other hand, the whole being dependent on pH and concentration. For a better understanding of the phenomenon, it is important to be able to study some parameters separately. In order to study the phenomenon of adsorption in aqueous medium on nanoporous carbons with structure and model morphology, nanoporous slit structures of turbostratic carbon type were generated numerically in C ++ language with thecalculation of the radial distribution function or pairs. The gas adsorption of a nonpolar or polar molecule and then oftwo polar molecules (H2O/CO2) and (H2O/C6H6O) was simulated by Grand Canonical Monte Carlo method on this model support (adsorption isotherm, adsorption heat, density of adsorbed molecules) as a function of temperature.The runtime has been drastically reduced by developing parallel codes optimized under MPI C ++. The influence of the shape and the pore size distribution was demonstrated by simulating the adsorption on the structure of an activated carbon already obtained by 3D reconstruction of the RMC type. Finally, from an experimental point of view, the intercalation of tetraalkylammonium ions electrochemically in slit carbons (HOPG and graphite) was explored in order to obtain nanoporous lamellar carbons ( ≈1 nm). The structure was characterized by X-ray diffraction.

In recent years, metal-organic frameworks (MOFs) have been identified as promising adsorbents in a number of industrially relevant, yet challenging, separations, including the removal of propane from propane/propylene mixtures and the separation of mixtures of xylene isomers. The highly tuneable nature of MOFs - wherein structures may be constructed from a variety of diverse building blocks – has resulted in the publication of a staggering number of frameworks incorporating a wide range of network topologies, pore shapes and pore diameters. As a result, there are a huge number of candidate adsorbents to consider for a given separation. Molecular simulation techniques allow the identification of those structural features and characteristics of a MOF which exert the greatest influence on the adsorption and separation of the compounds of interest, providing insights which can both guide the selection and accelerate the development of adsorbents for a specific application. The separation of propane/propylene mixtures via adsorption has typically focused on selective adsorption of the olefin, propylene, via specific olefin-adsorbent interactions. These propylene-selective MOFs result in processes which selectively remove the most abundant species in the process stream and are typically characterised by high heats of adsorption, resulting in large adsorption units and adsorbents which are difficult to regenerate. In this work, the capability of MOFs to selectively adsorb propane over propylene is explored, potentially allowing for the design of smaller and more energy-efficient adsorption units. By studying a range of different MOFs as well as carbon-based model pores, it was found that the low-pressure selectivity of the structure is determined by the strength of the electrostatic interaction between propylene and the framework, while the adsorptive preference at industrially-relevant pressures is dominated by the enhanced packing efficiency of propylene over propane. The confinement of C3 molecules, however, may be employed to negate this entropic advantage and guide the development of materials which selectively adsorb propane over propylene. It has recently been reported that the adsorptive preference of a MOF for one xylene isomer over another may be predicted based solely on the pore size distribution of the structure. In this work, the impact of pore size on selectivity was studied systematically in both one-dimensional model pore systems of varying geometries and analogous published MOF structures. The ability of the framework to discriminate between xylene molecules in these systems was found to be determined primarily by the different packing arrangements available to the different isomers – while small pores were found to favour the slimmest of the isomers, larger pores were found to favour the more compact ortho- isomer. Finally, the adsorption and diffusion of xylene isomers in a more complex MOF, UiO-66(Zr), was studied in depth. Simulations were able to correctly predict the previously-reported preference of the MOF for ortho-xylene (oX). The smaller volume of the oX molecule compared to the other isomers was found to be responsible both for an enhanced entropic contribution and higher guest-host interaction energies. The importance of framework flexibility in the diffusion of xylene isomers in UiO-66(Zr) was also explored, with distortion of the structure in response to interaction with adsorbed molecules found to be essential in allowing xylenes to diffuse through the pore space.

Confining walls induce qualitative changes in adsorbed fluids. Among the most intriguing phenomena is hysteresis, where a pore fills with fluid at a greater pressure than it empties. The causes and mechanisms by which this occurs are intensely investigated yet still poorly understood. Ordered mesoporous silicas, recently discovered materials with well-defined pore size distributions, provide an opportunity to deepen our understanding of the fundamental physics of the interaction of fluids with complex solids.In part of this computational investigation we examine idealized pores. In agreement with other recent studies, we find that in 'inkbottle'-shaped pores, where a large cavity is accessible to the bulk fluid only by constrictions, there is no evidence of the long-hypothesized phenomenon of `pore blocking', where the constrictions inhibit fluid desorption from the cavity. We find that even in these simple systems the mechanism of hysteresis depends on pore characteristics, fluid properties and external conditions.For silicas containing cylindrical holes of nearly uniform diameter, such as MCM-41, the state-of-the-art is to consider only a single pore, but the poor qualitative agreement of theoretical with experimental results has improved little as wall representations of increasing sophistication have been developed. Using only a one-dimensional potential, we reproduce features of isotherms, including in the hysteresis region, by averaging over a narrow distribution of pore sizes. The qualitative behavior is shown to be a collective phenomenon not representative of any individual pore. Adding surface roughness and a constriction to the pores yields results quantitatively nearly indistinguishable from experiments.For materials larger than MCM-41, a continuum simulation proves too computationally taxing. Thus, a lattice model with adjustable fineness of site spacing is developed. It is found that a surprisingly low level of fineness is needed for confined systems to closely approximate continuum results. This model is applied to cubically symmetric materials, such as MCM-48 and SBA-16, finding that simulations are able to reproduce much of the qualitative behavior seen experimentally, but the lack of existing knowledge of the nature of silica walls proves to be a limiting factor.

The detrimental effects of rising CO₂ levels on the global climate have made carbon abatement technologies one of the most widely researched areas of recent times. In this thesis, we first present a techno-economic analysis of a novel approach to directly capture CO₂ from air (Air Capture) using highly selective adsorbents. Our process modeling calculations suggest that the monetary cost of Air Capture can be reduced significantly by identifying adsorbents that have high capacities and optimum heats of adsorption. The search for the best performing material is not limited to Air Capture, but is generally applicable for any adsorption-based separation. Recently, a new class of nanoporous materials, Metal-Organic Frameworks (MOFs), have been widely studied using both experimental and computational techniques. In this thesis, we use a combined quantum chemistry and classical simulations approach to predict macroscopic properties of MOFs. Specifically, we describe a systematic procedure for developing classical force fields that accurately represent hydrocarbon interactions with the MIL-series of MOFs using Density Functional Theory (DFT) calculations. We show that this force field development technique is easily extended for screening a large number of complex open metal site MOFs for various olefin/paraffin separations. Finally, we demonstrate the capability of DFT for predicting MOF topologies by studying the effect of ligand functionalization during CuBTC synthesis. This thesis highlights the versatility and opportunities of using multiscale modeling approach that combines process modeling, classical simulations and quantum chemistry calculations to study nanoporous materials for adsorptive separations.

L'adsorption de CO dans la faujasite échangée au CuI et au Na+ a été modélisée à l'aide des approches quantiques (DFT) et classiques (Monte Carlo). Grâce à l'approche DFT, la surface d'énergie potentielle de la faujasite a été explorée. Différents types d'interactions de CO avec les cations ont été identifiés, pour chacune les effets induits par l'adsorption de CO aux niveaux structural et énergétique ont été analysés, et le calcul de la fréquence de vibration de CO a été réalisé. Grâce aux valeurs obtenues, une nouvelle attribution des spectres d'adsorption de CO dans CuY et NaY a été établie. D'un autre côté, grâce aux simulations Monte Carlo dans l'ensemble Grand Canonique, les propriétés d'adsorption (isothermes et enthalpies) de la faujasite vis-à-vis de CO ont été modélisées, et le mécanisme microscopique d'adsorption de CO a été établi. La mise en œuvre de ces simulations a nécessité de paramétrer un nouveau champ de force destiné à décrire les interactions CO/faujasite et CO/CO. / CO adsorption in CuI and Na+ exchanged faujasite has been modeled by mean of quantum (DFT) and classical (Monte Carlo) approaches. By mean of the DFT calculations, faujasite potential energy surface has been explored. Different types of CO interactions with the cations have been highlighted, for each one of them CO adsorption effects on the structural and energetic parameters have been analyzed, and calculations of the CO stretching frequency have been performed. Thanks to our calculated values, a new attribution of CO adsorption spectra in CuY and NaY has been established. On another side, by mean of Monte Carlo simulations in the Grand Canonical ensemble, faujasite adsorption properties regarding CO (isotherms and enthalpies) have been modeled, and the CO adsorption mechanism has been established at the microscopic level. The implementation of these simulations has required the derivation of a new force field describing the CO/faujasite and CO/CO interactions.

The present Doctoral Thesis focuses on the investigation, characterization and influence of polyacrylic materials in different scientific and technological disciplines via a detailed computer simulation using the Molecular Dynamics (MD) technique, in conjunction with the very accurate, all-atom Dreiding force-field. The main research concepts and objectives are discussed and analyzed in three separate parts. In the first part, atomistic configurations of two model pressure-sensitive acrylic adhesives (PSAs), the atactic homopolymer poly(n-BA) [poly(n-butyl acrylate)] and the atactic copolymer poly(n-BA-co-AA) [poly(n-butyl acrylate-co-acrylic acid)] in the bulk phase or confined between two selected substrates, glassy silica (SiO2) and metallic α-ferrite (α-Fe), were built and simulated by MD in the NPT statistical ensemble. First, an equilibration cycle consisting of temperature annealings and coolings was followed, in order to generate well-equilibrated configurations of the PSA systems. Detailed results from the atomistic simulations are presented concerning their volumetric behavior, glass transition temperature, conformational, structural, viscoelastic and dynamic properties. Particular emphasis was given to the analysis and characterization of the hydrogen bonds that form in the poly(n-BA-co-AA) system. By analyzing the MD trajectories, poly(n-BA-co-AA) was found to exhibit a higher density than poly(n-BA) by about 7% at all temperatures, to be characterized by smaller-size chains for a given molecular weight (MW), to exhibit significantly slower terminal and segmental dynamics properties, and to be characterized by a glass transition temperature that was approximately 40% higher than that of poly(n-BA). We also examined the type and degree of adsorption of the two acrylic systems on the selected substrates by analyzing the MD results for the local mass density as a function of distance from the solid plane and the distribution of adsorbed chain segments in train, loop, and tail conformations, and by computing the work of adhesion at the two substrates. The results revealed a stronger adsorption for both acrylics on the SiO2 surface due to highly attractive interactions between polymer molecules and substrate atoms, and as a consequence a higher value for the work of adhesion compared to that on the α-Fe surface. Furthermore, we have developed a generalized non-equilibrium molecular dynamics (NEMD) algorithm to simulate the mechanical response of the two adhesives under a uniaxial stretching deformation. In the second part of the Thesis, results have been obtained from a hierarchical simulation methodology that led to the prediction of the thermodynamic, conformational, structural, dynamic and mechanical properties of two polymer nanocomposites based on syndiotactic poly(methyl methacrylate) or sPMMA. The first was reinforced with uniformly dispersed graphene sheets and the second with fullerene particles. How graphene functionalization affects the elastic constants of the resulting nanocomposite has also been examined. The phase behavior of the nanocomposite (in particular as we varied the relative size between the sPMMA chains and the diameter of fullerene molecules) has also been studied as a function of fullerene volume fraction. The simulation strategy entailed three steps: 1) Generation of an initial structure, which was then subjected to potential energy minimization and detailed molecular dynamics (MD) simulations at T = 500K and P = 1atm to obtain well relaxed melt configurations of the nanocomposite. 2) Gradual cooling of selected configurations down to room temperature to obtain a good number of structures representative of the glassy phase of the polymer nanocomposite. 3) Molecular mechanics (MM) calculations of its mechanical properties following the method originally proposed by Theodorou and Suter. By analyzing the results under constant temperature and pressure, all nanocomposite systems were found to exhibit slower terminal and segmental relaxation dynamics than the pure polymer matrices. The addition of a small fraction of graphene sheets led in all cases to the enhancement of the elastic constants; this was significantly more pronounced in the case of functionalized graphene sheets. We further mention that, for all polymer/fullerene nanocomposites addressed here, no phase separation or variation of polymer chain dimensions was observed as a function of fullerene size and/or fullerene volume fraction. In the third part of the Thesis, and motivated by the use of acrylic polymers for the design of membranes with aligned carbon nanotubes (CNTs) for several separation technologies (such as water desalination and wastewater treatment), we report results from a detailed computer simulation study for the nano-sorption and mobility of four different small molecules (water, tyrosol, vanillic acid, and p-coumaric acid) inside smooth single-wall CNTs (SWCNTs). Most of the results have been obtained with the molecular dynamics (MD) method, but especially for the most narrow of the CNTs considered, the results for water molecule were further confirmed through an additional Grand Canonical (μVT) Monte Carlo (GCMC) simulation using a value for the water chemical potential μ pre-computed with the particle deletion method. Issues addressed in the Thesis include molecular packing and ordering inside the nanotube for the four molecules, average number of sorbed molecules per unit length of the tube, and mean residence time and effective axial diffusivities, all as a function of tube diameter and tube length. In all cases, a strong dependence of the results on carbon nanotube diameter was observed, especially in the way the different molecules are packed and organized inside the CNT. For water for which predictions of properties such as local structure and packing were computed with both methods (MD and GCMC), the two sets of results were found to be fully self-consistent for all types of SWCNTs considered. Water diffusivity inside the CNT (although, strongly dependent on the CNT diameter) was computed with two different methods, both of which gave identical results. For large enough CNT diameters (larger than about 13 Å), this was found to be higher than the corresponding experimental value in the bulk by about 55%. Surprisingly enough, for the rest of the (phenolic) molecules simulated in this Thesis, the simulations revealed no signs of mobility inside nanotubes with a diameter smaller than the (20, 20) tube. This has been attributed to strong phenyl-phenyl attractive interactions, also to favorable interactions of these molecules with the CNT walls, which cause them to form highly ordered, very stable structures inside the nanotube, especially under strong confinement. The interaction, in particular, of the methyl group (present in tyrosol, vanillic acid, and p-coumaric acid) with the CNT walls seems to play a key role in all these compounds causing them to remain practically immobile inside nanotubes characterized by diameters smaller than about 26 Å. It was only for larger-diameter CNTs that tyrosol, vanillic acid, and p-coumaric acid were observed to demonstrate appreciable mobility. / Η παρούσα Διδακτορική Διατριβή εστιάζει στη μελέτη της σχέσης μεταξύ δομής και μακροσκοπικών φυσικών ιδιοτήτων υλικών από πολυακρυλικά μέσω μίας λεπτομερούς προσομοίωσης στον υπολογιστή με τη μέθοδο της Μοριακής Δυναμικής (ΜΔ), σε συνδυασμό με ένα πολύ επακριβές πεδίο δυνάμεων (το Dreiding) σε ατομιστική λεπτομέρεια. Οι κύριες ερευνητικές έννοιες καθώς και οι στόχοι συζητιούνται και αναλύονται σε τρία ξεχωριστά μέρη. Στο πρώτο μέρος, ατομιστικές απεικονίσεις δύο προτύπων πίεσο-ευαίσθητων συγκολλητικών υλικών (acrylic pressure sensitive adhesives ή PSAs), του ατακτικού πολυ-βουτυλικού-ακρυλικού εστέρα (poly(n-BA)) και του συμπολυμερούς του με ακρυλικό οξύ (poly(n-BA-co-AA)), τόσο μακριά όσο και κοντά σε υποστρώματα σίλικας (SiO2) και α-φερρίτη (α-Fe), μελετήθηκαν στη βάση ενός φάσματος ιδιοτήτων (θερμοδυναμικές, δομικές, ιξωδοελαστικές, δυναμικές, και συγκολλητικές), όπως και η μηχανική τους απόκριση υπό συνθήκες μονοαξονικής εκτατικής παραμόρφωσης. Στο δεύτερο μέρος παρουσιάζονται τα αποτελέσματα που εξήχθησαν από μία ιεραρχική μεθοδολογία προσομοίωσης που οδήγησε στην πρόβλεψη της φασικής συμπεριφοράς και των μηχανικών ιδιοτήτων νανοσύνθετων πολυμερικών υλικών (polymer nanocomposites ή PNCs) βασισμένων στο συνδιοτατκτικό πολυ-μεθακρυλικό μεθυλεστέρα (syndiotactic poly(methyl methacrylate) ή sPMMA), ενισχυμένο με ομοιόμορφα διεσπαρμένα φύλλα γραφενίου (graphene sheets) ή σωματίδια φουλερενίου (fullerene particles). Στο τρίτο μέρος, υποκινούμενοι από τη χρήση των ακρυλικών πολυμερών στο σχεδιασμό μεμβρανών με ενσωματωμένους ευθυγραμμισμένους νανοσωλήνες άνθρακα (ΝΑ, carbon nanotubes ή CNTs) σε διάφορες τεχνολογίες διαχωρισμού μορίων (με έμφαση στον καθαρισμό του νερού), παρουσιάζουμε αποτελέσματα από προσομοιώσεις, για τη νανο-ρόφηση και την κινητικότητα τεσσάρων διαφορετικών μικρών μορίων (water, tyrosol, vanilic acid, και p-coumaric acid) στο εσωτερικό λείων μονο-στρωματικών ΝΑ (single-wall CNTs ή SWCNTs). Τα θέματα που εξετάζονται περιλαμβάνουν τη μοριακή διευθέτηση και τη διάταξη στο εσωτερικό Ν.Α. των τεσσάρων μορίων, το μέσο χρόνο παραμονής τους, καθώς και τους αξονικούς συντελεστές διάχυσής του, συναρτήσει της διαμέτρου και του μήκους των ΝΑ.