We numerically investigate the long-time behavior of the density-density auto-correlation function in driven lattice gases, with particle exclusion and periodic boundary conditions in one, two, and three dimensions using precise Monte Carlo simulations of larger system sizes than previous studies. In the one-dimensional asymmetric exclusion process on a ring with half the lattice sites occupied, we find that correlations induce extremely slow relaxation to the asymptotic power law decay We compare the crossover functions obtained from our simulations with various analytic results in the literature, and analyze the characteristic oscillations that occur in finite systems away from half-filling. As expected, correlations are weak in three dimensions and consequently the mean-field description is adequate. We also investigate the relaxation towards the non-equilibrium steady state in the two-time density-density auto-correlations, starting from strongly correlated initial conditions. We obtain simple aging scaling behavior in one, two, and three dimensions, with the expected power laws. We numerically investigate the behavior of driven lattice gases with nearest neighbor interactions at half-filling with periodic boundary conditions below and at the critical temperature using Monte Carlo simulations of very large lattices in two dimensions. This work is one of few that explores the relaxation to a non-equilibrium steady state. We obtain data collapse for the finite-size scaling form of density-density auto-correlation function at the critical point. We achieve data collapse using finite-size scaling of the time-dependent order parameter during the transient regime starting from strongly correlated initial conditions. We present simple aging scaling of the density-density auto-correlation function at the critical point starting from strongly correlated initial conditions using Monte Carlo simulations of two different lattice anisotropies. We thus unambiguously confirm the critical exponents determined by renormalization group methods using measurement of dynamic quantities in the transient regime. Measuring these dynamic quantities in the transient regime provides more conclusive measurements of the critical exponents than previous studies measuring static quantities in the stationary state. We provide qualitative arguments that the lattice anisotropy determines the steady-state for sub-critical quenches. / Ph. D.

Soumise à une sollicitation sismique importante, une structure de génie civil peut être endommagée. Cela se traduit par l'apparition de fissures et donc une réduction de ses propriétés élastiques de la structure. Cependant, dans le cas d'un événement sismique d'amplitude inférieure, la variation observée est seulement transitoire. Elle consiste en une perturbation brutale de la fréquence et de l'amortissement suivie d'un recouvrement lent vers leurs valeurs initiales. Ce phénomène non linéaire de relaxation est appelé Dynamique Lente. Il s'explique par la fermeture progressive des fissures initialement présentes dans la structure et qui ont été ouvertes par la sollicitation. Nous avons observé en laboratoire que l'analyse de la Dynamique Lente dans une poutre avant et après son endommagement permet de détecter cette augmentation de la densité de fissures. La sensibilité différentielle des modes à une dégradation locale indique, de plus, une piste pour le développement d'une méthode de localisation de l'endommagement. L'étude de la Dynamique Lente que nous avons menée dans deux structures de génie civil a montré que nous pouvions également y détecter l'apparition de l'endommagement. Le suivi de l'amortissement de manière continue a mis en évidence une relation linéaire entre l'amortissement et l'intensité des vibrations ambiantes à la fois dans les poutres et les structures de génie civil. Nous expliquons celle-ci via l'application du théorème Fluctuation-Dissipation à ces systèmes. Les résultats en laboratoire et la proximité de l'expression de la relation avec les équations utilisées dans le cas de la Dynamique Lente suggère une dépendance de ce phénomène vis à vis de la densité des hétérogénéités dans le matériau composant la structure. De plus amples recherches sont cependant nécessaires afin d'expliquer complètement nos observations et ainsi pouvoir les utiliser dans le cadre de la surveillance des structures uniquement à partir de données de très faible amplitude. / During strong seismic loadings, a structure may be damaged. This results in the appearing of cracks and then a reduction of the elastic properties of the structure. The degradation remains only transitory in the case of smaller seismic events. It consists in a sharp disruption of both the frequency and damping followed by their slow recovery to their initial values. This non linear phenomenon is called Slow Dynamics. It is explained by the gradual closing of the cracks which were initially present in the material and which were opened during the loading. We observed in the laboratory that the analysis of the Slow Dynamics in a beam before and after it is damaged allows to detect the increase of the crack density. The different sensitivties of the modes regarding a local damaging indicates a track to develop a method to locate the damages. The study of the Slow Dynamics in civil engineering structures demonstrated the possibility to detect the damaging also in this kind of system. The continuous monitoring of the damping highlighted a linear relationship between damping and the intensity of the ambiante vibration in the case of both the beams and real case structures. We explain this relation by applying the Fluctuation-Dissipation to these systems. Laboratory results and the proximity of the expression of the linear relationship wit hthe equations used in the theory of the Slow Dynamics suggest a dependency of this phenonmenon on the density of heterogeneities in the structure. Further research is however required in order to fully explain our observations and thus, to use them to monitor the state of health of structures.

