Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Energy separation of neutrons scattered at small angles from silicon using time-of-flight techniques

The time-of-flight technique is used on a small-angle neutron scattering instrument to separate the energies of the scattered neutrons, in order to determine the origin of the temperature-dependent scattering observed from silicon at Q > similar to 0.1 angstrom(-1). A quantitative analysis of the results in comparison with the phonon dispersion curves, determined by Dolling using a triple-axis neutron spectrometer, shows that the temperature-dependent scattering can be understood in terms of Umklapp processes whereby neutrons gain energy from phonons.

The vibrational spectrum and ultimate modulus of polyethylene

We have performed the first completely ab initio lattice dynamics calculation of the full orthorhombic cell of polyethylene using periodic density functional theory in the local density approximation (LDA) and the generalized gradient approximation (GGA). Contrary to current perceptions, we show that LDA accurately describes the structure whereas GGA fails. We emphasize that there is no parametrization of the results. We then rigorously tested our calculation by computing the phonon dispersion curves across the entire Brillouin zone and comparing them to the vibrational spectra, in particular the inelastic neutron scattering (INS) spectra, of polyethylene (both polycrystalline and aligned) and perdeuteriopolyethylene. The F-point frequencies (where the infrared and Raman active modes occur) are in good agreement with the latest low temperature data. The near-perfect reproduction of the INS spectra, gives confidence in the results and allows Lis to deduce a number of physical properties including the elastic moduli, parallel and perpendicular to the chain. We find that the Young's modulus for an infinitely long, perfectly crystalline polyethylene is 360.2 GPa at 0 K. The highest experimental value is 324 GPa, indicating that current high modulus fibers are similar to 90% of their maximum possible strength.

Lattice dynamics study of microstructures of electrorheological fluids

The lattice dynamics method is used to study the stability of the chain structures formed in electrorheological (ER) fluids. The appearance of the soft modes in the phonon dispersion of the structures indicates that the chains tend to distort and aggregate into thicker columns due to the electrostatic attractive forces and thermal generated forces between them. The results show that the stability of the chains relies on their width and the separation between them. The complete chain structures are more stable than the chains with defects. The results can be used to elucidate the densification phenomenon of the chains in the structuring process of ER fluids in the quiescent state.

Electronic, vibrational and related properties of group IV metal oxides by ab initio calculations

We present our theoretical results for the structural, electronic, vibrational and optical properties of MO(2) (M = Sn, Zr, Hf and Ti) obtained by first-principles calculations. Relativistic effects are demonstrated to be important for a realistic description of the detailed structure of the electronic frequency-dependent dielectric function, as well as of the carrier effective masses. Based on our results, we found that the main contribution of the high values calculated for the oxides dielectric constants arises from the vibrational properties of these oxides, and the vibrational static dielectric constant values diminish with increasing pressure. (c) 2008 Elsevier B.V. All rights reserved.

Confinamento de fónons ópticos em estruturas piezoelétricas periódicas e quasiperiódicas

We study the optical-phonon spectra in periodic and quasiperiodic (Fibonacci type) superlattices made up from III-V nitride materials (GaN and AlN) intercalated by a dielectric material (silica - SiO2). Due to the misalignments between the silica and the GaN, AlN layers that can lead to threading dislocation of densities as high as 1010 cm−1, and a significant lattice mismatch (_ 14%), the phonon dynamics is described by a coupled elastic and electromagnetic equations beyond the continuum dielectric model, stressing the importance of the piezoelectric polarization field in a strained condition. We use a transfer-matrix treatment to simplify the algebra, which would be otherwise quite complicated, allowing a neat analytical expressions for the phonon dispersion relation. Furthermore, a quantitative analysis of the localization and magnitude of the allowed band widths in the optical phonon s spectra, as well as their scale law are presented and discussed

Ashcroft model potential study of lattice dynamics of α-iron and barium

Ashcroft model potential has been used to compute phonon dispersion relations along the three principal symmetry directions, i.e. [k00], [kk0] and [kkk] for alpha-iron and barium. The computed phonons gave a reasonable agreement with the experimental ones in all the three principal summetry directions expect for the T-2 branch in [KK0] direction where the present study failed to reproduce the experimental findings.

LATTICE-DYNAMICS OF ALKALI-METALS IN A 3-BODY INTERACTION

The original model of Das et al. is modified in extending the electron-ion interaction on a three-body forces and including the crystal equilibrium condition to reduce one independent parameter. We studied the phonon dispersion relations along the three principal symmetry directions i.e. [xi, 0, 0], [xi, xi, 0] and [xi, xi, xi] and theta-T curves of alkali metals, Na, K, Rb, Cs and Li. There is close agreement between the computed results and the experimental observations.

