Το αντικείμενο της έρευνας της παρούσας εργασίας είναι η ατομιστική προσομοίωση της ρόφησης αερίων μέσα σε ζεολιθους, καθώς και η επίδραση της θερμοκρασίας στην δομή του και την ροφητική του ικανότητα. Από τους σχεδόν 200 συνθετικούς και φυσικούς ζεολίθους, επιλέχθηκε ο φωγιασίτης, ο οποίος αποτελεί ένα πολλά υποσχόμενο υλικό με πλήθος εφαρμογών, όπως η ρόφηση, η αποθήκευση και ο διαχωρισμός αερίων, αλλά και σε τομείς όπως η κατάλυση και η δημιουργία συσκευών ανίχνευσης αερίων. Στα πλαίσια της παρούσας εργασίας εξετάστηκε η αλληλεπίδραση των κατιόντων της δομής με μόρια διαφόρων αερίων, τα οποία ροφούνται στον φωγιασίτη, καθώς και ο βαθμός στον οποίο αυτή η αλληλεπίδραση επηρεάζει την ποσότητα του αερίου που μπορεί να ροφηθεί. Πιο συγκεκριμένα, αρχικά ανακατασκευάστηκε στον υπολογιστή η δομή της μοναδιαίας κυψελίδας του φωγιασίτη και έγινε μοντελοποίηση της ρόφησης μορίων CO2 και Η2 μέσα στην ανακατασκευασμένη κυψελίδα ως συνάρτηση της θερμοκρασίας, του λόγου Si/Al στον κρύσταλλο, και της σύστασής του σε μη πλεγματικά κατιόντα. Στο τμήμα αυτό έγινε η θεώρηση ότι η δομή μένει αμετάβλητη με την είσοδο των μορίων του αερίου. Στην συνέχεια μελετήθηκε η επίδραση της θερμοκρασίας στην θέση των ιόντων, των ατόμων του πλέγματος και στο μέγεθος της μοναδιαίας κυψελίδας, καθώς και η περαιτέρω επίδραση που φέρουν αυτές οι αλλαγές στη ρόφηση. Τα αποτελέσματα της ρόφησης συγκρίθηκαν με πειραματικά δεδομένα. Τέλος, μοντελοποιήθηκε η διάχυση του CO2 μέσα στον φωγιασίτη και υπολογίστηκε ο συντελεστής διαχύσεως σε διάφορες θερμοκρασίες. Τα αποτελέσματα της μοντελοποίησης συγκρίθηκαν με αντίστοιχα αποτελέσματα από πειραματική μέτρηση του συντελεστή διαχύσεως του αερίου στον φωγιασίτη. Η αναπαράσταση του υλικού και των φαινομένων της ρόφησης έγινε σε ατομιστικό επίπεδο, με την χρήση μεθόδων Monte Carlo, ενώ η μοντελοποίηση της συμπεριφοράς του φωγιασίτη με τις μεταβολές της θερμοκρασίας πραγματοποιήθηκε χρησιμοποιώντας τη μέθοδο της προσομοιωμένης ανόπτησης. Η διάχυση μοντελοποιήθηκε με την μέθοδο της μοριακής δυναμικής. / -

