Lattice dynamics and Raman spectra of isotopically mixed diamond

Isotope Effect in Thermal Conductivity of Polycrystalline CVD-Diamond: Experiment and Theory

We measured the thermal conductivity κ(T) of polycrystalline diamond with natural (natC) and isotopically enriched (12C content up to 99.96 at.%) compositions over a broad temperature T range, from 5 to 410 K. The high quality polycrystalline diamond wafers were produced by microwave plasma chemical vapor deposition in CH4-H2 mixtures. The thermal conductivity of 12C diamond along the wafer, as precisely determined using a steady-state longitudinal heat flow method, exceeds much that of the natC sample at T>60 K. The enriched sample demonstrates the value of κ(298K)=25.1±0.5 W cm−1 K−1 that is higher than the ever reported conductivity of natural and synthetic single crystalline diamonds with natural isotopic composition. A phenomenological theoretical model based on the full version of Callaway theory of thermal conductivity is developed which provides a good approximation of the experimental data. The role of different resistive scattering processes, including due to minor isotope 13C atoms, defects, and grain boundaries, is estimated from the data analysis. The model predicts about a 37% increase of thermal conductivity for impurity and dislocation free polycrystalline chemical vapor deposition (CVD)-diamond with the 12C-enriched isotopic composition at room temperature.

Isotopic Disorder: The Prevailing Mechanism in Limiting the Phonon Lifetime in Hexagonal BN

The phonon linewidth of isotopically controlled hexagonal boron nitride (h-BN) single crystals has been determined by Raman scattering. The scattering by isotopic mass disorder induces a phonon broadening that is largest for boron 11 fractions around 0.65. Lowest-order perturbation theory does not suffice to explain the dependence of the isotopic broadening on isotopic composition. A multiple-scattering theory based on the coherent potential approximation provides a good quantitative account of the phonon shift and broadening with isotopic composition observed in the experiments. Isotopic-disorder scattering is shown to have a prominent role in limiting the optical-phonon lifetime in h-BN.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Isotope engineering of van der Waals interactions in hexagonal boron nitride

Hexagonal boron nitride is a model lamellar compound where weak, non-local van der Waals interactions ensure the vertical stacking of two-dimensional honeycomb lattices made of strongly bound boron and nitrogen atoms. We study the isotope engineering of lamellar compounds by synthesizing hexagonal boron nitride crystals with nearly pure boron isotopes (¹⁰B and ¹¹B) compared to those with the natural distribution of boron (20 at% ¹⁰B and 80 at% ¹¹B). On the one hand, as with standard semiconductors, both the phonon energy and electronic bandgap varied with the boron isotope mass, the latter due to the quantum effect of zero-point renormalization. On the other hand, temperature-dependent experiments focusing on the shear and breathing motions of adjacent layers revealed the specificity of isotope engineering in a layered material, with a modification of the van der Waals interactions upon isotope purification. The electron density distribution is more diffuse between adjacent layers in ¹⁰BN than in ¹¹BN crystals. Our results open perspectives in understanding and controlling van der Waals bonding in layered materials.

Vibrational monitor of early demineralization in tooth enamel after in vitro exposure to phosphoridic liquid

The Raman spectroscopic method has been applied to quantitatively assess the in vitro degree of demineralization in healthy human teeth. Based on previous evaluations of Raman selection rules (empowered by an orientation distribution function (ODF) statistical algorithm) and on a newly proposed analysis of phonon density of states (PDOS) for selected vibrational modes of the hexagonal structure of hydroxyapatite, a molecular-scale evaluation of the demineralization process upon in vitro exposure to a highly acidic beverage (i.e., CocaCola™ Classic, pH=2.5) could be obtained. The Raman method proved quite sensitive and spectroscopic features could be directly related to an increase in off-stoichiometry of the enamel surface structure since the very early stage of the demineralization process (i.e., when yet invisible to other conventional analytical techniques). The proposed Raman spectroscopic algorithm might possess some generality for caries risk assessment, allowing a prompt non-contact diagnostic practice in dentistry.

Chemical dynamics simulations of energy transfer, surface-induced dissociation, soft-landing, and reactive-landing in collisions of protonated peptide ions with organic surfaces

Different simulation approaches like MM, QM + MM, and QM/MM, were used to study surface-induced dissociation, soft-landing, and reactive-landing for the peptide-H⁺ + surface collisions.

Importance of shattering fragmentation in the surface-induced dissociation of protonated octaglycine