Sensitive nonlinear acoustic methods are suitable for material characterization. This thesis describes three nonlinear acoustic methods that are proven useful for detection of defects like cracks and delaminations in solids. They offer the possibility to use relatively low frequencies which is advantageous because attenuation and diffraction effects are smaller for low frequencies. Therefore large and multi-layered complete objects can be investigated in about one second. Sometimes the position of the damage is required. But it is in general difficult to limit the geometrical extent of low-frequency acoustic waves. A technique is presented that constrains the wave field to a localized trapped mode so that damage can be located.

This dissertation seeks to merge the ideas from robust control theory such as H-Infinity control design and the Small Gain Theorem, L stability theory and Lyapunov stability from nonlinear control, and recent theoretical achievements in adaptive control. The fusion of frequency domain and linear time domain ideas allows the derivation of an H-Infinity Norm Minimization Approach (H-Infinity-NMA) for adaptive control architecture that permits a control designer to simplify the adaptive tuning process and tune the uncertainty compensation characteristics via linear control design techniques, band limit the adaptive control signal, efficiently handle redundant actuators, and handle unmatched uncertainty and matched uncertainty in a single design framework. The two stage design framework is similar to that used in robust control, but without sacrificing performance. The first stage of the design considers an ideal system with the system uncertainty completely known. For this system, a control law is designed using linear H-Infinity theory. Then in the second stage, an adaptive process is implemented that emulates the behavior of the ideal system. If the linear H-Infinity design is applied to control the emulated system, it then guarantees closed loop system stability of the actual system. All of this is accomplished while providing notions of transient performance bounds between the ideal system and the true system. Extensions to the theory include architectures for a class of output feedback systems, limiting the authority of an adaptive control system, and a method for improving the performance of an adaptive system with slow dynamics without any modification terms. Applications focus on using aerodynamic flow control for aircraft flight control and the Crew Launch Vehicle.

Ce travail est consacré à l'étude, par des techniques originales de diffusion dynamique de la lumière, de la dynamique lente de deux systèmes colloïdaux modèles. On a utilisée une technique de Corrélation Résolue dans le Temps (TRC) et des simulations numériques pour sonder, sur une grande plage de concentrations et sur sept décades en temps de relaxation, la dynamique d'un système de sphères dures colloïdales approchant de la transition vitreuse. A faible fraction volumique les temps de relaxation peuvent être représentés par une loi de puissance critique comme prévu par la Théorie de Couplage des Modes (MCT). A des concentrations élevées le temps de relaxation augmente selon une inattendue loi exponentielle critique avec une fraction volumique critique supérieure à celle prévue par MCT. L'étude des fluctuations de la dynamique confirme la présence d'un régime de dynamique activée: ce scénario est le même que l'on retrouve pour les systèmes vitreux moléculaires. Le deuxième sujet étudié est le comportement d'une suspension de colloïdes attractifs soumis à la contrainte gravitationnelle. Initialement, des agrégats de particules sédimentent en formant un gel sur le fond de la cellule. Les techniques spéciales employées nous ont permis de mesurer l'évolution des profils de concentration, des profils de vitesse et la dynamique locale dans la phase gel: la sédimentation peut être complètement décrite par le taux de déformation locale qui permet aussi une remise à l'échelle de la dynamique locale. Le rôle de la friction aux parois des cellules a été sondé par des expériences de polarimétrie et par la résolution numérique d'un modèle poroélastique pour la sédimentation du gel. / We investigated, using original light scattering techniques, the slow dynamics of two different model colloidal systems. We used a Time Resolved Correlation (TRC) technique and numerical simulations to study the dynamics of systems of colloidal hard spheres over a broad range of densities and seven decades in relaxation time. At low volume fraction the equilibrium relaxation time is described by the algebraic divergence predicted by mode-coupling theory (MCT). At higher density the relaxation time increases with an unexpected critical exponential form, whose the critical volume fraction is much larger then the MCT singularity. The fluctuations of the dynamics confirm the presence of an activated regime, as founded in molecular glass formers close to the glass transition. The second research subject concerns the behaviour of suspensions of attractive colloidal particles under gravitational stress. Initially, clusters fall forming a gelled deposit at the cell bottom. Our apparatus allows us to follow over time the concentration profile, the velocity profile and the local dynamics in the gel phase: the settling may be fully described by the local strain rate and dynamics exhibits remarkable scaling properties when time is normalized by strain rate. The role of the solid friction at the cell walls has been investigated by polarimetry experie nces and solving numerically a poroelastic model of gel sedimentation.

