No description available.

The present thesis reports computer simulation studies of several phase transition related phenomena in a range of soft-condensed matter systems. A coherent unifying theme of the thesis is the understanding of dynamics of phase transitions through free energy calculations using recently developed efficient non-Boltzmann sampling methods. Based on the system/phenomena of interest, the thesis has been classified into four major parts: I. Isotropic-nematic (IN) phase transition in liquid crystals. II. Nucleation phenomena in gas-liquid transition with particular emphasis on the systems close to the spinodal curve. III. Collapse transition in linear hydrocarbon (n-alkane) chains for a varying range of length, solvent and temperature. IV. Crystallization of unbranched polymer chains in dilute solution, with particular emphasis on the temperature dependent crossover between the rod-like crystalline state and spherical molten globule state. The thesis has been further divided into ten chapters running through the four parts mentioned before. In the following we provide a brief chapter-wise outline of the thesis. Part I deals with the power law relaxation and glassy dynamics in thermotropic liquid crystals close to the IN transition and consists of two chapters. To start with, Chapter I.1 provides an introduction to thermotropic liquid crystals. Here we briefly introduce various liquid crystalline phases, the order parameter used to characterize the IN transition, a few well established theoretical models, and we conclude with describing the recent experimental and computer simulation studies that have motivated the work described in the next chapter. In Chapter I.2, we present our molecular dynamics simulation studies on single particle and collective orientational dynamics across the IN transition for Lebwohl Lasher model, which is a well-known lattice model for thermotropic liquid crystals. Even this simplified model without any translational degrees of freedom successfully captures the short-tointermediate time power law decay recently observed in optical heterodyne detected optical Kerr effect (OHDOKE) measurements near the IN transition. The angular velocity time correlation function also exhibits a rather pronounced power law decay near the IN boundary. In the mean squared angular displacement at comparable time scales, we observe the emergence of a sub-diffusive regime which is followed by a super-diffusive regime before the onset of the longtime diffusive behavior. We observe signature of dynamical heterogeneity through pronounced non-Gaussian behavior in the orientational motion particularly at lower temperatures. Interestingly, this behavior closely resembles what is usually observed in supercooled liquids. We obtain the free energy as a function of orientational order parameter by the use of recently developed transition matrix Monte Carlo (TMMC) method. The free energy surface is flat for the system considered here and the barrier between isotropic and nematic phases is vanishingly small for this weakly first-order transition, hence allowing for large scale, collective, and correlated orientational density fluctuations. We attribute this large scale fluctuations as the reason for the observed power law decay of the orientational time correlation functions. Part II consists of three chapters, where we focus on the age old problem of nucleation and growth, both from the perspective of thermodynamics and kinetics. We account for the rich history of the problem in the introductory Chapter II.1. In this chapter we describe various types and examples of the nucleation phenomena, and a brief account of the major theoretical approaches used so far. We begin with the most successful Classical Nucleation Theory (CNT), and then move on to more recent applications of Density Functional Theory (DFT) and other mean-field types of models. We conclude with a comparison between the experiments, theories and computational studies. In the next chapter (Chapter II.2) we attempt to elucidate the mechanism of nucleation near the gas-liquid spinodal from a microscopic point of view. Here we construct a multidimensional free energy surface of nucleation of the liquid phase from the parent supercooled and supersaturated vapor phase near the gas-liquid spinodal. In particular, we remove the Becker-Doring constraint of having only one growing cluster in the system. The free energy, as a function of the size of the largest cluster, develops a pronounced minimum at a subcritical cluster size close to the spinodal. This signifies a two step nature of the process of nucleation, where the rapid formation of subcritical nuclei is followed by further growth by slower density fluctuations on an uphill free energy surface. An alternative free energy pathway involving the participation of many subcritical clusters is envisaged near the spinodal where the growth of the nucleus is found to be promoted by a coalescence mechanism in contrast to the single particle addition assumption within CNT. The growth of the stable phase becomes progressively collective and spatially diffuse, and the significance of a “critical nucleus” is lost for deeper quenches. In this chapter we present our studies both in 3dimensional Lennard-Jones (LJ) system and Ising model (both 2and 3dimensions). Our general findings seem to be independent of the model chosen. While the previous chapter focuses on relatively well-studied 3-dimensional (3D) LJ system, in Chapter II.3 we present our studies on the characteristics of the nucleation phenomena in 2dimensional (2D) Lennard-Jones fluid. To the best of our knowledge this is the first extensive computer simulation study to check the accuracy of CNT in 2D. Using various Monte Carlo methods, we calculate the free energy barrier for nucleation, line tension, and bulk densities of equilibrium liquid and vapor phases, and also investigate the size and shape of the critical nucleus. The study is carried out at an intermediate level of supersaturation (away from the spinoidal limit). In 2D, a surprisingly large cutoff (rc ≥ 7.0σ where σ is the diameter of LJ particles) in the truncation of the LJ potential is required to obtain converged results. A lower cutoff leads to a substantial error in the values of the line tension, nucleation barrier, and characteristics of the critical cluster. Note that typically 2.5σ is sufficient for 3D LJ fluids. We observe that in 2D system CNT fails to provide a reliable estimate of the free energy barrier. While it is known to slightly overestimate the nucleation barrier in 3D, it underestimates the barrier by as much as 50% at the saturation ratio S = 1.1(defined as S = P/Pc, where Pc is the coexistence pressure) and at the reduced temperature T* = 0.427(defined as T* = KBT/ ε, where ε is the depth of the potential well). The reason for the marked inadequacy of the CNT in 2D can be attributed to the non-circular nature of the critical clusters. Although the shape becomes increasingly circular and the clusters become more compact with increase in cutoff radius, an appreciable non-circular nature remains even without any cutoff to make the simple CNT inaccurate. Part III again consists of three chapters and focuses on the conformational equilibria. Collapse transition and self-organized structures of n-alkanes in solution. In Chapter III.1 we carry out a brief survey of the existing theories of polymer in solution, with particular emphasis on the collapse process in poor solvents. We also introduce the concept of “hydrophobicity” and “hydrophobic collapse”, which is now a subject enormous interest, partly because it my help in understanding the initial processes involved in protein folding. We briefly discuss the subject of formation of beautiful self-organized structures by block copolymers, and also simple homopolymers which is essentially the focus of the work embodied in the next two chapters. In Chapter III.2 we demonstrated a chain length dependent crossover in the structural properties of linear hydrocarbon (n-alkane) chains using detailed atomistic simulations in explicit water. We identify a number of exotic structures o the polymer chain through energy minimization of representative snapshots collected from molecular dynamics trajectory. While the collapsed state is ring-like(circular) for small chains(CnH2n+2; n ≤ 20) and spherical for very long ones( n = 100), we find the emergence of ordered helical structures at intermediate lengths (n ~ 40). We find different types of disordered helices and toroid-like structures at n = 60. We also report a sharp transition in the stability of the collapsed state as a function of the chain length through relevant free energy calculations. While the collapsed state is only marginally metastable for C20H42, a clear bistable free energy surface emerges only when the chain is about 30 monomers long. For n = 30, the polymer exhibits an intermittent oscillation(characterized by well-developed 1/f noise, where f is the frequency ) between the collapsed and the coil structures, characteristic of two stable states separated by a small barrier. This appears to support a weakly first order phase transition between the extended and the collapsed states. Chapter III.3 extends the study of previous chapter to much longer chains (n ≥ 100), which irreversibly collapse in water into globular forms. Even though the collapsed form has a nearly spherical shape, close inspection shows a propensity towards local ordering in the alignment of the polymer segments. This tendency to maintain alignment in order to maximize the number of contacts leads to a core-shell like structure, where the shell is often characterized by a bent rod-like shape consisting of two adjacent segments running in parallel. A key event associated with the initial stage of collapse seems to be the formation of a skewed ring (or loop) that serves as a “nucleation center” for rest of the chain to collapse into. Time evolution of the radial distribution function of water surrounding the polymer, shows that the density of neighboring water decreases by only about 15-20% from that of bulk water. Even though interior of the ting-like structures is fully devoid of water, solvent accessible surface representation shows that these regions are geometrically/spatially inaccessible to water molecules. We suggest that the role of water is to stabilize such ring-like structures once formed by natural conformational fluctuations of the polymer chain. This view is confirmed by observation of spontaneous formation and melting away of such ring-like entities in a polar aprotic solvent(DMSO). We also comment on the role of the flexibility of polymer chains in determining the collapse kinetics. The last part(Part IV) of the thesis consists of two chapters that deal with the crystallization of linear polymer chains from dilute solution. The way long chain polymers crystallize is drastically different from their small molecule counterparts due to their topological connectivity. Linear polymers often crystallize from dilute solution in the form of thin lamellae with well-defined crystallographic features. In Chapter IV.1 we briefly survey the current theoretical understanding and confusions associated with the highly debated field of polymer crystallization. While the last few decades have seen the development of many successful phenomenological theories, the molecular mechanism of formation of such self-organized lamellae is extremely complex and very poorly understood. There are clearly two distinct steps in polymer crystallization. Firstly, the individual linear polymers must self-organize into bundles of somewhat regular structures. These structures then further aggregate to lamellar form and crystallize into a lattice. In this respect , it has marked similarity to the problem of protein crystallization. In chapter IV.2 we present Brownian dynamics simulation studies of a single polythelene chain of length 500. Such systems can reasonably mimic the process of crystallization from dilute solutions. Our simulations could successfully reproduce some of the interesting phenomena observed in experiments and very recent computer simulation studies, including multi-center nucleation of rod-like structures within a single polymer chain, an inverse relation between lamellar thickness and temperature etc. But our primary focus has been to understand the nature of the phase transition as one traverses along the melting temperature and the underlying free energy surface. Near the melting temperature we observe a very intriguing fluctuation between the disordered molten globule state and the ordered rod-like crystalline, where these two forms have highly different shape and structure. These fluctuations have strong signature of 1/f noise or intermittency. This clearly indicates the existence of a weakly first order transition, where two widely different states with large difference in values of order parameter are separated by a rather small free energy barrier. This can be related to the experimentally observed density fluctuations that resemble spinodal decomposition. It is important to note that very similar fluctuations have been observed in our previous studies on liquid crystals (Chapter 1.2) and intermediate sized alkalines in water(Chapter III.2) that signifies a universal underlying energy landscape for these systems. We have discussed the scope of future work at the end of each chapter whenever appropriate.

No description available.

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.