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.

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.

Spectroscopic studies of gas-phase molecular clusters

Spectroscopic investigations of hydrogen-bonding and van der Waals' interactions m molecular clusters were studied by the techniques of infrared predissociation and resonance-enhanced multiphoton ionization spectroscopies (REMPI). Ab initio calculations were applied in conjunction for data interpretation.

The infrared predissociation spectroscopy of CN^-•(H_2O)_n (n = 2 - 6) clusters was reported in the region of 2950 - 3850 cm^(-1). The hydrogen bondings for the C-site and N-site binding, and among the water molecules were identified for n = 2 to 4. A spectral transition was observed for n = 5 and 6, implying that the anion was surface-bound onto the water aggregates in larger clusters.

The infrared predissociation spectroscopy of Br^-•(NH_3) and I^-•(NH_3)_n (n =1-3) clusters was reported in the region of 3050-3450 cm^(-1). For the Br^-•(NH_3) complex, a dominating ionic NH stretch appeared at 3175 cm^(-1), and the weaker free NH stretch appeared at 3348 cm^(-1). The observed spectrum was consistent to the structure in which there was one nearly linear hydrogen bond between Br^- and the NH_3 moiety. For the I^- •(NH_3) complex, five distinct IR absorption bands were observed in the spectrum. The spectrum was not consistent with basic frequency patterns of three geometries considered in the ab initio calculations - complex with one, two and three hydrogen bondings between I^- and the NH_3 moiety. Substantial inhomogenous broadening were displayed in the spectra for I^-•(NH_3)_n (n =2-3), suggesting the presence of multiple isomers.

The REMPI spectroscopy of the bound 4p ^2П 1/2 and ^2П 3/2 states, and the dissociative 3d ^2Σ^+ 1/2 state in the Al•Ar complex was reported. The dissociative spectrum at Al^+ channel suggested the coupling of the 4p ^2П 1/2,3/2 states to the repulsive 3d ^2Σ^+ 1/2 state. The spin-electronic coupling was further manifested in the dissociative Al^+ spectrum of the 3d ^2Σ^+ 1/2 state. Using the potential energy curves obtained from ab initio calculations, a bound → continuum Franck-Condon-intensity simulation was performed and compared with the one-photon 3d ^2Σ^+ 1/2 profile. The agreement provided evidence for the petturbation above the Al(3d)Ar dissociation limit, and the repulsive character of the 3d ^2Σ^+ 1/2 state.

Microwave absorption and magnetization tunneling in Mn12-acetate molecular clusters

In this paper we study the effect of microwave absorption on the quantum relaxation rate of Mn12 molecular clusters. We have determined first the resonant frequencies of a microwave resonator containing a single crystal of Mn12-Acetate and measured initial isothermal magnetization curves while microwave power was put into the resonator. We have found that the tunneling rate changes one order of magnitude for certain frequencies. This suggests that the microwave shaking of the nuclear spin and molecular vibrational degrees of freedom is responsible for the huge increasing of the tunneling rate.

Following the creation of active gold nanocatalysts from phosphine-stabilized molecular clusters

The phosphine-stabilised gold cluster [Au6(Ph2P-o-tolyl)6](NO3)2 is converted into an active nanocatalyst for the oxidation of benzyl alcohol through low-temperature peroxide-assisted removal of the phosphines, avoiding the high-temperature calcination process. The process was monitored using in-situ X-ray absorption spectroscopy, which revealed that after a certain period of the reaction with tertiary butyl hydrogen peroxide, the phosphine ligands are removed to form nanoparticles of gold which matches with the induction period seen in the catalytic reaction. Density functional theory calculations show that the energies required to remove the ligands from the [Au6Ln]2+ increase significantly with successive removal steps, suggesting that the process does not occur at once but sequentially. The calculations also reveal that ligand removal is accompanied by dramatic re-arrangements in the topology of the cluster core.

Molecular-Based Magnetism in High-Spin Molecular Clusters and Three-Dimensional Networks Based on Cyanometalate Building Blocks

The field of molecule-based magnets is a relatively new branch of chemistry, which involves the design and study of molecular compounds that exhibit a spontaneous magnetic ordering below a critical temperature, Tc. One major goal involves the design of materials with tuneable Tc's for specific applications in memory storage devices. Molecule-based magnets with high magnetic ordering temperatures have recently been obtained from bimetallic and mixed-valence transition metal μ-cyanide complexes of the Prussian blue family. Since the μ-cyanide linkages permit an interaction between paramagnetic metal ions, cyanometalate building blocks have found useful applications in the field of molecule-based magnets. Our work involves the use of octacyanometalate building blocks for the self-assembly of two new classes of magnetic materials namely, high-spin molecular clusters which exhibit both ferromagnetic intra- and intercluster coupling, and specific extended network topologies which show long-range ferromagnetic ordering.

Energetic protons emitted from coulomb explosion dynamics of large-sized hydrogen clusters driven by an ultrashort intense laser pulse