A QM + MM direct chemical dynamics simulation was performed to study collisions of protonated octaglycine, gly(8)-H(+), with the diamond {111} surface at an initial collision energy E(i) of 100 eV and incident angle theta(i) of 0 degrees and 45 degrees. The semiempirical model AM1 was used for the gly(8)-H(+) intramolecular potential, so that its fragmentation could be studied. Shattering dominates gly(8)-H(+) fragmentation at theta(i) = 0 degrees, with 78% of the ions dissociating in this way. At theta(i) = 45 degrees shattering is much less important. For theta(i) = 0 degrees there are 304 different pathways, many related by their backbone cleavage patterns. For the theta(i) = 0 degrees fragmentations, 59% resulted from both a-x and b-y cleavages, while for theta(i) = 45 degrees 70% of the fragmentations occurred with only a-x cleavage. For theta(i) = 0 degrees, the average percentage energy transfers to the internal degrees of freedom of the ion and the surface, and the energy remaining in ion translation are 45%, 26%, and 29%. For 45 degrees these percentages are 26%, 12%, and 62%. The percentage energy-transfer to DeltaE(int) for theta(i) = 0 degrees is larger than that reported in previous experiments for collisions of des-Arg(1)-bradykinin with a diamond surface at the same theta(i). This difference is discussed in terms of differences between the model diamond surface used in the simulations and the diamond surface prepared for the experiments.

Energy calibration of a high-resolution inelastic x-ray scattering spectrometer

The energy scale of a triple-axis x-ray spectrometer with meV energy resolution based on perfect silicon crystal optics is calibrated, utilizing the most recent determination of the silicon lattice parameter and its thermal expansion coefficient and recording the dispersion of longitudinal acoustic and optical phonons in a diamond single crystal and the molecular vibration mode in liquid nitrogen. Comparison of the x-ray results with previous inelastic neutron and Raman scattering results as well as with ab initio phonon dispersion calculations yields an overall agreement better than 2%.

Role of Projectile and Surface Temperatures in the Energy Transfer Dynamics of Protonated Peptide Ion Collisions with the Diamond {111} Surface

The effects of temperature on energy transfer during collisions of protonated diglycine ions, Gly(2)-H(+), with a diamond {111} surface were investigated by chemical dynamics simulations. The simulations were performed for a collision energy of 70 eV and angle of 0 degrees with respect to the surface normal. In one set of simulations the initial surface temperature, T(surf), was varied from 300 to 2000 K, while the Gly(2)-H(+) vibrational and rotational temperatures were maintained at 300 K. For the second set of simulations the Gly(2)-H(+) vibrational temperature, T(vib), was varied from 300 to 2000 K, keeping both the Gly(2)-H(+) rotational and surface temperatures at 300 K. Increasing either the surface temperature or Gly(2)-H(+) vibrational temperature to values as high as 2000 K has, at most, only a negligible effect on the partitioning of the incident collision energy to the surface and to the vibrational and rotational modes of Gly(2)-H(+). To a good approximation, the initial surface and peptide ion energies are nearly adiabatic during the collisional energy transfer. This adiabaticity of the initial peptide ion energy agrees with experiments (J. Phys. Chem. A 2004, 108, 1). A more quantitative analysis of the effects of T(vib) and T(surf) shows there are small, but noticeable, effects on the energy transfer efficiencies. Namely, increasing the vibrational or surface temperature results in a near-linear decrease in the energy transfer to the degrees of freedom associated with this temperature.

Dynamics and Kinetics of Heat Transfer at the Interface of Model Diamond {111} Nanosurfaces

A molecular dynamics simulation was performed to study the effect of an applied force on heat transfer at the interface of model diamond [111] nanosurfaces. The force was applied to a small, hot nanosurface at 800, 1000, or 1200 K brought into contact with a larger, colder nanosurface at 300 K. The relaxation of the initial nonequilibrium interfacial force occurs on a subpicosecond time scale, much shorter than that required for heat transfer. Heat transfer occurs with exponential kinetics and a rate constant that increases linearly with the interfacial force according to 7 x 10(-4) ps(-1)/nN. This rate constant only increases by at most 10% as the temperature of the hot surface is increased from 800 to 1200 K. Replacing the interfacial H-atoms on both surfaces by D atoms also has a very small effect on the heat transfer. However, if one nanosurface has H atoms on its interface and the other nanosurface's interface has D atoms, then there is a marked 25% decrease in the rate constant for heat transfer. Increasing the size of the hot surface, and, thus, the interfacial contact area, increases the rate of heat transfer but not the rate constant. For the same interfacial force, different anharmonic models for the nanosurfaces' potential energy function give the same heat transfer rate constant. The possibility of quantum effects for heat transfer across the diamond interface is considered.

