Local numerical modelling of magnetoconvection and turbulence - implications for mean-field theories

During the last decades mean-field models, in which large-scale magnetic fields and differential rotation arise due to the interaction of rotation and small-scale turbulence, have been enormously successful in reproducing many of the observed features of the Sun. In the meantime, new observational techniques, most prominently helioseismology, have yielded invaluable information about the interior of the Sun. This new information, however, imposes strict conditions on mean-field models. Moreover, most of the present mean-field models depend on knowledge of the small-scale turbulent effects that give rise to the large-scale phenomena. In many mean-field models these effects are prescribed in ad hoc fashion due to the lack of this knowledge. With large enough computers it would be possible to solve the MHD equations numerically under stellar conditions. However, the problem is too large by several orders of magnitude for the present day and any foreseeable computers. In our view, a combination of mean-field modelling and local 3D calculations is a more fruitful approach. The large-scale structures are well described by global mean-field models, provided that the small-scale turbulent effects are adequately parameterized. The latter can be achieved by performing local calculations which allow a much higher spatial resolution than what can be achieved in direct global calculations. In the present dissertation three aspects of mean-field theories and models of stars are studied. Firstly, the basic assumptions of different mean-field theories are tested with calculations of isotropic turbulence and hydrodynamic, as well as magnetohydrodynamic, convection. Secondly, even if the mean-field theory is unable to give the required transport coefficients from first principles, it is in some cases possible to compute these coefficients from 3D numerical models in a parameter range that can be considered to describe the main physical effects in an adequately realistic manner. In the present study, the Reynolds stresses and turbulent heat transport, responsible for the generation of differential rotation, were determined along the mixing length relations describing convection in stellar structure models. Furthermore, the alpha-effect and magnetic pumping due to turbulent convection in the rapid rotation regime were studied. The third area of the present study is to apply the local results in mean-field models, which task we start to undertake by applying the results concerning the alpha-effect and turbulent pumping in mean-field models describing the solar dynamo.

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some of the current issues of the solar simulations are put forward.

Exotic superfluid states of ultra-cold Fermi gases

This is a study of ultra-cold Fermi gases in different systems. This thesis is focused on exotic superfluid states, for an example on the three component Fermi gas and the FFLO phase in optical lattices. In the two-components case, superfluidity is studied mainly in the case of the spin population imbalanced Fermi gases and the phase diagrams are calculated from the mean-field theory. Different methods to detect different phases in optical lattices are suggested. In the three-component case, we studied also the uniform gas and harmonically trapped system. In this case, the BCS theory is generalized to three-component gases. It is also discussed how to achieve the conditions to get an SU(3)-symmetric Hamiltonian in optical lattices. The thesis is divided in chapters as follows: Chapter 1 is an introduction to the field of cold quantum gases. In chapter 2 optical lattices and their experimental characteristics are discussed. Chapter 3 deals with two-components Fermi gases in optical lattices and the paired states in lattices. In chapter 4 three-component Fermi gases with and without a harmonic trap are explored, and the pairing mechanisms are studied. In this chapter, we also discuss three-component Fermi gases in optical lattices. Chapter 5 devoted to the higher order correlations, and what they can tell about the paired states. Chapter 6 concludes the thesis.

Non-Gaussianities and Preheating from N-flation and Magnetogenesis via rotating cosmic string loops

In this thesis we examine multi-field inflationary models of the early Universe. Since non-Gaussianities may allow for the possibility to discriminate between models of inflation, we compute deviations from a Gaussian spectrum of primordial perturbations by extending the delta-N formalism. We use N-flation as a concrete model; our findings show that these models are generically indistinguishable as long as the slow roll approximation is still valid. Besides computing non-Guassinities, we also investigate Preheating after multi-field inflation. Within the framework of N-flation, we find that preheating via parametric resonance is suppressed, an indication that it is the old theory of preheating that is applicable. In addition to studying non-Gaussianities and preheatng in multi-field inflationary models, we study magnetogenesis in the early universe. To this aim, we propose a mechanism to generate primordial magnetic fields via rotating cosmic string loops. Magnetic fields in the micro-Gauss range have been observed in galaxies and clusters, but their origin has remained elusive. We consider a network of strings and find that rotating cosmic string loops, which are continuously produced in such networks, are viable candidates for magnetogenesis with relevant strength and length scales, provided we use a high string tension and an efficient dynamo.

On cluster approximations of cellular automata

