Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

Light-matter Interactions in Semiconductors and Metals: From Nitride Optoelectronics to Quantum Plasmonics

This thesis puts forth a theory-directed approach coupled with spectroscopy aimed at the discovery and understanding of light-matter interactions in semiconductors and metals.

The first part of the thesis presents the discovery and development of Zn-IV nitride materials.The commercial prominence in the optoelectronics industry of tunable semiconductor alloy materials based on nitride semiconductor devices, specifically InGaN, motivates the search for earth-abundant alternatives for use in efficient, high-quality optoelectronic devices. II-IV-N2 compounds, which are closely related to the wurtzite-structured III-N semiconductors, have similar electronic and optical properties to InGaN namely direct band gaps, high quantum efficiencies and large optical absorption coefficients. The choice of different group II and group IV elements provides chemical diversity that can be exploited to tune the structural and electronic properties through the series of alloys. The first theoretical and experimental investigation of the ZnSnxGe1−xN2 series as a replacement for III-nitrides is discussed here.

The second half of the thesis shows ab−initio calculations for surface plasmons and plasmonic hot carrier dynamics. Surface plasmons, electromagnetic modes confined to the surface of a conductor-dielectric interface, have sparked renewed interest because of their quantum nature and their broad range of applications. The decay of surface plasmons is usually a detriment in the field of plasmonics, but the possibility to capture the energy normally lost to heat would open new opportunities in photon sensors, energy conversion devices and switching. A theoretical understanding of plasmon-driven hot carrier generation and relaxation dynamics in the ultrafast regime is presented here. Additionally calculations for plasmon-mediated upconversion as well as an energy-dependent transport model for these non-equilibrium carriers are shown.

Finally, this thesis gives an outlook on the potential of non-equilibrium phenomena in metals and semiconductors for future light-based technologies.

Classical and Quantum Effects in Plasmonic Metals

The field of plasmonics exploits the unique optical properties of metallic nanostructures to concentrate and manipulate light at subwavelength length scales. Metallic nanostructures get their unique properties from their ability to support surface plasmons– coherent wave-like oscillations of the free electrons at the interface between a conductive and dielectric medium. Recent advancements in the ability to fabricate metallic nanostructures with subwavelength length scales have created new possibilities in technology and research in a broad range of applications.

In the first part of this thesis, we present two investigations of the relationship between the charge state and optical state of plasmonic metal nanoparticles. Using experimental bias-dependent extinction measurements, we derive a potential- dependent dielectric function for Au nanoparticles that accounts for changes in the physical properties due to an applied bias that contribute to the optical extinction. We also present theory and experiment for the reverse effect– the manipulation of the carrier density of Au nanoparticles via controlled optical excitation. This plasmoelectric effect takes advantage of the strong resonant properties of plasmonic materials and the relationship between charge state and optical properties to eluci- date a new avenue for conversion of optical power to electrical potential.

The second topic of this thesis is the non-radiative decay of plasmons to a hot-carrier distribution, and the distribution’s subsequent relaxation. We present first-principles calculations that capture all of the significant microscopic mechanisms underlying surface plasmon decay and predict the initial excited carrier distributions so generated. We also preform ab initio calculations of the electron-temperature dependent heat capacities and electron-phonon coupling coefficients of plasmonic metals. We extend these first-principle methods to calculate the electron-temperature dependent dielectric response of hot electrons in plasmonic metals, including direct interband and phonon-assisted intraband transitions. Finally, we combine these first-principles calculations of carrier dynamics and optical response to produce a complete theoretical description of ultrafast pump-probe measurements, free of any fitting parameters that are typical in previous analyses.

Hole-Phonon Relaxation and Photocatalytic Properties of Titanium Dioxide and Zinc Oxide: First-Principles Approach

First-principles calculations for the temporal characteristics of hole-phonon relaxation in the valence band of titanium dioxide and zinc oxide have been performed. A first-principles method for the calculations of the quasistationary distribution function of holes has been developed. The results show that the quasistationary distribution of the holes in TiO2 extends to an energy level approximately 1eV below the top of the valence band. This conclusion in turn helps to elucidate the origin of the spectral dependence of the photocatalytic activity of TiO2. Analysis of the analogous data for ZnO shows that in this material spectral dependence of photocatalytic activity in the oxidative reactions is unlikely.

Triel Bonds, pi-Hole-pi-Electrons Interactions in Complexes of Boron and Aluminium Trihalides and Trihydrides with Acetylene and Ethylene

