Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Theoretical and computational approaches on heterogeneous nucleation

Nucleation is the first step of a first order phase transition. A new phase is always sprung up in nucleation phenomena. The two main categories of nucleation are homogeneous nucleation, where the new phase is formed in a uniform substance, and heterogeneous nucleation, when nucleation occurs on a pre-existing surface. In this thesis the main attention is paid on heterogeneous nucleation. This thesis wields the nucleation phenomena from two theoretical perspectives: the classical nucleation theory and the statistical mechanical approach. The formulation of the classical nucleation theory relies on equilibrium thermodynamics and use of macroscopically determined quantities to describe the properties of small nuclei, sometimes consisting of just a few molecules. The statistical mechanical approach is based on interactions between single molecules, and does not bear the same assumptions as the classical theory. This work gathers up the present theoretical knowledge of heterogeneous nucleation and utilizes it in computational model studies. A new exact molecular approach on heterogeneous nucleation was introduced and tested by Monte Carlo simulations. The results obtained from the molecular simulations were interpreted by means of the concepts of the classical nucleation theory. Numerical calculations were carried out for a variety of substances nucleating on different substances. The classical theory of heterogeneous nucleation was employed in calculations of one-component nucleation of water on newsprint paper, Teflon and cellulose film, and binary nucleation of water-n-propanol and water-sulphuric acid mixtures on silver nanoparticles. The results were compared with experimental results. The molecular simulation studies involved homogeneous nucleation of argon and heterogeneous nucleation of argon on a planar platinum surface. It was found out that the use of a microscopical contact angle as a fitting parameter in calculations based on the classical theory of heterogeneous nucleation leads to a fair agreement between the theoretical predictions and experimental results. In the presented cases the microscopical angle was found to be always smaller than the contact angle obtained from macroscopical measurements. Furthermore, molecular Monte Carlo simulations revealed that the concept of the geometrical contact parameter in heterogeneous nucleation calculations can work surprisingly well even for very small clusters.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Late Transition Metal and Aluminum Complexes for the Polymerization of Ethene and Acrylates

Polyethene, polyacrylates and polymethyl acrylates are versatile materials that find wide variety of applications in several areas. Therefore, polymerization of ethene, acrylates and methacrylates has achieved a lot attention during past years. Numbers of metal catalysts have been introduced in order to control the polymerization and to produce tailored polymer structures. Herein an overview on the possible polymerization pathways for ethene, acrylates and methacrylates is presented. In this thesis iron(II) and cobalt(II) complexes bearing tri- and tetradentate nitrogen ligands were synthesized and studied in the polymerization of tertbutyl acrylate (tBA) and methyl methacrylate (MMA). Complexes are activated with methylaluminoxane (MAO) before they form active combinations for polymerization reactions. The effect of reaction conditions, i.e. monomer concentration, reaction time, temperature, MAO to metal ratio, on activity and polymer properties were investigated. The described polymerization system enables mild reaction conditions, the possibility to tailor molar mass of the produced polymers and provides good control over the polymerization. Moreover, the polymerization of MMA in the presence of iron(II) complex with tetradentate nitrogen ligands under conditions of atom transfer radical polymerization (ATRP) was studied. Several manganese(II) complexes were studied in the ethene polymerization with combinatorial methods and new active catalysts were found. These complexes were also studied in acrylate and methacrylate polymerizations after MAO activation and converted into the corresponding alkyl (methyl or benzyl) derivatives. Combinatorial methods were introduced to discover aluminum alkyl complexes for the polymerization of acrylates and methacrylates. Various combinations of aluminum alkyls and ligands, including phosphines, salicylaldimines and nitrogen donor ligands, were prepared in situ and utilized to initiate the polymerization of tBA. Phosphine ligands were found to be the most active and the polymerization MMA was studied with these active combinations. In addition, a plausible polymerization mechanism for MMA based on ESI-MS, 1H and 13C NMR is proposed.

Genetic background of bipolar disorder and related cognitive impairments