In this thesis I examine one commonly used class of methods for the analytic approximation of cellular automata, the so-called local cluster approximations. This class subsumes the well known mean-field and pair approximations, as well as higher order generalizations of these. While a straightforward method known as Bayesian extension exists for constructing cluster approximations of arbitrary order on one-dimensional lattices (and certain other cases), for higher-dimensional systems the construction of approximations beyond the pair level becomes more complicated due to the presence of loops. In this thesis I describe the one-dimensional construction as well as a number of approximations suggested for higher-dimensional lattices, comparing them against a number of consistency criteria that such approximations could be expected to satisfy. I also outline a general variational principle for constructing consistent cluster approximations of arbitrary order with minimal bias, and show that the one-dimensional construction indeed satisfies this principle. Finally, I apply this variational principle to derive a novel consistent expression for symmetric three cell cluster frequencies as estimated from pair frequencies, and use this expression to construct a quantitatively improved pair approximation of the well-known lattice contact process on a hexagonal lattice.

Soluautomaattien approksimaatioista

Tässä tutkielmassa tarkastelen soluautomaatteja ja niiden yleisesti käytettyj ä keskivertokenttä- ("mean field") ja pariapproksimaatioita. Määrittelen soluautomaatin käsitteen yleisellä tasolla ja kuvailen kolme eri tyypillisesti käytettyä soluautomaattityyppiä, joista kahdesta esimerkkeinä mainitsen John Conwayn tunnetun "Game of Life" -soluautomaatin sekä epidemiologian alalla yleisesti käytetyn hilakontaktiprosessin. Tutkielman toisessa osassa esittelen eri tapoja, joilla soluautomaattien käyttäytymistä voidaan analyyttisesti approksimoida ja karakterisoida, mukaanlukien niin kutsuttu soluautomaatin aktiviteetti (λ) sekä yllämainitut keskivertokenttä- ja pariapproksimaatiot. Kahta jälkimmäistä sovellan hilakontaktiprosessiin, ja vertailen näin saatuja tuloksia.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

Effects of seed origin latitude on the timing of height growth cessation and field performance of silver birch

The aim of this thesis was to increase our knowledge about the effects of seed origin on the timing of height growth cessation and field performance of silver birch from different latitudes, with special attention paid to the browsing damage by moose in young birch plantations. The effect of seed origin latitude and sowing time on timing of height growth cessation of first-year seedlings was studied in a greenhouse experiment with seven seed origins (lat. 58º - 67ºN). Variation in critical night length (CNL) for 50 % bud set within two latitudinally distant stands (60º and 67ºN) was studied in three phytotron experiments. Browsing by moose on 5-11 -year-old silver birch saplings from latitudinally different seed origins (53º - 67ºN) was studied in a field experiment in southern Finland. Yield and stem quality of 22-year-old silver birch trees of Baltic, Finnish and Russian origin (54º - 63ºN) and the effect of latitudinal seed transfers were studied in two provenance trials at Tuusula, southern and Viitasaari, central Finland. The timing of height growth cessation depended systematically on latitude of seed origin and sowing date. The more northern the seed origin, the earlier the growth cessation and the shorter the growth period. Later sowing dates delayed growth cessation but also shortened the growth period. The mean CNL of the southern ecotype was longer, 6.3 ± 0.2 h (95 % confidence interval), than that of the northern ecotype, 3.1 ± 0.3 h. Within-ecotype variance of the CNL was higher in the northern ecotype (0.484 h2) than in the southern ecotype (0.150 h2). Browsing by moose decreased with increasing latitude of seed origin and sapling height. Origins transferred from more southern latitudes were more heavily browsed than the more northern native ones. Southern Finnish seed origins produced the highest volume per unit area in central Finland (lat. 63º11'N). Estonian and north Latvian stand seed origins, and the southern Finnish plus tree origins, were the most productive ones in southern Finland (lat. 60º21'N). Latitudinal seed transfer distance had a significant effect on survival, stem volume/ha and proportion of trees with a stem defect. The relationship of both survival and stem volume/ha to the latitudinal seed transfer distance was curvilinear. Volume was increased by transferring seed from ca. 2 degrees of latitude from the south. A longer transfer from the south, and transfer from the north, decreased the yield. The proportion of trees with a stem defect increased linearly in relation to the latitudinal seed transfer distance from the south.

Noncommutative Quantum Field Theories and UV/IR Mixing