Orientador: Maurice de Koning / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Fisica Gleb Wataghin / Made available in DSpace on 2018-08-08T20:27:57Z (GMT). No. of bitstreams: 1 PomaBernaola_AdolfoMaximo_M.pdf: 3697892 bytes, checksum: a07c1ad647a61d9862283f697732410e (MD5) Previous issue date: 2007 / Resumo: Nesta dissertação abordamos o estudo de uma das limitações da simulação atomística denominada o evento raro, quem é responsável pela limitação temporal, exemplos de problemas que envolvem os eventos raros são, o enovelamento de proteínas, mudanças conformacionais de moléculas, reações químicas (em solução), difusão de sólidos e os processos de nucleação numa transição de fase de 1a ordem, entre outros. Métodos convencionais como Dinâmica Molecular (MD) ou Monte Carlo (MC) são úteis para explorar a paisagem de energia potencial de sistemas muito complexos, mas em presença de eventos raros se tornam muito ineficientes, devido à falta de estatística na amostragem do evento. Estes métodos gastam muito tempo computacional amostrando as configurações irrelevantes e não as transições de interesse. Neste sentido o método Transition Path Sampling (TPS), desenvolvido por D. Chandler e seus colaboradores, consegue explorar a paisagem de energia potencial e obter um conjunto de verdadeiras trajetórias dinâmicas que conectam os estados metaestáveis em presença de evento raros. A partir do ensemble de caminhos a constante de reação e o mecanismo de reação podem ser extraídos com muito sucesso. Neste trabalho de mestrado implementamos com muito sucesso o método TPS e realizamos uma comparação quantitativa em relação ao método MC configuracional num problema padrão da isomerização de uma molécula diatômica imersa num líquido repulsivo tipo Weeks-Chandler-Andersen (WCA). A aplicação destes métodos mostrou como o ambiente, na forma de solvente, pode afetar a cinética de um evento raro / Abstract: In this dissertation we aproach the study of one of the limitations of the atomistic simulation called the rare event, which is responsible for the temporal limitation. Examples of problems that involve the rare event are the folding protein, conformational changes in molecules, chemical reactions (in solution), solid diffusion, and the processes of nucleation in a first-order phase transition, among other. Conventional methods as Molecular Dynamics (MD) or Monte Carlo (MC) are useful to explore the potencial energy landscape of very complex systems, but in presence of rare events they become very inefficient, due to lack of statistics in the sampling of the event. These methods spend much computational time sampling the irrelevant configurations and not the transition of interest. In this sense, the Transition Path Sampling (TPS) method, developed by D. Chandler and his collaborators, can explore the potential energy landscape and get a set of true dynamical trajectories that connect the metastable states in presence of the rare events. From this ensemble of trajectories the rate constant and the mechanism of reaction can be extracted with great success. In this work we implemented the TPS method and carried out a quantitative comparison in relation to the configurational MC method in a standard problem of the isomerization of a diatomic molecule immersed in a Weeks-Chandler-Andersen (WCA) repulsive fluid. The application of these methods showed as the environment, in the form of solvent, can affect the kinetic of a rare event / Mestrado / Física Estatistica e Termodinamica / Mestre em Física

Orientador: Charlles Rubber de Almeida Abreu / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-18T21:35:53Z (GMT). No. of bitstreams: 1 Martins_TiagoDias_M.pdf: 4166409 bytes, checksum: 78baf81dcc9cf9d746302c21e813f9e2 (MD5) Previous issue date: 2011 / Resumo: O estudo a nível microscópico do Equilíbrio Sólido-Líquido (ESL) só se desenvolveu a partir do século XX e ainda existem aspectos de tal fenômeno que não são bem compreendidos, mesmo para os sistemas mais simples. Por esse motivo, é necessário que métodos eficientes para simular diretamente o ESL sejam formulados. O principal objetivo deste trabalho foi desenvolver e implementar Métodos Monte Carlo Multicanônicos para simular o ESL de sistemas de Lennard-Jones e, então, avaliar a aplicabilidade de tais métodos com base na sua eficiência. Primeiramente, buscou-se melhorar os resultados do método desenvolvido por Muguruma e Okamoto (2008), que visava amostrar a energia do sistema uniformemente e o volume com probabilidade de Boltzmann, enquanto mantinha-se fixo o número de partículas da caixa de simulação cúbica. Em seguida, uma nova abordagem foi proposta, cuja principal característica é não restringir o formato da caixa, embora a amostragem fosse realizada da mesma forma. Em ambos os métodos, a densidade de estados do sistema foi estimada usando a técnica da matriz de transições e uma abordagem de paralelização das simulações foi empregada a fim de se obter uma amostragem mais efetiva. Os resultados obtidos com o primeiro método mostraram que uma simulação como essa depende significativamente do número de partículas. Além disso, a amostragem dos estados de energia ficou presa na fase sólida. Já os resultados obtidos com a nova metodologia foram ligeiramente melhores que os obtidos com o método anterior. A amostragem para o sistema com o menor número de partículas foi equivalente para as duas fases em questão. No entanto, as simulações com sistemas maiores apresentaram o mesmo problema apontado no método anterior. A princípio, sugeriu-se a existência de uma transição polimórfica. No entanto, simulações na fase sólida com a nova metodologia apontam o contrário. Apesar dos resultados insatisfatórios, este trabalho se constitui em um passo importante no estudo de um tema complexo, que é a previsão de Equilíbrio Sólido-Líquido utilizando métodos de simulação direta / Abstract: The study at the microscopic level of Solid-Liquid Equilibrium (SLE) was only developed in the 20th century and still exist some aspects of this phenomenon that are not well understood, even for the simplest systems. Therefore, it is necessary the development of efficient methods to directly simulate the SLE. The main objective of this work was to develop and implement Multicanonical Monte Carlo methods to simulate the SLE of the Lennard-Jones system and then evaluate their applicability based on efficiency. First, we sought to improve the performance of a method used by Muguruma and Okamoto (2008), which aimed to sample the system energy uniformly and the volume with the Boltzmann probability, while keeping fixed the number of particles in a cubic simulation box. Then a new approach was proposed, whose main feature was let the shape of the box free, although the sampling was performed in the same way. In both methods, the density of states of the system was estimated using the transition matrix technique and an approach for parallelization of the simulations was employed to obtain a more effective sampling. The results obtained with the first method showed that such a simulation depends on the number of particles present in the system. Moreover, the sampling of energy states got stuck in the solid phase. The results obtained with the new method were slightly better than those obtained with the previous method. The sampling for the system with the lowest number of particles was equivalent for the two phases concerned. However, simulations with larger systems presented the same problem of the previous method. At first, it was suggested the existence of a polymorphic transition. However, simulations in the solid phase with the new methodology show the opposite. Despite the unsatisfactory results, this work constitutes na important step in studying a complex subject, which is the prediction of Solid-Liquid Equilibrium using direct simulation methods / Mestrado / Engenharia de Processos / Mestre em Engenharia Química