Protons with very high kinetic energy of about 10keV and the saturation effect of proton energy for laser intensity have been observed in the interaction of an ultrashort intense laser pulse with large-sized hydrogen clusters. Including the cluster-size distribution as well as the laser-intensity distribution on the focus spot, the theoretical calculations based on a simplified Coulomb explosion model have been compared with our experimental measurements, which are in good agreement with each other.

Nuclear fusion from Coulomb explosions of deuterated methane clusters subjected to ultraintense femtosecond laser pulses

This paper reports that Coulomb explosions taken place in the experiment of heteronuclear deuterated methane clusters ((CD4)(n)) in a gas jet subjected to intense femtosecond laser pulses (170 mJ, 70 fs) have led to table-top laser driven DD nuclear fusion. The clusters produced in supersonic expansion had an average energies of deuterons produced in the laser-cluster interaction were 60 and 1.5 KeV, respectively. From DD collisons of energetic deuterons, a yield of 2.5(+/-0.4)x10(4) fusion neutrons of 2.45 MeV per shot was realized, giving rise to a neutron production efficiency of about 1.5 x 10(5) per joule of incident laser pulse energy. Theoretical calculations were performed and a fairly good agreement of the calculated neutron yield with that obtained from the present experiment was found.

Coulomb explosion of hydrogen clusters irradiated by an ultrashort intense laser pulse

The explosion dynamics of hydrogen clusters driven by an ultrashort intense laser pulse has been analyzed analytically and numerically by employing a simplified Coulomb explosion model. The dependence of average and maximum proton kinetic energy on cluster size, pulse duration, and laser intensity has been investigated respectively. The existence of an optimum cluster size allows the proton energy to reach the maximum when the cluster size matches with the intensity and the duration of the laser pulse. In order to explain our experimental results such as the measured proton energy spectrum and the saturation effect of proton energy, the effects of cluster size distribution as well as the laser intensity distribution on the focus spot should be considered. A good agreement between them is obtained.

Preparation of I-2 clusters and their absorption spectra

Samples have been prepared at different temperatures by loading It molecules into the cages of zeolite 5A, and the measurements of the absorption spectra have been carried out for the prepared samples. It is shown that 12 molecular clusters are formed in the cages of zeolite 5A, and it is also found that molecular clusters which are bonded with intermolecular forces have an important feature, namely, the intermolecular distance in molecular clusters can be changed on different preparing conditions and the blue shift of absorption edges can not be as the criterion of forming molecular clusters.

Estudo da arquitetura molecular de filmes layer-by-layer de ftalocianina tetrasulfonada de cobre

Pós-graduação em Ciência e Tecnologia de Materiais - FC

Theoretical study of the XP3 (X = Al, B, Ga) clusters

The lowest singlet and triplet states of AlP3, GaP3 and BP3 molecules with C-s, C-2v and C-3v symmetries were characterized using the B3LYP functional and the aug-cc-pVTZ and aug-cc-pVQZ correlated consistent basis sets. Geometrical parameters and vibrational frequencies were calculated and compared to existent experimental and theoretical data. Relative energies were obtained with single point CCSD(T) calculations using the aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets, and then extrapolating to the complete basis set (CBS) limit. (C) 2011 Elsevier B.V. All rights reserved.

Tris(trimethylsilyl)silyl stabilized phosphorus and lead clusters

