Skip to content
6CCVD Research

Fully analytic valence force field model for the elastic and inner elastic properties of diamond and zincblende crystals

At a Glance

Section titled “At a Glance”
MetadataDetails
Publication Date2019-09-30
JournalPhysical review. B./Physical review. B
AuthorsDaniel S. P. Tanner, A. Miguel, Stefan Schulz, Eoin P. O’Reilly
InstitutionsAalto University, University College Cork
Citations11
AnalysisFull AI Review Included

Technical Documentation and Analysis: Valence Force Field Modeling for Diamond and Zincblende Materials

Section titled “Technical Documentation and Analysis: Valence Force Field Modeling for Diamond and Zincblende Materials”

Executive Summary

Section titled “Executive Summary”

This paper presents a robust, fully analytic Valence Force Field (VFF) model that provides superior accuracy for predicting the mechanical and vibrational properties of diamond and zincblende (ZB) crystals. This advanced modeling capability is essential for engineers designing high-precision, strained diamond-based heterostructures and quantum devices.

  • Advanced Modeling: Developed a novel Coulombic VFF model (Model C) capable of fully describing the elastic energy density of diamond and ZB crystals.
  • Precision Elastic Constants: The VFF approach provides an exact description of the three cubic elastic constants ($C_{11}, C_{12}, C_{44}$) and the Kleinman internal strain parameter ($\zeta$).
  • Targeted Stability: A critical stabilization criterion was met by analytically fitting the VFF effective charge parameter ($S$) directly to the inner elastic constant ($E_{11}$), ensuring crystal stability against internal strain.
  • Validation: Model C demonstrated superior performance in benchmarking against LDA Density Functional Theory (DFT) calculations, achieving maximum bond length differences ($\Delta r_{ij}$) of only 0.21% in strained GaAs/InAs supercells.
  • Improved Spectra: The new free parameterization produces improved acoustic and optical phonon spectra compared to conventional models, vital for thermal and optical device design.
  • Application Focus: Provides computationally efficient and high-accuracy interatomic potentials suitable for calculating strain and relaxation in large, non-homogeneous systems (e.g., diamond quantum dots or strained PCD interfaces).

Technical Specifications

Section titled “Technical Specifications”

The following key material parameters and performance metrics were derived or targeted using the VFF models, illustrating the level of precision required for advanced crystal engineering.

ParameterValue (GaAs)UnitContext
LDA DFT $C_{11}$115GPaLongitudinal Elastic Constant Input
LDA DFT $C_{44}$58GPaShear Elastic Constant Input
LDA DFT $E_{11}$34GPa Å-2Inner Elastic Constant (Crucial for VFF Model C Fitting)
Kleinman Parameter $\zeta$0.547DimensionlessMeasures atomic sub-lattice response to shear strain
Anisotropy Factor $A$1.83Dimensionless$A = 2C_{44} / (C_{11} - C_{12})$. Must be A < 2 for covalent VFF stability
Model C $\Delta r_{ij}$ Error0.21%Average difference in calculated bond length vs. DFT (All supercells)
Model C $\Delta \theta$ Error0.32%Average difference in calculated bond angle vs. DFT (All supercells)
SCD Stability Criterion$E_{11} > 0$GPa Å-2Required condition for crystal stability against internal strain

Key Methodologies

Section titled “Key Methodologies”

The core of the research focused on deriving explicit analytic expressions for the VFF force constants ($k_{r}, k_{\theta}, k_{rr}, k_{r\theta}$) and the effective charge parameter ($S$) based on known macroscopic elastic properties ($C_{ij}$) and internal strain properties ($E_{11}, \zeta$).

  1. VFF Potential Base: Martin’s VFF model was utilized, incorporating short-range bond-stretching ($k_{r}$), bond-bending ($k_{\theta}$), bond-bond stretching ($k_{rr}$), and bond-stretching angle-bending ($k_{r\theta}$) terms, augmented by long-range Coulombic terms.
  2. Input Data Acquisition: Baseline material properties (elastic constants $C_{ij}$, inner elastic constants $E_{11}$, and Kleinman parameter $\zeta$) were obtained from high-accuracy LDA Density Functional Theory (DFT) calculations.
  3. Covalent VFF Derivation (Model A): Analytic expressions for force constants were solved by inverting the relations between VFF terms and elastic constants, neglecting Coulomb interaction ($S=0$). This model was shown to be unstable for materials where the anisotropy factor $A > 2$ (e.g., GaN, AlN).
  4. Conventional Coulombic VFF Derivation (Model B): The effective charge parameter ($S$) was calculated using the zone-center transverse and longitudinal optical phonon mode splitting ($\omega_{LO}, \omega_{TO}$). This stabilized mildly ionic materials but failed for highly ionic III-N compounds.
  5. Free Parameterization of Effective Charge (Model C - Optimal): The effective charge parameter ($S$) was determined analytically by explicitly fitting the model to ensure exact reproduction of the internal elastic constant $E_{11}$.
  6. Benchmarking: Structural relaxations were performed on large, strained supercells (e.g., [001] and [111] GaAs/InAs interfaces) using the developed VFF potentials and compared against the DFT reference data to quantify accuracy.