Direct dynamics study of N-protonated diglycine surface-induced dissociation. Influence of collision energy

A quantum mechanical and molecular mechanical (QM + MM) direct dynamics classical trajectory simulation is used to study energy transfer and fragmentation in the surface-induced dissociation (SID) of N-protonated diglycine, (gly)2H+. The peptide ion collides with the hydrogenated diamond [111] surface. The Austin Model 1 (AM1) semiempirical electronic structure theory is used for the (gly)2H+ intramolecular potential and molecular mechanical functions are used for the diamond surface potential and peptide/surface intermolecular potential. The simulations are performed at collision energies Ei of 30, 50, 70, and 100 eV and collision angle of 0 degrees (perpendicular to the surface). The percent energy transfer to the peptide ion is nearly independent of Ei, while energy transfer to the surface increases with increase in Ei. A smaller percent of the energy remains in peptide translation as Ei is increased. These trends in energy transfer are consistent with previous trajectory simulations of SID. At each Ei the most likely initial pathway leading to fragmentation is rupture of the +H3NCH2-CONHCH2COOH bond. Fragmentation occurs by two general mechanisms. One is the traditional Rice-Ramsperger-Kassel-Marcus (RRKM) model in which the peptide ion is activated by its collision with the surface, "bounces off", and then dissociates after undergoing intramolecular vibrational energy redistribution (IVR). The other mechanism is shattering in which the ion fragments as it collides with the surface. Shattering is the origin of the large increase in number of product channels with increase in Ei, i.e., 6 at 30 eV, but 59 at 100 eV. Shattering becomes the dominant dissociation mechanism at high Ei.

Dynamics of energy transfer in peptide-surface collisions

Classical trajectory simulations are performed to study energy transfer in collisions of protonated triglycine (Gly)(3) and pentaglycine (Gly)(5) ions with n-hexyl thiolate self-assembled monolayer (SAM) and diamond [111] surfaces, for a collision energy E(i) in the range of 10-110 eV and a collision angle of 45 degrees. Energy transfer to the peptide ions' internal degrees of freedom is more efficient for collision with the diamond surface; i.e., 20% transfer to peptide vibration/rotation at E(i) = 30 eV. For collision with diamond, the majority of E(i) remains in peptide translation, while the majority of the energy transfer is to surface vibrations for collision with the softer SAM surface. The energy-transfer efficiencies are very similar for (Gly)(3) and (Gly)(5). Constraining various modes of (Gly)(3) shows that the peptide torsional modes absorb approximately 80% of the energy transfer to the peptide's internal modes. The energy-transfer efficiencies depend on E(i). These simulations are compared with recent experiments of peptide SID and simulations of energy transfer in Cr(CO)(6)(+) collisions with the SAM and diamond surfaces.

Photoluminescence of isotopically purified silicon: how sharp are bound exciton transitions?

We report the first high resolution photoluminescence studies of isotopically pure Si (99.896% (28)Si). New information is obtained on isotopic effects on the indirect band gap energy, phonon energies, and phonon broadenings, which is in good agreement with calculations and previous results obtained in Ge and diamond. Remarkably, the linewidths of the no-phonon boron and phosphorus bound exciton transitions in the (28)Si sample are much narrower than in natural Si and are not well resolved at our maximum instrumental resolution of approximately 0.014 cm(-1). The removal of the dominant broadening resulting from isotopic randomness in natural Si reveals new fine structure in the boron bound exciton luminescence.

Isotope-disorder-induced line broadening of phonons in the Raman spectra of SiC

The width of phonon lines in the Raman spectra of ideal isotopically pure solids is determined by inelastic scattering processes. In solids that contain a mixture of different isotopes of one atomic constituent, elastic scattering due to isotopic mass disorder opens up decay channels that result in additional line broadening. We use different polytypes of SiC with an associated number of Raman active modes in order to experimentally validate the proportionality between linewidth and phonon density of states predicted by a simple elastic scattering theory.

The thermal conductivity of CVD diamond films

The thermal conductivity of chemical vapour deposition diamond films is controlled by the microstructure, impurity content and carbon double bonds in the films. In high conductivity films, dislocation scattering is dominant at low temperatures, while phonon-phonon scattering limits the conductivity at room temperature. In lower quality films, hydrogen and metal impurities as well as carbon double bonds constrain the conductivity up to room temperature. Significant anisotropies and gradients in the thermal conductivity exist in some films because of their micro structure.