Our present-day understanding of fundamental constituents of matter and their interactions is based on the Standard Model of particle physics, which relies on quantum gauge field theories. On the other hand, the large scale dynamical behaviour of spacetime is understood via the general theory of relativity of Einstein. The merging of these two complementary aspects of nature, quantum and gravity, is one of the greatest goals of modern fundamental physics, the achievement of which would help us understand the short-distance structure of spacetime, thus shedding light on the events in the singular states of general relativity, such as black holes and the Big Bang, where our current models of nature break down. The formulation of quantum field theories in noncommutative spacetime is an attempt to realize the idea of nonlocality at short distances, which our present understanding of these different aspects of Nature suggests, and consequently to find testable hints of the underlying quantum behaviour of spacetime. The formulation of noncommutative theories encounters various unprecedented problems, which derive from their peculiar inherent nonlocality. Arguably the most serious of these is the so-called UV/IR mixing, which makes the derivation of observable predictions especially hard by causing new tedious divergencies, to which our previous well-developed renormalization methods for quantum field theories do not apply. In the thesis I review the basic mathematical concepts of noncommutative spacetime, different formulations of quantum field theories in the context, and the theoretical understanding of UV/IR mixing. In particular, I put forward new results to be published, which show that also the theory of quantum electrodynamics in noncommutative spacetime defined via Seiberg-Witten map suffers from UV/IR mixing. Finally, I review some of the most promising ways to overcome the problem. The final solution remains a challenge for the future.

Conformal Field Theory Methods for Variants of Schramm-Loewner Evolutions

This thesis consists of an introduction, four research articles and an appendix. The thesis studies relations between two different approaches to continuum limit of models of two dimensional statistical mechanics at criticality. The approach of conformal field theory (CFT) could be thought of as the algebraic classification of some basic objects in these models. It has been succesfully used by physicists since 1980's. The other approach, Schramm-Loewner evolutions (SLEs), is a recently introduced set of mathematical methods to study random curves or interfaces occurring in the continuum limit of the models. The first and second included articles argue on basis of statistical mechanics what would be a plausible relation between SLEs and conformal field theory. The first article studies multiple SLEs, several random curves simultaneously in a domain. The proposed definition is compatible with a natural commutation requirement suggested by Dubédat. The curves of multiple SLE may form different topological configurations, ``pure geometries''. We conjecture a relation between the topological configurations and CFT concepts of conformal blocks and operator product expansions. Example applications of multiple SLEs include crossing probabilities for percolation and Ising model. The second article studies SLE variants that represent models with boundary conditions implemented by primary fields. The most well known of these, SLE(kappa, rho), is shown to be simple in terms of the Coulomb gas formalism of CFT. In the third article the space of local martingales for variants of SLE is shown to carry a representation of Virasoro algebra. Finding this structure is guided by the relation of SLEs and CFTs in general, but the result is established in a straightforward fashion. This article, too, emphasizes multiple SLEs and proposes a possible way of treating pure geometries in terms of Coulomb gas. The fourth article states results of applications of the Virasoro structure to the open questions of SLE reversibility and duality. Proofs of the stated results are provided in the appendix. The objective is an indirect computation of certain polynomial expected values. Provided that these expected values exist, in generic cases they are shown to possess the desired properties, thus giving support for both reversibility and duality.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Aspects of Inflationary Models at Low Energy Scales

Cosmological inflation is the dominant paradigm in explaining the origin of structure in the universe. According to the inflationary scenario, there has been a period of nearly exponential expansion in the very early universe, long before the nucleosynthesis. Inflation is commonly considered as a consequence of some scalar field or fields whose energy density starts to dominate the universe. The inflationary expansion converts the quantum fluctuations of the fields into classical perturbations on superhorizon scales and these primordial perturbations are the seeds of the structure in the universe. Moreover, inflation also naturally explains the high degree of homogeneity and spatial flatness of the early universe. The real challenge of the inflationary cosmology lies in trying to establish a connection between the fields driving inflation and theories of particle physics. In this thesis we concentrate on inflationary models at scales well below the Planck scale. The low scale allows us to seek for candidates for the inflationary matter within extensions of the Standard Model but typically also implies fine-tuning problems. We discuss a low scale model where inflation is driven by a flat direction of the Minimally Supersymmetric Standard Model. The relation between the potential along the flat direction and the underlying supergravity model is studied. The low inflationary scale requires an extremely flat potential but we find that in this particular model the associated fine-tuning problems can be solved in a rather natural fashion in a class of supergravity models. For this class of models, the flatness is a consequence of the structure of the supergravity model and is insensitive to the vacuum expectation values of the fields that break supersymmetry. Another low scale model considered in the thesis is the curvaton scenario where the primordial perturbations originate from quantum fluctuations of a curvaton field, which is different from the fields driving inflation. The curvaton gives a negligible contribution to the total energy density during inflation but its perturbations become significant in the post-inflationary epoch. The separation between the fields driving inflation and the fields giving rise to primordial perturbations opens up new possibilities to lower the inflationary scale without introducing fine-tuning problems. The curvaton model typically gives rise to relatively large level of non-gaussian features in the statistics of primordial perturbations. We find that the level of non-gaussian effects is heavily dependent on the form of the curvaton potential. Future observations that provide more accurate information of the non-gaussian statistics can therefore place constraining bounds on the curvaton interactions.

Näkökulmia eurooppalaiseen sosiaalipolitiikkaan : malli, väestö, resurssit ja kommunikaatio

