Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Effects of environmental variation on ecological and evolutionary dynamics

Environmental variation is a fact of life for all the species on earth: for any population of any particular species, the local environmental conditions are liable to vary in both time and space. In today's world, anthropogenic activity is causing habitat loss and fragmentation for many species, which may profoundly alter the characteristics of environmental variation in remaining habitat. Previous research indicates that, as habitat is lost, the spatial configuration of remaining habitat will increasingly affect the dynamics by which populations are governed. Through the use of mathematical models, this thesis asks how environmental variation interacts with species properties to influence population dynamics, local adaptation, and dispersal evolution. More specifically, we couple continuous-time continuous-space stochastic population dynamic models to landscape models. We manipulate environmental variation via parameters such as mean patch size, patch density, and patch longevity. Among other findings, we show that a mixture of high and low quality habitat is commonly better for a population than uniformly mediocre habitat. This conclusion is justified by purely ecological arguments, yet the positive effects of landscape heterogeneity may be enhanced further by local adaptation, and by the evolution of short-ranged dispersal. The predicted evolutionary responses to environmental variation are complex, however, since they involve numerous conflicting factors. We discuss why the species that have high levels of local adaptation within their ranges may not be the same species that benefit from local adaptation during range expansion. We show how habitat loss can lead to either increased or decreased selection for dispersal depending on the type of habitat and the manner in which it is lost. To study the models, we develop a recent analytical method, Perturbation expansion, to enable the incorporation of environmental variation. Within this context, we use two methods to address evolutionary dynamics: Adaptive dynamics, which assumes mutations occur infrequently so that the ecological and evolutionary timescales can be separated, and via Genotype distributions, which assume mutations are more frequent. The two approaches generally lead to similar predictions yet, exceptionally, we show how the evolutionary response of dispersal behaviour to habitat turnover may qualitatively depend on the mutation rate.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Social dynamics for sustainable food systems : Actors' orientations towards sustainability in primary production and public consumption

The modern food system and sustainable development form a conceptual combination that suggests sustainability deficits in environmental impacts and nutritional status of western populations. This study explores actors orientations towards sustainability by probing into social dynamics for sustainability within primary production and public consumption. If actors within these two worlds were to express converging orientations for sustainability, the system dynamics of the market would enable more sustainable growth in terms of production dictated by consumption. The study is based on a constructivist research approach with qualitative text analyses. The findings were validated by internal and external food system actors and are suggested to represent current social dynamics within Finnish food system. The key findings included primary producers social skilfulness, which enabled networking with other actors in very different paths of life, learning in order to promote one s trade, and trusting reflectively in partners in order to expand business. These activities extended the supply chain in a spiral fashion by horizontal and vertical forward integration, until large retailers were met for negotiations on a more equal basis. This mode of chain level coordination, typically building around the core of social and partnership relations, was coined as a socially overlaid network, and seen as sustainable coordination mode for endogenous growth. The caterers exhibited more or less committed professional identity for sustainability within their reach. The facilitating approaches for professional identities dealt successfully with local and organic food in addition to domestic food, and also imported food. The co-operation with supply chains created innovative solutions and savings for the business parties to be shared. There were also more complicated identities as juggling, critical and delimited approaches for sustainability, with less productive efforts due to restrictions such as absence of organisational sustainability strategy, weak presence of local and organic suppliers, limited understanding about sustainability and no organisational resources for informed choices for sustainability. The convergence between producers and caterers existed to an extent allowing suggestion that increased clarity about sustainable consumption and production by actors could be constructed using advanced tools. The study looks for introduction of more profound environmental and socio-economic knowledge through participatory research with supply chain actors. Learning in the workplace about food system reality in terms of supply chain co-operation may prove to be a change engine that leads to advanced network operations and a more sustainable food system.

Social dynamics for sustainable food systems. Actors’ orientations towards sustainability in primary production and public consumption