MP2/aug-cc-pVTZ calculations were performed on complexes of aluminium and boron trihydrides and trihalides with acetylene and ethylene. These complexes are linked through triel bonds where the triel center (B or Al) is characterized by the Lewis acid properties through its -hole region while -electrons of C2H2 or C2H4 molecule play the role of the Lewis base. Some of these interactions possess characteristics of covalent bonds, i.e., the Al--electrons links as well as the interaction in the BH3-C2H2 complex. The triel--electrons interactions are classified sometimes as the 3c-2e bonds. In the case of boron trihydrides, these interactions are often the preliminary stages of the hydroboration reaction. The Quantum Theory of Atoms in Molecules as well as the Natural Bond Orbitals approach are applied here to characterize the -hole--electrons interactions.

A influência dos elementos da primeira série de transição nas propriedades eletrônicas de aglomerados de paládio: um estudo baseado na Teoria do Funcional de Densidade (TFD)

Na área da catálise, interações de moléculas pequenas como H2,N2,O2, CO, NO etc, com aglomerados metálicos, têm nas suas adsorções químicas etapa fundamental nas atividades catalíticas. O paládio está entre os metais mais empregados em catalisadores usados em muitos processos catalíticos heterogêneos. Neste presente trabalho foi realizada uma série de cálculos baseado na Teoria do Funcional de Densidade (TFD ou DFT) empregando o método BP86/LANL2DZ/6-311+G(d,p) para estudar a adsorção da molécula NO que foi utilizada como protótipo molecular para interagir sobre aglomerados de paládio puros e dopados com metais de transição da primeira série de transição (Pd3M e Pd9M). Neste contexto, o escopo deste trabalho foi: (i) obter informações a respeito das possíveis alterações estruturais e eletrônicas ocorridas em aglomerados de paládio quando dopados com metais da primeira série de transição (MT), (ii) estudar as alterações no comportamento do átomo de paládio, nos seus respectivos aglomerados puros, frente aos aglomerados dopados com metais de transição da primeira série,(iii) buscar possíveis padrões de influência dos MT nas propriedades de aglomerados de paládio (ex: energias dos orbitais de fronteira, transferências eletrônicas de doação e retrodoação, densidades de spin, hibridizações dos orbitais responsáveis pelas ligações com moléculas adsorventes etc). Após a realização dos cálculos baseados na DFT, os resultados mostram que, para o estudo sobre os aglomerados contendo quatro átomos, foi obtida correlação direta entre a carga adquirida pelo metal M dentro do aglomerado metálico e os orbitais moleculares de fronteira LUMO calculados para Pd3M. Esta correlação direta não se mantém quando o tamanho do aglomerado de paládio é aumento de quatro para dez átomos. A energia de adsorção apresentada pela molécula de NO apresenta boa correlação com a energia do LUMO, independentemente do número de átomos e da geometria do aglomerado e da natureza do metal dopante. A molécula de NO adsorve mais favoravelmente no modo Bridge, independentemente de qual metal está dopando o aglomerado de Pd9M. Entretanto, para o aglomerado de Pd3M, o modo de adsorção dependerá da natureza do metal dopante. A mudança na geometria e no número de átomos de paládio existentes no aglomerado provoca mudanças no modo e na energia de adsorção da molécula de NO adsorvente. Uma investigação mais aprofundada deve encontrar outras possíveis correlações entre propriedades dos metais da primeira série e aglomerados de paládio dopados com os mesmos

A cohesive approach to thin-shell fracture and fragmentation

We develop a finite-element method for the simulation of dynamic fracture and fragmentation of thin-shells. The shell is spatially discretized with subdivision shell elements and the fracture along the element edges is modeled with a cohesive law. In order to follow the propagation and branching of cracks, subdivision shell elements are pre-fractured ab initio and the crack opening is constrained prior to crack nucleation. This approach allows for shell fracture in an in-plane tearing mode, a shearing mode, or a bending of hinge mode. The good performance of the method is demonstrated through the simulation of petalling failure experiments in aluminum plates. © 2005 Elsevier B.V. All rights reserved.

Comparison of quantum chemistry calculation methods of nucleic acid base-quartets

Nine base-quartets were calculated by six semi-empirical methods and ab initio Hartree-Fork method using STO-3G basis set. The results showed that PM3 method can be use to calculate base quartets, the results of PM3 calculations are close to the ab initio

Ballasted & electrically steerable carbon nanotube field emitters

Here we present our on-going efforts toward the development of stable ballasted carbon nanotube-based field emitters employing hydrothermally synthesized zinc oxide nanowires and thin film silicon-on-insulator substrates. The semiconducting channel in each controllably limits the emission current thereby preventing detrimental burn-out of individual emitters that occurs due to unavoidable statistical variability in emitter characteristics, particularly in their length. Fabrication details and emitter characterization are discussed in addition to their field emission performance. The development of a beam steerable triode electron emitter formed from hexagonal carbon nanotube arrays with central focusing nanotube electrodes, is also described. Numerical ab-initio simulations are presented to account for the empirical emission characteristics. Our engineered ballasted emitters have shown some of the lowest reported lifetime variations (< 0.7%) with on-times of < 1 ms, making them ideally-suited for next-generation displays, environmental lighting and portable x-rays sources. © 2012 SPIE.

