Skip to content

Thermal Transport Properties of Diamond Phonons by Electric Field

MetadataDetails
Publication Date2022-09-28
JournalNanomaterials
AuthorsYongsheng Zhao, Fengyun Yan, Xue Liu, Hongfeng Ma, Zhenyu Zhang
InstitutionsLanzhou University of Technology
Citations3
AnalysisFull AI Review Included

Technical Documentation: Electric Field Modulation of Diamond Thermal Transport

Section titled “Technical Documentation: Electric Field Modulation of Diamond Thermal Transport”

This document analyzes the research on controlling diamond’s thermal transport properties via external electric fields, providing technical specifications and outlining how 6CCVD’s advanced MPCVD diamond materials and customization services enable the replication and extension of this groundbreaking work for next-generation thermal management applications.


The research demonstrates the feasibility of actively tuning the lattice thermal conductivity ($\kappa$) of diamond using external electric fields, a critical breakthrough for high-power electronics and quantum technology.

  • Active Thermal Modulation: External electric fields (E-fields) can significantly modulate diamond’s thermal conductivity, increasing it by up to 36.1% or decreasing it by 34% at 300 K.
  • Anisotropic Response: The thermal response is highly dependent on crystal orientation, with the [111] direction showing the most pronounced modulation effect.
  • Mechanism Confirmed: The E-field breaks the diamond lattice symmetry, altering the electron density distribution, which subsequently changes the phonon mean free path (MFP) and phonon-phonon interaction strength, primarily affecting phonons below 30 THz.
  • High-Power Relevance: This tunability is essential for developing advanced diamond heat sinks capable of dynamic thermal management in high-power Insulated Gate Bipolar Transistor (IGBT) equipment and integrated circuits.
  • Material Requirement: Achieving these theoretical thermal properties requires ultra-high purity, low-defect Single Crystal Diamond (SCD) substrates, which 6CCVD specializes in producing via MPCVD.

The following table summarizes the key quantitative findings regarding diamond’s thermal response to electric fields, extracted from the first-principles calculations (T = 300 K).

ParameterValueUnitContext
Benchmark Lattice Thermal Conductivity ($\kappa$)1960.80W·m-1K-1Calculated at 0 a.u. E-field
Maximum $\kappa$ Achieved2654W·m-1K-1[111] orientation, E = +0.004 a.u.
Minimum $\kappa$ Achieved1283W·m-1K-1[111] orientation, E = -0.004 a.u.
Maximum $\kappa$ Increase Rate ($\Delta \kappa / \kappa$)36.1%[111] orientation, E = +0.004 a.u.
Maximum $\kappa$ Decrease Rate ($\Delta \kappa / \kappa$)-34%[111] orientation, E = -0.004 a.u.
Critical E-Field Strength (Microscopic)±0.004a.u.Used for maximum modulation
Critical E-Field Strength (Macroscopic)$\approx \pm 2.06 \times 10^{9}$V·cm-1Applied along [111] direction
Temperature Range Studied240 - 500KThermal conductivity analysis range
Phonon Frequency Range Affected< 30THzRange showing primary E-field response
C-C Bond Length (Fixed)1.547ÅUnchanged by E-field application

The study utilized advanced first-principles calculations to model the electron-phonon coupling and thermal transport response of diamond under a uniform electric field.

  1. Computational Framework: First-principles calculations were performed using the ABINIT software package, leveraging Density-Functional Perturbation Theory (DFPT) to analyze phonon properties.
  2. Thermal Transport Calculation: Lattice thermal conductivity ($\kappa$) was calculated using the Phono3py package, which employs the finite displacement method to determine second-order and third-order force constants.
  3. Transport Equation Solution: The Boltzmann phonon transport equations (BTE) were solved using the single-mode relaxation time approximation (RTA) method.
  4. Electric Field Application: A uniform electric field was applied directly to the diamond primitive cell, simulating two parallel electrode plates sandwiching the bulk material.
  5. Crystallographic Analysis: Calculations were performed across three major crystal orientation groups: <100>, <110>, and <111>, to investigate anisotropic effects.
  6. Convergence Parameters: High-precision convergence was ensured using a $12 \times 12 \times 12$ grid for the Brillouin zone and a $56 \times 56 \times 56$ q grid for thermal conductivity calculations.

This research confirms that SCD diamond is not merely a passive heat spreader but an actively tunable thermal material, provided the material quality is sufficient to support the theoretical thermal properties. 6CCVD is uniquely positioned to supply the high-purity, customized SCD required for experimental validation and device integration of this technology.

