MODE RESOLVED THERMAL TRANSPORT IN AMORPHOUS SILICON: A MOLECULAR DYNAMICS STUDY
|
|
- Laura Mills
- 6 years ago
- Views:
Transcription
1 MODE RESOLVED THERMAL TRANSPORT IN AMORPHOUS SILICON: A MOLECULAR DYNAMICS STUDY Christopher H. Baker and Pamela M. Norris Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, Virginia 2294 Amorphous materials have found application in all manner of consumer products and continue to show promise for next generation technologies. Yet thermal transport in this class of materials remains a point of dispute. We studied the mechanisms of thermal transport in amorphous silicon using molecular dynamics simulation and wavelet transform techniques. Perturbation of a single atom in the equilibrated structure introduced broadband vibrational energy, which propagated outward. Analyzing the system dynamics with the wavelet transform provided information on the spectral wave velocity and diffusivity. We compare initial results to the assumptions underlying two current theories of thermal transport in amorphous materials. Introduction Amorphous materials have applications in many technologies where heat transfer is of concern. Examples include common window glass, semiconductor alloys of silicon, 1 and phase-change memory, 2 among others. Just as for crystalline materials, the small length scales found in developing technologies requires a microscopic theory of thermal transport in amorphous materials. A successful theory must not only explain the bulk conductivity based on the microscopic structure, but it must also correctly describe the length dependent conductivity at small length scales, when device dimensions are on the order of the mean-free-path of energy carriers. Despite the theoretical and practical interest in amorphous materials over the last century or so, there remains controversy over how to model thermal transport through the highly disordered atomic structure found in amorphous materials. Theories of thermal transport in solid systems revolve around the vibrational modes of the system. In crystals, the normal modes of vibration arising from interatomic forces and the symmetry of the crystalline lattice are specifically referred to as phonons. Vibrations also no doubt play a role in the amorphous phase, where there is no lattice, and we will examine two theories based on this principle, those pioneered by 1) Cahill and Pohl 3,4 and 2) Feldman, Allen, and Bickham. 5 It should be noted that the models presented by these authors are by no means isolated works and rely on concepts developed by the physics and microscale thermal transport communities over decades and that this does not constitute an exhaustive list of amorphous transport models, e.g. the fracton model of transport. 6,7 The Cahill-Pohl model 3,4 approaches the issue from the crystalline state. The thermal conductivity is calculated by assuming a crystalline dispersion, but with the key modification that these phonons have a mean-free-path equal to one half their wavelength, instead of on the order of tens of nanometers to microns as they would in the crystalline phase. 8 Feldman et al. begin with the eigenstates of a simulated amorphous system itself, which give the allowed vibrational states. 5 By comparing the eigenvectors with the wavevectors of a crystal, the authors define three types of vibrons (analogous to phonons) that transmit energy through different mechanisms: propagons, diffusons, and locons, which transmit energy ballistically, diffusively, and not at all. We present a simulation methodology which enables the direct observation of microscopic heat conduction in amorphous materials, where the properties of heat carriers can be discriminated by frequency. A pure 28 Si system is used for this study because of the abundance of experimental and computational literature on the same system. In this way, the assumptions behind current theories of thermal transport in amorphous materials can be directly tested. Methods The methodology consists of two stages. In the first stage, an equilibrium amorphous structure at zero kelvin is generated. In the second stage, kinetic energy is injected into the system by giving one atom an instantaneous non-zero velocity. The deposited energy then radiates outward and its spatial evolution measured with respect to time. Since the Fourier spectrum of a step function contains energy at all frequencies, this method introduces kinetic energy at all frequencies into the system (but not with equal weight). The time and frequency component of the energy at a certain distance from the epicenter can be discerned using the wavelet transform. 9,1 Baker 1
2 The open source molecular dynamics program LAMMPS 11 was used for the simulation. To begin, a 15 3 unit cell cube of 28 Si was constructed. The Stillinger-Weber potential 12 was used to describe the atomic interactions with the parameter set formulated by Vink et al.. 13 We chose this parameter set since it is optimized for amorphous silicon (a- Si) specifically. A time step of.5 femtoseconds was used throughout the simulation. To produce the amorphous structure, the atomic velocities were set to a distribution corresponding to 5 K and run in the NPT ensemble (constant number, pressure, and temperature) for 4 picoseconds at 1 bar and 3 K, melting the crystal. The thermostat was then modified to cool the system down to about 1 K over the course of 7.5 ps, which is a cooling rate of about 4 x 1 14 K/s, a little faster than rates used in the literature. 14 At this point, an NPH ensemble (constant enthalpy) was employed in conjunction with a fictitious damping coefficient of.1 ev ps/å 2 to drain the remaining thermal energy over the course of 2 ps, leaving the system at a temperature on the order of 1 8 K. The atomic coordinates of the equilibrated structure were printed. Finally, all atomic velocities were set to zero in preparation for the dynamic portion of the simulation. Broadband vibrational energy was introduced to the system by setting the velocity of a single atom to 1. Å/ps (kinetic energy of about 3 mev) with a random directional orientation. The system was then allowed to evolve for 3 ps, printing the velocity components of each atom every 3 time steps. This generated 21 values for velocity as a function of time, where the velocity oscillates as a function of time. The time series of each atom was analyzed using the wavelet transform 9,1 to extract the kinetic energy content of the atom as a function of time and frequency. The resulting description of the energy was summed over atoms with the same radial distance (bin size of 1.95 Å) from the atom at the epicenter. The wavelet transform is essentially a compromise between analysis of a signal in the time domain versus the frequency domain; it contains both temporal and spectral information. For analysis of the data, a mother wavelet frequency of THz was used to give equal resolution in time and frequency at about 5 THz. Other parameters used in the transform were the same as used in a previous work, which should be consulted for details of the transform. 1 g(r) amorph. cryst r(å) Figure 1: Radial distribution function of the equilibrated structure. Triangles denote neighbor shells in a crystal of equivalent density (46.8 nm 3 ). The third nearest neighbor spike has disappeared which is characteristic of amorphous silicon. 15,16 Over a length scale on the order of about 8 Å, the atoms are essentially randomly distributed. n(r) (atoms/å) Sim. Expt r(å) Figure 2: Number distribution function of the equilibrated structure and experimental measurements of pure a-si. 17,18 The simulated system captures the amorphous structure with slight quantitative differences. Baker 2
3 g(θ) (atoms/θ) θ (deg.) Figure 3: Angular distribution function of the equilibrated structure. The peak corresponds to the crystalline tetrahedral angle, but has some finite spread because of the disorder in the system. Results A large amount of information can be derived from the equilibrated structure. The most important fingerprint of an amorphous material is its radial distribution function (RDF), which gives the number of atoms as a function of distance from an arbitrary central atom. The RDF is visualized in two ways: normalized with respect to the RDF of an ideal gas of equivalent density, or normalized such that integration yields the total number of atoms within a sphere of radius r. Both versions have been referred to as the RDF in the literature; here we use RDF for the former and number distribution function (NDF) for the latter. The RDF is plotted in Fig. 1. A short range order is apparent due to the spikes in g(r) within the first several angstroms which quickly decay in magnitude, yielding a completely random structure on the length scale of about 8 angstroms and greater. The third nearest neighbor shell in the crystal structure has disappeared, which is characteristic of a-si. The NDF, n(r), is compared to the experimentally determined function for a-si 17,18 in Fig. 2. The simulated a-si system shows good agreement with experiment. One notable difference is in the first neighbor peak, where the experimental distribution is slightly narrower. The angular distribution function, g(θ), also displays properties expected for an amorphous system, the distribution of angles between nearest neighbors is centered at the tetrahedral value of 19.5 and the width indicates the disorder in the system (Fig. 3). A notable feature is the tail at larger angles which is also present in another simulation of a-si. 14 Despite the good agreement with experimental probability (%) % 5.2% neighbors (atoms) Figure 4: Distribution of coordination number. Despite the close correspondence between n(r) with experimental data, the coordination number is overestimated in the simulated system; it is predominantly 4 in real a-si with a small percentage of 3 and 5 coordinated atoms. results for the number distribution function and the expected shape of g(θ), other metrics indicate that the obtained amorphous system does not completely reproduce experimental a-si. Counting the number of nearest neighbors for each atom, defined as those that lie within the first minimum in n(r) results in a histogram of the coordination number of the system (Fig. 4). Whereas atoms are predominantly 5- fold coordinated in the simulated system, the experimental value is an average of 3.88 coordination, where most atoms are four-fold coordinated. 17 Even though amorphous systems are completely disordered, defects are typically defined in this context as atoms that are under- or over-coordinated with respect to the crystalline system. Hence the simulated system is heavily defected (over-abundance of 5-fold coordinated atoms) with respect to experimental a-si. This over coordination, perhaps arising from the interatomic potential overestimating the Si-Si bond strength, may also be the cause of a difference in the vibrational properties between the simulated system and the experimental system. The vibrational density of states was calculated using the harmonic lattice dynamics program GULP 19 and compared with experimental values found using inelastic neutron scattering 2,21 (Fig. 5). Agreement with the simulated DOS is obtained only after scaling the experimental DOS by a factor of 1.55 in frequency. After shifting, there is good qualitative agreement between the two curves, indicating that, while not exactly reproducing the amorphous state of a-si, the simulated system does show all the characteristics of an amorphous system. Baker 3
4 DOS (arb. units) Sim. Expt f (THz) Figure 5: Density of states of the simulated and experimental systems. 2,21 The experimental DOS has been shifted by a factor of 1.55 in frequency and both distributions are normalized to have unit area. The shift causes good qualitative agreement between the simulated and experimental systems. Figure 6: Energy distribution versus time for 2.9 THz frequency. The color is proportional to the wavelet energy density (ev/ps THz); absolute magnitudes should not be compared between figures. The analysis of the dynamic portion of the simulation is shown in Figs The figures present the radial distribution of energy as a function of time for increasing frequencies. The wavelet transform introduces some uncertainty into the results, having a lower temporal resolution at low frequencies compared to higher frequencies. This manifests itself in the width of the energy ridge in Fig. 6. At higher frequencies, the temporal resolution increases (with the tradeoff occurring in the frequency resolution), leading to a narrower stem near the origin (e.g. Figs. 7 and 8). As the frequency increases, the energy goes from pure wave transport (Fig. 6), to diffusive-wave transport (Fig. 7), to almost purely diffusive transport (Fig. 8). Linearly fitting the wave propagation in Fig. 6 results in a velocity of about 2.4 km/s. Energy distribution versus time for 8.21 THz fre- Figure 7: quency. Discussion We now discuss the implications of our data with respect to two theories of thermal transport in amorphous materials. In the Cahill-Pohl model, 3,4 which is often used to estimate the thermal conductivity of amorphous solids, the key assumptions are that heat carriers travel at the crystalline sound speed (for longitudinal and transverse acoustic branches) and that the mean-free-path of these modes is one-half their wavelength. This results in an expression for the conductivity: Energy distribution versus time for THz fre- Figure 8: quency. Baker 4
5 f (THz) LO LA TO TA k (A 1 ) Figure 9: Dispersion of crystalline phase in the (1) direction. Acoustic, optical, longitudinal, and transverse modes are noted. Dashed line is the zone center extrapolation, equivalent to Debye type dispersion. π Λ = 6 1/3 T 2 Θi kb n 2/3 /T v i i Θ i x 3 e x (e x 1) 2 dx, hf (1) where k B is Boltzmann s constant, n is the atomic number density, v i is the Debye velocity of mode i, T is temperature, Θ i is the mode specific Debye temperature (proportional to the frequency at maximum wavevector), and x is the reduced phonon energy: k B T. The strength of this theory is that it has no adjustable parameters and requires only basic knowledge of the crystalline phase of the material in question. Yet, from Fig. 6, the measured velocity in the amorphous phase is only 2.4 km/s compared to Debye velocities of 11. km/s (longitudinal) and 7.1 km/s (transverse) in the crystalline phase at the same density (Fig. 9). This suggests two possible explanations for why the Cahill-Pohl model still yields accurate estimates of Λ. The first is that it overestimates group velocities while underestimating meanfree-paths. Referring to Figs. 9 and 6 the halfwavelength of Debye modes with a frequency of 2.9 THz is about 26.4 Å (LA) or 17. Å (TA). The energy in Fig. 6 clearly travels farther than this distance without showing diffusive energy loss. Furthermore, experimental evidence has shown that the sound speed of a-si is less than that of the crystalline phase, 22 although their value is about 8% of the weighted crystalline sound speed versus our value of 29%. Our much lower value of sound speed may be due to the over-coordinated, highly defected structure. It is unclear how this would affect comparisons made to the Cahill-Pohl model as our highly defected structure would surely lead to increased scattering of phonons. Further study with simulated amorphous structures with more accurate coordination will be necessary. In the model of Feldman et al., 5 energy carriers in amorphous structures are divided into three characteristic types. Propagons travel ballistically with long mean-free-paths, diffusons transmit energy through an intrinsic diffusivity, and locons do not contribute to transport. The authors estimate the transitions between carrier types to occur at about 4.8 THz for the propagon-diffuson transition and about 17.4 THz for the diffuson-locon transition in a-si. These frequencies are in reasonable agreement with our findings, that the transitions occur at about 3.95 THz and THz, respectively. The agreement is not entirely surprising since Feldman et al. used the Stillinger-Weber potential in their calculations. The diffuson-locon transition appears to correlate with the acoustic to optical transition in the dispersion of about 17.4 THz for the (1) direction (Fig. 9). Disagreement between our results and those of experiment, especially with regard to the amorphous structure, result from two sources. The first is the preparation procedure used. It is possible that reducing the cooling rate will allow the structure to equilibrate to a lower energy state with four-fold coordination. The second is in the use of the interatomic potential. The Stillinger-Weber potential may not be able to recreate experimental a-si with the desired degree of fidelity, no matter what parameter set it used, and in this case, a different potential like the Tersoff 23,24 or density functional theory based Car-Parinello molecular dynamics 25,26,27 would more closely reproduce the desired system. The accuracy of the results may also be improved by adjusting the parameters used in the wavelet transform and performing the same study on independently generated structures or on the same structure for different orientations of the initial kinetic energy. Ultimately, fitting the results to a model of ballistic diffusive transport will allow the frequency resolved velocity and diffusivity (related to the mean-free-path) to be measured quantitatively. Summary In summary, we have presented a method by which the velocity and diffusivity of different vi- Baker 5
6 brational modes can be directly measured in simulated amorphous systems. Initial results indicate that the Cahill-Pohl model may achieve its accuracy by overestimating mode velocities and underestimating mean-free-paths. The data agreed qualitatively with the theory and findings of Feldman et al.. Further improvements to the study, through a more refined structure generation procedure and a more accurate description of interatomic forces, will enable quantitative comparison to theoretical treatment and experimental findings for transport in amorphous structures. Acknowledgements The authors would like to acknowledge financial support from the Air Force Office of Scientific Research (Grant No. FA ). C. H. B. is grateful for support from the Virginia Space Grant Consortium through a Graduate STEM Research Fellowship and the Department of Defense through a NDSEG Fellowship. References [1] P. G. Le Comber, Journal of Non-Crystalline Solids 115, 1 (1989). [2] M. Wuttig and N. Yamada, Naure Materials 6, 824 (27). [3] D. G. Cahill and R. O. Pohl, Solid State Communications 7, 927 (1989). [4] D. G. Cahill, S. K. Watson, and R. O. Pohl, Physical Review B 46, 6131 (1992). [5] J. L. Feldman, P. B. Allen, and S. R. Bickham, Physical Review B 59, 3551 (1999). [6] A. Aharony, S. Alexander, O. Entin-Wohlman, and R. Orbach, Physical Review Letters 58, 132 (1987). [7] A. Jagannathan, R. Orbach, and O. Entin- Wohlman, Physical Review B 39, (1989). [8] W. Kim, R. Wang, and A. Majumdar, Nano Today 2, 4 (27). [9] G. Kaiser, A Friendly Guide to Wavelets (Birkhauser, Boston, 1994). [1] C. H. Baker, D. A. Jordan, and P. M. Norris, Physical Review B 86, 1436 (212). [11] S. Plimpton, Journal of Computational Physics 117, 1 (1995), lammps.sandia.gov. [12] F. H. Stillinger and T. A. Weber, Physical Review B 31, 5262 (1985). [13] R. L. C. Vink, G. T. Barkema, W. F. van der Weg, and N. Mousseau, Journal of Non- Crystalline Solids 282, 248 (21). [14] M. Ishimaru, S. Munetoh, and T. Motooka, Physical Review B 56, (1997). [15] J. Fortner and J. S. Lannin, Physical Review B 39, 5527 (1989). [16] D. A. Drabold, P. A. Fedders, O. F. Sankey, and J. D. Dow, Physical Review B 42, 5135 (199). [17] K. Laaziri, S. Krycia, S. Roorda, M. Chicoine, J. L. Robertson, J. Wang, and S. C. Moss, Physical Review Letters 82, 346 (1999). [18] K. Laaziri, S. Krycia, S. Roorda, M. Chicoine, J. L. Robertson, J. Wang, and S. C. Moss, Physical Review B 6, 1352 (1999). [19] J. D. Gale and A. L. Rohl, Molecular Simulation 29, 291 (23). [2] W. A. Kamitakahara, H. R. Shanks, J. F. McClelland, U. Buchenau, F. Gompf, and L. Pintschovius, Physical Review Letters 52, 644 (1984). [21] W. A. Kamitakahara, C. M. Soukoulis, H. R. Shanks, U. Buchenau, and G. S. Grest, Physical Review B 36, 6539 (1987). [22] I. R. Cox-Smith, H. C. Liang, and R. O. Dillon, Journal of Vacuum Science and Technology A 3, 674 (1985). [23] J. Tersoff, Physical Review B 37, 6991 (1988). [24] J. Tersoff, Physical Review B 39, 5566 (1989). [25] R. Car and M. Parrinello, Physical Review Letters 55, 2471 (1985). [26] R. Car and M. Parrinello, Physical Review Letters 6, 24 (1988). [27] I. Štich, R. Car, and M. Parrinello, Physical Review B 44, 1192 (1991). Baker 6
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationRaman Spectra of Amorphous Silicon
Chapter 6 Raman Spectra of Amorphous Silicon In 1985, Beeman, Tsu and Thorpe established an almost linear relation between the Raman transverse-optic (TO) peak width Ɣ and the spread in mean bond angle
More informationResolving the Vibrational and Electronic Contributions to Thermal Conductivity of Silicon Near the Solid-Liquid Transition: Molecular Dynamics Study
I. J. Trans. Phenomena, Vol. 13, pp. 143 150 Reprints available directly from the publisher Photocopying permitted by license only 013 Old City Publishing, Inc. Published by license under the OCP Science
More informationBasic mechanisms of structural relaxation and diffusion in amorphous silicon
Basic mechanisms of structural relaxation and diffusion in amorphous silicon G. T. Barkema, Normand Mousseau, R.L.C. Vink and Parthapratim Biswas (1) Theoretical Physics, Utrecht University, Utrecht, The
More informationProceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition
Proceedings of the ASME 211 International Mechanical Engineering Congress & Exposition IMECE211 Proceedings of the ASME 211 International Mechanical November Engineering 11-17, 211, Congress Denver, &
More informationMolecular dynamics simulations of EXAFS in germanium
Cent. Eur. J. Phys. 93 2011 710-715 DOI: 10.2478/s11534-010-0074-0 Central European Journal of Physics Molecular dynamics simulations of EXAFS in germanium Research Article Janis Timoshenko Alexei Kuzmin
More informationSupporting Information. Influence of Vapor Deposition on Structural. and Charge Transport Properties of. Ethylbenzene Films
Supporting Information Influence of Vapor Deposition on Structural and Charge Transport Properties of Ethylbenzene Films Lucas W. Antony, Nicholas E. Jackson,, Ivan Lyubimov, Venkatram Vishwanath, Mark
More informationMicron School of Materials Science and Engineering. Problem Set 9 Solutions
Problem Set 9 Solutions 1. Mobility in extrinsic semiconductors is affected by phonon scattering and impurity scattering. Thoroughly explain the mobility plots for the following figures from your textbook
More informationThermal characterization of Au-Si multilayer using 3- omega method
Thermal characterization of Au-Si multilayer using 3- omega method Sunmi Shin Materials Science and Engineering Program Abstract As thermal management becomes a serious issue in applications of thermoelectrics,
More informationCoherent THz Noise Sources. T.M.Loftus Dr R.Donnan Dr T.Kreouzis Dr R.Dubrovka
Coherent THz Noise Sources T.M.Loftus Dr R.Donnan Dr T.Kreouzis Dr R.Dubrovka 1 Noise Source An unusual source Broadband Incoherent Lambertian emission Why consider it? 2 Power from various devices in
More informationComplex superlattice unit cell designs for reduced thermal conductivity. Abstract
Complex superlattice unit cell designs for reduced thermal conductivity E. S. Landry Department of Mechanical Engineering Carnegie Mellon University Pittsburgh, PA 15213 M. I. Hussein Department of Aerospace
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 8 Oct 1996
December 21, 2013 arxiv:cond-mat/9610066v1 [cond-mat.stat-mech] 8 Oct 1996 Some Finite Size Effects in Simulations of Glass Dynamics Jürgen Horbach, Walter Kob, Kurt Binder Institut für Physik, Johannes
More informationAtomic Scale Modeling of the Effect of Irradiation on Silica Optical Fibers
Atomic Scale Modeling of the Effect of Irradiation on Silica Optical Fibers Harish Govindarajan *, Materials Science Engineering, The Ohio State University Introduction Optical fibers and optically-based
More informationEnergy Density and Thermal Diffusivity of Ioffe-Regel Confined Vibrations
(carolinesgorham.com) Energy Density and Thermal Diffusivity of Ioffe-Regel Confined Vibrations Caroline S. Gorham MRS Fall Meeting Symposium ii8: Phonons in Nano and Bulk Materials Session Chair(s): J.
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 5: Specific Heat of Lattice Waves Outline Review Lecture 4 3-D Elastic Continuum 3-D Lattice Waves Lattice Density of Modes Specific Heat of Lattice Specific
More informationThe samples used in these calculations were arranged as perfect diamond crystal of
Chapter 5 Results 5.1 Hydrogen Diffusion 5.1.1 Computational Details The samples used in these calculations were arranged as perfect diamond crystal of a2 2 2 unit cells, i.e. 64 carbon atoms. The effect
More informationChapter 3 Properties of Nanostructures
Chapter 3 Properties of Nanostructures In Chapter 2, the reduction of the extent of a solid in one or more dimensions was shown to lead to a dramatic alteration of the overall behavior of the solids. Generally,
More informationComparisons of DFT-MD, TB- MD and classical MD calculations of radiation damage and plasmawallinteractions
CMS Comparisons of DFT-MD, TB- MD and classical MD calculations of radiation damage and plasmawallinteractions Kai Nordlund Department of Physics and Helsinki Institute of Physics University of Helsinki,
More informationElectro-Thermal Transport in Silicon and Carbon Nanotube Devices E. Pop, D. Mann, J. Rowlette, K. Goodson and H. Dai
Electro-Thermal Transport in Silicon and Carbon Nanotube Devices E. Pop, D. Mann, J. Rowlette, K. Goodson and H. Dai E. Pop, 1,2 D. Mann, 1 J. Rowlette, 2 K. Goodson 2 and H. Dai 1 Dept. of 1 Chemistry
More informationSize-dependent model for thin film and nanowire thermal conductivity
AIP/23-QED Size-dependent model for thin film and nanowire thermal conductivity Alan J. H. McGaughey,, a) Eric S. Landry,, 2 Daniel P. Sellan, 3 and Cristina H. Amon, 3 ) Department of Mechanical Engineering,
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationFig. 1: Raman spectra of graphite and graphene. N indicates the number of layers of graphene. Ref. [1]
Vibrational Properties of Graphene and Nanotubes: The Radial Breathing and High Energy Modes Presented for the Selected Topics Seminar by Pierce Munnelly 09/06/11 Supervised by Sebastian Heeg Abstract
More informationFYS Vår 2015 (Kondenserte fasers fysikk)
FYS410 - Vår 015 (Kondenserte fasers fysikk) http://www.uio.no/studier/emner/matnat/fys/fys410/v15/index.html Pensum: Introduction to Solid State Physics by Charles Kittel (Chapters 1-9 and 17, 18, 0)
More informationComplex superlattice unit cell designs for reduced thermal conductivity
Complex superlattice unit cell designs for reduced thermal conductivity E. S. Landry Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA M. I. Hussein
More informationMolecular dynamics simulations of coherent optical photon emission from shock waves in crystals
Molecular dynamics simulations of coherent optical photon emission from shock waves in crystals Evan J. Reed, 1,2, * Marin Soljačić, 2 Richard Gee, 1 and J. D. Joannopoulos 2 1 Lawrence Livermore National
More informationLecture 1. OUTLINE Basic Semiconductor Physics. Reading: Chapter 2.1. Semiconductors Intrinsic (undoped) silicon Doping Carrier concentrations
Lecture 1 OUTLINE Basic Semiconductor Physics Semiconductors Intrinsic (undoped) silicon Doping Carrier concentrations Reading: Chapter 2.1 EE105 Fall 2007 Lecture 1, Slide 1 What is a Semiconductor? Low
More informationMONTE CARLO MODELING OF HEAT GENERATION IN ELECTRONIC NANOSTRUCTURES
Proceedings of IMECE 2 22 ASME International Mechanical Engineering Congress and Exposition November 17-22, 22, New Orleans, Louisiana, USA IMECE2/HT-32124 MONTE CARLO MODELING OF HEAT GENERATION IN ELECTRONIC
More informationSECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS
2753 SECOND PUBLIC EXAMINATION Honour School of Physics Part C: 4 Year Course Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS TRINITY TERM 2011 Wednesday, 22 June, 9.30 am 12.30
More informationNanophysics: Main trends
Nano-opto-electronics Nanophysics: Main trends Nanomechanics Main issues Light interaction with small structures Molecules Nanoparticles (semiconductor and metallic) Microparticles Photonic crystals Nanoplasmonics
More informationQuantum Condensed Matter Physics Lecture 5
Quantum Condensed Matter Physics Lecture 5 detector sample X-ray source monochromator David Ritchie http://www.sp.phy.cam.ac.uk/drp2/home QCMP Lent/Easter 2019 5.1 Quantum Condensed Matter Physics 1. Classical
More informationClassical Theory of Harmonic Crystals
Classical Theory of Harmonic Crystals HARMONIC APPROXIMATION The Hamiltonian of the crystal is expressed in terms of the kinetic energies of atoms and the potential energy. In calculating the potential
More informationLaser Diodes. Revised: 3/14/14 14: , Henry Zmuda Set 6a Laser Diodes 1
Laser Diodes Revised: 3/14/14 14:03 2014, Henry Zmuda Set 6a Laser Diodes 1 Semiconductor Lasers The simplest laser of all. 2014, Henry Zmuda Set 6a Laser Diodes 2 Semiconductor Lasers 1. Homojunction
More informationSupporting Information
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2015 Supporting Information Single Layer Lead Iodide: Computational Exploration of Structural, Electronic
More informationDoctor of Philosophy
FEMTOSECOND TIME-DOMAIN SPECTROSCOPY AND NONLINEAR OPTICAL PROPERTIES OF IRON-PNICTIDE SUPERCONDUCTORS AND NANOSYSTEMS A Thesis Submitted for the degree of Doctor of Philosophy IN THE FACULTY OF SCIENCE
More informationNon-Fourier heat conduction in a single-walled carbon nanotube: Classical molecular dynamics simulations
Non-Fourier heat conduction in a single-walled carbon nanotube: Classical molecular dynamics simulations Junichiro Shiomi and Shigeo Maruyama * Department of Mechanical Engineering, The University of Tokyo
More informationX-ray, Neutron and e-beam scattering
X-ray, Neutron and e-beam scattering Introduction Why scattering? Diffraction basics Neutrons and x-rays Techniques Direct and reciprocal space Single crystals Powders CaFe 2 As 2 an example What is the
More informationHomework 2: Structure and Defects in Crystalline and Amorphous Materials
Homework 2: Structure and Defects in Crystalline and Amorphous Materials Muhammad Ashraful Alam Network of Computational Nanotechnology Discovery Park, Purdue University. In Lectures 4-6, we have discussed
More informationLecture 12: Phonon heat capacity
Lecture 12: Phonon heat capacity Review o Phonon dispersion relations o Quantum nature of waves in solids Phonon heat capacity o Normal mode enumeration o Density of states o Debye model Review By considering
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationPhysics with Neutrons I, WS 2015/2016. Lecture 11, MLZ is a cooperation between:
Physics with Neutrons I, WS 2015/2016 Lecture 11, 11.1.2016 MLZ is a cooperation between: Organization Exam (after winter term) Registration: via TUM-Online between 16.11.2015 15.1.2015 Email: sebastian.muehlbauer@frm2.tum.de
More informationAtomic Motion via Inelastic X-Ray Scattering
Atomic Motion via Inelastic X-Ray Scattering Cheiron School Beamline Practical - Tuesday ONLY at BL43LXU Alfred Q.R. Baron with H. Uchiyama We will introduce students to the use of inelastic x-ray scattering,
More informationPhonons I - Crystal Vibrations (Kittel Ch. 4)
Phonons I - Crystal Vibrations (Kittel Ch. 4) Displacements of Atoms Positions of atoms in their perfect lattice positions are given by: R 0 (n 1, n 2, n 3 ) = n 10 x + n 20 y + n 30 z For simplicity here
More informationNanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons
Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons Gang Chen Massachusetts Institute of Technology OXFORD UNIVERSITY PRESS 2005 Contents Foreword,
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2009 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 Spring 2009 POP QUIZ
More informationAtomic Motion via Inelastic X-Ray Scattering
Atomic Motion via Inelastic X-Ray Scattering Cheiron School Beamline Practical - Monday ONLY at BL35 Alfred Q.R. Baron & Satoshi Tsutsui We will introduce students to the use of inelastic x-ray scattering,
More informationSelf Formation of Porous Silicon Structure: Primary Microscopic Mechanism of Pore Separation
Solid State Phenomena Vols. 97-98 (2004) pp 181-184 (2004) Trans Tech Publications, Switzerland Journal doi:10.4028/www.scientific.net/ssp.97-98.181 Citation (to be inserted by the publisher) Copyright
More informationLecture 11: Periodic systems and Phonons
Lecture 11: Periodic systems and Phonons Aims: Mainly: Vibrations in a periodic solid Complete the discussion of the electron-gas Astrophysical electrons Degeneracy pressure White dwarf stars Compressibility/bulk
More informationOlivier Bourgeois Institut Néel
Olivier Bourgeois Institut Néel Outline Introduction: necessary concepts: phonons in low dimension, characteristic length Part 1: Transport and heat storage via phonons Specific heat and kinetic equation
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2014 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 POP QUIZ Phonon dispersion relation:
More informationMatter Waves. Chapter 5
Matter Waves Chapter 5 De Broglie pilot waves Electromagnetic waves are associated with quanta - particles called photons. Turning this fact on its head, Louis de Broglie guessed : Matter particles have
More informationIntroduction to Semiconductor Physics. Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India
Introduction to Semiconductor Physics 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India http://folk.uio.no/ravi/cmp2013 Review of Semiconductor Physics Semiconductor fundamentals
More informationPoint defects in pure amorphous silicon and their role in structural relaxation: A tight-binding molecular-dynamics study
Point defects in pure amorphous silicon and their role in structural relaxation: A tight-binding molecular-dynamics study Xavier Urli, Cristiano L. Dias,* Laurent J. Lewis, and Sjoerd Roorda Département
More informationAP* Chapter 10. Liquids and Solids. Friday, November 22, 13
AP* Chapter 10 Liquids and Solids AP Learning Objectives LO 1.11 The student can analyze data, based on periodicity and the properties of binary compounds, to identify patterns and generate hypotheses
More informationThe Dielectric Function of a Metal ( Jellium )
The Dielectric Function of a Metal ( Jellium ) Total reflection Plasma frequency p (10 15 Hz range) Why are Metals Shiny? An electric field cannot exist inside a metal, because metal electrons follow the
More informationCourse Name: AP Physics. Team Names: Jon Collins. Velocity Acceleration Displacement
Course Name: AP Physics Team Names: Jon Collins 1 st 9 weeks Objectives Vocabulary 1. NEWTONIAN MECHANICS and lab skills: Kinematics (including vectors, vector algebra, components of vectors, coordinate
More informationNon-Continuum Energy Transfer: Phonons
Non-Continuum Energy Transfer: Phonons D. B. Go Slide 1 The Crystal Lattice The crystal lattice is the organization of atoms and/or molecules in a solid simple cubic body-centered cubic hexagonal a NaCl
More informationAb initio molecular dynamics simulation on temperature-dependent properties of Al Si liquid alloy
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 16 (4) 57 514 PII: S953-8984(4)7691-8 Ab initio molecular dynamics simulation on temperature-dependent properties
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationChapter 37 Early Quantum Theory and Models of the Atom. Copyright 2009 Pearson Education, Inc.
Chapter 37 Early Quantum Theory and Models of the Atom Planck s Quantum Hypothesis; Blackbody Radiation Photon Theory of Light and the Photoelectric Effect Energy, Mass, and Momentum of a Photon Compton
More informationPhonon II Thermal Properties
Phonon II Thermal Properties Physics, UCF OUTLINES Phonon heat capacity Planck distribution Normal mode enumeration Density of states in one dimension Density of states in three dimension Debye Model for
More informationNon-equilibrium Green s functions: Rough interfaces in THz quantum cascade lasers
Non-equilibrium Green s functions: Rough interfaces in THz quantum cascade lasers Tillmann Kubis, Gerhard Klimeck Department of Electrical and Computer Engineering Purdue University, West Lafayette, Indiana
More informationdynamics computer code, whiccan simulate the motion of atoms at a surface and We coembininewton's equation for the nuclei,
DUI FILE COpy T DTIC ANNUAL LETTER REPORT BY 0. F. ELEC33 20KE AD-A224 069 P.1. FOR N00014-85-K-0442 U Theory of electronic states and formations energies of.defect comp exes, interstitial defects, and
More informationThe bulk modulus of covalent random networks
J. Phys.: Condens. Matter 9 (1997) 1983 1994. Printed in the UK PII: S0953-8984(97)77754-3 The bulk modulus of covalent random networks B R Djordjević and M F Thorpe Department of Physics and Astronomy
More informationAcoustic metamaterials in nanoscale
Acoustic metamaterials in nanoscale Dr. Ari Salmi www.helsinki.fi/yliopisto 12.2.2014 1 Revisit to resonances Matemaattis-luonnontieteellinen tiedekunta / Henkilön nimi / Esityksen nimi www.helsinki.fi/yliopisto
More informationUnderstanding Phonon Dynamics via 1D Atomic Chains
Understanding Phonon Dynamics via 1D Atomic Chains Timothy S. Fisher Purdue University School of Mechanical Engineering, and Birck Nanotechnology Center tsfisher@purdue.edu Nanotechnology 501 Lecture Series
More informationHomework Week 3: Nanoscale and macroscale characterization Thermoelectricity: From Atoms to Systems
Homework Week 3: Nanoscale and macroscale characterization Thermoelectricity: From Atoms to Systems Je-Hyeong Bahk and Ali Shakouri nanohub-u Fall 2013 Answer the thirteen questions including all the sub-questions
More informationTransient lattice dynamics in fs-laser-excited semiconductors probed by ultrafast x-ray diffraction
Transient lattice dynamics in fs-laser-excited semiconductors probed by ultrafast x-ray diffraction K. Sokolowski-Tinten, M. Horn von Hoegen, D. von der Linde Inst. for Laser- and Plasmaphysics, University
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationChapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice
Chapter 5 Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice In chapter 3 and 4, we have demonstrated that the deformed rods, rotational rods and perturbation
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationON THE ASSUMPTION OF DETAILED BALANCE IN PREDICTION OF DIFFUSIVE TRANSMISSION PROBABILITY DURING INTERFACIAL TRANSPORT
Nanoscale and Microscale Thermophysical Engineering, 14: 21 33, 2010 Copyright Ó Taylor & Francis Group, LLC ISSN: 1556-7265 print / 1556-7273 online DOI: 10.1080/15567260903530379 ON THE ASSUMPTION OF
More informationSelf-compensating incorporation of Mn in Ga 1 x Mn x As
Self-compensating incorporation of Mn in Ga 1 x Mn x As arxiv:cond-mat/0201131v1 [cond-mat.mtrl-sci] 9 Jan 2002 J. Mašek and F. Máca Institute of Physics, Academy of Sciences of the CR CZ-182 21 Praha
More informationSupporting Information: Heat-transport mechanisms in molecular building blocks of inorganic/organic hybrid superlattices
1 2 3 Supporting Information: Heat-transport mechanisms in molecular building blocks of inorganic/organic hybrid superlattices 4 5 6 7 8 9 Ashutosh Giri, 1 Janne-Petteri Niemelä, 2 Tommi Tynell, 2 John
More informationCalculation of thermal expansion coefficient of Fe3Al with the addition of transition metal elements Abstract Introduction Methodology
Calculation of thermal expansion coefficient of Fe Al with the addition of transition metal elements Tatiana Seletskaia, Leonid Muratov, Bernard Cooper, West Virginia University, Physics Departement, Hodges
More informationDevice-size atomistic models of amorphous silicon
PHYSICAL REVIEW B, VOLUME 64, 245214 Device-size atomistic models of amorphous silicon R. L. C. Vink* and G. T. Barkema Institute for Theoretical Physics, Utrecht University, Leuvenlaan 4, 3584 CE Utrecht,
More informationSIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS
Journal of Optoelectronics and Advanced Materials Vol. 7, No. 4, August 2005, p. 1915-1922 SIMULATIONAL ANALYSIS OF GLASSES PREPARED VIA DIFFERENT INTERATOMIC POTENTIALS Y. Sano, J. K ga, F. Yonezawa Department
More information5.62 Physical Chemistry II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 5.62 Physical Chemistry II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 5.62 Spring 2008 Lecture
More informationBranislav K. Nikolić
First-principles quantum transport modeling of thermoelectricity in nanowires and single-molecule nanojunctions Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark,
More informationDepartment of Physics, University of Maryland, College Park MIDTERM TEST
PHYSICS 731 Nov. 5, 2002 Department of Physics, University of Maryland, College Park Name: MIDTERM TEST Budget your time. Look at all 5 pages. Do the problems you find easiest first. 1. Consider a D-dimensional
More informationMonte Carlo Study of Thermal Transport of Direction and Frequency Dependent Boundaries in High Kn Systems
Monte Carlo Study of Thermal Transport of Direction and Frequency Dependent Boundaries in High Kn Systems N.A. Roberts and D.G. Walker Department of Mechanical Engineering Vanderbilt University May 30,
More informationSilica glass structure generation for ab initio calculations using small samples of amorphous silica
Silica glass structure generation for ab initio calculations using small samples of amorphous silica Renée M. Van Ginhoven* University of Washington, Seattle, Washington 98512, USA Hannes Jónsson University
More informationDominant structural defects in amorphous silicon
Journal of Physics: Condensed Matter PAPER Dominant structural defects in amorphous silicon To cite this article: Paule Dagenais et al 2015 J. Phys.: Condens. Matter 27 345004 View the article online for
More informationarxiv: v1 [cond-mat.mes-hall] 29 Apr 2017
Theoretical treatment of anharmonicity of vibrational modes of single-walled carbon nanotubes Hengjia Wang, 1 Doyl Dickel, 2 and Murray S. Daw 1 1 Department of Physics and Astronomy, arxiv:1705.00201v1
More informationESE 372 / Spring 2013 / Lecture 5 Metal Oxide Semiconductor Field Effect Transistor
Metal Oxide Semiconductor Field Effect Transistor V G V G 1 Metal Oxide Semiconductor Field Effect Transistor We will need to understand how this current flows through Si What is electric current? 2 Back
More informationB 2 P 2, which implies that g B should be
Enhanced Summary of G.P. Agrawal Nonlinear Fiber Optics (3rd ed) Chapter 9 on SBS Stimulated Brillouin scattering is a nonlinear three-wave interaction between a forward-going laser pump beam P, a forward-going
More informationPhonons (Classical theory)
Phonons (Classical theory) (Read Kittel ch. 4) Classical theory. Consider propagation of elastic waves in cubic crystal, along [00], [0], or [] directions. Entire plane vibrates in phase in these directions
More informationNonlinear Electrodynamics and Optics of Graphene
Nonlinear Electrodynamics and Optics of Graphene S. A. Mikhailov and N. A. Savostianova University of Augsburg, Institute of Physics, Universitätsstr. 1, 86159 Augsburg, Germany E-mail: sergey.mikhailov@physik.uni-augsburg.de
More informationStructural properties of low-density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. 65, No. 6 journal of December 2005 physics pp. 1085 1096 Structural properties of low-density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2
More informationPart II. Fundamentals of X-ray Absorption Fine Structure: data analysis
Part II Fundamentals of X-ray Absorption Fine Structure: data analysis Sakura Pascarelli European Synchrotron Radiation Facility, Grenoble, France Page 1 S. Pascarelli HERCULES 2016 Data Analysis: EXAFS
More informationTopic 5-1: Introduction to Phonons Kittel pages: 91, 92
Topic 5-1: Introduction to Phonons Kittel pages: 91, 92 Summary: In this video we introduce the concept that atoms are not rigid, fixed points within the lattice. Instead we treat them as quantum harmonic
More informationChapter 6. Transport in Ag 2 Se
86 Chapter 6 Transport in Ag 2 Se This chapter is principally about the thermoelectric performance of Ag 2 Se near its phase transition. In the sample studied here the Hall carrier concentration does not
More informationPlasmonics. The long wavelength of light ( μm) creates a problem for extending optoelectronics into the nanometer regime.
Plasmonics The long wavelength of light ( μm) creates a problem for extending optoelectronics into the nanometer regime. A possible way out is the conversion of light into plasmons. They have much shorter
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationAn Introduction to Lattice Vibrations
An Introduction to Lattice Vibrations Andreas Wacker 1 Mathematical Physics, Lund University November 3, 2015 1 Introduction Ideally, the atoms in a crystal are positioned in a regular manner following
More informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationSpatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media Karim G. Sabra, School of Mechanical Engineering,
More informationReview of Semiconductor Physics
Solid-state physics Review of Semiconductor Physics The daunting task of solid state physics Quantum mechanics gives us the fundamental equation The equation is only analytically solvable for a handful
More informationStructural properties of low density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. xx, No. x journal of xxxxxx 2005 physics pp. 1 12 Structural properties of low density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2 1 The
More information