Bipolar disorder (BP) is a complex psychiatric disorder characterized by episodes of mania and depression. BP affects approximately 1% of the world’s population and shows no difference in lifetime prevalence between males and females. BP arises from complex interactions among genetic, developmental and environmental factors, and it is likely that several predisposing genes are involved in BP. The genetic background of BP is still poorly understood, although intensive and long-lasting research has identified several chromosomal regions and genes involved in susceptibility to BP. This thesis work aims to identify the genetic variants that influence bipolar disorder in the Finnish population by candidate gene and genome-wide linkage analyses in families with many BP cases. In addition to diagnosis-based phenotypes, neuropsychological traits that can be seen as potential endophenotypes or intermediate traits for BP were analyzed. In the first part of the thesis, we examined the role of the allelic variants of the TSNAX/DISC1 gene cluster to psychotic and bipolar spectrum disorders and found association of distinct allelic haplotypes with these two groups of disorders. The haplotype at the 5’ end of the Disrupted-in-Schizophrenia-1 gene (DISC1) was over-transmitted to males with psychotic disorder (p = 0.008; for an extended haplotype p = 0.0007 with both genders), whereas haplotypes at the 3’ end of DISC1 associated with bipolar spectrum disorder (p = 0.0002; for an extended haplotype p = 0.0001). The variants of these haplotypes also showed association with different cognitive traits. The haplotypes at the 5’ end associated with perseverations and auditory attention, while the variants at the 3’ end associated with several cognitive traits including verbal fluency and psychomotor processing speed. Second, in our complete set of BP families with 723 individuals we studied six functional candidate genes from three distinct signalling systems: serotonin-related genes (SLC6A4 and TPH2), BDNF -related genes (BDNF, CREB1 and NTRK2) and one gene related to the inflammation and cytokine system (P2RX7). We replicated association of the functional variant Val66Met of BDNF with BP and better performance in retention. The variants at the 5’ end of SLC6A4 also showed some evidence of association among males (p = 0.004), but the widely studied functional variants did not yield any significant results. A protective four-variant haplotype on P2RX7 showed evidence of association with BP and executive functions: semantic and phonemic fluency (p = 0.006 and p = 0.0003, respectively). Third, we analyzed 23 bipolar families originating from the North-Eastern region of Finland. A genome-wide scan was performed using the 6K single nucleotide polymorphism (SNP) array. We identified susceptibility loci at chromosomes 7q31 with a LOD score of 3.20 and at 9p13.1 with a LOD score of 4.02. We followed up both linkage findings in the complete set of 179 Finnish bipolar families. The finding on chromosome 9p13 was supported (maximum LOD score of 3.02), but the susceptibility gene itself remains unclarified. In the fourth part of the thesis, we wanted to test the role of the allelic variants that have associated with bipolar disorder in recent genome-wide association studies (GWAS). We could confirm findings for the DFNB31, SORCS2, SCL39A3, and DGKH genes. The best signal in this study comes from DFNB31, which remained significant after multiple testing corrections. Two variants of SORCS2 were allelic replications and presented the same signal as the haplotype analysis. However, no association was detected with the PALB2 gene, which was the most significantly associated region in the previous GWAS. Our results indicate that BP is heterogeneous and its genetic background may accordingly vary in different populations. In order to fully understand the allelic heterogeneity that underlies common diseases such as BP, complete genome sequencing for many individuals with and without the disease is required. Identification of the specific risk variants will help us better understand the pathophysiology underlying BP and will lead to the development of treatments with specific biochemical targets. In addition, it will further facilitate the identification of environmental factors that alter risk, which will potentially provide improved occupational, social and psychological advice for individuals with high risk of BP.

The Finnish Leadership Style in Transition - A Study of Leadership Criteria in the Insurance Business, 1997-2004 (summary section only)