The study analyses European social policy as a political project that proceeds under the guidance of the European Commission. In the name of modernisation, the project aims to build a new idea for the welfare state. To understand the project, it is necessary to distance oneself from both the juridical competence of the European Union and the traditional national welfare state models. The question is about sharing problems, as well as solutions to them: it is the creation and sharing of common views, concepts and images that play a key role in European integration. Drawing on texts and speeches produced by the European Commission, the study throws light on the development of European social policy during the first years of the 2000s. The study "freeze-frames" the welfare debate having its starting points in the nation states in the name of the entity of Europe. The first article approaches the European social model as a story in itself, a preparatory, persuasive narrative that concerns the management of change. The article shows how the audience can be motivated to work towards a set target by using discursive elements in a persuasive manner: the function of a persuasive story is to convince the target audience of the appropriateness of the chosen direction and to shape their identity so that they are favourably disposed to the desired political targets. This is a kind of "intermediate state" where the story, despite its inner contradictions and inaccuracies, succeeds in appearing as an almost self-evident path towards a modern social policy that Europe is currently seen to be in need of. The second article outlines the European social model as a question of governance. Health as a sector of social policy is detached from the old political order, which was based on the welfare state, and is closely linked to economy. At the same time the population is primarily seen as an economic resource. The Commission is working towards a "Europe of Health" that grapples with the problem of governance with the help of the "healthisation" of society, healthy citizenship and health economics. The way the Commission speaks is guided by the Union's powerful interest to act as "Europe" in the field of welfare policy. At the same time, the traditional separateness of health policy is effaced in order to be able to make health policy reforms a part of the Union's wider modernisation targets. The third article then shows the European social policy as its own area of governance. The article uses an approach based on critical discourse analysis in examining the classification systems and presentation styles adopted by Commission communications, as well as the identities that they help build. In analysing the "new start" of the Lisbon strategy from the perspective of social policy, the article shows how the emphasis has shifted from the persuasive arguments for change with necessary common European targets in the early stages of the strategy towards the implementation of reforms: from a narrative to a vision and from a diagnosis to healing. The phase of global competition represents "the modern" with which European society with its culture and ways of life now has to be matched. The Lisbon strategy is a way to direct this societal change, thus building a modern European social policy. The fourth article describes how the Commission uses its communications policy to build practices and techniques of governance and how it persuades citizens to participate in the creation of a European project of change. This also requires a new kind of agency: agents for whom accountability and responsibilities mean integration into and commitment to European society. Accountability is shaped into a decisive factor in implementing the European Union's strategy of change. As such it will displace hierarchical confrontations and emphasise common action with a view to modernising Europe. However, the Union's discourse cannot be described as being a political language that would genuinely rouse and convince the audience at the level of everyday life. Keywords: European social policy, EU policy, European social model, European Commission, modernisation of welfare, welfare state, communications, discoursiveness.

Conformal Field Theory Methods for Variants of Schramm-Loewner Evolutions

This thesis consists of an introduction, four research articles and an appendix. The thesis studies relations between two different approaches to continuum limit of models of two dimensional statistical mechanics at criticality. The approach of conformal field theory (CFT) could be thought of as the algebraic classification of some basic objects in these models. It has been succesfully used by physicists since 1980's. The other approach, Schramm-Loewner evolutions (SLEs), is a recently introduced set of mathematical methods to study random curves or interfaces occurring in the continuum limit of the models. The first and second included articles argue on basis of statistical mechanics what would be a plausible relation between SLEs and conformal field theory. The first article studies multiple SLEs, several random curves simultaneously in a domain. The proposed definition is compatible with a natural commutation requirement suggested by Dubédat. The curves of multiple SLE may form different topological configurations, ``pure geometries''. We conjecture a relation between the topological configurations and CFT concepts of conformal blocks and operator product expansions. Example applications of multiple SLEs include crossing probabilities for percolation and Ising model. The second article studies SLE variants that represent models with boundary conditions implemented by primary fields. The most well known of these, SLE(kappa, rho), is shown to be simple in terms of the Coulomb gas formalism of CFT. In the third article the space of local martingales for variants of SLE is shown to carry a representation of Virasoro algebra. Finding this structure is guided by the relation of SLEs and CFTs in general, but the result is established in a straightforward fashion. This article, too, emphasizes multiple SLEs and proposes a possible way of treating pure geometries in terms of Coulomb gas. The fourth article states results of applications of the Virasoro structure to the open questions of SLE reversibility and duality. Proofs of the stated results are provided in the appendix. The objective is an indirect computation of certain polynomial expected values. Provided that these expected values exist, in generic cases they are shown to possess the desired properties, thus giving support for both reversibility and duality.