Influence of Magnesium Source on the Mechanochemical Synthesis of Magnesium-Substituted Hydroxyapatite.

Magnesium, as one of the most abundant cations in the human body, plays an important role in both physiological and pathological processes. In this study, it was shown that a promising biomedical material, Mg-substituted hydroxyapatite (Mg-HA), can be synthesized via a fast mechanochemical method. For this method, the nature of magnesium-containing carriers was shown to be important. When using magnesium oxide as a source of magnesium, the partial insertion of magnesium cations into the apatite structure occurs. In contrast, when magnesium hydroxide or monomagnesium phosphate is used, single-phase Mg-HA is formed. Both experimental and theoretical investigations showed that an increase in the Mg content leads to a decrease in the lattice parameters and unit cell volume of Mg-HA. Density functional theory calculations showed the high sensitivity of the lattice parameters to the crystallographic position of the calcium site substituted by magnesium. It was shown experimentally that the insertion of magnesium cations decreases the thermal stability of hydroxyapatite. The thermal decomposition of Mg-HA leads to the formation of a mixture of stoichiometric HA, magnesium oxide, and Mg-substituted tricalcium phosphate phases.

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.

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.

Molecular modeling of the piezoelectric properties of ferroelectric composites containing polyvinylidene fluoride (PVDF) and either graphene or graphene oxide.

Molecular modeling of ferroelectric composites containing polyvinylidene fluoride (PVDF) and either graphene (G) or graphene oxide (GO) were performed using the semi-empirical quantum approximation PM3 in HyperChem. The piezo properties of the composites were analyzed and compared with experimental data obtained for P(VDF-TrFE)-GO films. Qualitative agreement was obtained between the results of the modeling and the experimental results in terms of the properties of the measured effective piezoelectric coefficient d 33eff and its decrease in the presence of G/GO in comparison with the average computed piezoelectric coefficient . When models incorporating one or several G layers with 54 carbon atoms were investigated, the average piezoelectric coefficient was found to decrease to -9.8 pm/V for the one-sided model PVDF/G and to -18.98 pm/V for the sandwich model G/PVDF/G as compared with the calculated piezoelectric coefficient for pure PVDF ( = -42.2 pm/V computed in present work, and = -38.5 pm/V, obtained from J Mol Model 35 (2013) 19:3591-3602). When models incorporating one or several GO layers with 98 carbon atoms were considered, the piezoelectric coefficient was found to decrease to -14.6 pm/V for the one-sided PVDF/GO model and to -29.8 pm/V for the sandwich GO/PVDF/GO model as compared with the same calculated piezoelectric coefficient for pure PVDF.

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.

Micromechanics of base pair unzipping in the DNA duplex.

All-atom molecular dynamics (MD) simulations of DNA duplex unzipping in a water environment were performed. The investigated DNA double helix consists of a Drew-Dickerson dodecamer sequence and a hairpin (AAG) attached to the end of the double-helix chain. The considered system is used to examine the process of DNA strand separation under the action of an external force. This process occurs in vivo and now is being intensively investigated in experiments with single molecules. The DNA dodecamer duplex is consequently unzipped pair by pair by means of the steered MD. The unzipping trajectories turn out to be similar for the duplex parts with G·C content and rather distinct for the parts with A·T content. It is shown that during the unzipping each pair experiences two types of motion: relatively quick rotation together with all the duplex and slower motion in the frame of the unzipping fork. In the course of opening, the complementary pair passes through several distinct states: (i) the closed state in the double helix, (ii) the metastable preopened state in the unzipping fork and (iii) the unbound state. The performed simulations show that water molecules participate in the stabilization of the metastable states of the preopened base pairs in the DNA unzipping fork.