Kapitza conductance and phonon scattering at grain boundaries by simulation

We use a nonequilibrium molecular-dynamics method to compute the Kapitza resistance of three twist grain boundaries in silicon, which we find to increase significantly with increasing grain boundary energy, i.e., with increasing structural disorder at the grain boundary. The origin of this Kapitza resistance is analyzed directly by studying the scattering of packets of lattice vibrations of well-defined polarization and frequency from the grain boundaries. We find that scattering depends strongly on the wavelength of the incident wave packet. In the case of a high-energy grain boundary, the scattering approaches the prediction of the diffuse mismatch theory at high frequencies, i.e., as the wavelength becomes comparable to the lattice parameter of the bulk crystal. We discuss the implications of our results in terms of developing a general model of scattering probabilities that can be applied to mesoscale models of heat transport in polycrystalline systems

Scattering of g-process longitudinal optical phonons at hotspots in silicon

Transistors with gate lengths below 100 nm generate phonon hotspots with dimensions on the order of 10 nm and peak power densities of about 50 W/mum(3). This work employs molecular dynamics to investigate the impact of lattice energy density on phonon scattering at the hotspot. The hotspot studied in this work consists of longitudinal optical phonons involved in the g-type intervalley scattering of conduction electrons in silicon. A comparison of the decay modes in hotspots with high and moderate energy densities reveals that the decay mechanisms are the same but the relaxation rates differ. Scattering occurs through a three phonon process of the form LO-- \u3e LA+TA, involving the zone-edge transverse acoustic modes. An increase in the energy density from a moderate value of 5 to 125 W/mum(3) changes the relaxation time from 79 to 16 ps, approximately proportional to the the maximum initial amplitude of the phonons. This work improves the accuracy of the scattering rates of optical phonons and helps in advancing the electro-thermal modeling of nanotransistors

Phonon-defect scattering in doped silicon by molecular dynamics simulation

Molecular dynamics simulations are used to study the scattering of phonon wave packets of well-defined frequency and polarization from individual point defects and from a field of point defects in Si. The relative amounts of energy in the transmitted and reflected phonon fields are calculated and the parameters that influence the phonon scattering process are determined. The results show that the fractions of transmitted and reflected energies strongly depend on the frequency of the incident phonons and on the mass and concentration of the defects. These results are compared with the classic formula for the scattering strength for point defects derived by Klemens, which we find to be valid when each phonon-defect scattering event is independent. The Klemens formula fails when coupled multiple scattering dominates. The phonon density of states is used to characterize the effects of point defects on mode mixing

Atomistic simulation of nanocrystalline materials

Atomistic simulations show that high-energy grain boundaries in nanocrystalline copper and nanocrystalline silicon are highly disordered. In the case of silicon the structures of the grain boundaries are essentially indistinguishable from that of bulk amorphous silicon. Based on a free-energy argument, we suggest that below a critical grain size nanocrystalline materials should be unstable with respect to the amorphous phase

Impact of Homogeneous Strain On Uranium Vacancy Diffusion In Uranium Dioxide

We present a detailed mechanism of, and the effect of homogeneous strains on, the migration of uranium vacancies in UO2. Vacancy migration pathways and barriers are identified using density functional theory and the effect of uniform strain fields are accounted for using the dipole tensor approach. We report complex migration pathways and noncubic symmetry associated with the uranium vacancy in UO2 and show that these complexities need to be carefully accounted for to predict the correct diffusion behavior of uranium vacancies. We show that under homogeneous strain fields, only the dipole tensor of the saddle with respect to the minimum is required to correctly predict the change in the energy barrier between the strained and the unstrained case. Diffusivities are computed using kinetic Monte Carlo simulations for both neutral and fully charged state of uranium single and divacancies. We calculate the effect of strain on migration barriers in the temperature range 800–1800 K for both vacancy types. Homogeneous strains as small as 2% have a considerable effect on diffusivity of both single and divacancies of uranium, with the effect of strain being more pronounced for single vacancies than divacancies. In contrast, the response of a given defect to strain is less sensitive to changes in the charge state of the defect. Further, strain leads to anisotropies in the mobility of the vacancy and the degree of anisotropy is very sensitive to the nature of the applied strain field for strain of equal magnitude. Our results suggest that the influence of strain on vacancy diffusivity will be significantly greater when single vacancies dominate the defect structure, such as sintering, while the effects will be much less substantial under irradiation conditions where divacancies dominate

Efficiency of free energy calculations of spin lattices by spectral quantum algorithms

Quantum algorithms are well-suited to calculate estimates of the energy spectra for spin lattice systems. These algorithms are based on the efficient calculation of the discrete Fourier components of the density of states. The efficiency of these algorithms in calculating the free energy per spin of general spin lattices to bounded error is examined. We find that the number of Fourier components required to bound the error in the free energy due to the broadening of the density of states scales polynomially with the number of spins in the lattice. However, the precision with which the Fourier components must be calculated is found to be an exponential function of the system size.Comment: 9 pages, 4 figures; corrected typographical and minor mathematical error

Comparison of the structure of grain boundaries in silicon and diamond by molecular-dynamics simulations

Molecular-dynamics simulations were used to synthesize nanocrystalline silicon with a grain size of up to 75 {angstrom} by crystallization of randomly misoriented crystalline seeds from the melt. The structures of the highly-constrained interfaces in the nanocrystal were found to be essentially indistinguishable from those of high-energy bicrystalline grain boundaries (GBs) and similar to the structure of amorphous silicon. Despite disorder, these GBs exhibit predominantly four-coordinated (sp{sup 3}-like) atoms and therefore have very few dangling bonds. By contrast, the majority of the atoms in high-energy bicrystalline GBs in diamond are three-coordinated (sp{sup 2}-like). Despite the large fraction of three-coordinated GB carbon atoms, they are rather poorly connected amongst themselves, thus likely preventing any type of graphite-like electrical conduction through the GBs

Comparison of the Structure of Grain Boundaries in Silicon and Diamond by Molecular-Dynamics Simulations

Molecular-dynamics simulations were used to synthesize nanocrystalline silicon with a grain size of up to 75 {angstrom} by crystallization of randomly misoriented crystalline seeds from the melt. The structures of the highly-constrained interfaces in the nanocrystal were found to be essentially indistinguishable from those of high-energy bicrystalline grain boundaries (GBs) and similar to the structure of amorphous silicon. Despite disorder, these GBs exhibit predominantly four-coordinated (sp{sup 3}-like) atoms and therefore have very few dangling bonds. By contrast, the majority of the atoms in high-energy bicrystalline GBs in diamond are three-coordinated (sp{sup 2}-like). Despite the large fraction of three-coordinated GB carbon atoms, they are rather poorly connected amongst themselves, thus likely preventing any type of graphite-like electrical conduction through the GBs