La surveillance de la réponse structurale est fondamentale pour estimer la performance des bâtiments et réduire les pertes lors de futurs séismes. Un moyen pratique de détecter les changements de comportement structural consiste à analyser les variations des propriétés élastiques lors d'excitations dynamiques. Dans ce travail, on montre que les variations de la fréquence fondamentale des bâtiments lors de tremblements de terre (faibles à forts) pourraient être expliquées par des processus élastiques non linéaires qui se produisent à l'intérieur du matériau, et qui finalement affectent le comportement macroscopique global des bâtiments. Ces processus élastiques non linéaires sont responsables de la diminution temporaire ou permanente de la rigidité structurale, pouvant expliquer les processus de récupération des propriétés élastiques observés à la suite d'événements sismiques. Cette étude comble le fossé entre des expériences de laboratoire à l'échelle microscopique et des observations sismologiques à l'échelle macroscopique, où l’élasticité non linéaire est également observée. Dans un premier temps, une base de données sismiques établie dans le cadre de cette thèse est présentée, incluant des réponses de bâtiments instrumentés de façon permanente dans le monde: des milliers d’enregistrements de mouvements sismiques et plusieurs bâtiments du Japon et des États-Unis ont été traités, apportant des connaissances utiles pour le domaine du génie parasismique, notamment pour la prédiction empirique de la réponse structurale en fonction de mesures d'intensité du mouvement au sol. Les incertitudes associées à la prédiction d’endommagement sont présentées, ainsi que l'évaluation de la vulnérabilité d'un bâtiment sous forme de courbes de fragilité. Ensuite, la base de données est utilisée pour analyser les signatures élastiques non linéaires dans les bâtiments, en particulier les effets de la dynamique lente (ou relaxation). Les variations des fréquences de résonance sont étudiées à court et à long terme, en estimant la contribution du sol à la réponse du système sol-structure. Différents états structuraux sont déduits en fonction des amplitudes de chargement et propriétés observées via les enregistrements. Des modèles de relaxation développés en laboratoire sont ensuite adaptés aux données des bâtiments afin de caractériser la densité de fissuration et les hétérogénéités, en effectuant des comparaisons entre les états structuraux avant et après de fortes excitations telles que le séisme de 2011 (Mw=9) de Tohoku (Japon). Les effets des chargements sont observés lors de la récupération des séquences de répliques. Les résultats sont étendus à différentes typologies de bâtiments, en analysant l'influence du matériau et des caractéristiques de chargement, notamment les taux de déformation. Enfin, quelques conclusions générales sont présentées, ainsi qu'une perspective de travail utilisant des outils de machine learning pour prédire la réponse de bâtiments en fonction de signatures élastiques non linéaires observées. / Monitoring structural response is fundamental for evaluating the performance of buildings and reducing losses during future earthquakes. One practical way to detect changes in structural behavior is analyzing variations of elastic properties during dynamic excitations. Here we show that variations in the fundamental frequency of buildings during (weak -to- strong) earthquakes might be explained by nonlinear elastic processes carried out within the structural material, which affect the global macroscopic structural behavior. These nonlinear elastic processes are responsible for both transitory and permanent structural softening, and might explain the intriguing recovery effects observed in the fundamental frequency of buildings following seismic events. This study bridges the gap between microscale laboratory experiments and macroscale seismological observations, where nonlinear elasticity is also observed. In the first part of this study, a new seismic database of building responses is presented: thousands strong motion recordings and several buildings from Japan and US were processed, providing useful tools for the earthquake engineering community, notably for the empirical prediction of structural response as a function of several ground motion intensity measures. Examples of uncertainties associated to damage prediction are presented, as well as the vulnerability assessment of a building throughout fragility curves. Next, the seismic database is used to analyze nonlinear elastic signatures in buildings, particularly the slow dynamics or relaxation effects. Variations of resonant frequencies are monitored at both short and long-term, estimating the contribution of soil in the response of the system soil-structure. Different levels of damage are inferred according to loading amplitudes and structural states. Some laboratory-based models of relaxation are adapted to the building data in order to infer crack-density and heterogeneities over time, making comparisons between structural states before and after large excitations such as the Mw 9 Tohoku earthquake. Conditioning effects are observed during the backbone recovery of aftershocks sequences. The results are extended to different building typologies, analyzing the influence of structural material and loading features, notably strain-rates. Finally, some general conclusions are presented, together with a perspective work using machine learning to predict building response based on nonlinear elastic signatures.