Les excellentes propriétés électroniques, thermiques et mécaniques du graphène confèrent à ce matériau planaire (bi-dimensionnel) un énorme potentiel applicatif, notamment en électronique. Néanmoins, ce matériau présente de sérieux inconvénients qui pourraient limiter son champ d'applications. Par exemple, sa structure de bandes électronique sans bande interdite rend difficile le blocage du courant dans un dispositif. De plus, pour les applications thermoélectriques, sa forte conductance thermique est aussi une forte limitation. Il y a donc beaucoup de défis à relever pour rendre ce matériau vraiment utile pour des applications. Cette thèse porte sur l'étude des propriétés électroniques et thermoélectriques dans les hétérostructures planaires constituées de graphène et de nitrure de bore hexagonal (BN). Différentes configuration de ce nouveau matériau hybride permettent de moduler la bande interdite, la conductance thermique et le coefficient Seebeck. Cette étude a été menée au moyen de calculs atomistiques basés sur les approches des liaisons fortes (TB) et du modèle à constantes de force (FC). Le transport d'électrons et de phonons a été simulé dans le formalisme des fonctions de Green hors équilibre. Les résultats montrent que, grâce à la modulation de la bande interdite, des transistors à base d'hétérostructures de BN et de graphène peuvent présenter un très bon rapport courant passant / bloqué d'environ 10⁴ à 10⁵. En outre, nous montrons l'existence d'états quantiques hybrides à l'interface zigzag entre le graphène et le BN donnant lieu à des propriétés de transport électronique très intéressantes. Enfin, ce travail montre qu'en agençant correctement des nano-flocons de BN sur les côtés d'un nanoruban de graphène, la conductance des phonons peut être fortement réduite alors que l'ouverture de bande interdite conduit à un accroissement important du coefficient Seebeck. Il en résulte qu'un facteur de mérite thermoélectrique ZT plus grand que l'unité peut être réalisé à température ambiante. / Graphene is a fascinating 2-dimensional material exhibiting outstanding electronic, thermal and mechanical properties. Is this expected to have a huge potential for a wide range of applications, in particular in electronics. However, this material also suffers from a strong drawback for most electronic devices due to the gapless character of its band structure, which makes it difficult to switch off the current. For thermoelectric applications, the high thermal conductance of this material is also a strong limitation. Hence, many challenges have to be taken up to make it useful for actual applications. This thesis work focuses on the theoretical investigation of a new strategy to modulate and control the properties of graphene that consists in assembling in-plane heterostructures of graphene and Boron Nitride (BN). It allows us to tune on a wide range the bandgap, the thermal conductance and the Seebeck coefficient of the resulting hybrid nanomaterial. The work is performed using atomistic simulations based on tight binding (TB), force constant (FC) models for electrons and phonons, respectively, coupled with the Green's function formalism for transport calculation. The results show that thanks to the tunable bandgap, it is possible to design graphene/BN based transistors exhibiting high on/off current ratio in the range 10⁴-10⁵. We also predict the existence hybrid quantum states at the zigzag interface between graphene and BN with appealing electron transport. Finally this work shows that by designing properly a graphene ribbon decorated with BN nanoflakes, the phonon conductance is strongly reduced while the bandgap opening leads to significant enhancement of Seebeck coefficient. It results in a thermoelectric figure of merit ZT larger than one at room temperature.