The aim of the study was to explore the importance of evaluating leadership criteria in Finland at leader/subordinate levels of the insurance industry. The overall purpose of the thesis is tackled and analyzed from two different perspectives: - by examining the importance of the leadership criteria and style of Finnish insurance business leaders and their subordinates - by examining the opinions of insurance business leaders regarding leadership criteria in two culturally different countries: the US and Finland. This thesis consists of three published articles that scrutinise the focal phenomena both theoretically and empirically. The main results of the study do not lend support to the existence of a universal model of leadership criteria in the insurance business. As a matter of fact, the possible model seems to be based more on the special organizational and cultural circumstances of the country in question. The leadership criteria seem to be quite stable irrespective of the comparatively short research time period (3–5 years) and hierarchical level (subordinate/leader). Leaders have major difficulties in changing their leadership style. In fact, in order to bring about an efficient organizational change in the company you have to alternate the leader. The cultural dimensions (cooperation and monitoring) identified by Finnish subordinates were mostly in line with those of their managers, whilst emphasizing more the aspect of monitoring employees, which could be seen from their point of view as another element of managers’ optimizing/efficiency requirements. In Finnish surveys the strong emphasis on cooperation and mutual trust become apparent by both subordinates and managers. The basic problem is still how to emphasize and balance them in real life in such a way that both parties are happy to work together on a common basis. The American surveys suggests hypothetically that in a soft market period (buyer’s market) managers employ a more relationship-oriented leadership style and correspondingly adapt their leadership style to a more task-oriented approach in a hard market phase (seller’s market). In making business better Finnish insurance managers could probably concentrate more on task-oriented items such as reviewing, budgeting, monitoring and goal-orientation. The study also suggests that the social safety net of the European welfare state ideology has so far shielded the culture-specific sense of social responsibility of Finnish managers from the hazards of free competition and globalization.

Cognitive functioning and its heritability in bipolar I disorder

Bipolar I disorder is a severe psychiatric disorder characterized by episodic mood alterations that can be manic, depressive or mixed. Bipolar disorder seems to be highly genetic, but the etiology of this complex disorder has remained elusive. In recent years, studies have found that euthymic patients with bipolar disorder may have impairments particularly in executive functioning, verbal learning and memory. These impairments may be present also among some of the relatives of these patients, who may be vulnerable to the disorder. Using neuropsychological variables as endophenotypes, i.e. intermediate phenotypes between genes and the phenotypes, has been suggested to aid search for the etiological background of the disorder, but evidence is sparse on whether these variables fulfill the criteria for endophenotypes. The present thesis is part of the Genetic Epidemiology and Molecular Genetics of Severe Mental Disorders in Finland project. The specific aim was to investigate whether neuropsychological test variables would indicate genetic liability to the disorder and could therefore be regarded as endophenotypes. Thus, cognitive functions and their heritability were studied in bipolar I disorder patients and in their unaffected first-degree relatives from a population-based sample of families, comparing them to a population-based control group. In order to add homogeneity to the subgroups of bipolar disorder patients and their relatives, cognitive functions and their heritability were further studied in a group of families affected by bipolar I disorder only (bipolar families) and another group of families affected by both bipolar I disorder and schizophrenia or schizoaffective disorders (mixed families). Finally, the effect of processing speed on other cognitive functions was investigated. The study showed that especially executive functioning and processing speed fulfilled the endophenotype criteria. Impairments in these functions were found in bipolar patients and in their relatives irrespective of other severe psychopathology in the family. These functions were highly heritable in these families. Study also showed that generalized impairment in verbal memory may associate more with bipolar disorder than to vulnerability to other psychotic disorders, and be more related to fully developed disease; impairments in verbal learning and memory were found only in patients, and they were not found to be highly heritable. Finally, the most potential endophenotype, i.e. processing speed, seemed to contribute to a range of other cognitive dysfunctions seen in bipolar disorder patients. Processing speed, in particular, has also been shown to be a valid endophenotype in subsequent association analyses in psychiatric genetics in Finland and internationally. Information concerning cognitive impairments and their association with the psychosocial consequences of bipolar disorder is important in planning treatment. It is also important to understand and acknowledge that patients may have cognitive impairments that affect their everyday life. Psychosocial interventions and neuropsychological rehabilitation may supplement other conventional treatments for bipolar patients.