Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Ab initio Molecular Dynamics Simulations of the Structural Response of Solids to Ultrashort Laser and XUV Pulses

The theoretical model and underlying physics described in this thesis are about the interaction of femtosecond-laser and XUV pulses with solids. The key to understand the basics of such interaction is to study the structural response of the materials after laser interaction. Depending on the laser characteristics, laser-solid interaction can result in a wide range of structural responses such as solid-solid phase transitions, vacuum phonon squeezing, ultrafast melting, generation of coherent phonons, etc. During my research work, I have modeled the systems irradiated by low-, medium- and high-laser intensities, and studied different types of structural dynamics of solids at various laser fluences.

On the transition from localized to delocalized electronic states in divalent-metal clusters

The transition from van der Waals to covalent bonding, which is expected to occur in divalent-metal clusters with increasing cluster size, is discussed. We propose a model which takes into account, within the same electronic theory, the three main competing contributions, namely the kinetic energy of the electrons, the Coulomb interactions between electrons, and the s \gdw p intraatomic transitions responsible for van der Waals like bonding. The model is solved by taking into account electron correlations using a generalized Gutzwiller approximation (slave boson method). The occurrence of electron localization is studied as a function of the interaction parameters and cluster size.

Theoretical study of the physicochemical properties of the light transactinides

A theoretical study of the physicochemical properties of elements 104, 105, and 106 and their compounds in the gas phase and aqueous solutions has been undertaken using relativistic atomic and molecular codes. Trends in properties such as bonding, ionization potentials, electron affinities, energies of electronic transitions, stabilities of oxidation states etc. have been defined within the corresponding chemical groups and within the transactinides. These trends are shown to be determined by increasing relativistic effects within the groups. The behaviour of some gas phase compounds and complexes in solutions is predicted for the gas chromatography and solvent extraction experiments. Redox potentials in aqueous solutions of these elements are estimated.

Electronic structure and properties of the group 4, 5 and 6 highest chlorides including elements 104, 105, and 106

Results of the Dirac-Slater discrete variational calculations for the group 4, 5, and 6 highest chlorides including elements 104, 105, and 106 have shown that the groups are not identical with respect to trends in the electronic structure and bonding. The charge density distribution data show that notwithstanding the basic increase in covalency within the groups this increase diminishes in going from group 4 to group 6. As a result, E106Cl_6 will be less stable toward thermal decomposition than WCl_6, which is confirmed by an estimated low E106-Cl bond energy. \delta H_form equal to -90.3 ± 6 kcal/rnol is obtained for E106Cl_6 in the gas phase, which is indicative of a very low stability of this compound. The stability of the maximum oxidation state is shown to decrease in the direction E104(+4) > E105(+5) > E106(+6).

Volume isotope shifts in low lying transitions of Au I

Electronic factors in the volume isotope shift have been calculated in an ab initio way with the relativistic Dirac-Fock method for a number of different optical single/and two-photon transitions in Au I. The agreement with a semi-empirical method is within 10% for the resonance transition. For this one and a few other transitions the effect of core excitation has been analyzed with the Multi-configuration Dirac-Fock method as well, and it was found to reduce the electronic factor in the order of 5 %.

The electronic structure of the group 6 oxyanions [MO_4]^2-, where M = Cr, Mo, W, and element 106

The electronic structure of the group 6 oxyanions [MO_4]^2-, where M = Cr, Mo, W, and element 106 have been calculated using the Dirac-Slater Discrete Variational method. Results of the calculations show a relative decrease in the metal-oxygen bond strengths for the [E106O_4]^2- ion in the solid state compared to that for the [WO_4]^2- anion. Calculated energies of the electronic charge-transfer transitions are indicative of a strong possible luminescence of [El06O_4]^2- in the blue-violet area. In solutions [El06O_4]^2- will be the most stable ion out of the entire series. Estimated reduction potential E^0 (El06O^2-_4/E106O^3-_4) equal to -1.60V shows only a slightly increased stability of the +6 oxidation state for element 106 in comparison with W.

Synthesis and characterisation of novel nitrogen and fluorine containing spirobifluorene derivatives for electronic applications and crystal engineering

Among organic materials, spirobifluorene derivatives represent a very attractive class of materials for electronic devices. These compounds have high melting points, glass transitions temperatures and morphological stability, which makes these materials suitable for organic electronic applications. In addition, some of spirobifluorenes can form porous supramolecular associations with significant volumes available for the inclusion of guests. These molecular associations based on the spirobifluorenes are noteworthy because they are purely molecular analogues of zeolites and other microporous solids, with potential applications in separation, catalysis, sensing and other areas.

First-principles electronic theory of non-collinear magnetic order in transition-metal nanowires

The structural, electronic and magnetic properties of one-dimensional 3d transition-metal (TM) monoatomic chains having linear, zigzag and ladder geometries are investigated in the frame-work of first-principles density-functional theory. The stability of long-range magnetic order along the nanowires is determined by computing the corresponding frozen-magnon dispersion relations as a function of the 'spin-wave' vector q. First, we show that the ground-state magnetic orders of V, Mn and Fe linear chains at the equilibrium interatomic distances are non-collinear (NC) spin-density waves (SDWs) with characteristic equilibrium wave vectors q that depend on the composition and interatomic distance. The electronic and magnetic properties of these novel spin-spiral structures are discussed from a local perspective by analyzing the spin-polarized electronic densities of states, the local magnetic moments and the spin-density distributions for representative values q. Second, we investigate the stability of NC spin arrangements in Fe zigzag chains and ladders. We find that the non-collinear SDWs are remarkably stable in the biatomic chains (square ladder), whereas ferromagnetic order (q =0) dominates in zigzag chains (triangular ladders). The different magnetic structures are interpreted in terms of the corresponding effective exchange interactions J(ij) between the local magnetic moments μ(i) and μ(j) at atoms i and j. The effective couplings are derived by fitting a classical Heisenberg model to the ab initio magnon dispersion relations. In addition they are analyzed in the framework of general magnetic phase diagrams having arbitrary first, second, and third nearest-neighbor (NN) interactions J(ij). The effect of external electric fields (EFs) on the stability of NC magnetic order has been quantified for representative monoatomic free-standing and deposited chains. We find that an external EF, which is applied perpendicular to the chains, favors non-collinear order in V chains, whereas it stabilizes the ferromagnetic (FM) order in Fe chains. Moreover, our calculations reveal a change in the magnetic order of V chains deposited on the Cu(110) surface in the presence of external EFs. In this case the NC spiral order, which was unstable in the absence of EF, becomes the most favorable one when perpendicular fields of the order of 0.1 V/Å are applied. As a final application of the theory we study the magnetic interactions within monoatomic TM chains deposited on graphene sheets. One observes that even weak chain substrate hybridizations can modify the magnetic order. Mn and Fe chains show incommensurable NC spin configurations. Remarkably, V chains show a transition from a spiral magnetic order in the freestanding geometry to FM order when they are deposited on a graphene sheet. Some TM-terminated zigzag graphene-nanoribbons, for example V and Fe terminated nanoribbons, also show NC spin configurations. Finally, the magnetic anisotropy energies (MAEs) of TM chains on graphene are investigated. It is shown that Co and Fe chains exhibit significant MAEs and orbital magnetic moments with in-plane easy magnetization axis. The remarkable changes in the magnetic properties of chains on graphene are correlated to charge transfers from the TMs to NN carbon atoms. Goals and limitations of this study and the resulting perspectives of future investigations are discussed.