Des techniques combinées d'analyse/synthèse d'images et de simulation atomistique ont permis d’étudier la nanostructure/-texture de matériaux carbonés anisotropes et denses de type pyrocarbone (PyC) laminaire hautement texturé. Des représentations atomiques d’un PyC laminaire rugueux tel que préparé (AP) ainsi que d’un PyC laminaire régénéré AP et après plusieurs traitements thermiques (HT) ont été reconstruites pour mieux caractériser ces matériaux. Ces modèles comportent des domaines graphéniques de quelques nanomètres, joints entre eux par des lignes de défauts formées de paires de cycles à 5 et 7 carbones dans le plan et par des dislocations vis et des atomes tétravalents entre les plans. Les modèles les plus ordonnés ont des domaines plus étendus et un plus faible taux de connexions inter-plan. Les propriétés mécaniques et thermiques prédites à partir de ces modèles sont proches de celles du graphite et augmentent avec la cohérence intra-plan et la densité de connexions inter-plans. Des modèles de graphène polycristallins ont aussi été générés. Ils sont apparus, du point de vue structural et des propriétés mécaniques, très proches des feuillets de carbones des PyCs. Ils ont permis d'étudier la réorganisation structurale se produisant au cours du HT : formation de lignes de défauts, réparation de lacunes, … Il s'agit d'un premier pas vers l'étude de la graphitation des PyCs. La méthode de reconstruction a enfin été adaptée à l'étude de l'évolution structurale d'un graphite au cours de son irradiation par les électrons. Cela a permis d'observer à l'échelle atomique la création et la propagation des défauts au cours de l'irradiation. / Combined images analysis/synthesis techniques and atomistic simulation methods have allowed studying the nanostructure/-texture of anisotropic dense carbons of the highly textured laminar pyrocarbon (PyC) type.Atomic representations of an as-prepared (AP) rough laminar PyC as well as a regenerative laminar PyC AP and after several heat treatments (HT) were reconstructed to better characterize these materials.The models contain nanosized graphene domains connected between them by line defects formed by pairs of rings with 5 and 7 carbons inside layers and by screw dislocations and fourfold atoms between layers. The most ordered models have larger domains and a lower percentage of connections between the layers.Mechanical and thermal properties predicted from these models are close to those of graphite and increase with the coherence inside layers and the density of connections between layers.Models of polycrystalline graphene were also generated, showing structure and mechanical properties very close to those of the carbon layers extracted from PyCs. The structural reorganization occurring during the HT of such materials was studied: thinning of line defects and vacancy healing were observed. This represents a first step towards the study of the graphitization of PyCs.The reconstruction method was eventually adapted to study the structural evolution of a nuclear-grade graphite during its irradiation by electrons, allowing us to observe how defects are created and propagate during irradiation.

The presence of the ordered $\gamma^{\prime}$ phase ($\text{Ni}_{3}\text{Al}$) in Ni-base superalloys is fundamental to the performance of engineering components such as turbine disks and blades which operate at high temperatures and loads. Hence for these alloys it is important to optimize their microstructure and phase composition. This is typically done by varying their chemistry and heat treatment to achieve an appropriate balance between $\gamma^{\prime}$ content and other constituents such as carbides, borides, oxides and topologically close packed phases. In this work we have set out to investigate the onset of $\gamma^{\prime}$ ordering in Ni-Al single crystals and in Ni-Al bicrystals containing coincidence site lattice grain boundaries (GBs) and we do this at high temperatures, which are representative of typical heat treatment schedules including quenching and annealing. For this we use the atomistic simulation methods of molecular dynamics (MD) and density functional theory (DFT). In the first part of this work we develop robust Bayesian classifiers to identify the $\gamma^{\prime}$ phase in large scale simulation boxes at high temperatures around 1500 K. We observe significant \gamma^{\prime} ordering in the simulations in the form of clusters of $\gamma^{\prime}$-like ordered atoms embedded in a $\gamma$ host solid solution and this happens within 100 ns. Single crystals are found to exhibit the expected homogeneous ordering with slight indications of chemical composition change and a positive correlation between the Al concentration and the concentration of $\gamma^{\prime}$ phase. In general, the ordering is found to take place faster in systems with GBs and preferentially adjacent to the GBs. The sole exception to this is the $\Sigma3 \left(111\right)$ tilt GB, which is a coherent twin. An analysis of the ensemble and time lag average displacements of the GBs reveals mostly `anomalous diffusion' behaviour. Increasing the Al content from pure Ni to Ni 20 at.% Al was found to either consistently increase or decrease the mobility of the GB as seen from the changing slope of the time lag displacement average. The movement of the GB can then be characterized as either `super' or `sub-diffusive' and is interpreted in terms of diffusion induced grain boundary migration, which is posited as a possible precursor to the appearance of serrated edge grain boundaries. In the second part of this work we develop a method for the training of empirical interatomic potentials to capture more elements in the alloy system. We focus on the embedded atom method (EAM) and use the Ni-Al system as a test case. Recently, empirical potentials have been developed based on results from DFT which utilize energies and forces, but neglect the electron densities, which are also available. Noting the importance of electron densities, we propose a route to include them into the training of EAM-type potentials via Bayesian linear regression. Electron density models obtained for structures with a range of bonding types are shown to accurately reproduce the electron densities from DFT. Also, the resulting empirical potentials accurately reproduce DFT energies and forces of all the phases considered within the Ni-Al system. Properties not included in the training process, such as stacking fault energies, are sometimes not reproduced with the desired accuracy and the reasons for this are discussed. General regression issues, known to the machine learning community, are identified as the main difficulty facing further development of empirical potentials using this approach.