To replicate and extend this research, the highest quality, low-defect material is essential to minimize intrinsic scattering and achieve the benchmark thermal conductivity ($\approx 2000 \text{ W}\cdot\text{m}^{-1}\text{K}^{-1}$).

  • Optical Grade Single Crystal Diamond (SCD): Required for achieving ultra-high purity and low nitrogen/defect concentrations, ensuring the intrinsic phonon properties necessary for E-field modulation are preserved.
  • Custom Orientation SCD: Since the [111] direction showed the maximum thermal modulation (36.1% increase), 6CCVD offers custom SCD growth specifically oriented along the [111] axis to maximize device performance.

Implementing E-field modulation requires precise material engineering and integration of electrodes capable of handling high voltages (up to $2 \times 10^{9} \text{ V}\cdot\text{cm}^{-1}$ across the lattice).

Requirement from Research6CCVD Customization CapabilityBenefit to Researcher/Engineer
High-Purity SubstratesSCD plates/wafers up to 500 ”m thick, Ra < 1 nm polished.Ensures low intrinsic scattering, maximizing the baseline thermal conductivity and E-field effect.
Specific Crystal OrientationCustom growth of [111] oriented SCD.Enables maximum thermal tunability (36.1% modulation) as predicted by the calculations.
Electrode IntegrationIn-house metalization services: Au, Pt, Pd, Ti, W, Cu.Allows for direct fabrication of parallel plate electrodes necessary to apply the uniform E-field for experimental validation.
Custom DimensionsPlates/wafers up to 125 mm (PCD) or custom SCD dimensions.Supports prototyping and scaling of E-field tunable diamond heat sinks for high-power devices (e.g., IGBTs).
Substrate ThicknessSCD/PCD substrates up to 10 mm thick.Provides robust platforms for high-voltage testing and integration into complex thermal stacks.

6CCVD’s in-house PhD team specializes in the fundamental physics and engineering of MPCVD diamond. We offer comprehensive support for projects focused on:

  • E-field Tunable Thermal Management: Assisting researchers in selecting the optimal SCD orientation and purity grade to experimentally validate the predicted thermal modulation effects.
  • High-Frequency Phonon Engineering: Consulting on material specifications (e.g., isotopic purity) to further manipulate phonon scattering mechanisms in conjunction with external fields.
  • Advanced Semiconductor Integration: Providing material solutions for high-power IGBT and quantum sensor applications where dynamic thermal control is paramount.

For custom specifications or material consultation, visit 6ccvd.com or contact our engineering team directly.

View Original Abstract

For the preparation of diamond heat sinks with ultra-high thermal conductivity by Chemical Vapor Deposition (CVD) technology, the influence of diamond growth direction and electric field on thermal conductivity is worth exploring. In this work, the phonon and thermal transport properties of diamond in three crystal orientation groups (<100>, <110>, and <111>) were investigated using first-principles calculations by electric field. The results show that the response of the diamond in the three-crystal orientation groups presented an obvious anisotropy under positive and negative electric fields. The electric field can break the symmetry of the diamond lattice, causing the electron density around the C atoms to be segregated with the direction of the electric field. Then the phonon spectrum and the thermodynamic properties of diamond were changed. At the same time, due to the coupling relationship between electrons and phonons, the electric field can affect the phonon group velocity, phonon mean free path, phonon-phonon interaction strength and phonon lifetime of the diamond. In the crystal orientation [111], when the electric field strength is ±0.004 a.u., the thermal conductivity is 2654 and 1283 W·m−1K−1, respectively. The main reason for the change in the thermal conductivity of the diamond lattice caused by the electric field is that the electric field has an acceleration effect on the extranuclear electrons of the C atoms in the diamond. Due to the coupling relationship between the electrons and the phonons, the thermodynamic and phonon properties of the diamond change.

  1. 2013 - First Principle Study of Valence-Band Offsets at AlN/Diamond Heterojunctions [Crossref]
  2. 2020 - Investigation on the Interface Characteristic between TiN and Diamond by First-Principles Calculation [Crossref]
  3. 1993 - Thermal Conductivity of Isotopically Modified Single Crystal Diamond [Crossref]
  4. 2001 - Thermal Conductivity Measurements on CVD Diamond [Crossref]
  5. 1997 - Measuring Local Thermal Conductivity in Polycrystalline Diamond with a High Resolution Photothermal Microscope [Crossref]
  6. 2018 - Lattice Thermal Transport in Superhard Hexagonal Diamond and Wurtzite Boron Nitride: A Comparative Study with Cubic Diamond and Cubic Boron Nitride [Crossref]
  7. 2017 - Design of Interfacial Cr3C2 Carbide Layer via Optimization of Sintering Parameters Used to Fabricate Copper/Diamond Composites for Thermal Management Applications [Crossref]
  8. 2019 - Tunable Coef Fi Cient of Thermal Expansion of Cu-B/Diamond Composites Prepared by Gas Pressure in Fi Ltration [Crossref]
  9. 2020 - Fabrication of High Thermal Conductivity Copper/Diamond Composites by Electrodeposition under Potentiostatic Conditions [Crossref]