A systematic and versatile method to derive accurate and efficient Circuit Oriented Large Signal Average Models (COLSAMs) that approximate the slow dynamics manifold of the moving average values of the relevant state variables for Pulse-Width Modulated (PWM) dc to dc and three-phase to dc power converters is developed. These COLSAMs can cover continuous conduction mode (CCM) as well as discontinuous conduction mode (DCM) of operation and they are over one order of magnitude cheaper, computation wise, than the switching models. This method leads primarily to simple and effective input-output oriented models that represent transfer as well as loading characteristics of the converter. Sine these models consist of time invariant continuous functions they can be linearized at an operating point in order to obtain small-signal transfer functions that approximate the dynamics of the original PWM system around an orbit. The models are primarily intended for software circuit simulators (i.e. Spice derived types, Saber, Simplorer, etc), to take advantage of intrinsic features such as transient response, linearization, transfer function, harmonic distortion calculations, without having to change simulation environment. Nevertheless, any mathematics simulator for ordinary differential equations can be used with the set of equations obtained through application of Kirchoff's laws to the COLSAMs. Furthermore, the COLSAMs provide physical insight to help with power stage and control design, and they allow easy interconnection among themselves, as well as with switching models, for complete analysis at different scales (time, signal level, complexity; interconnectivity). A new average model for the Zero-Voltage Switched Full-Bridge (ZVS-FB) PWM Converter is developed with the above method and its high accuracy is verified with simulations from a switching behavioral model for several circuit component values for both CCM and DCM. Intrinsic positive damping effects and special delay characteristics created by an energy holding element in a saturable reactor-based Zero-Voltage Zero-Current Switched Full-Bridge (ZVZCS-FB) PWM converter are explained for the first time by a new average model. Its large signal predictions match very well those from switch model simulations whereas its small-signal predictions are verified with experimental results from 3.5 kW prototype modules. The latter are used in a multi-module converter to supply the DC power bus in and aircraft. The design of control loops for the converter is based on the new model and its linearization. The ZVZCS-FB PWM converter's average model above is extended to deal with interconnection issues and constraints in a Quasi-Single Stage (QSS) Zero-Voltage Zero-Current Switched (ZVZCS) Three-Phase Buck Rectifier. The new model reveals strong nonlinear transfer characteristics for standard Space Vector Modulation (SVM), which lead to high input current distortion and output voltage ripple inadmissible in telecommunications applications. Physical insight provided by this average model led to the development of a combined modified SVM and feed-forward duty-cycle compensation scheme to reliably minimize the output voltage ripple. Experimental results from a 6 kW prototype validate large signal model for standard and modified SVM, with and without duty-cycle compensation scheme. / Ph. D.