Abstract The modern food system and sustainable development form a conceptual combination that suggests sustainability deficits in the ways we deal with food consumption and production - in terms of economic relations, environmental impacts and nutritional status of western population. This study explores actors’ orientations towards sustainability by taking into account actors’ embedded positions within structures of the food system, actors’ economic relations and views about sustainability as well as their possibilities for progressive activities. The study looks particularly at social dynamics for sustainability within primary production and public consumption. If actors within these two worlds were to express converging orientations for sustainability, the system dynamics of the market would enable more sustainable growth in terms of production dictated by consumption. The study is based on a constructivist research approach with qualitative text analyses. The data consisted of three text corpora, the ‘local food corpus’, the ‘catering corpus’ and the ‘mixed corpus’. The local food actors were interviewed about their economic exchange relations. The caterers’ interviews dealt with their professional identity for sustainability. Finally, the mixed corpus assembled a dialogue as a participatory research approach, which was applied in order to enable researcher and caterer learning about the use of organic milk in public catering. The data were analysed for theoretically conceptualised relations, expressing behavioural patterns in actors’ everyday work as interpreted by the researcher. The findings were corroborated by the internal and external communities of food system actors. The interpretations have some validity, although they only present abstractions of everyday life and its rich, even opaque, fabric of meanings and aims. The key findings included primary producers’ social skilfulness, which enabled networking with other actors in very different paths of life, learning in order to promote one’s trade, and trusting reflectively in partners in order to extend business. These activities expanded the supply chain in a spiral fashion by horizontal and vertical forward integration, until large retailers were met for negotiations on a more equal or ‘other regarding’ basis. This kind of chain level coordination, typically building around the core of social and partnership relations, was coined as a socially overlaid network. It supported market access of local farmers, rooted in their farms, who were able to draw on local capital and labour in promotion of competitive business; the growth was endogenous. These kinds of chains – one conventional and one organic – were different from the strategic chain, which was more profit based and while highly competitive, presented exogenous growth as it depended on imported capital and local employees. However, the strategic chain offered learning opportunities and support for the local economy. The caterers exhibited more or less committed professional identity for sustainability within their reach. The facilitating and balanced approaches for professional identities dealt successfully with local and organic food in addition to domestic food, and also imported food. The co-operation with supply chains created innovative solutions and savings for the business parties to be shared. The rule-abiding approach for sustainability only made choices among organic supply chains without extending into co-operation with actors. There were also more complicated and troubled identities as juggling, critical and delimited approaches for sustainability, with less productive efforts due to restrictions such as absence of organisational sustainability strategy, weak presence of local and organic suppliers, limited understanding about sustainability and no organisational resources to develop changes towards a sustainable food system. Learning in the workplace about food system reality in terms of supply chain co-operation may prove to be a change engine that leads to advanced network operations and a more sustainable food system. The convergence between primary producers and caterers existed to an extent allowing suggestion that increased clarity about sustainable consumption and production by actors could be approached using advanced tools. The study looks for introduction of more profound environmental and socio-economic knowledge through participatory research with supply chain actors in order to promote more sustainable food systems. Summary of original publications and the authors’ contribution I Mikkola, M. & Seppänen, L. 2006. Farmers’ new participation in food chains: making horizontal and vertical progress by networking. In: Langeveld, H. & Röling N. (Eds.). Changing European farming systems for a better future. New visions for rural areas. Wageningen, The Netherlands. Wageningen Academic Publishers: 267–271. II Mikkola, M. 2008. Coordinative structures and development of food supply chains. British Food Journal 110 (2): 189–205. III Mikkola, M. 2009. Shaping professional identity for sustainability. Evidence in Finnish public catering. Appetite 53 (1): 56–65. IV Mikkola, M. 2009. Catering for sustainability: building a dialogue on organic milk. Agronomy Research 7 (Special issue 2): 668–676. Minna Mikkola has been responsible for developing the generic research frame, particular research questions, the planning and collection of the data, their qualitative analysis and writing the articles I, II, III and IV. Dr Laura Seppänen has contributed to the development of the generic research frame and article I by introducing the author to the basic concepts of economic sociology and by supporting the writing of article II with her critical comments. Articles are printed with permission from the publishers.

Modelling vacuum arcs : from plasma initiation to surface interactions

A better understanding of vacuum arcs is desirable in many of today's 'big science' projects including linear colliders, fusion devices, and satellite systems. For the Compact Linear Collider (CLIC) design, radio-frequency (RF) breakdowns occurring in accelerating cavities influence efficiency optimisation and cost reduction issues. Studying vacuum arcs both theoretically as well as experimentally under well-defined and reproducible direct-current (DC) conditions is the first step towards exploring RF breakdowns. In this thesis, we have studied Cu DC vacuum arcs with a combination of experiments, a particle-in-cell (PIC) model of the arc plasma, and molecular dynamics (MD) simulations of the subsequent surface damaging mechanism. We have also developed the 2D Arc-PIC code and the physics model incorporated in it, especially for the purpose of modelling the plasma initiation in vacuum arcs. Assuming the presence of a field emitter at the cathode initially, we have identified the conditions for plasma formation and have studied the transitions from field emission stage to a fully developed arc. The 'footing' of the plasma is the cathode spot that supplies the arc continuously with particles; the high-density core of the plasma is located above this cathode spot. Our results have shown that once an arc plasma is initiated, and as long as energy is available, the arc is self-maintaining due to the plasma sheath that ensures enhanced field emission and sputtering. The plasma model can already give an estimate on how the time-to-breakdown changes with the neutral evaporation rate, which is yet to be determined by atomistic simulations. Due to the non-linearity of the problem, we have also performed a code-to-code comparison. The reproducibility of plasma behaviour and time-to-breakdown with independent codes increased confidence in the results presented here. Our MD simulations identified high-flux, high-energy ion bombardment as a possible mechanism forming the early-stage surface damage in vacuum arcs. In this mechanism, sputtering occurs mostly in clusters, as a consequence of overlapping heat spikes. Different-sized experimental and simulated craters were found to be self-similar with a crater depth-to-width ratio of about 0.23 (sim) - 0.26 (exp). Experiments, which we carried out to investigate the energy dependence of DC breakdown properties, point at an intrinsic connection between DC and RF scaling laws and suggest the possibility of accumulative effects influencing the field enhancement factor.