Optoelectronics and defect levels in hydroxyapatite by first-principles.

Hydroxyapatite (HAp) is an important component of mammal bones and teeth, being widely used in prosthetic implants. Despite the importance of HAp in medicine, several promising applications involving this material (e.g., in photo-catalysis) depend on how well we understand its fundamental properties. Among the ones that are either unknown or not known accurately, we have the electronic band structure and all that relates to it, including the bandgap width. We employ state-of-the-art methodologies, including density hybrid-functional theory and many-body perturbation theory within the dynamically screened single-particle Green's function approximation, to look at the optoelectronic properties of HAp. These methods are also applied to the calculation of defect levels. We find that the use of a mix of (semi-)local and exact exchange in the exchange-correlation functional brings a drastic improvement to the band structure. Important side effects include improvements in the description of dielectric and optical properties not only involving conduction band (excited) states but also the valence. We find that the highly dispersive conduction band bottom of HAp originates from anti-bonding σ* states along the ⋯OH-OH-⋯ infinite chain, suggesting the formation of a conductive 1D-ice phase. The choice of the exchange-correlation treatment to the calculation of defect levels was also investigated by using the OH-vacancy as a testing model. We find that donor and acceptor transitions obtained within semi-local density functional theory (DFT) differ from those of hybrid-DFT by almost 2 eV. Such a large discrepancy emphasizes the importance of using a high-quality description of the electron-electron interactions in the calculation of electronic and optical transitions of defects in HAp.

Effect of Magnesium Substitution on Structural Features and Properties of Hydroxyapatite.

Hydroxyapatite (HAP) is the main mineral component of bones and teeth. It is widely used in medicine as a bone filler and coating for implants to promote new bone growth. Ion substitutions into the HAP structure highly affect its properties. One of the most important substituents is magnesium. This paper presents new results obtained using high-precision hybrid density functional theory calculations for Mg/Ca substitutions in HAP in a wide magnesium concentration range within a 2 × 2 × 2 supercell model. Experimental data on the mechanochemical synthesis of HAP-Mg samples with different Mg concentrations are also presented. A comparison between the experiment and the theory showed good agreement: the HAP-Mg unit cell parameters and volume decreased with increasing degree of Mg/Ca substitution. The changes in the distances between the Ca and O, Ca and H, and Mg and O ions upon Mg/Ca substitution in different calcium positions was analyzed. The resulting asymmetry and distortion of the cell parameters were evaluated. It was shown that bulk modulus, energy levels, and band gap depend on the degree of Mg substitutions in the Ca1 and Ca2 positions. The formation energies of Mg/Ca substitutions showed non-monotonic behavior that was different for Ca1 and Ca2 positions. The Ca2 position had a slightly higher probability (~5 meV/f.u.) of substitution than Ca1 position at a Mg concentration x = 0.5. At x = 1, substitution in both positions can coexist. The simulated IR spectra for different Mg/Ca substitutions showed that Mg in the Ca2 position changes the IR spectrum more significantly than Mg in the Ca1 position. Similar changes were recorded in the IR spectra of the synthesized samples. The electronic structure is shown to be sensitive to the number and position of substitutions, which may be used to tweak the optical properties of the HAP-Mg material.

Molecular modeling of the piezoelectric effect in the ferroelectric polymer poly(vinylidene fluoride) (PVDF).

In this work, computational molecular modeling and exploration was applied to study the nature of the negative piezoelectric effect in the ferroelectric polymer polyvinylidene fluoride (PVDF), and the results confirmed by actual nanoscale measurements. First principle calculations were employed, using various quantum-chemical methods (QM), including semi-empirical (PM3) and various density functional theory (DFT) approaches, and in addition combined with molecular mechanics (MM) methods in complex joint approaches (QM/MM). Both PVDF molecular chains and a unit cell of crystalline ß-phase PVDF were modeled. This computational molecular exploration clearly shows that the nature of the so-called negative piezo-electric effect in the ferroelectric PVDF polymer has a self-consistent quantum nature, and is related to the redistribution of the electron molecular orbitals (wave functions), leading to the shifting of atomic nuclei and reorganization of all total charges to the new, energetically optimal positions, under an applied electrical field. Molecular modeling and first principles calculations show that the piezoelectric coefficient d 33 has a negative sign, and its average values lies in the range of d 33 ~ -16.6 to -19.2 pC/N (or pm/V) (for dielectric permittivity ε = 5) and in the range of d 33 ~ -33.5 to -38.5 pC/N (or pm/V) (for ε = 10), corresponding to known data, and allowing us to explain the reasons for the negative sign of the piezo-response. We found that when a field is applied perpendicular to the PVDF chain length, as polarization increases the chain also stretches, increasing its length and reducing its height. For computed value of ε ~ 5 we obtained a value of d31 ~ +15.5 pC/N with a positive sign. This computational study is corroborated by measured nanoscale data obtained by atomic force and piezo-response force microscopy (AFM/PFM). This study could be useful as a basis for further insights into other organic and molecular ferroelectrics.