Lattice dynamics of noble metals on effective three-body interaction

Quite recently we modified the original model of Sarkar et al. for cubic metals in extending the ion-ion interaction, ion-electron interaction and the introduction of crystal equilibrium condition. We applied our scheme to alkali metals. We studied here the lattice dynamics of noble metals on our approach by calculating phonon dispersion relations along the three principal symmetry directions, [ξ00], [ξξ00] and [ξξξ] and the (θ-T) curves of three noble metals: copper, silver and gold. We obtained reasonable agreement with the experimental findings.

Simulative Investigation on the Electronic, Vibrational and Optical Properties of the Ge2Sb2Te5 Chalcogenide

Chalcogenides are chemical compounds with at least one of the following three chemical elements: Sulfur (S), Selenium (Sn), and Tellurium (Te). As opposed to other materials, chalcogenide atomic arrangement can quickly and reversibly inter-change between crystalline, amorphous and liquid phases. Therefore they are also called phase change materials. As a results, chalcogenide thermal, optical, structural, electronic, electrical properties change pronouncedly and significantly with the phase they are in, leading to a host of different applications in different areas. The noticeable optical reflectivity difference between crystalline and amorphous phases has allowed optical storage devices to be made. Their very high thermal conductivity and heat fusion provided remarkable benefits in the frame of thermal energy storage for heating and cooling in residential and commercial buildings. The outstanding resistivity difference between crystalline and amorphous phases led to a significant improvement of solid state storage devices from the power consumption to the re-writability to say nothing of the shrinkability. This work focuses on a better understanding from a simulative stand point of the electronic, vibrational and optical properties for the crystalline phases (hexagonal and faced-centered cubic). The electronic properties are calculated implementing the density functional theory combined with pseudo-potentials, plane waves and the local density approximation. The phonon properties are computed using the density functional perturbation theory. The phonon dispersion and spectrum are calculated using the density functional perturbation theory. As it relates to the optical constants, the real part dielectric function is calculated through the Drude-Lorentz expression. The imaginary part results from the real part through the Kramers-Kronig transformation. The refractive index, the extinctive and absorption coefficients are analytically calculated from the dielectric function. The transmission and reflection coefficients are calculated using the Fresnel equations. All calculated optical constants compare well the experimental ones.

New path integral simulation algorithms and their application to creep in the quantum sine-Gordon chain

A path integral simulation algorithm which includes a higher-order Trotter approximation (HOA)is analyzed and compared to an approach which includes the correct quantum mechanical pair interaction (effective Propagator (EPr)). It is found that the HOA algorithmconverges to the quantum limit with increasing Trotter number P as P^{-4}, while the EPr algorithm converges as P^{-2}.The convergence rate of the HOA algorithm is analyzed for various physical systemssuch as a harmonic chain,a particle in a double-well potential, gaseous argon, gaseous helium and crystalline argon. A new expression for the estimator for the pair correlation function in the HOA algorithm is derived. A new path integral algorithm, the hybrid algorithm, is developed.It combines an exact treatment of the quadratic part of the Hamiltonian and thehigher-order Trotter expansion techniques.For the discrete quantum sine-Gordon chain (DQSGC), it is shown that this algorithm works more efficiently than all other improved path integral algorithms discussed in this work. The new simulation techniques developed in this work allow the analysis of theDQSGC and disordered model systems in the highly quantum mechanical regime using path integral molecular dynamics (PIMD)and adiabatic centroid path integral molecular dynamics (ACPIMD).The ground state phonon dispersion relation is calculated for the DQSGC by the ACPIMD method.It is found that the excitation gap at zero wave vector is reduced by quantum fluctuations. Two different phases exist: One phase with a finite excitation gap at zero wave vector, and a gapless phase where the excitation gap vanishes.The reaction of the DQSGC to an external driving force is analyzed at T=0.In the gapless phase the system creeps if a small force is applied, and in the phase with a gap the system is pinned. At a critical force, the systems undergo a depinning transition in both phases and flow is induced. The analysis of the DQSGC is extended to models with disordered substrate potentials. Three different cases are analyzed: Disordered substrate potentials with roughness exponent H=0, H=1/2,and a model with disordered bond length. For all models, the ground state phonon dispersion relation is calculated.

Complex materials via colloidal crystallization

