Phase field crystal models and fast dynamics

Abstract

Emergent microscale structural properties of crystalline materials have been puzzling scientist since the dawn of modern solid state physics. The microstructure of crystalline materials is affected by disorder: topological defects such as vacancies and grain boundaries can fundamentally change both the electronic and structural properties of materials. The microstructure is determined by the manufacturing process, the most famous being solidification from melt. Understanding the solidification process can be used to produce materials with specific properties in a controlled way. Classical density functional theory (CDFT) and its computationally more efficient counterpart, the phase field crystal (PFC) model, have emerged as efficient tools for studying microscopic crystalline properties of materials at diffusive time scales, avoiding the unnecessary short time scale of thermal vibrations. In this thesis we further develop the PFC methodology. First, we use the DFT and PFC methods to calculate the free energy of the solid–liquid boundary of the well-known Yukawa system. In the second part – the main body of the work – we develop fast dynamics for the PFC framework. The last part is dedicated to computational work: we propose an efficient numerical scheme for solving the dynamics of the PFC amplitude system. Developing fast dynamics tackles an old problem of diffusive dynamical systems. The presence of collective vibrations of atoms is characteristic of metallic materials. In principle all the dynamical theories that describe elasticity of materials should also include these vibrations. However, diffusive dynamics are unable to describe such vibrations. To this end, we propose two approaches. In the first approach the elastic excitations that create these vibrations are explicitly equilibrated. In the second approach we develop a hydrodynamic theory for the PFC framework based on physical conservation laws. In both cases we demonstrate the viability and the importance of the approach. The latter theory is more general and can be used to describe a larger spectrum of fast phenomena such as advective mass transport. The theories developed in this thesis allow studying previously inaccessible problems such as fast solidification. In addition they should produce better understanding of many problems where fast dynamics are present. These include problems where relaxation of elastic excitations is important, such as polycrystalline coarsening, and problems where mass transport is important e.g. certain type of dendritic solidification. Another advantage of the methods presented here is that the theoretical framework is general and can be extended to different systems.

Aineen emergentit mikrotason rakenteelliset ominaisuudet ovat kutkuttaneet tieteilijöitä nykyaikaisen kiinteän aineen fysiikan alkupäivistä alkaen. Epäjärjestys vaikuttaa aineen mikroskooppiseen rakenteeseen: topologiset virheet, kuten vakanssit ja kiderajapinnat voivat muuttaa kokonaisvaltaisesti sekä aineen elektroniset että rakenteelliset ominaisuudet. Valmistusprosessi – joista kuuluisin on kenties sulasta aineesta jähmettäminen – määrittää aineen mikrorakenteen. Jähmettymisprosessin syvä ymmärtäminen mahdollistaa syntyvän materiaalin erityislaatuisten ominaisuuksien säätelyn. Klassinen tiheysfunktionaaliteoria (KTFT) ja siitä johdettu laskennallisesti tehokkaampi faasikenttäkidemalli (FKK-malli) ovat osoittautuneet tehokkaiksi keinoiksi kiteiden mikroskooppisten rakenteiden tutkimiseen diffusiivisilla aikaskaaloilla. Nämä metodit välttävät diffusiivisten ilmiöiden kannalta tarpeettoman atomien värähtelyiden laskemisen. Tässä väitöskirjassa kehitämme FKK-metodologiaa. Ensimmäisessä osassa laskemme tunnetun Yukawa-systeemin neste–kiinteä -rajapinnan vapaaenergian käyttäen KTFT ja FKK -metodeja. Toisessa osassa, joka kattaa merkittävimmän osan tästä väitöstyöstä, kehitämme nopean dynamiikan FKK-metodologialle. Väitöskirjan viimeisessä osassa esittelemme tehokkaan numeerisen menetelmän FKK-mallin perinteisen dynamiikan laskemiseksi. Nopea dynamiikka auttaa diffusiivisten dynaamisten systeemien kanssa ilmenevän vanhan ongelman ratkaisemisessa. Atomien kollektiivinen värähtely on yksi metallien tärkeimmistä ominaisuuksista. Periaatteessa kaikkien aineen elastisia ominaisuuksia kuvaavien dynaamisten teorioiden pitäisi sisällyttää kyseinen värähtely, mutta diffusiivinen dynamiikka ei tähän pysty. Esitämme kaksi tapaa tämän ongelman ratkaisemiseksi. Ensimmäisessä lähestymistavassa tasapainotamme elastiset jännitykset, jotka synnyttävät edellä mainitut värähtelyt. Toisessa lähestymistavassa kehitämme hydrodynaamisen teorian FKK-malleille. Osoitamme lähestymistapojen käyttökelpoisuuden ja tärkeyden. Jälkimmäinen lähestymistapa on yleisempi ja pystyy kuvaamaan useampia nopeaan dynamiikkaan liittyviä ilmiöitä, kuten massan kuljetus nopeuskentän avulla. Tässä työssä kehitetyt teoriat mahdollistavat uusien ongelmien, kuten nopean jähmettymisen, tutkimisen. Teoriat mahdollistavat myös monien jo tutkittujen ongelmien paremman ymmärtämisen. Tällaisista ongelmista mainittakoon esimerkiksi tapaukset, joissa elastiset jännitykset ovat olennaisia, kuten monikiteisen aineen karkeistuminen sekä ongelmat, joissa pitää huomioida massan kuljetus, kuten tietynlainen dendriittinen jähmettyminen. Yksi olennainen tässä työssä esitettyjen teorioiden ominaisuus on, että ne ovat yleisyytensä vuoksi helposti sovellettavissa moniin erilaisiin ongelmiin.