THE VALIDITY OF CLASSICAL NUCLEATION THEORY AND ITS APPLICATION TO DISLOCATION NUCLEATION
A DISSERTATION SUBMITTED TO THE DEPARTMENT OF PHYSICS AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY
Seunghwa Ryu August 2011
© 2011 by Seunghwa Ryu. All Rights Reserved. Re-distributed by Stanford University under license with the author.
This work is licensed under a Creative Commons Attribution- 3.0 United States License. http://creativecommons.org/licenses/by/3.0/us/
This dissertation is online at: http://purl.stanford.edu/rx036ms4124
ii I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Wei Cai, Primary Adviser
I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Douglas Osheroff, Co-Adviser
I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Paul McIntyre
I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
William Nix
Approved for the Stanford University Committee on Graduate Studies. Patricia J. Gumport, Vice Provost Graduate Education
This signature page was generated electronically upon submission of this dissertation in electronic format. An original signed hard copy of the signature page is on file in University Archives.
iii Abstract
Nucleation has been the subject of intense research because it plays an important role in the dynamics of most first-order phase transitions. The standard theory to describe the nucleation phenomena is the classical nucleation theory (CNT) because it cor- rectly captures the qualitative features of the nucleation process. However potential problems with CNT have been suggested by previous studies. We systematically test the individual components of CNT by computer simulations of the Ising model and find that it accurately predicts the nucleation rate if the correct droplet free energy computed by umbrella sampling is provided as input. This validates the fundamental assumption of CNT that the system can be coarse grained into a one dimensional Markov chain with the largest droplet size as the reaction coordinate. Employing similar simulation techniques, we study the dislocation nucleation which is essential to our understanding of plastic deformation, ductility, and me- chanical strength of crystalline materials. We show that dislocation nucleation rates can be accurately predicted over a wide range of conditions using CNT with the ac- tivation free energy determined by umbrella sampling. Our data reveal very large activation entropies, which contribute a multiplicative factor of many orders of mag- nitude to the nucleation rate. The activation entropy at constant strain is caused by thermal expansion, with negligible contribution from the vibrational entropy. The ac- tivation entropy at constant stress is significantly larger than that at constant strain, as a result of thermal softening. The large activation entropies are caused by anhar- monic effects, showing the limitations of the harmonic approximation widely used for rate estimation in solids. Similar behaviors are expected to occur in other nucleation processes in solids.
iv Acknowledgements
First of all, I am very much indebted to my principal adviser, Professor Wei Cai. It is of great fortune for me to work with such a bright and gentle person who I want to follow as a role model both as a scientist and as a gentleman. I have learned how to approach a scientific problem and how to tackle it: under his guidance, my random ideas transformed into a well defined research project, and a seemingly formidable problem into a series of small problems that can be handled systematically. Pro- fessor Cai has also helped me every single step that I need to grow as a scientist, such as writing a concise paper and delivering an insightful presentation. In addition to academic advices, I have also learned the virtues of a gentleman: he has consis- tently shown positive attitude on life, humility in the quest of knowledge, respect on other people, and dedication to his family. I am especially grateful for having many discussions with him on non-academic subjects regarding various aspects of life and being able to listen to his advices. Past four years that I worked with Professor Cai have been one of the most important periods in my life in which I have grown both intellectually and mentally. I owe very special thanks to my co-adviser, Professor Douglas Osheroff, under whom I had worked during the first three years of my graduate study. I had no experience on the experimental physics when I arrived at Stanford, and joined his group in the hope that I could learn a completely new subject under the guidance of a famous Nobel Laureate. Indeed, I have learned a lot from his wizardly expertise and rich experiences on the low temperature physics experiment. I still remember the first moment when I saw the signature of He-3 superfluid transition with full of joy, after several months of struggles to fix the dilution fridge. My experience in
v Professor Osheroff’s group has aided and will continue to aid me to communicate and cooperate with experimentalists. I deeply appreciate the generosity he showed me when I decided to leave his group after I realized my natural preference toward theoretical studies. Since then, gratefully, he has been my co-adviser and given me precious advices on research, career, and life in general. His dedication to science education for general public and humble attitude have shown me the way I want to follow when becoming a senior scientist in the far future. I would like to thank Professor William Nix and Paul McIntyre to serve on my thesis committee. The dislocation course that I took from Professor Nix and the kinetic process course from Professor McIntyre have provided me the theoretical basis for the dissertation project. I appreciate Professor Nix for the discussion and valuable assessments on the dislocation nucleation study. He exemplifies the ideal life as a senior professor: he still actively works and enjoys the life at the same time, and willingly shares his valuable time for helping students and young faculties. I would like to thank Professor McIntyre for the invitation to his nanowire group meeting and insightful advices on the nanowire growth simulation project, another branch of my doctoral study. He exemplifies the quality of a real professional: he gives critical assessments on students working in various subjects with his deep and broad understanding in materials science research and organizes collaborations with several groups very efficiently. I also want to thank Professor Evan Reed for serving as the chair for my thesis defense meeting. I am happy to thank two special seniors in our group, Dr. Keonwook Kang and Dr. Eunseok Lee. I was going through a difficult time when I moved to Cai group due to the anxiety from starting a completely new field in the midst of graduate study and the ignorance in computational work. Without their moral support and help on the technical skills on computer simulations, I would have not succeeded in changing my research field so smoothly. I would like to thank all Cai group mem- bers who shared valuable discussions on my project. And many thanks to former lab mates in Osheroff group for the training on the low temperature physics exper- iments and to McIntyre group members for sharing interesting experimental results on semiconductor nanowire growth.
vi Besides spending time in the lab doing research, I have been nourished by having good friends and sharing unforgettable memories with them. I would like to thank my fellow KAIST alumni at Stanford, friends in Cornerstone Community Church (special thanks to Dr. Jungjoon Lee), Professor Lew group members, fellow Korean students in physics and mechanical engineering departments, and all other close friends not included in these groups. I have been so comfortable and relaxed with my friends having many trips and parties, playing sports and games, tasting delicious foods and liquors, and going museums and concerts together. Their encouragements and advices on many aspects of life are also priceless. I want to thank my home university, Korea Advanced Institute of Science and Technology (KAIST) and the department of physics where I acquired a solid foun- dation in physics as an undergraduate. Special thanks to my undergraduate adviser, Professor Mahn Won Kim, and Professor Hawoong Jeong for their encouragement and invaluable advices in my career. I appreciate the financial supports from the Stanford Graduate Fellowship and the Korea Science and Engineering Foundation Fellowship, which allowed me to choose research projects more freely. Lastly, I would like to thank my family for all their love, support, and encourage- ment. I do not know how I can repay what I have received from my parents for the rest of my life. The devotion, patience, and responsibility that they have shown in their life have been and will be the source of my strength. I thank my younger brother who kept encouraging me for the past years. I would like to thank my grandmother in heaven who had dedicated her life to her family and showed me what the true love is. During seven years of graduate study, I have gradually realized that I owe every- thing I accomplished to people around me and how important it is to interact well with others. This precious lesson on the relationship, as well as the expertise I gained in my doctoral research, is the true gem that I will cherish for my life.
vii Preface
During my doctoral study with Professor Wei Cai, I have worked on a diverse spec- trum of research projects in computational physics and computational materials sci- ence. Since it was impossible to assemble all the published works in a coherent manner, I have opted for presenting a tightly knit story with the theme of the com- putational investigation of nucleation phenomena, out of some portion of them. For that purpose, I have written a long introduction to well establish a niche, and a care- ful review on existing theories and experiments, and computational methods, which will help readers better understand the materials presented in this dissertation. As noted in the abstract, the main text of this dissertation addresses two corre- lated projects sharing common theory and numerical algorithms: (1) the test of the classical nucleation theory using the Ising model and (2) the prediction of disloca- tion nucleation rate from the classical nucleation theory using atomistic simulations. Readers who prefer a more distilled presentation are refered to following journal ar- ticles: for the Ising model,
Seunghwa Ryu and Wei Cai,“The Validity of Classical Nucleation Theory for • Ising Models”, Phys. Rev. E (Rapid Communications) 81, 030601 (R) (2010).
Seunghwa Ryu and Wei Cai, “Numerical Tests of Nucleation Theories for the • Ising Models”, Phys. Rev. E 82, 011603 (2010).
The first one is a concise letter containing the crux of the work and the second one is a full paper including more extensive discussions which contains most contents of Chapter 4.
viii for dislocation nucleation,
Seunghwa Ryu, Keonwook Kang and Wei Cai, “Entropic Effect on the Rate of • Dislocation Nucleation”, Proc. Natl. Acad. Sci. USA 108, 5174 (2011).
Seunghwa Ryu, Keonwook Kang, and Wei Cai, “Predicting the Dislocation Nu- • cleation Rate as a Function of Temperature and Stress”, J. Mater. Res. (2011), in press.
Sylvie Aubry, Keonwook Kang, Seunghwa Ryu and Wei Cai, “Energy Barrier • for Homogeneous Dislocation Nucleation: Comparing Atomistic and Continuum Models”, Scripta Mater. 64, 1043 (2011).
The first one is a condensed article highlighting the entropic effect and the second one is a full paper including more in-depth discussions which contains most contents of Chapter 5. Third one considers the comparison between atomistic and contin- uum model. This work is not included since I am not a main contributor, but it is an interesting work that may attracts readers engrossed in dislocation nucleation research. Another major branch of my doctoral study is the simulation of the gold-catalyzed growth of silicon nanowires via vapor-liquid-solid (VLS) mechanism. For this project, thermal properties obtained from existing atomistic models are investigated for gold, silicon, and various other materials. We have devised efficient methods for computing the free energies of solid and liquid alloys, to improve and develop a gold-silicon poten- tial that is fitted to the experimental binary phase diagram. These works are included as appendices that are referred in Chapter 3 where various interatomic potential mod- els and simulation methods are reviewed. Because the series of studies constitute a complete set of story by themselves, readers interested solely in the nanowire growth mechanism can skip the main text of the dissertation and read through Appendices D, E, and F. Each of Appendices D, E, F contains the contents of following three journal articles, with one-to-one correspondence:
Seunghwa Ryu and Wei Cai, “Comparison of Thermal Properties Predicted by • Interatomic Potential Models”, Modell. Simul. Mater. Sci. Eng. 16, 085005
ix (2008).
Seunghwa Ryu, Christopher R. Weinberger, Michael I. Baskes, and Wei Cai, • “Improved Modified Embedded-Atom Method Potentials for Gold and Silicon”, Modell. Simul. Mater. Sci. Eng. 17, 075008 (2009).
Seunghwa Ryu and Wei Cai, “A Gold-Silicon Potential Fitted to the Binary • Phase Diagram”, J. Phys.: Condens. Matter. 22, 055401 (2010).
We have published a pioneering simulation of the silicon nanowire growth using the potential developed in above studies, which is not included in the dissertation due to weaker link with the main text. We investigated the origins of the orientation dependence of nanowire growth rate and the shift of phase diagram in nanoscale. Readers interested in this study are invited to read the following journal article.
Seunghwa Ryu and Wei Cai, “Molecular Dynamics Simulations of Gold-Catalyzed • Growth of Silicon Bulk Crystals and Nanowires”, J. Mater. Res. (2011), in press.
We have also done an interesting piece of work on quantum entanglement for intellectual amusement. A fast algorithm to calculate the entanglement of formation of a mixed state was developed with which we obtain the statistics of the entanglement of formation on ensembles of random density matrices of higher dimensions than possible before. The correlations between the entanglement of formation and other quantities that are easier to compute, such as participation ratio and negativity are studied. The details of this work can be found in the following journal article.
Seunghwa Ryu, Wei Cai and Alfredo Caro, “Quantum Entanglement of Forma- • tion between Qudits”, Phys. Rev. A 77, 052312 (2008).
x Contents
Abstract iv
Acknowledgements v
Preface viii
1 Introduction 1 1.1 Computational Investigation of Nucleation ...... 1 1.2 ScopeoftheDissertation...... 4
2 Background and Motivation 7 2.1 BasicThermodynamicsofNucleation ...... 7 2.2 Nucleation Rate Predictions From Nucleation Theories ...... 11 2.3 NucleationExperiments ...... 14 2.4 Dislocation Nucleation and Materials Strength at Small Scale . . . . 20
3 Computational Methods 25 3.1 IsingModel ...... 26 3.2 InteratomicPotential ...... 29 3.2.1 List of Empirical Potentials ...... 30 3.2.2 The Benchmarks of EAM Copper Potential ...... 34 3.3 Molecular Dynamics Simulation ...... 38 3.4 Monte Carlo Simulation ...... 39 3.5 AdvancedSamplingMethods ...... 41
xi 3.5.1 Forward Flux Sampling ...... 42 3.5.2 Umbrella Sampling ...... 45 3.5.3 Computing Ratefrom Becker-D¨oring Theory ...... 48
4 Numerical Tests of Nucleation Theories 50 4.1 Introduction...... 50 4.2 Nucleation Theories Applied to the Ising Model ...... 53 4.2.1 Becker-D¨oringTheory ...... 53 4.2.2 Langer’s Field Theory ...... 57 4.3 ComputationalMethods ...... 57 4.4 Results...... 59 4.4.1 NucleationRate...... 59 4.4.2 Critical Droplet Size and Shape ...... 61 4.4.3 DropletFreeEnergyof2DIsingModel ...... 63 4.4.4 DropletFreeEnergyof3DIsingmodel ...... 66 4.4.5 EffectiveEntropyofNucleation ...... 70 4.5 SummaryandDiscussion...... 72
5 Predicting the Dislocation Nucleation Rate 75 5.1 Introduction...... 75 5.2 ThermodynamicsofNucleation ...... 78 5.2.1 Activation Free Energies ...... 78 5.2.2 Activation Entropies ...... 80 5.2.3 Difference between the Two Activation Entropies ...... 83 5.2.4 Previous Estimates of Activation Entropy ...... 86 5.3 ComputationalMethods ...... 89 5.3.1 Simulation Cell ...... 89 5.3.2 Nucleation Rate Calculation ...... 92 5.4 Results...... 95 5.4.1 Benchmark with MD Simulations ...... 95 5.4.2 Homogeneous Dislocation Nucleation in Bulk Cu ...... 96 5.4.3 Heterogeneous Dislocation Nucleation in Cu Nano-Rod . . . 99
xii 5.5 Discussion...... 100 5.5.1 Testing the “Thermodynamic Compensation Law” ...... 100 5.5.2 Entropic Effect on Nucleation Rate and Yield Strength . . . . 103 5.6 Summary ...... 106
6 Summary and Outlook 107 6.1 Conclusion...... 107 6.2 FutureWorks ...... 108
A Derivations 111 A.1 Nucleation Rate Prediction from Classical Nucleation Theory . . . . . 112 A.2 NucleationTheorems ...... 117
B More Data on the Nucleation in the Ising Model 121 B.1 AttachmentRate ...... 121 B.2 DropletShape...... 122 B.3 TheConstantTerminDropletFreeEnergy ...... 123 B.4 Free Energy Curves F (n)fortheIsingModel...... 126
C More Discussion on the Dislocation Nucleation 131 σ γ C.1 Equality of Critical Sizes nc and nc ...... 131 C.2 Equality of Activation Gibbs and Helmholtz Free Energies ...... 132
C.3 Physical Interpretation of Activation Entropy Difference ∆Sc ..... 133
C.4 Approximation of Sc(σ) ...... 135 C.5 ActivationFreeEnergyData...... 136 C.6 Activation Volume and Critical Loop Size ...... 139
D Thermal Properties from Interatomic Potentials 142 D.1 Introduction...... 142 D.2 Comparison between Model Predictions and Experiments ...... 144 D.2.1 Semiconductors: SiandGe...... 145 D.2.2 FCCMetals:Au,Cu,AgandPb ...... 147 D.2.3 BCCMetals:Mo,TaandW...... 147
xiii D.3 Free Energy Method for Melting Point Calculation ...... 148 D.3.1 SolidFreeEnergy...... 150 D.3.2 LiquidFreeEnergy ...... 153 D.3.3 Melting Point and Error Estimate ...... 155 D.4 Summary ...... 158 D.5 Error Estimates in Free Energy Calculations ...... 158
E MEAM Potentials for Pure Au and Pure Si 160 E.1 Introduction...... 160 E.2 ProblemStatement ...... 162 E.2.1 Limitations of the MEAM Gold Potential ...... 163 E.2.2 Limitations of MEAM Silicon Potentials ...... 165 E.3 MethodsandResults ...... 167 E.3.1 Multi-body Screening Function ...... 168 E.3.2 PairPotentialandEquationofState ...... 170 E.4 Summary ...... 174 E.5 Multi-bodyScreeningFunction ...... 174 E.6 Further Benchmarks of the MEAM† Potentials...... 177
F Gold-Silicon Binary Potential 180 F.1 Introduction...... 180 F.2 MEAMModelforGoldandSilicon ...... 181 F.2.1 FunctionalForm ...... 181 F.2.2 DeterminingtheParameters ...... 182 F.3 ConstructionofBinaryPhaseDiagram ...... 185 F.3.1 Free Energy of Solid with Impurities ...... 186 F.3.2 Free Energy of Liquid Alloy ...... 188 F.3.3 ConstructionofBinaryPhaseDiagram ...... 191 F.4 Summary ...... 192 F.5 FurtherBenchmarks ...... 194
Bibliography 196
xiv List of Tables
3.1 Lattice properties of Cu predicted by the EAM potential [111]. . . . 34
3.2 The intrinsic stacking fault energy γSF and the unstable stacking fault
energy γUSF of Cu from experiments, EAM potential [111], and ab ini- tio calculation [114]. ab initio data depend on the exchange-correlation functionalsusedinthestudy...... 36
C.1 Data for homogeneous nucleation: σxy in GPa, Ec, E˜c and Fc in eV, + 14 −1 ˜ fc in 10 s . γxy and Γ are dimensionless. The error in Ec is about 0.003 eV, due to the small errors in equilibrating the simulation cell to
achieve the pure shear stress state. The error in Fc is about 0.5 kBT , i.e. approximately 0.01 eV, due to the statistical error in umbrella sam- pling. The error in Zeldovich factor Γ is within 0.01. The attachment ± rate f + has relative error of 50%...... 138 c ± + C.2 Data for heterogeneous nucleation: σzz in GPa, Ec, and Fc in eV, fc 14 −1 in 10 s . γxy and Γ are dimensionless. The error in Fc is about
0.5 kBT , i.e. approximately 0.01 eV, due to the statistical error in umbrella sampling. The error in Zeldovich factor Γ is within 0.01. ± The attachment rate f + has relative error of 50%. Notice that, due c ± to the existence of thermal strain, the elastic strain values are slightly differentatdifferenttemperatures...... 139
xv D.1 Thermal properties of various elements as predicted by several empir- ical potentials and compared with experiments [183, 207, 208]. The
properties include the melting point Tm (in K), latent heat of fusion
L (in J/g), solid and liquid entropy at melting point, SS and SL (in J/mol K), and thermal expansion coefficient α (in 10−6K−1) at 300 K. The MEAM∗-Au and MEAM∗-Cu entries correspond to a modifica-
tion of the original MEAM model by changing cmin from 2.0 to 0.8. The MEAM† entries of BCC metals are computed by the new MEAM model that includes second nearest neighbor interactions [209, 210]. . 146
D.2 The estimated free energy difference ∆Fi and its standard deviation in the 5 different adiabatic switching steps for the melting point calcula- tionsofSW-SiandMEAM-Simodels...... 159
E.1 Model predictions and experimental data [207, 230, 231] on thermal and mechanical properties of gold. The computation models include original MEAM [206], EAM [204], 2nn-MEAM [227] and two modifi- cations made in this study (2nn-MEAM∗ and 2nn-MEAM†), and first- principles calculation with DFT/LDA [232]. The properties include the −1 melting point Tm (in K), latent heat of fusion L (in kJ mol ), solid and −1 −1 liquid entropy at the melting point, SS and SL (in J mol K ), diffu- −9 2 −1 sion constant of the liquid D at Tm (in 10 m s ), thermal expansion coefficient α of the solid (in 10−6K−1) at 300 K, ideal shear strength
τc (in GPa) of the solid at zero temperature, and volume change on
melting ∆Vm/Vsolid (in %). Statistical errors are on the order of the lastdigitofthepresentedvalues...... 163
xvi E.2 Model predictions and experimental data [207, 230, 196] of thermal and mechanical properties of silicon. The ideal shear strength [233] obtained from first principle calculation is also presented for compari- son. The computational models include the original MEAM [206] and MEAM second nearest-neighbor (2nn) [228]. Superscripts ∗ and † rep- resent modifications to these models introduced in this work. SW [105] and Tersoff [106] models are also included for comparison. The vari- ables and their units are identical to those in Table E.1...... 166
E.3 The elastic constants, C11, C12 and C44 (in GPa), vacancy formation en- 2 ergy Ev (in eV) and surface energies E100, E110, and E111 (in ergs/cm ) from experimental measurements [183, 206, 241, 242], first principle calculations [241, 243, 244, 245] (except silicon elastic constants com- puted here using DFT/LDA) and various MEAM models considered in this study. The unrelaxed energies are given in parenthesis. The DFT results are from several different pseudopotentials. Surface reconstruc- tion such as dimer structure is not considered in calculation of relaxed surfaceenergyofsilicon...... 178
F.1 Equilibrium lattice constant a, bulk modulus B, cohesive energy Ec,
and cubic elastic constants C11 and C44 for DC structure of Si, FCC structure of Au, and B1 structure of Au-Si. For pure Si and Au, the differences between the experimental and ab initio values are listed in the column labelled “offset”. Their average is the expected “offset” value for the hypothetical B1 structure. The values marked with ∗ are the ab initio values plus the correction terms given in the “offset” column. The last column is what the MEAM model is fitted to or predicts...... 183
F.2 MEAM and ab initio (DFT/LDA) predictions of impurity energies. E1 is the energy needed to substitute an atom in an FCC Au crystal by
a Si atom. E2 is the energy needed to substitute an atom in a DC Si crystalwithbyaAuatom...... 184
xvii F.3 Parameters for the Au-Si MEAM cross-potential using B1 as the ref-
erence structure. Cmin(i, j, k) are cut-off parameters in the multi-body screening function. They describe the screening effect on the interac- tion between atoms of type i and j by their common neighbor of type
k, where i, j, k = 1 (Au) or 2 (Si). The same Cmax(i, j, k) is used for every combination of i, j, k ...... 185 F.4 Comparison between MEAM and ab initio predictions on the energy
and elastic properties of B1, B2 and L12 structures of Au-Si. a is the
equilibrium lattice constant, B is bulk modulus, Ec is the cohesive en-
ergy, and C11 and C44 are cubic elastic constants. MEAM data should be compared with ab initio (DFT/LDA) results that have been ad- justed for known differences from experimental values for pure elements.195
xviii List of Figures
2.1 ∆G versus n showing volume and surface contributions resulting in a
peak at n = nc...... 9 2.2 Schematic of a solid spherical cap with radius of curvature r and con- tact angle θ formingawall...... 10 2.3 Comparison of experimental nucleation rates of water from vapor (cir- cles) with the predictions of the classical theory [72]. The full lines belong to the classical Becker-D¨oring (BD) theory. Logarithmic scales areusedforbothaxes...... 17 2.4 Schematic of dislocations [80]. (a) An edge dislocation is a defect where an extra half-plane of atoms is introduced mid way through the crystal, distorting nearby plane of atoms. (b) A screw dislocation is a shear rippleextendingfromsidetoside...... 20 2.5 Schematic of edge dislocation motion that induces plastic strain [80]. 21 2.6 (a) and (b) Two consecutive in situ TEM compression tests on a FIB microfacbricated 160-nm-top-diameter Ni pillar with 111 orien- tation [89]. (a) Dark-field TEM image of the pillar before the tests; note the high initial dislocation density. (b) Dark-field TEM image of the sam pillar after the first test; the pillar is now free of dislocations. (c) MD simulation of nanoindentation process [95]. Snapshot of dis- location nucleation at the first plastic yield point on Au(111), Lower figure shows a two-layer-thick cross section of a (111) plane containing the partial dislocation loop on the right in the upper figure...... 22
xix 1 3.1 (a) Magnetization M = N si of the 2D square lattice Ising model as function of temperature in the absence of external field, i.e. h = 0. At T < T = 2/ ln(1 + √2)J/k 2.269J/k , majority of spins in c B ∼ B the system are spontaneously aligned in either up or down. However,
T > Tc, spontaneous ordering disappears because of thermal fluctua-
tion. (b) schematic of ferromagnetic state below Tc. (c) schematic of
paramagnetic state above Tc...... 27 3.2 (a) Analogy between the Ising model demagnetization and nucleation phenomena. (b) The snapshot of Monte Carlo simulation of demagne- tization in the presence of an external field h> 0. The transition from the down spin dominated meta-stable state to the up spin dominated stable state occurs via nucleation of small island of up spins. . . . . 28 3.3 (a) The stacking pattern of the 111 planes in FCC crystals. (b) Perfect
Burgers vectors b1, b2, b3 and partial Burgers vectors bp1, bp2, bp3 on the 111 plane. Figures are taken from the literature [113]...... 35 3.4 Generalized stacking fault energy on (111) plane along a/6[211] direc- tion calculated with EAM potential [111]...... 36 3.5 Linear thermal expansion of Cu calculated in the quasi-harmonic ap- proximation (QHA) and by the Monte Carlo method using the EAM
potential [111]. The melting point of Cu (Tm) is indicated...... 37
xx 3.6 Schematic of forward flux sampling. (a) To compute the nucleation
rate of small droplet with size n0, we count the number q that droplets
larger than n0 forms for time duration t. Then, nucleation rate I0 becomes q/t. Here, it is not desirable to count small fluctuation in n
around n0 as separate events. The counter reset only after n comes
back to the original basin n < nA and becomes ready to count another
event. nA can be considered as an error margin in the digital signal processing. (b) At each interface i, we have an ensemble of configu-
rations having largest droplet size ni. N independent MC simulations are performed, starting from a randomly chosen configuration from the ensemble. Then, number M of reaching next interface i +1 before re-
turning back to nA is counted as successful forward flux. Then, the probability of reaching next interface P (n n ) becomes N/M. ... 42 i+1| i 3.7 (color online) The probability P (λ λ ) (solid line) of reaching interface i| 0 λi from λ0 and average committor probability PB(λi) (circles) over
interface λi at (kBT, h)=(1.5, 0.05) for the 2D Ising model. The 50% committorpointismarkedby*...... 44 3.8 (a) For a given free energy landscape, we can sample a limited region
within the kBT range. (b) We can sample other regions with higher free energy with umbrella sampling that limits the Monte Carlo move withintherangeofbiasfunction...... 45 3.9 (From the top left cornet, in clock wide direction). Processing the raw histogram taken from umbrella sampling simulations. (1) We first ob- ′′ tain the biased relative distribution P (n) at each window from series of ′ umbrella sampling. (2) Unbiased relative distribution P (n) can be ob- 1 2 tained by multiplying the inverse of weighing function, exp[ 2 k(n−ni) ], kBT to the data of each window i. (3) From the overlapping histogram methods [125], we can merge the distributions into a single curve. Af- ter normalization, we obtain the absolute probability P (n). (4) The free energy curve ∆G(n) can be obtained from k T ln(∆G(n)). .. 46 − B
xxi 3.10 (color online) (a) Droplet free energy F (n) obtained by US at kBT = 1.5 and h =0.05 in the 2D Ising model. (b) Fluctuation of droplet size ∆n2(t) asafunctionoftime...... 48
4.1 (color online) Effective surface free energy σeff as a function of temper- ature for the 2D Ising model from analytic expression [47]. The free
energy of the surface parallel to the sides of the squares, σ(10), is also plottedforcomparison...... 54 4.2 Schematic of numerical test on the CNT nucleation rate...... 58 4.3 The nucleation rate I computed by FFS (open symbols) and Becker- D¨oring theory with US free energies (filled symbols) in the (a) 2D and (c) 3D Ising models. The ratio between nucleation rates obtained by FFS and Becker-D¨oring theory at different temperatures in the (b) 2D and (d) 3D Ising models. The symbols in (b) and (d) match those definedin(a)and(c),respectively...... 60 4.4 (a) For the 2D Ising model, the critical droplet size n obtained from
FFS (filled symbols) and umbrella sampling (open symbols). nc pre- dicted by Becker-D¨oring theory (dotted line) and by field theoretic equation (solid line) are plotted for comparison. (b) For the 3D Ising model, the critical droplet size n obtained from FFS (filled symbols) andumbrellasampling(opensymbols)...... 61 4.5 (a) Histogram of committor probability in an ensemble of spin configu-
rations with n = 496 for the 2D Ising model at kBT =1.5 and h =0.05. Representative droplets are also shown, with black and white squares corresponding to +1 and 1 spins, respectively. (b) Histogram of com- − mittor probability in an ensemble of spin configurations with n = 524
for the 3D Ising model at kBT =2.20 and h =0.40...... 62
xxii 4.6 (a) Droplet free energy curve F (n) of the 2D Ising model at kBT =1.5 and h = 0.05 obtained by US (circles) is compared with Eq. (4.10) (solid line) and Eq. (4.3) (dashed line). Logarithmic correction term 5 4 kBT ln n (dot-dashed line) and the constant term d (dotted line) are also drawn for comparison. (b) Magnified view of (a) near n = 0, together with the results from analytic expressions (squares) available for n 17(seeAppendixB.3)...... 64 ≤ 4.7 (a) Droplet free energy F (n) of the 3D Ising model at kBT = 2.40 and h = 0 obtained by US (circles) is compared with Eq. (4.22) (solid
line) and Eq. (4.11) (dots). Logarithmic term τkBT ln n is also plotted (dot-dashed line). The difference in predictions by classical expres- sion Eq. (4.11) and field theory Eq. (4.22) are very small compared to F (n) itself and cannot be observed at this scale. (b) Magnified view of (a) near n = 0, together with the analytic solution of small droplets (squares, see Appendix B.3) and the exponential correction term(dashedline)...... 66 4.8 (a) Surface free energies of the 3D Ising model as functions of tempera-
ture. Circles are fitted values of σeff from Eq. (4.22), dashed line is the
expected behavior of σeff over a wider range of temperature, and solid line is the free energy of the (100) surface [141]. Numerically fitted
values of σeff from Heermann et al. [137] are plotted as +. (b) τ values that give the best fit to the free energy data from US. τ can be roughly described by a linear function of T shown as a straight line. No abrupt
change is observed near the roughening temperature TR...... 69
xxiii 4.9 (a) Droplet free energy as a function of droplet size n at h = 0.1 and
different kBT for the 2D Ising model. The critical droplet free energy is marked by circles. (b) Critical droplet free energy (circles) from (a)
as a function of kBT for the 2D Ising model. The solid line is a linear fit of the data, and the dashed line is the prediction of Eq. (4.5). (c) Droplet free energy as a function of droplet size n at h = 0.45 and
different kBT for the 3D Ising model. (d) Critical droplet free energy
from (c) as a function of kBT for the 3D Ising model. The solid line is a linear fit of the data, and the dashed line is the prediction of Eq. (4.13). 71
5.1 Homogeneous dislocation nucleation rate per lattice site in Cu under pure shear stress σ = 2.0 GPa on the (111) plane along the [112] direction as a function of T −1, predicted by Becker-D¨oring theory using free energy barrier computed from umbrella sampling (See Section 5.3). The solid line is a fit to the predicted data (in circles). The slope of the
line is Hc/kB, while the intersection point of the extrapolated line with
the vertical axis is ν0 exp(Sc/kB). Dashed line presents the nucleation rate predicted by ν exp( H /k T ), in which the activation entropy is 0 − c B completely ignored, leading to an underestimate of the nucleation rate by 20ordersofmagnitude...... 81 ∼ 5.2 Schematics of simulation cells designed for studying (a) homogeneous and (c) heterogeneous nucleation. In (a), the spheres represent atoms enclosed by the critical nucleus of a Shockley partial dislocation loop. In (c), atoms on the surface are colored by gray and atoms enclosed by the dislocation loop are colored by magenta. Shear stress-strain curves of the Cu perfect crystal (before dislocation nucleation) at different temperatures for (b) homogeneous and (d) heterogeneous nucleation simulation cells...... 90
xxiv 5.3 (a) The Helmholtz free energy of the dislocation loop as a function of its
size n during homogeneous nucleation at T = 300 K, σxy = 2.16 GPa
(γxy = 0.135) obtained from umbrella sampling. (b) Size fluctuation of critical nuclei from MD simulations...... 92 5.4 Atomistic configurations of dislocation loops at (a) 0 K and (b) 300 K. 93 5.5 The fraction of 192 MD simulations in which dislocation nucleation
has not occurred at time t, Ps(t), at T = 300 K and σxy = 2.16 GPa (γ =0.135). Dotted curve presents the fitted curve exp( IMDt) with xy − IMD =2.5 108s−1...... 95 × 5.6 Activation Helmholtz free energy for homogeneous dislocation nucle-
ation in Cu. (a) Fc as a function of shear strain γ at different T . (b)
Gc as a function of shear strain σ at different T . Squares represent um- brella sampling data and dots represent zero temperature MEP search results using simulation cells equilibrated at different temperatures.
(c) Fc as a function of T at γ = 0.092. (d) Gc as a function of T at σ = 2.0 GPa. Circles represent umbrella sampling data and dashed linesrepresentapolynomialfit...... 97 5.7 Activation free energy for heterogeneous dislocation nucleation from
the surface of a Cu nanorod. (a) Fc as a function of compressive
strain ǫzz at different T . (b) Gc as a function of compressive stress
σzz at different T . Squares represent umbrella sampling data and dots represent zero temperature MEP search results using simulation cells equilibratedatdifferenttemperatures...... 99
5.8 The relation between Ec and Sc in the temperature range of zero to 300 K for (a) homogeneous and (b) heterogeneous nucleation. The
relation between Hc and Sc for (c) homogeneous and (d) heterogeneous nucleation. The solid lines represent simulation data and the dashed ∗ ∗ lines are empirical fits of the form Sc = Ec/T or Sc = Hc/T . .... 102
xxv 5.9 Contour lines of (a) homogeneous and (b) heterogeneous dislocation nucleation rate per site I as a function of T and σ. The predictions with
and without accounting for the activation entropy Sc(σ) are plotted in thick and thin lines, respectively. The nucleation rate of I 106 s−1 ∼ per site is accessible in typical MD timescales whereas the nucleation rate of I 10−4 10−9 is accessible in typical experimental timescales, ∼ − depending on the number of nucleation sites...... 103 5.10 (a) Nucleation stress of our bulk sample (containing 14,976 atoms) un- der constant shear strain loading rateγ ˙ = 10−3 and (b) nucleation stress of the nanorod under constant compressive strain loading rate ǫ˙ = 10−3. The strain rate 10−3 is experimentally accessible loading rate. The solid lines are the prediction based on the activation free energy computed by umbrella sampling. The dashed lines are the nu- cleation stress prediction when the activation entropy is neglected. The dotted line in (b) is the prediction based on the approximation by Zhu etal.[153]...... 104
+ B.1 (color online) (a) The pre-exponential factor fc Γ in 2D computed from + Monte Carlo and US. (b) The ratio between the attachment rate fc in 2D computed by Monte Carlo and that predicted by Eq.(4.7). (a) + The pre-exponential factor fc Γ in 3D computed from Monte Carlo and + US. (b) The ratio between the attachment rate fc in 3D computed by MonteCarloandthatpredictedbyEq.(4.15)...... 122 B.2 Droplets in (a) 2D and (b) 3D Ising models randomly chosen from FFS simulations at different (T, h) conditions. n is the size of the droplet. 124
B.3 The free energy curve F (n) of 2D Ising system at kBT = 1.0 and (a) h = 0.05, (b) h = 0.06, (c) h = 0.07, (d) h = 0.08, (e) h = 0.09, (f) h = 0.10 obtained by US (circles) is compared with Eq. (6) (solid 5 line) and Eq. (8) (dashed line). Logarithmic correction term 4 kBT ln n (dot-dashed line) and the constant term d (dotted line) are also drawn forcomparison...... 127
xxvi B.4 The free energy curve F (n) of 2D Ising system at kBT = 1.5 and (a) h = 0.04, (b) h = 0.05, (c) h = 0.06, (d) h = 0.07, (e) h = 0.08, (f) h =0.09, (g) h =0.10, (h) h =0.11, (i) h =0.12, (j) h =0.13 obtained by US (circles) is compared with Eq. (6) (solid line) and Eq. (8) (dashed 5 line). Logarithmic correction term 4 kBT ln n (dot-dashed line) and the constant term d (dotted line) are also drawn for comparison. . . . . 129
B.5 The free energy curve F (n) of 2D Ising system at kBT = 1.9 and (a) h = 0.015, (b) h = 0.017, (c) h = 0.022, (d) h = 0.025, (e) h = 0.03, (f) h =0.035 obtained by US (circles) is compared with Eq. (6) (solid 5 line) and Eq. (8) (dashed line). Logarithmic correction term 4 kBT ln n (dot-dashed line) and the constant term d (dotted line) are also drawn forcomparison...... 130
C.1 The relation between critical dislocation size nc and the activation volume Ω ∂Gc .(a) homogeneous nucleation (b) heterogeneous nu- c ≡ − ∂σ cleation. Circles represent the activation volume obtained from the
derivative of Gc with respect to σ. Squares represent the activation volume data multiplied by 1/S where S is the Schmid factor. Dashed linesarelinearfitstothedata...... 140
D.1 Gibb’s free energy per atom for both the solid phase (solid line) and liquid phase (dashed line). The symbols represent data points in Broughton and Li [198] with squares for the solid phase and circles fortheliquidphase...... 157
E.1 (Color Online) Generalized stacking fault energy of different poten- tial models for gold: (a) 2nn-MEAM [227] (dashed line), 2nn-MEAM∗ (dotted line), and 2nn-MEAM† (solid line); (b) EAM [204] (dashed line), MEAM [206] (dotted line) and DFT/LDA (solid line)...... 165
xxvii E.2 (Color Online) Generalized stacking fault energy of different potential models for silicon. (a) MEAM [206](dashed curve), MEAM∗ (dotted curve), MEAM† (solid curve). (b) 2nn-MEAM [228](dashed curve). (c) Tersoff [106](solid curve), SW [105] (dashed curve) (d) DFT/LDA [233] (solidcurve) ,DFT/GGA[233](dashedcurve) ...... 167 E.3 (Color Online) (a) Pair-correlation functions of the solid (solid line) and liquid (dotted line) phases of gold described by the 2nn-MEAM [227] potential at its melting point. (b) The equation of state function in the 2nn-MEAM (dotted line) potential and the new 2nn-MEAM† potential (solid line). (c) The Gibbs free energy of the 2nn-MEAM (thick lines) and 2nn-MEAM† (thin lines) potentials for gold. Solid lines for the solid phase and dashed lines for the liquid phase...... 171 E.4 (Color Online) (a) Pair-correlation functions of the solid (solid line) and liquid (dotted line) phases of silicon described by the MEAM [206] potential at its melting point. (b) The equation of state function in the original MEAM (dotted line) potential and the new MEAM† potential (solid line). (c) The Gibbs free energy of the MEAM (thick lines) and MEAM† (thin lines) potentials for silicon. Solid lines for the solid phase and dashed lines for the liquid phase...... 173 E.5 (Color Online) Ellipses defined in Eq. (E.7) for different values of C (0.8,2.0,2.8). The line segments represent nearest-neighbor bonds i-j and j-k in BCC (square), FCC (circle), and Diamond-Cubic (aster- isk) crystal structures, scaled by the second nearest-neighbor distance
rik. In these three crystal structures, the atom j lies on the ellipses (not shown) corresponding to C =0.5, 1.0 and 2.0, respectively. . . . 176 E.6 (Color Online) Bond angle distribution functions of liquid phase of sili- con described by the MEAM† (solid line) ,DFT/LDA (dashed line) [246] andSW(dottedline)[246]...... 179
xxviii F.1 Binary phase diagram of Au-Si. MEAM prediction is plotted in thick line and experimental phase diagram is plotted in thin line. L corre- sponds to the liquid phase. Au(s) and Si(s) correspond to the Au-rich and Si-rich solid phases, respectively...... 186
F.2 Gibbs free energy ∆gimp(T ) and enthalpy ∆himp(T ) (a) a Si impurity within Au crystal. (b) a Au impurity within Si crystal...... 189
F.3 (a) Liquid free energy Gliq(x, T ) at T = 1250 K. Circles are simula- tion results, which are fitted to a spline (solid line). A straight line connecting the liquid free energy of pure Au and pure Si is drawn for
comparison. (b) The free energy of mixing Gmix(x, T ) for the liquid phase at T = 1250 K. Predictions from the MEAM potential is plotted
in thick line, which is the difference between Gliq(x, T ) and the straight line shown in (a). Free energy obtained from CALPHAD method [256] are plotted in thin line...... 190
F.4 (a) Enthalpy of mixing ∆Hliq(x, T ) at T = 1373 K from experiments (circles), MEAM (thick line), and CALPHAD (thin line). (b) Excess XS free energy of mixing ∆Gliq (x, T ) at T = 1685 K from experiments (circles), MEAM (thick line), and CALPHAD (thin line)...... 192 F.5 Common tangent method to construct binary phase diagram from free
energy curves. (a) Gibbs free energy of the three phases, GFCC(x, T ),
GDC(x, T ) and Gliq(x, T ), as a function of composition x at T = 700 K. All of them are referenced to free energies of pure Au liquid and pure
Si liquid. Common tangent lines are drawn between GFCC(x, T ) and
Gliq(x, T ) (from x1 =0.011 to x2 =0.225), and between Gliq(x, T ) and G (x, T ) (from x = 0.251 to x 1). (b) Binary phase diagram of DC 3 4 ≈ the MEAM Au-Si potential. The phase boundaries at T = 700 K are determinedfromthedatain(a)...... 193
xxix Chapter 1
Introduction
1.1 Computational Investigation of Nucleation
Nucleation refers to the formation of small region of a new phase in the background of a meta-stable phase during a first order phase transition such as melting and freezing. Such phase transitions are ubiquitous in the nature and hence have been investigated in a wide range of scientific disciplines, including physics, materials sci- ence, meteorology, medical science, biology, and nano technology. Phase transitions that proceed via nucleation include supercooled fluids [1, 2, 3, 4], cloud formation [5], kidney stone formation [6], polymerization [7], electro-weak phase transitions [8], and nano materials [9]. To understand and have control over the nucleation processes, it is important to predict the nucleation rate I as a function of ∆ and T . Here, ∆ refers to the chemical potential difference between the meta-stable phase and the new phase, and T is the absolute temperature. It is found that the nucleation rate is exponentially sensitive to ∆ and T , which makes the prediction of nucleation rate very difficult. While ∆ depends solely on T in most single component systems (ignoring stress effects), it depends on both T and composition in multi-component systems, making the understanding and prediction of the nucleation process even more challenging. The nucleation phenomena has been investigated by three different approaches:
1 CHAPTER 1. INTRODUCTION 2
experiment, theory, and computer simulation [10, 11]. If appropriate tools are avail- able, direct experiments on the system of interest would be an ideal way to understand the nucleation. However, it is likely that nucleation experiments are affected by the existence of impurities and surfaces that are hard to eliminate. Even if impurities are controlled and nucleation rate can be measured with a reasonable accuracy, it is still difficult to reveal the microscopic detail of the nucleation process from experiments. Deeper understanding of the microscopic process is crucial to apply the insight gained from an experiment on a simple system to predict phase transitions of more complex systems where experimental investigation is very difficult. Theoretical study can complement experiments by revealing the microscopic detail based on known physics, and provides rate predictions that can be compared with the experimental results. In many cases, the qualitative trend obtained from the theory can be well compared with experiments, but the quantitative matching is difficult to achieve [10, 11, 12, 13]. Because most theories rely on multiple assumptions that are difficult to be validated, it is hard to pinpoint which part of the theory causes the discrepancy. Besides, the theoretical expression on the nucleation rate requires several pre-determined paramters as input, such as surface energy σ, characteristic vibration frequency ν, and the chemical potential difference ∆ at given temperature, composition, external stress and so on. While the rate prediction is sensitive to these parameters, it is difficult to accurately measure them in many cases. As an alternative approach, computer simulations have also played an important role in probing nucleation processes [14, 15]. The advantage of simulations lies in creating a model system which is difficult to prepare in a laboratory. For example, we can model the solidification of liquid that has no impurity, which is an ideal testbed of homogeneous nucleation theories. We can trace positions and velocities of all atoms in the model system, which allows a very accurate prediction of the nucleation rate as well as the detailed microscopic processes. Of course, the interactions among particles cannot be perfectly modeled and it is impossible to perform a virtual experiment reproducing what happens in the nature exactly. Still, if we have model systems that mimick real systems reasonably well, i.e. when the computational models capture important characteristics of chemical bonding and predict phase diagrams close to CHAPTER 1. INTRODUCTION 3
experiments, computer simulations often leads to new discoveries that have not been predicted from existing theories. For instance, it has been found that multiple order parameters other than the size of nucleus affect nucleation rates of crystallization in many circumstances [16, 17]. It was a surprising discovery that the critical nuclei in a binary suspension of oppositely charged colloid is not the one with the lowest free energy barrier for nucleation, but the one with the fastest growth rate [18]. While computer simulations have been proved as a powerful tool to study the nucleation phenomena, it is important to devise smart algorithms that can overcome the limitations arising from the limited computing power, length scale limitation and time scale limitation. Length scale limitation refers to the limited number of atoms that can be modeled by computer simulations. For example, while a drop (about 0.3 cc) of water consists of 1022 molecules, it takes a few hours of modern CPU (central ∼ processor unit) time to simulate the dynamics of 10, 000 water molecules for a few hundreds picoseconds even though computationally cheap empirical potential is used to model water molecules. Fortunately, the length scale problem is relatively easy to solve and is not a critical obstacle. Because atoms involved in nucleation event can be as small as a few hundred atoms in many substances, computer simulations of a few tens of thousands particles are often big enough to capture the formation of the critical nuclei and to describe the nucleation pathway. Even larger critical nucleus can be handled by using multiple processors simultaneously. For example, a molecular dynamics simulation has recently reached the 1011 atoms simulating a 1.5 m cubed box, employing state-of-the-art parallel computing algorithms [19]. The major limiting factor in the simulation of nucleation process is the time scale problem. The time step of molecular dynamics simulations must be on the order of femtoseconds to stablize numerical integrators used to trace the motion of particles, because the characteristic vibration frequency is around 1013s−1 in most condensed matter systems. It takes a few days to proceed a few million time steps of simulations which correspond to a few nano-seconds of simulations time. However, typical time scale of nucleation events is a few millisecond to a few seconds which is many orders of magnitude larger than the time scale of conventional molecular dynamics simulations. It is recognized that a system spends most of time fluctuating around a meta-stable CHAPTER 1. INTRODUCTION 4
phase, while a successful nucleation event is extremely rare. When it occurs, the formation of a stable nucleus can happen within picoseconds. Because it is such a rare event that controls the onset of phase transitions, many versions of advanced sampling methods have been developed that captures such rare events selectively [20], which allows us to estimate the nucleation rate. This dissertation is devoted mainly to the study of nucleation processes via com- puter simulations. Using the Ising model [21], the simplest and well-investigated model of phase transitions (as well as of ferromagnetisms), we systematically test the validity of the classical nucleation theory (CNT) [22, 23] which is a standard theory that has been used to describe the nucleation phenomena for almost a century [24]. The validated part of the classical nucleation theory, in combination with computer simulations, has been applied to predict the rate of dislocation nucleation [25] which is essential to our understanding of plastic deformation, ductility, and mechanical strength of crystalline materials. We have employed advanced sampling methods to overcome the time scale problems when studying both the Ising model and dislocation nucleation.
1.2 Scope of the Dissertation
The dissertation is organized as follows. Chapter 2 introduces nucleation theories and experiments relevant to this work. We begin with a short description of the basic thermodynamics of nucleation such as the chemical potential difference, nucleation barrier, and population of droplets of new phase. Nucleation rate predictions from three different nucleation theories will be presented with special focus on the classical nucleation theory (CNT) and its two fundamental assumptions. With brief historic review of nucleation experi- ments, discrepancies between the CNT prediction and experimental results will be highlighted. In the last part of the chapter, we briefly introduce the concepts of dislo- cation and explain why dislocation nucleation plays an important role in determining the mechanical behavior of materials at small scale. CHAPTER 1. INTRODUCTION 5
Chapter 3 summarizes the computational methods used in this work. We in- troduce the Ising model whose demagnetization process has a close analogy to the nucleation dynamics. We review interatomic potentials that were developed to de- scribe different bonding mechanisms, with a special focus on the Cu embedded-atom- method (EAM) potential which is used in dislocation nucleation study. An overview of molecular dynamics (MD) and Monte Carlo (MC) methods is presented with pros and cons of each method. We also describe two advanced sampling methods that can overcome the timescale limit of conventional MD and MC simulations. In chapter 4, we test the validity of the classical nucleation theory (CNT) by calculating the individual components of CNT via computer simulations of the Ising models. We open this chapter with a brief description on how nucleation theories are applied to the Ising model. Using two independent simulation techniques, we confirm the fundamental assumption that nucleation process can be described by 1D Markov chain, under a wide range of conditions in both 2D and 3D Ising models. However, it is found that the free energy predicted by CNT does not match with numerical results, unless appropriate correction term is added. Our analysis confirms that the nucleation rate by CNT can be predicted accurately if a correct free energy barrier obtained by umbrella sampling is used as an input. Chapter 5 provides an in-depth description on the prediction of dislocation nucle- ation rate based on the classical nucleation theory in combination with the umbrella sampling technique. The results reveal very large activation entropies, originated from the anharmonic effects, which can alter the nucleation rate by many orders of magnitude. Here we discuss the thermodynamics and algorithms underlying these calculations in great detail. In particular, we prove that the activation Helmholtz free energy equals the activation Gibbs free energy in the thermodynamic limit, and explain the large difference in the activation entropies in the constant stress and con- stant strain ensembles. We also discuss the origin of the large activation entropies for dislocation nucleation, along with previous theoretical estimates of the activation entropy. Finally, chapter 6 reviews the results presented in the dissertation and discuss future research opportunities. One possibility is the investigation of the gold catalyzed CHAPTER 1. INTRODUCTION 6
growth of silicon nanowire via vapor-liquid-solid (VLS) mechanism which involves silicon crystal nucleation inside gold-silicon eutectic liquid alloy. As a preparation for the project, we have developed a Au-Si potential that is fitted to the binary phase diagram and efficient free energy calculation method for solid and liquid alloy. These contributions are presented in Appendices D, E, and F. Chapter 2
Background and Motivation
This chapter reviews theoretical and experimental studies of nucleation phenomena. We begin with a brief explanation on the thermodynamic origin of the nucleation barrier in the first order transitions. We present nucleation rate predictions from three different nucleation theories, which is followed by a concise review of experiments for testing the nucleation theories. The relation between dislocation nucleation and materials strength at small scale will be discussed in the last section of the chapter.
2.1 Basic Thermodynamics of Nucleation
Most first order phase transitions require the appearance of small nuclei of the new phase as a prerequisite. Nucleation refers to such localized budding of new phases in the background of the ambient phases, i.e. meta-stable phases [10, 11]. Some examples of the emergent phases include gaseous bubbles, small crystallites, and liquid droplets in the volume of supersaturated solution of gas, undercooled liquid, and undercooled vapor, respectively. It is statistical fluctuations that create nuclei that undergo the transient appearance and disappearance. Only when a “critical” size is exceeded, the dissolution probability of nuclei becomes small enough and the new phase evolves into a macroscopic size. The work of formation of the critical nucleus, so called “nucleation barrier”, is supplied by thermal fluctuation. To quantitatively describe the process in terms of physics, we consider a volume
7 CHAPTER 2. BACKGROUND AND MOTIVATION 8
containing a original phase with chemical potential 1 (i.e. the Gibbs free energy per particle) which is a function of temperature T . We will consider only a single component system and ignore stress effects for simplicity. Formation of a droplet
of new phase with chemical potential 2 costs the surface free energy S σ where S is the surface area of the droplet and σ is the interface free energy between two phases. Hence, the change of the Gibbs free energy upon the formation of the droplet containing n particle is
∆G(n)= n( )+ S(n)σ. (2.1) − 1 − 2
The chemical potential difference ∆ = is the thermodynamics driving force 1 − 2 that induces the phase transition and becomes positive at conditions where the new
phase is thermodynamically favored, i.e. 1 > 2. For example, the chemical potential