The work presented in this thesis deals with complex materials, which were obtained by self-assembly of monodisperse colloidal particles, also called colloidal crystallization. Two main fields of interest were investigated, the first dealing with the fabrication of colloidal monolayers and nanostructures, which derive there from. The second turned the focus on the phononic properties of colloidal particles, crystals, and glasses. For the fabrication of colloidal monolayers a method is introduced, which is based on the sparse distribution of dry colloidal particles on a parent substrate. In the ensuing floating step the colloidal monolayer assembles readily at the three-phase-contact line, giving a 2D hexagonally ordered film under the right conditions. The unique feature of this fabrication process is an anisotropic shrinkage, which occurs alongside with the floating step. This phenomenon is exploited for the tailored structuring of colloidal monolayers, leading to designed hetero-monolayers by inkjet printing. Furthermore, the mechanical stability of the floating monolayers allows the deposition on hydrophobic substrates, which enables the fabrication of ultraflat nanostructured surfaces. Densely packed arrays of crescent shaped nanoparticles have also been synthesized. It is possible to stack those arrays in a 3D manner allowing to mutually orientate the individual layers. In a step towards 3D mesoporous materials a methodology to synthesize hierarchically structured inverse opals is introduced. The deposition of colloidal particles in the free voids of a host inverse opal allows for the fabrication of composite inverse opals on two length scales. The phononic properties of colloidal crystals and films are characterized by Brillouin light scattering (BLS). At first the resonant modes of colloidal particles consisting of polystyrene, a copolymer of methylmethacrylate and butylacrylate, or of a silica core-PMMA shell topography are investigated, giving insight into their individual mechanical properties. The infiltration of colloidal films with an index matching liquid allows measuring the phonon dispersion relation. This leads to the assignment of band gaps to the material under investigation. Here, two band gaps could be found, one originating from the fcc order in the colloidal crystal (Bragg gap), the other stemming from the vibrational eigenmodes of the colloidal particles (hybridization gap).

Coupling light into graphene plasmons with the help of surface acoustic waves

We propose a scheme for coupling laser light into graphene plasmons with the help of electrically generated surface acoustic waves. The surface acoustic wave forms a diffraction grating which allows us to excite the long lived phononlike branch of the hybridized graphene plasmon-phonon dispersion with infrared laser light. Our approach avoids patterning the graphene sheet, does not rely on complicated optical near-field techniques, and allows us to electrically switch the coupling between far-field radiation and propagating graphene plasmons.

Prediction of a superhard carbon-rich C-N compound comparable to diamond

Until now, it has been a challenge both in experiment and in theory to design new superhard materials with high hardness values that are comparable to that of diamond. Here, by using first-principles calculations, we have introduced two new phases for a carbon-rich C-N compound with stoichiometry C3N, which is predicted to be energetically stable or metastable with respect to graphite and solid N2 at ambient pressure. It is found that C3N has a layered structure containing graphitic layers sandwiched with freely rotated N2 molecules. The layer-structured C3N is calculated to transform into a three-dimensional C2221 structure at 9 GPa with sp3-hybridized C atoms and sp2-hybridized N atoms. Phonon dispersion and elastic constant calculations reveal the dynamical and mechanical stability of the C2221 phase of C3N at ambient pressure. Significantly, first-principles ideal strength calculations indicate that the C2221 phase of C3N is a superhard material with an estimated Vickers hardness (∼76 GPa) comparable to that of diamond (60-120 GPa). The present results shed strong light on designing new superhard materials in the C-N system.

Single layer lead iodide : computational exploration of structural, electronic and optical properties, strain induced band modulation and the role of spin-orbital-coupling

Graphitic like layered materials exhibit intriguing electronic structures and thus the search for new types of two-dimensional (2D) monolayer materials is of great interest for developing novel nano-devices. By using density functional theory (DFT) method, here we for the first time investigate the structure, stability, electronic and optical properties of monolayer lead iodide (PbI2). The stability of PbI2 monolayer is first confirmed by phonon dispersion calculation. Compared to the calculation using generalized gradient approximation, screened hybrid functional and spin-orbit coupling effects can not only predicts an accurate bandgap (2.63 eV), but also the correct position of valence and conduction band edges. The biaxial strain can tune its bandgap size in a wide range from 1 eV to 3 eV, which can be understood by the strain induced uniformly change of electric field between Pb and I atomic layer. The calculated imaginary part of the dielectric function of 2D graphene/PbI2 van der Waals type hetero-structure shows significant red shift of absorption edge compared to that of a pure monolayer PbI2. Our findings highlight a new interesting 2D material with potential applications in nanoelectronics and optoelectronics.