Aromaticity on the fly: cyclic transition state stabilization at finite temperature.

We study the transition state of pericyclic reactions at elevated temperature with unbiased ab initio molecular dynamics. We find that the transition state of the intramolecular rearrangements for barbaralane and bullvalene remains aromatic at high temperature despite the significant thermal atomic motions. Structural, magnetic, and electronic properties of the dynamical transition state show the concertedness and aromatic character. Free-energy calculations also support the validity of the transition state theory for the present rearrangement reactions. The calculations demonstrate that cyclic delocalization represents a strong force to synchronize the thermal atomic motions even at high temperatures.

Does the orbital degeneracy play any role in the high thermopower of lamellar cobaltites?

This article investigates the role of the CoO6 octahedron distortion on the electronic properties and more particularly on the high value of the Seebeck coefficient in the BiCaCoO lamellar cobaltites. Our measurements provide clues indicating that the t2g orbital degeneracy lifting has to be considered to account for the observed high temperature limit of the thermopower. They also provide experimental arguments for locating the a1g and eg′ orbitals levels on the energy scale, through the compression of the octahedron. These results are in agreement with recent ab initio calculation including the electronic correlations and concluding for the inversion of these levels as compared to the expectation from the crystal field theory. © 2007 American Institute of Physics.

Near-field characterization of chemical vapor deposition graphene in the microwave regime

Near-field measurements were performed at X-band frequencies for graphene on copper microstrip transmission lines. An improvement in radiation of 0.88 dB at 10.2 GHz is exhibited from the monolayer graphene antenna which has dc sheet resistivity of 985 Ω/sq. Emission characteristics were validated via ab initio simulations and compared to empirical findings of geometrically comparable copper patches. This study contributes to the current knowledge of the electronic properties of graphene. © 2013 AIP Publishing LLC.

Tensile and compressive mechanical behavior of twinned silicon carbide nanowires

Molecular dynamics simulations with the Tersoff potential were used to study the response of twinned SiC nanowires under tensile and compressive strain. The critical strain of the twinned nanowires can be enhanced by twin stacking faults, and their critical strains are larger than those of perfect nanowires with the same diameters. Under axial tensile strain, the bonds of the nanowires are stretched just before failure. The failure behavior is found to depend on the twin segment thickness and the diameter of the nanowires. An atomic chain is observed for thin nanowires with small twin segment thickness under tension strain. Under axial compressive strain, the collapse of twinned SiC nanowires exhibits two different failure modes, depending on the length and diameter of the nanowires, i.e., shell buckling for short nanowires and columnar buckling for longer nanowires.

Mechanical and chemical bonding properties of ground state BeH2

The crystal structure, mechanical properties and electronic structure of ground state BeH2 are calculated employing the first-principles methods based on the density functional theory. Our calculated structural parameters at equilibrium volume are well consistent with experimental results. Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. The bulk modulus B, Shear modulus G, Young's modulus E, and Poisson's ratio upsilon are deduced from the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics in band structure, density of states, and charge distribution. The ionicity in the Be-H bond is mainly featured by charge transfer from Be 2s to H 1s atomic orbitals while its covalency is dominated by the hybridization of H 1s and Be 2p states. The Bader analysis of BeH2 and MgH2 are performed to describe the ionic/covalent character quantitatively and we find that about 1.61 (1.6) electrons transfer from each Be (Mg) atom to H atoms.

Defects in gallium nitride nanowires: First principles calculations

Atomic configurations and formation energies of native defects in an unsaturated GaN nanowire grown along the [001] direction and with (100) lateral facets are studied using large-scale ab initio calculation. Cation and anion vacancies, antisites, and interstitials in the neutral charge state are all considered. The configurations of these defects in the core region and outermost surface region of the nanowire are different. The atomic configurations of the defects in the core region are same as those in the bulk GaN, and the formation energy is large. The defects at the surface show different atomic configurations with low formation energy. Starting from a Ga vacancy at the edge of the side plane of the nanowire, a N-N split interstitial is formed after relaxation. As a N site is replaced by a Ga atom in the suboutermost layer, the Ga atom will be expelled out of the outermost layers and leaves a vacancy at the original N site. The Ga interstitial at the outmost surface will diffuse out by interstitialcy mechanism. For all the tested cases N-N split interstitials are easily formed with low formation energy in the nanowires, indicating N-2 molecular will appear in the GaN nanowire, which agrees well with experimental findings.