Die Hypersilylgruppe (Me3Si)3Si stellt einen sehr sperrigen, Elektronen liefernden Substituenten dar und kann zur Stabilisierung niedriger Oxidationsstufen sowie ungewöhnlicher Strukturelemente dienen. Durch Reaktionen der base-freien Hypersilanide der Alkalimetalle sowie des Dihypersilylplumbandiyls mit unterschiedlichsten phosphorhaltigen Reagenzien konnten eine Reihe hypersilyl-stabilisierter Phosphor- und Bleicluster-Verbindungen erhalten werden. Kaliumhypersilanid reagiert in Toluol glatt mit weißem Phosphor bei Raumtemperatur in Toluol unter quantitativer Bildung von rotem Kalium-bis(hypersilyl)tetraphosphenid [(Me3Si)3Si]2P4K2 (1), einem Kaliumsalz des Tetraphosphens (Me3Si)3Si-PH-P=P-PH-Si(SiMe3)3. In Benzol oder Toluol steht 1 im Gleichgewicht mit dem dimeren Octaphosphanid [(Me3Si)3Si]4P8K4 (2). Bei längerem Stehen der toluolischen Lösungen zerfällt 1 langsam vermutlich in Folge einer Protolyse zum gelben Pentaphosphanid [(Me3Si)3Si]3P5K2 (4). Aus benzolischer Lösung konnte hingegen ein weiteres Oktaphosphanid, [(Me3Si)3Si]3P8K3 (5), isoliert werden. Führt man die Reaktion Kaliumhypersilanid mit P4 in stärker koordinierenden Lösungsmitteln wie Diethylether durch, so entstehen neben 1 größere Mengen des Triphosphenids [(Me3Si)3Si]2P3K (3); dieses enthält ein Triphosphaallyl-Anion mit partieller P-P-Doppelbindung. Setzt man Lithiumhypersilanid mit weißem Phosphor um, so beobachtet man eine vollständig andere Produktpallette. Als Hauptprodukte lassen Polyphosphane wie beispielsweise [(Me3Si)3Si]2P4 (6) nachweisen, das zu 1 analoge [(Me3Si)3Si]2P4Li2 (7) entsteht nur in vergleichsweise kleinen Mengen. In der Gegenwart von Hexahydro-1,3,5-trimethyl-S-triazin, entsteht aus Lithiumhypersilanid und P4 hingegen im wesentlichen [(Me3Si)3Si]2P3Li (8) neben beträchtlichen Mengen von (Me3Si)4Si. Dessen Bildung erfordert eine Si-Si-Bindungsspaltung im Verlauf der Reaktion. Die Reaktion von Natriumhypersilanid mit P4 verläuft sehr unübersichtlich, das Pentaphosphanid [(Me3Si)3Si]3P5Na2 (9) ist das einzige isolierbare Produkt. Setzt man 1 mit [(Me3Si)2Si]2Sn um, so bilden sich überraschenderweise, je nach verwendetem Solvens [(Me3Si)3Si]3P4SnK (10) oder [(Me3Si)3Si]2[(Me3Si)2N]P4SnK (11). Alle neuen Verbindungen wurden NMR-spektroskopisch charakterisiert, die Phosphenide 1, 7, 8 sowie die Phosphanide 2, 4, 5, 9, 10 darüber hinaus durch Kristallstrukturanalysen. Dihypersilylplumbandiyl und -stannandiyl reagieren bei tiefer Temperatur mit P4, MPH2 (M=Li, K), PMe3, and PH3 zu formalen Lewis-Säure-Base-Addukten. Die Addukte {[(Me3Si)3Si]2PbPH2}M [M = Li (15), K (18)], {{[(Me3Si)3Si]2Pb}2PH2}M [M = Li (19), K (20)], und [(Me3Si)3Si]2EPMe3 [E = Pb (21), Sn (22)] wurden als kristalline Feststoffe erhalten und konnten vollständig charakterisiert werden. Die metastabilen Addukte {[(Me3Si)3Si]2E}4P4 (E = Pb, Sn) und [(Me3Si)3Si]2PbPH3 konnten lediglich NMR-spektroskopisch nachgewiesen werden. Bei Raumtemperatur entstehen in Folge von Ligandenaustausch-Prozessen die kristallographisch charakterisierten Heterokubane [(Me3Si)3Si]4P4E4 [E = Pb (12), Sn (14)], das Diphosphen (Me3Si)3SiP=PSi(SiMe3)3 (13) sowie der Pb2P2-Heterocyclus [(Me3Si)3SiPbP(H)Si(SiMe3)3]2 (17). Bei tiefer Temperatur wird aus einer sehr langsamen Reaktion von Dihypersilylplumbandiyl und PH3 in sehr kleinen Ausbeuten ein weiteres, völlig unerwartetes Produkt gebildet: der Bleicluster [(Me3Si)3Si]6Pb12 (23). Er weist ein verzerrt ikosaedrisches, zentrosymmetrisches Pb12-Gerüst auf. Nach jetzigen Erkenntnissen läuft seine Bildung über das nicht fassbare Hydridoplumbandiyl HPbSi(SiMe3)3, das intermediär durch Substituentenaustausch zwischen Pb[Si(SiMe3)3]2 and PH3 entsteht. Der Ersatz des Phosphans durch andere Hydridquellen wie (Ph3PCuH)6, (iBu)2AlH, and Me3NAlH3 führt ebenfalls zur Bildung von Bleiclustern, allerdings ist jetzt der Cluster [(Me3Si)3Si]6Pb10 (24) das Hauptprodukt. Beide Cluster, 23 und 24, gehorchen den Wade-Regeln.

Comparison of CBS-QB3, CBS-APNO, G2, and G3 thermochemical predictions with experiment for formation of ionic clusters of hydronium and hydroxide ions complexed with water

The GAUSSIAN 2, GAUSSIAN 3, complete basis set-QB3, and complete basis set-APNO methods have been used to calculate ΔH∘ and ΔG∘ values for ionic clusters of hydronium and hydroxide ions complexed with water. Results for the clusters H3O+(H2O)n andOH−(H2O)n, where n=1–4 are reported in this paper, and compared against experimental values contained in the National Institutes of Standards and Technology (NIST) database. Agreement with experiment is excellent for the three ab initio methods for formation of these clusters. The high accuracy of these methods makes them reliable for calculating energetics for the formation of ionic clusters containing water. In addition this allows them to serve as a valuable check on the accuracy of experimental data reported in the NIST database, and makes them useful tools for addressing unresolved issues in atmospheric chemistry.