This thesis, which contains fourteen chapters in two parts, presents theoretical and computer simulation studies of dynamics in supercooled liquids and thermotropic liquid crystals. These two apparently diverse physical systems are unified by a startling similarity in their complex slow dynamics. Part I consists of six chapters on supercooled liquids while Part II comprises seven chapters on thermotropic liquid crystals. The fourteenth chapter provides a concluding note. Part I starts with an introduction to supercooled liquids given in chapter 1. This chapter discusses basic features of supercooled liquids and the glass transition and portrays some of the theoretical frameworks and formalisms that are widely recognized to have contributed to our present understanding. Chapter 2 introduces a new model of binary mixture in order to study dynamics across the supercooled regime. The system consists of an equimolar mixture of the Lennard-Jones spheres and the Gay-Berne ellipsoids of revolution, and thus one of its components has orientational degrees of freedom (ODOF). A decoupling between trans-lational diffusion and rotational diffusion is found to occur below a temperature where the second rank orientational correlation time starts showing a steady deviation from the Arrhenius temperature behavior. At low temperatures, the optical Kerr effect (OKE) signal derived from the system shows a short-to-intermediate time power law decay with a very weak dependence on temperature, if at all, of the power law exponent as has been observed experimentally. At the lowest temperature investigated, jump motion is found to occur in both the translational and orientational degrees of freedom. Chapter 3 studies how the binary mixture, introduced in the previous chapter, explores its underlying potential energy landscape. The study reveals correlations between the decoupling phenomena, observed almost universally in supercooled molecular liquids, and the manner of exploration of the energy landscape of the system. A significant deviation from the Debye model of rotational diffusion in the dynamics of ODOF is found to begin at a temperature at which the average inherent structure energy of the system starts falling as the temperature decreases. Further, the coupling between rotational diffusion and translational diffusion breaks down at a still lower temperature, where a change occurs in the temperature dependence of the average inherent structure energy. Chapters 4-6 describe analytical and numerical approaches to solve kinetic models of glassy dynamics for various observables. The β process is modeled as a thermally activated event in a two-level system and the a process is described as a β relaxation mediated cooperative transition in a double-well. The model resembles a landscape picture, conceived by Stillinger [Science 267, 1935 (1995)], where the a process is assumed to involve a concerted series of the β processes, the latter being identified as elementary relaxations involving transitions between contiguous basins. For suitable choice of parameter values, the model could reproduce many of the experimentally observed features of anomalous heat capacity behavior during a temperature cycle through the glass transition as described in chapter 4. The overshoot of the heat capacity during the heating scan that marks the glass transition is found to be caused by a delayed energy relaxation. Chapter 5 shows that the model can also predict a frequency dependent heat capacity that reflects the two-step relaxation behavior. The high-frequency peak in the heat capacity spectra appears with considerably larger amplitude than the low-frequency peak, the latter being due to the a relaxation. The model, when simplified with a modified description of the a process that involves an irreversible escape from a metabasin, can be solved analytically for the relaxation time. This version of the model captures salient features of the structural relaxation in glassy systems as described in chapter 6. In Part II, thermotropic liquid crystals are studied in molecular dynamics simulations using primarily the family of the Gay-Berne model systems. To start with, chapter 7 provides a brief introduction to thermotropic liquid crystals, especially from the perspective of the issues discussed in the following chapters. This chapter ends up with a detail description of the family of the Gay-Berne models. Chapter 8 demonstrates that a model system for calamitic liquid crystal (comprising rod-like molecules) could capture the short-to-intermediate time power law decay in the OKE signal near the isotropic-nematic (I-N) phase transition as observed experimentally. The single-particle second rank orientational time correlation function (OTCF) for the model liquid crystalline system is also found to sustain a power law decay regime in the isotropic phase near the I-N transition. On transit across the I-N phase boundary, two power law decay regimes, separated by a plateau, emerge giving rise to a step-like feature in the single-particle second rank OTCF. When the time evolution of the rotational non-Gaussian parameter is monitored as a diagnostic of spatially heterogeneous dynamics, a dominant peak is found to appear following a shoulder at short times, signaling the growth of pseudonematic domains. These observations are compared with those relevant ones obtained for the supercooled binary mixture, as discussed in chapter 2, in the spirit of the analogy suggested recently by Fayer and coworkers [J. Chem. Phys. 118, 9303 (2003)]. In chapter 9, orientational dynamics across the I-N transition are investigated in a variety of model systems of thermotropic liquid crystals. A model discotic system that consists of disc-like molecules as well as a lattice system have been considered in the quest of a universal short-to-intermediate time power law decay in orientational relaxation, if any. A surprisingly general power law decay at short to intermediate times in orientational relaxation is observed in all these systems. While the power law decay of the OKE signal has been recently observed experimentally in calamitic systems near the I-N phase boundary and in the nematic phase by Fayer and coworkers [J. Chem. Phys. 116, 6339 (2002), J. Phys. Chem. B 109, 6514 (2005)], the prediction for the discotic system can be tested in experiments. Chapter 10 presents the energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals by determining the inherent structures of a family of one-component Gay-Berne model systems. This study throws light on the interplay between the orientational order and the translational order in the mesophases the systems exhibit. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures. The inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the I-N transition and this breakdown is found to occur at a temperature below which the system explores increasingly deeper potential energy minima. There exists a remarkable similarity in the manner of exploration of the potential energy landscape between the Gay-Berne systems studied here and the well known Kob-Andersen binary mixture reported previously [Nature, 393, 554 (1998)]. In search of a dynamical signature of the coupling between orientational order and translational order, anisotropic translational diffusion in the nematic phase has been investigated in the Gay-Berne model systems as described in chapter 11. The translational diffusion coefficient parallel to the director D// is found to first increase and then decrease as the temperature drops through the nematic phase. This reversal occurs where the smectic order parameter of the underlying inherent structures becomes significant for the first time. The non-monotonic temperature behavior of D// can thus be viewed from an energy landscape analysis as a dynamical signature of the coupling between orientational and translational order at the microscopic level. Such a view is likely to form the foundation of a theoretical framework to explain the anisotropic translation diffusion. Chapter 12 investigates the validity of the Debye model of rotational diffusion near the I-N phase boundary with a molecular dynamics simulation study of a Gay-Berne model system for calamitic liquid crystals. The Debye model is found to break down near the I-N phase transition. The breakdown, unlike the one observed in supercooled molecular liquids where a jump diffusion model is often invoked, is attributed to the growth of orientational pair correlation. A mode-coupling theory analysis is provided in support of the explanation. Chapter 13 presents a molecular dynamics study of a binary mixture of prolate ellipsoids of revolution with different aspect ratios interacting with each other through a generalized Gay-Berne potential. Such a study allows to investigate directly the aspect ratio dependence of the dynamical behavior. In the concluding note, chapter 14 starts with a brief summary of the outcome of the thesis and ends up with suggestion of a few relevant problems that may prove worthwhile to be addressed in future.