<p align='justify'>L’adaptation des surfaces pour des fonctions prédéterminées par le choix des matériaux métalliques ou des couches minces ayant des propriétés mécaniques avancées peut potentiellement permettre de réaliser des nouvelles applications à petites échelles. Concevoir de telles applications utilisant des nouveaux matériaux nécessite en premier lieu la connaissance des propriétés mécaniques des matériaux ciblés à l’échelle microscopique et nanoscopique. Une méthode souvent appliquée pour caractériser les matériaux à petites échelles est la nanoindentation, qui peut être vue comme une mesure de dureté à l’échelle nanoscopique.</p><p><p align='justify'>Ce travail présente une contribution relative à l'interprétation des résultats de la nanoindentation, qui fait intervenir un grand nombre de phénomènes physiques couplés à l'aide de simulations numériques. A cette fin une approche interdisciplinaire, adaptée aux phénomènes apparaissant à petites échelles, et située à l’intersection entre la physique, la mécanique et la science des matériaux a été utilisée. Des modèles numériques de la nanoindentation ont été conçus à l'échelle atomique (modèle discret) et à l'échelle des milieux continus (méthode des éléments finis), pour étudier le comportement du nickel pur. Ce matériau a été choisi pour ses propriétés mécaniques avancées, sa résistance à l'usure et sa bio-compatibilité, qui peuvent permettre des applications futures intéressantes à l'échelle nanoscopique, particulièrement dans le domaine biomédical. Des méthodes avancées de mécanique du solide ont été utilisées pour prendre en compte les grandes déformations locales du matériau (par la formulation corotationelle), et pour décrire les conditions de contact qui évoluent au cours de l'analyse dans le modèle à l'échelle des milieux continus (traitement des conditions de contact unilatérales et tangentielles par une forme de Lagrangien augmenté).</p><p><p align='justify'>L’application des modèles numériques a permis de contribuer à l’identification des phénomènes qui gouvernent la nanoindentation du nickel pur. Le comportement viscoplastique du nickel pur pendant nanoindentation a été identifié dans une étude expérimentale-numérique couplée, et l'effet cumulatif de la rugosité et du frottement sur la dispersion des résultats de la nanoindentation a été montré par une étude numérique (dont les résultats sont en accord avec des tendances expérimentales).</p> <p><p align='justify'>Par ailleurs, l’utilisation de l’outil numérique pour une autre application à petites échelles, la manipulation des objets par contact, a contribué à la compréhension de la variation de l’adhésion électrostatique pendant micromanipulation. La déformation plastique des aspérités de surface sur le bras de manipulateur (en nickel pur) a été identifiée comme une source potentielle d’augmentation importante de l'adhésion pendant la micromanipulation, qui peut potentiellement causer des problèmes de relâche et de précision de positionnement, observés expérimentalement.</p><p><p align='justify'>Les résultats présentés dans cette thèse montrent que des simulations numériques basées sur la physique du problème traité peuvent expliquer des tendances expérimentales et contribuer à la compréhension et l'interprétation d'essais couramment utilisé pour la caractérisation aux petites échelles. Le travail réalisé dans cette thèse s’inscrit dans un projet de recherche appelé "mini-micro-nano" (mµn), financé par la Communauté Française de Belgique dans le cadre de "l'Action de Recherche Concertée", convention 04/09-310.</p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

No description available.

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). Τα θέματα που εξετάζονται περιλαμβάνουν τη μοριακή διευθέτηση και τη διάταξη στο εσωτερικό Ν.Α. των τεσσάρων μορίων, το μέσο χρόνο παραμονής τους, καθώς και τους αξονικούς συντελεστές διάχυσής του, συναρτήσει της διαμέτρου και του μήκους των ΝΑ.