Τα πρωτότυπα αποτελέσματα της διατριβής περιέχονται στα κεφ. 2, 3 και 4. Στο κεφ. 2 μελετούμε με αριθμητικές μεθόδους ένα 3-διάστατο διακριτό μοντέλο μακροοικονομίας με σταθερές ισοτιμίες. Χρησιμοποιώντας μια μέθοδο πλέγματος, βρίσκουμε την περιοχή ευστάθειας στον παραμετρικό χώρο, προσδιορίζουμε την καμπύλη διακλαδώσεων Hopf-Neimark και θεωρούμε σύντομα την εμφάνιση “γλωσσών” Arnold. Υπολογίζονται διαγράμματα διακλαδώσεων και εκθέτη Lyapunov που δίνουν πληροφορίες για τους επιχειρηματικούς κύκλους και την πολύπλοκη δυναμική του μοντέλου και. παρουσιάζουμε παραδείγματα κυκλικών και χαοτικών ελκυστών. Στο κεφ. 3 μελετούμε με τις ίδιες μεθόδους ένα διακριτό μοντέλο αλληλεπίδρασης περιοχών με σταθερές ισοτιμίες, επέκταση του προηγούμενου μοντέλου σε 5 διαστάσεις. Στόχος ήταν να δείξουμε πόσο εφικτή και αποτελεσματική είναι μία αριθμητική μελέτη για ηπίως πολυδιάστατα διακριτά δυναμικά συστήματα με πολλές παραμέτρους. Βρήκαμε ότι η κίνηση κεφαλαίων δεν αρκεί για τη δημιουργία κύκλων όταν είναι χαμηλή η εμπορική αλληλεπίδραση. Το κατώφλι εμπορίου προβλέπεται περίπου στο 15% των εμπορικών συναλλαγών. Αντίθετα, το μοντέλο δεν προβλέπει αναγκαίο ελάχιστο επίπεδο κίνησης κεφαλαίων για την εμφάνιση των κύκλων. Δίνουμε παραδείγματα διαγραμμάτων διακλάδωσης και εκθέτη Lyapunov που δείχνουν την εμφάνιση κύκλων ή ακολουθίας διπλασιασμού περιόδου, και παραδείγματα της ανάπτυξής τους. Το κεφ. 4 περιέχει σύντομη περιγραφή συμπληρωματικών αποτελεσμάτων στα παραπάνω μοντέλα, και στα αντίστοιχα μοντέλα μεταβλητής ισοτιμίας συναλλάγματος, καθώς και κατευθύνσεις μελλοντικής έρευνας. Στο κεφ. 5 περιγράφονται σύντομα οι υπολογιστικές τεχνικές που χρησιμοποιήσαμε. Η διατριβή δείχνει την αποτελεσματικότητα της αριθμητικής προσέγγισης για πολυδιάστατα διακριτά μοντέλα. / The original results of the dissertation are contained in ch. 2, 3 and 4, and concern mainly the problem of business cycles. In ch. 2 we explore numerically a 3D discrete Kaldorian macrodynamic model of open economy with fixed exchange rates. Using a grid search method we determine the stability region in parameter space, and the Hopf-Neimark bifurcation curve, and discuss briefly the occurrence of Arnold tongues. Bifurcation and Lyapunov exponent diagrams are computed providing information on the business cycles and illustrating the complex dynamics involved. Examples of cycles and chaotic attractors are presented. In ch. 3 we explore a 5D extension of the previous model using the same methods. The aim was to demonstrate the feasibility and effectiveness of the numerical approach for discrete dynamical systems of moderately high dimensionality and several parameters. We found that capital movement is not sufficient to generate interregional business cycles when trade interaction is low. The trade threshold is predicted at about 15% of trade transactions. By contrast, no minimum level of capital mobility exists as a requirement for the emergence of business cycles. Examples of bifurcation and Lyapunov exponent diagrams illustrating the occurrence of cycles or period doubling, and examples of their development, are given. Ch. 4 contains a short description of complementary results on the above models, and on two other models which extend the previous models to the case of flexible exchange rates, as well as some lines of future research. In ch. 5, the computational techniques employed in the present study are briefly described. The dissertation indicates the effectiveness of the numerical approach for high dimensional discrete models.

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