An increase in the co-operativity in the motion of particles and a growth of a suitably defined dynamical correlation length seem to be generic features exhibited by all liquids upon supercooling. These features have been observed both in experiments and in numerical simulations of glass-forming liquids. Specially designed NMR experiments have estimated that the rough magnitude of this correlation length is of the order of a few nanometers near the glass transition. Simulations also predict that there are regions in the system which are more liquid-like than other regions. A complete theoretical understanding of this behaviour is not available at present. In recent calculations, Berthier, Biroli and coworkers [1, 2] extended the simple mode coupling theory (MCT) to incorporate the effects of dynamic heterogeneity and predicted the existence of a growing dynamical correlation length associated with the cooperativity of the dynamics. MCT also predicts a power law divergence of different dynamical quantities at the mode coupling temperature and at temperatures somewhat higher than the mode coupling temperature, these predictions are found to be consistent with experimental and simulation results. The system size dependence of these quantities should exhibit finite size scaling (FSS) similar to that observed near a continuous phase transition in the temperature range where they show power law growth. Hence we have used the method of finite size scaling in the context of the dynamics of supercooled liquids. In chapter 2, we present the results of extensive molecular dynamics simulations of a model glass forming liquid and extract a dynamical correlation length ξ associated with dynamic heterogeneity by performing a detailed finite size scaling analysis of a four-point dynamic susceptibility χ4(t) [3] and the associated Binder cumulant. We find that although these quantities show the “normal” finite size scaling behaviour expected for a system with a growing correlation length, the relaxation time τ does not. Thus glassy dynamics can not be fully understood in terms of “standard” critical phenomena. Inspired by the success of the empirical Adam-Gibbs relation [4] which relates dynamics with the configurational entropy, we have calculated the configurational entropy for different system sizes and temperatures to explain the nontrivial scaling behaviour of the relaxation time. We find that the behaviour of the relaxation time τ can be explained in terms of the Adam-Gibbs relation [4] for all temperatures and system sizes. This observation raises serious questions about the validity of the mode coupling theory which does not include the effects of the potential energy (or free energy) landscape on the dynamics. On the other hand, in the “random first order transition” theory (RFOT), introduced by Wolynes and coworkers [5], the configurational entropy plays a central role in determining the dynamics. So we also tried to explain our simulation results in terms of RFOT. However, this interpretation has the drawback that the value of one of the exponents of this theory extracted from our numerical results does not satisfy an expected physical bound, and there is no clear explanation for the obtained values of other exponents. Thus we find puzzling values for the exponents relevant to the applicability of RFOT, which are in need of explanation. This can be due to the fact that RFOT focuses only near the glass transition, while all our simulation results are for temperatures far above the glass transition temperature (actually, above the mode coupling temperature). Interestingly, results similar to ours were obtained in a recent analysis [6] of experimental data near the laboratory glass transition, on a large class of glass-forming materials. Thus right now we do not have any theory which can explain our simulation data consistently from all perspectives. There have been some attempts to extend the RFOT analysis to temperatures above the mode coupling temperature [7, 8] and to estimate a length scale associated with the configurational entropy at such temperatures. We compare our results with the predictions arising from these analyses. In chapter 3, we present simulation results that suggest that finite size scaling analysis is probably the only feasible method for obtaining reliable estimates of the dynamical correlation length for supercooled liquids. As mentioned before, although there exists a growing correlation length, the behaviour of all measured quantities (specifically, the relaxation time) is not in accordance with the behaviour expected in “standard” critical phenomena. So one might suspect the results for the correlation length extracted from the scaling analysis. To find out whether the results obtained by doing finite size scaling are correct, we have done simulations of very large system sizes for the same model glass forming liquid. In earlier studies, the correlation length has been extracted from the wave vector dependence of the dynamic susceptibility in the limit of zero wave vector, but to estimate the correlation length with reasonable accuracy one needs data in the small wave vector range. This implies that one needs to simulate very large systems. But as far as we know, in all previous studies typical system sizes of the order of 10, 000 particles have been used to do this analysis. In this chapter we show by comparing results for systems of 28, 000 and 350, 000 particles that these previous estimates are not reliable. We also show that one needs to simulate systems with at least a million particles to estimate the correlation length correctly near the mode coupling temperature and this size increases with decreasing temperature. We compare the correlation length obtained by analyzing the wave vector dependence of the dynamic susceptibility for a 350, 000particle system with the results obtained from the finite size scaling analysis. We were only able to compare the results in the high temperature range due to obvious reasons. However the agreement in the high temperature range shows that the finite size scaling analysis is robust and also establishes the fact that finite size scaling is the only practical method to extract reliable correlation lengths in supercooled liquids. In chapter 4, we present a free energy landscape analysis of dynamic heterogeneity for a monodisperse hard sphere system. The importance of the potential energy landscape for particles interacting with soft potentials is well known in the glass community from the work of Sastry et al. [9] and others, but the hard sphere system which does not have any well defined potential energy landscape also exhibits similar slow dynamics in the high density limit. Thus it is not clear how to treat the hard sphere systems within the same energy landscape formalism. Dasgupta et al. [10, 11, 12, 13, 14, 15] showed that one can explain the slow dynamics of these hard core systems in term of a free energy landscape picture. They and other researchers showed that these system have many aperiodic local minima in its free energy landscape, with free energy lower than that of the liquid. Using the Ramkrishnan-Yussouff free energy functional, we have performed multi parameter variational minimizations to map out the detailed density distribution of glassy free energy minima. We found that the distribution of the widths of local density peaks at glassy minima is spatially heterogeneous. By performing hard sphere event driven molecular dynamics simulation, we show that there exists strong correlation between these density inhomogeneity and the local Debye-Waller factor which provides a measure of the dynamic heterogeneity observed in simulations. This result unifies the system of hard core particles with the other soft core particles in terms of a landscapebased description of dynamic heterogeneity. In chapter 5, we extend the same free energy analysis to a polydisperse system and show that there is a critical polydispersity beyond which the crystal state is not stable and glassy states are thermodynamically stable. We also found a reentrant behaviour in the liquid-solid phase transition within this free-energy based formalism. These results are in qualitative agreement with experimental observations for colloidal systems.