6CCVD Solutions & Capabilities

Section titled “6CCVD Solutions & Capabilities”

This research validates the critical importance of precisely controlled internal strain parameters ($E_{11}, \zeta$) for the successful engineering of semiconductor heterostructures, including diamond-based systems. 6CCVD’s expertise in customized MPCVD growth and post-processing directly addresses the material requirements defined by these advanced theoretical models.

Applicable Materials

Section titled “Applicable Materials”

To replicate and extend this research for high-performance devices, researchers require materials that are intrinsically stable and processed to exacting dimensional tolerances.

  • Optical Grade Single Crystal Diamond (SCD): Required for applications relying on minimal strain, high purity, and superior thermal management. Essential for stable quantum emitter hosting (e.g., NV centers) where strain must be minimized or precisely managed.
    • Recommendation: High-Purity SCD plates, specified for extremely low residual strain (stress $<$ 100 GPa).
  • Polycrystalline Diamond (PCD) / Thin Film SCD: Applicable for substrates or heterostructure layers where specific elastic constants are designed for strain engineering, as modeled by the VFF.
    • Capability Match: We offer SCD/PCD films with thicknesses ranging from 0.1 ”m up to 500 ”m, allowing precise control over layer dimensions critical for superlattice modeling.
  • Boron-Doped Diamond (BDD): Used as a conductive layer in electronic heterostructures.
    • Capability Match: We supply BDD films with tunable doping concentration, enabling specific control over electronic and electro-mechanical properties relevant to inner elastic constants.

Customization Potential

Section titled “Customization Potential”

The VFF methodology validates the need for atomic-scale precision in heterostructures and strained alloys. 6CCVD provides the necessary customization to translate theoretical stability into physical reality.

Theoretical Requirement6CCVD Customization Solution
Precise Layer Thickness/DimensionPlates/wafers available up to 125mm (PCD); custom SCD thickness control from 0.1 ”m to 500 ”m.
Strain Management & Interface QualityPolishing capabilities achieve ultra-low roughness: Ra < 1 nm (SCD) and Ra < 5 nm (Inch-size PCD), crucial for high-fidelity, low-defect interfaces in strained systems.
Electronic/Contact IntegrationCustom Metalization Services (Au, Pt, Pd, Ti, W, Cu) available internally for electrode formation or bonding layers on diamond, supporting complex device architectures.
Unique Material GeometryPrecision laser cutting and shaping services to meet specific device geometry requirements (e.g., mesas, waveguides) derived from VFF stress analysis.

Engineering Support

Section titled “Engineering Support”

The complexity of fitting elastic constants and internal strain parameters, as demonstrated by the three VFF models in this paper, underscores the need for expert material consultation.

  • 6CCVD maintains an in-house PhD-level engineering team specialized in MPCVD physics and diamond properties.
  • Our experts can assist researchers in selecting the optimal MPCVD diamond material grade and processing specifications (thickness, polishing, metalization) necessary to implement device designs based on VFF or DFT modeling of strained heterostructures and quantum dots.
  • We offer global shipping (DDU default, DDP available) to ensure timely delivery of custom materials worldwide.

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

View Original Abstract

Using a valence force field model based on that introduced by Martin, we present three related methods through which we analytically determine valence force field parameters. The methods introduced allow easy derivation of valence force field parameters in terms of the Kleinman parameter ζ and bulk properties of zincblende and diamond crystals. We start with a model suited for covalent and weakly ionic materials, where the valence force field parameters are derived in terms of ζ and the bulk elastic constants C11, C12, and C44. We show that this model breaks down as the material becomes more ionic and specifically when the elastic anisotropy factor A=2C44/(C11−C12)>2. The analytic model can be stabilized for ionic materials by including Martin’s electrostatic terms with effective cation and anion charges in the valence force field model. Inclusion of effective charges determined via the optical phonon mode splitting provides a stable model for all but two of the materials considered (zincblende GaN and AlN). A stable model is obtained for all materials considered by also utilizing the inner elastic constant E11 to determine the magnitude of the effective charges used in the Coulomb interaction. Test calculations show that the models describe well structural relaxation in superlattices and alloys and reproduce key phonon band structure features.

Tech Support

Section titled “Tech Support”

Original Source

Section titled “Original Source”

References

Section titled “References”
  1. 1954 - Dynamical Theory of Crystal Lattices

Reads Similar to This

Age Estimation from Teeth in Longman’s Beaked Whales (Indopacetus Pacificus) Stranded in New Caledonia (South Pacific) Facet-Dependent Temporal and Spatial Changes in Boron-Doped Diamond Film Electrodes due to Anodic Corrosion Towards the scale‐up of electrolysis with diamond anodes - effect of stacking on the electrochemical oxidation of 2,4 D