Surface-enhanced Raman spectroscopy studies of organophosphorous model molecules and pesticides

This work reports the analytical application of surface-enhanced Raman spectroscopy (SERS) in the trace analysis of organophosphorous pesticides (trichlorfon and glyphosate) and model organophosphorous compounds (dimethyl methylphosphonate and o-ethyl methylphosphonothioate) bearing different functional groups. SERS measurements were carried out using Ag nanocubes with an edge square dimension of ca. 100 nm as substrates. Density functional theory (DFT) with the B3LYP functional was used for the optimization of ground state geometries and simulation of Raman spectra of the organophosphorous compounds and their silver complexes. Adsorption geometries and marker bands were identified for each of the investigated compound. Results indicate the usefulness of SERS methodology for the sensitive analyses of organophosphorous compounds through the use of vibrational spectroscopy.

Quantum Monte Carlo study of small aluminum clusters Al-n (n=2-13)

Using fixed node diffusion quantum Monte Carlo (FN-DMC) simulations and density functional theory (DFT) within the generalized gradient approximations, we calculate the total energies of the relaxed and unrelaxed neutral, cationic, and anionic aluminum clusters, Al-n (n = 1-13). From the obtained total energies, we extract the ionization potential and electron detachment energy and compare with previous theoretical and experimental results. Our results for the electronic properties from both the FN-DMC and DFT calculations are in reasonably good agreement with the available experimental data. A comparison between the FN-DMC and DFT results reveals that their differences are a few tenths of electron volt for both the ionization potential and the electron detachment energy. We also observe two distinct behaviors in the electron correlation contribution to the total energies from smaller to larger clusters, which could be assigned to the structural transition of the clusters from planar to three-dimensional occurring at n = 4 to 5.

Spectroscopic Study on the Structural Differences of Thermally Induced Cross-Linking Segments in Emeraldine Salt and Base Forms of Polyaniline

This paper reports the spectroscopic study on the structural differences of thermally induced cross-linking segments in polyaniline in its emeraldine salt (PANI-ES) and base (PANI-EB) forms. Casting films of PANI-ES (ES-film) and PANI-EB (EB-film) were prepared and heated at 150 degrees C under atmospheric air for 30 min. Raman spectra excited at 632.8 nm of heated ES-film presented the characteristic bands of phenazine-like structures at 1638, 1392, and 575 cm(-1), whereas EB-film showed lower relative intensities for these bands. The lower content of phenazine-like segments in heated EB-film is related to residual polaronic segments from preparation procedures, as revealed by Raman. This statement was confirmed by a sequence of thermal and doping experiments in both films. Quantum-chemical calculations by density functional theory (DFT) and time-dependent density functional theory (TD-DFT) showed that the phenazine-like structure presents the intense Raman band at 1350 cm(-1) due to heterocycle breathing mode, and the non-phenazine-like structure (substituted hydrophenazine-type) presents higher energy for HOMO-LUMO transition, indicating the lack of conjugation in the heterocycle compared with the phenazine-like structure. According to experimental and theoretical data reported here, it is proposed that only thermally treated PANI-ES presents phenazine-like rings, whereas PANI-EB presents heterocyclic non-aromatic structures.

VCD to determine absolute configuration of natural product molecules: secolignans from Peperomia blanda

The absolute configuration and solution-state conformers of three peperomin-type secolignans isolated from Peperomia blanda (Piperaceae) are unambiguously determined by using vibrational circular dichroism (VCD) spectroscopy associated with density functional theory (DFT) calculations. Advantages of VCD over the electronic form of CD for the analysis of diastereomers are also discussed. This work extends our growing knowledge about secondary metabolites within the Piperaceae family species while providing a definitive and straightforward method to assess the absolute stereochemistry of secolignans.

The 2:1 cycloadducts from [3+2] 1,3-dipolar cycloaddition of nitrile oxide and vinylacetic acid. Synthesis and liquid crystal behaviour

Four liquid crystals (LC) 3,7a-bis(4-alkyloxyphenyl)-7,7a-dihydro-6H-isoxazolo[2,3-d][1,2,4]oxadiazol-6-yl)acetic acid (7a-d) were synthesised and the mesomorphic behaviour reported. The LCs were characterised as 2: 1 bisadducts, which were obtained from a double [3+2] 1,3-dipolar cycloaddition. In the first step, the cycloaddition of 4-alkyloxyphenylnitrile oxide (4a-d) and vinylacetic acid (5) gave the initial unobserved 1:1 cycloadducts 2-[3-(4-alkyloxyphenyl)-4,5-dihydroisoxazol-5-yl]acetic acid (6a-d). In the second step, the addition of a second equivalent of 4 to 6 yielded the 2: 1 bisadducts 7a-d without any traces of 6. All compounds 7a-d were unstable during the transition from the mesophase to the isotropic state upon first heating as evidenced by the large peaks in the differential scanning calorimetry traces. Due to the chemical instability of the compounds upon heating, the transition temperature related to the smectic C to smectic A transitions was acquired by means of an image processing method. X-Ray diffraction experiments were also used to analyse the liquid-crystalline phases. A theoretical calculation was performed using density functional theory (DFT) methods at the PBE1PBE/6-311+G(2d,p) level (with solvent effect) in order to get information about the energetic profile of the 2: 1 cycloaddition. DFT studies revealed that the cycloaddition process is controlled by the HOMO(dipolarophile) - LUMO(1,3-dipole), and that the double [3+2] 1,3-dipolar cycloaddition reaction is quite possible.

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

First principles calculations for vacancies and antisites in PbSe and PbTe: bulk and nanowire.

The energetic stability and the electronic properties of vacancies (VX) and antisites (XY) in PbSe and PbTe are investigated. PbSe and PbTe are narrow band gap semiconductors and have the potential to be used in infrared detectors, laser, and diodes. They are also of special interest for thermoelectric devices (TE). The calculations are based in the Density Functional Theory (DFT) and the General Gradient Approximation (GGA) for the exchange-correlation term, as implemented in the VASP code. The core and valence electrons are described by the Projected Augmented Wave (PAW) and the Plane Wave (PW) methods, respectively. The defects are studied in the bulk and nanowire (NW) system. Our results show that intrinsec defects (vacancies and antisites) in PbTe have lower formation energies in the NW as compared to the bulk and present a trend in migrate to the surface of the NW. For the PbSe we obtain similar results when compare the formation energy for the bulk and NW. However, the Pb vacancy and the antisites are more stable in the core of the NW. The intrinsec defects are shallow defects for the bulk system. For both PbSe and PbTe VPb is a shallow acceptor defect and VSe and VT e are shallow donor defects for the PbSe and PbTe, respectively. Similar electronic properties are observed for the antisites. For the Pb in the anion site we obtain an n-type semiconductor for both PbSe and PbTe, SeP b is a p-type for the PbSe, and T eP b is a n-type for PbTe. Due the quantum con¯nement effects present in the NW (the band gap open), these defects have different electronic properties for the NW as compared to the bulk. Now these defects give rise to electronic levels in the band gap of the PbTe NW and the VT e present a metallic character. For the PbSe NW a p-type and a n-type semiconductor is obtained for the VP b and P bSe, respectively. On the other hand, deep electronic levels are present in the band gap for the VSe and SePb. These results show that due an enhanced in the electronic density of states (DOS) near the Fermi energy, the defective PbSe and PbTe are candidates for efficient TE devices.

Reaction intermediates of ethanol electro-oxidation on platinum investigated by SFG spectroscopy

Although electrochemical oxidation of simple organic molecules on metal catalysts is the basic ingredient of fuel cells, which have great technological potential as a renewable source of electrical energy, the detailed reaction mechanisms are in most cases not completely understood. Here, we investigate the ethanol-platinum interface in acidic aqueous solution using infrared-visible sum frequency generation (SFG) spectroscopy and theoretical calculations of vibrational spectra in order to identify the intermediates present during the electro-oxidation of ethanol. The complex vibrational spectrum in the fingerprint region imply on the coexistence of several adsorbates. Based on spectra in ultra-high-vacuum (UHV) and electrochemical environment from the literature and our density functional theory (DFT) calculations of vibrational spectra, new adsorbed intermediates, never before observed with conventional infrared (IR) spectroscopy, are proposed here: g2-acetaldehyde, g2-acetyl, ethylidyne, monodentate acetate, methoxy, tertiary methanol derivative, COH residue, g2-formaldehyde, mono and bidentate formate, CH3 and CH2 residues. In addition, we present new evidences for an ethoxy intermediate, a secondary ethanol derivative and an acetyl species, and we confirm the presence of previously observed adsorbates: a tertiary ethanol derivative, bidentate acetate, and COad. These results indicate that the platinum surface is much more reactive, and the reaction mechanism for ethanol electro-oxidation is considerably more complex than previously considered. This might be also true for many other molecule-catalyst systems.

Chemical Bonding of Hydrocarbons to Metal Surfaces

Using x-ray absorption spectroscopy (XAS), x-ray emission spectroscopy (XES) and x-ray photoelectron spectroscopy (XPS) in combination with density functional theory (DFT) the changes in electronic and geometric structure of hydrocarbons upon adsorption are determined. The chemical bonding is analyzed and the results provide new insights in the mechanisms responsible for dehydrogenation in heterogeneous catalysis. In the case of alkanes, n-octane and methane are studied. XAS and XES show significant changes in the electronic structure upon adsorption. XES shows new adsorption induced occupied states and XAS shows quenching of CH*/Rydberg states in n-octane. In methane the symmetry forbidden gas phase lowest unoccupied molecular orbital becomes allowed due to broken symmetry. New adsorption induced unoccupied features with mainly metal character appear just above the Fermi level in XA spectra of both adsorbed methane and n-octane. These changes are not observed in DFT total energy geometry optimizations. Comparison between experimental and computed spectra for different adsorbate geometries reveals that the molecular structures are significantly changed in both molecules. The C-C bonds in n-octane are shortened upon adsorption and the C-H bonds are elongated in both n-octane and methane. In addition ethylene and acetylene are studied as model systems for unsaturated hydrocarbons. The validity of both the Dewar-Chatt-Duncanson chemisorption model and the alternative spin-uncoupling picture is confirmed, as well as C-C bond elongation and upward bending of the C-H bonds. The bonding of ethylene to Cu(110) and Ni(110) are compared and the results show that the main difference is the amount of back-donation into the molecular π* orbital, which allows the molecule to desorb molecularly from the Cu(110) surface, whereas it is dehydrogenated upon heating on the Ni(110) surface. Acetylene is found to adsorb in two different adsorption sites on the Cu(110) surface at liquid nitrogen temperature. Upon heating the molecules move into one of these sites due to attractive adsorbate-adsorbate interaction and only one adsorbed species is present at room temperature, at which point the molecules start reacting to form benzene. The bonding of the two species is very similar in both sites and the carbon atoms are rehybridized essentially to sp2.

Local Structure of Hydrogen-Bonded Liquids

Ordinary yet unique, water is the substance on which life is based. Water seems, at first sight, to be a very simple molecule, consisting of two hydrogen atoms attached to one oxygen. Its small size belies the complexity of its action and its numerous anomalies, central to a broad class of important phenomena, ranging from global current circulation, terrestrial water and CO2 cycles to corrosion and wetting. The explanation of this complex behavior comes from water's unique ability to form extensive three-dimensional networks of hydrogen-bonds, whose nature and structures, in spite of a great deal of efforts involving a plethora of experimental and theoretical techniques, still lacks a complete scientific understanding. This thesis is devoted to the study of the local structure of hydrogen-bonded liquids, with a particular emphasis on water, taking advantage of a combination of core-level spectroscopies and density functional theory spectra calculations. X-ray absorption, in particular, is found to be sensitive to the local hydrogen-bond environment, thus offering a very promising tool for spectroscopic identification of specific structural configurations in water, alcohols and aqueous solutions. More specifically, the characteristic spectroscopic signature of the broken hydrogen-bond at the hydrogen side is used to analyze the structure of bulk water, leading to the finding that most molecules are arranged in two hydrogen-bond configurations, in contrast to the picture provided by molecular dynamics simulations. At the liquid-vapor interface, an interplay of surface sensitive measurements and theoretical calculations enables us to distinguish a new interfacial species in equilibrium with the gas. In a similar approach the cluster form of the excess proton in highly concentrated acid solutions and the different coordination of methanol at the vacuum interface and in the bulk can also be clearly identified. Finally the ability of core-level spectroscopies, aided by sophisticated density functional theory calculations, to directly probe the valence electronic structure of a system is used to observe the nature of the interaction between water molecules and solvated ions in solution. Water around transition metal ions is found to interact with the solute via orbital mixing with the metal d-orbitals. The hydrogen-bond between water molecules is explained in terms of electrostatic interactions enhanced by charge rehybridization in which charge transfer between connecting molecules is shown to be fundamental.

Spectroscopic signatures from first principles calculations: from surface adsorbates to liquids and polymers

In dieser Arbeit werden drei wasserstoffverbrückte Systeme in der kondensierten Phase mit Hilfe von first-principles-Elektronenstruktur-Rechnungen untersucht, die auf der Dichtefunktionaltheorie (DFT) unter periodischen Randbedingungen basieren. Ihre lokalen Konformationen und Wasserstoffbrückenbindungen werden mittels ab-initio Molekulardynamiksimulationen berechnet und weiterhin durch die Bestimmung ihrer spektroskopischen Parameter charakterisiert. Der Schwerpunkt liegt dabei auf lokalen Strukturen und auf schnellen Fluktuationen der Wasserstoffbrückenbindungen, welche von zentraler Bedeutung für die physikalischen und chemischen Eigenschaften der betrachteten Systeme sind. Die für die lokalen, instantanen Konformationen berechneten Spektren werden verwendet, um die physikalischen Prozesse, die hinter den untersuchten Phänomenen stehen, zu erklären: die Wasseradsorption auf metallischen Oberflächen, die Ionensolvatisierung in wässrigen Lösungen und der Protonentransport in protonleitenden Polymeren, welche Prototypen von Membranen für Brennstoffzellen sind. Die Möglichkeit der Vorhersage spektroskopischer Parameter eröffnet vielfältige Möglichkeiten des Dialogs zwischen Experimenten und numerischen Simulationen. Die in dieser Arbeit vorgestellten Ergebnisse zeigen, dass die Zuverlässigkeit dieser theoretischen Berechnungen inzwischen für viele experimentell relevante Systeme ein quantitatives Niveau erreicht hat.

Computational studies on the structure and stabilities of magnetic inter-metallic compounds

The purpose of this thesis is to further the understanding of the structural, electronic and magnetic properties of ternary inter-metallic compounds using density functional theory (DFT). Four main problems are addressed. First, a detailed analysis on the ternary Heusler compounds is made. It has long been known that many Heusler compounds ($X_2YZ$; $X$ and $Y$ transition elements, $Z$ main group element) exhibit interesting half-metallic and ferromagnetic properties. In order to understand these, the dependence of magnetic and electronic properties on the structural parameters, the type of exchange-correlation functional and electron-electron correlation was examined. It was found that almost all Co$_2YZ$ Heusler compounds exhibit half-metallic ferromagnetism. It is also observed that $X$ and $Y$ atoms mainly contribute to the total magnetic moment. The magnitude of the total magnetic moment is determined only indirectly by the nature of $Z$ atoms, and shows a trend consistent with Slater-Pauling behaviour in several classes of these compounds. In contrast to experiments, calculations give a non-integer value of the magnetic moment in certain Co$_2$-based Heusler compounds. To explain deviations of the calculated magnetic moment, the LDA+$U$ scheme was applied and it was found that the inclusion of electron-electron correlation beyond the LSDA and GGA is necessary to obtain theoretical description of some Heusler compounds that are half-metallic ferromagnets. The electronic structure and magnetic properties of substitutional series of the quaternary Heusler compound Co$_2$Mn$_{1-x}$Fe$_x$Si were investigated under LDA+$U$. The calculated band structure suggest that the most stable compound in a half-metallic state will occur at an intermediate Fe concentration. These calculated findings are qualitatively confirmed by experimental studies. Second, the effect of antisite disordering in the Co$_2$TiSn system was investigated theoretically as well as experimentally. Preservation of half-metallicity for Co$_2$TiSn was observed with moderate antisite disordering and experimental findings suggest that the Co and Ti antisites disorder amounts to approximately 10~% in the compound. Third, a systematic examination was carried out for band gaps and the nature (covalent or ionic) of bonding in semiconducting 8- and 18-electron or half-metallic ferromagnet half-Heusler compounds. It was found that the most appropriate description of these compounds from the viewpoint of electronic structures is one of a $YZ$ zinc blende lattice stuffed by the $X$ ion. Simple valence rules are obeyed for bonding in the 8- and 18-electron compounds. Fourth, hexagonal analogues of half-Heusler compounds have been searched. Three series of compounds were investigated: GdPdSb, GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP. GdPdSb is suggested as a possible half-metallic weak ferromagnet at low temperature. GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP were investigated because they exhibit interesting bonding, structural and magnetic properties. The results qualitatively confirm experimental studies on magnetic and structural behaviour in GdPdSb, GdAutextit{X} (textit{X} = Mn, Cd and In) and EuNiP compounds. ~

Numerische und analytische Untersuchungen stark korrelierter fermionischer Mehrbandsysteme

In dieser Arbeit werden vier unterschiedliche, stark korrelierte, fermionische Mehrbandsysteme untersucht. Es handelt sich dabei um ein Mehrstörstellen-Anderson-Modell, zwei Hubbard-Modelle sowie ein Mehrbandsystem, wie es sich aus einer ab initio-Beschreibung für ein korreliertes Halbmetall ergibt.rnrnDie Betrachtung des Mehrstörstellen-Anderson-Modells konzentriert sich auf die Untersuchung des Einflusses der Austauschwechselwirkung und der nicht-lokalen Korrelationen zwischen zwei Störstellen in einem einfach-kubischen Gitter. Das zentrale Resultat ist die Abstandsabhängigkeit der Korrelationen der Störstellenelektronen, welche stark von der Gitterdimension und der relativen Position der Störstellen abhängen. Bemerkenswert ist hier die lange Reichweite der Korrelationen in der Diagonalrichtung des Gitters. Außerdem ergibt sich, dass eine antiferromagnetische Austauschwechselwirkung ein Singulett zwischen den Störstellenelektronen gegenüber den Kondo-Singuletts der einzelnen Störstellen favorisiert und so den Kondo-Effekt der einzelnen Störstellen behindert.rnrnEin Zweiband-Hubbard-Modell, das Jz-Modell, wird im Hinblick auf seine Mott-Phasen in Abhängigkeit von Dotierung und Kristallfeldaufspaltung auf dem Bethe-Gitter untersucht. Die Entartung der Bänder ist durch eine unterschiedliche Bandbreite aufgehoben. Wichtigstes Ergebnis sind die Phasendiagramme in Bezug auf Wechselwirkung, Gesamtfüllung und Kristallfeldparameter. Im Vergleich zu Einbandmodellen kommen im Jz-Modell sogenannte orbital-selektive Mott-Phasen hinzu, die, abhängig von Wechselwirkung, Gesamtfüllung und Kristallfeldparameter, einerseits metallischen und andererseits isolierenden Charakter haben. Ein neuer Aspekt ergibt sich durch den Kristallfeldparameter, der die ionischen Einteilchenniveaus relativ zueinander verschiebt, und für bestimmte Werte eine orbital-selektive Mott-Phase des breiten Bands ermöglicht. Im Vergleich mit analytischen Näherungslösungen und Einbandmodellen lassen sich generische Vielteilchen- und Korrelationseffekte von typischen Mehrband- und Einteilcheneffekten differenzieren.rnrnDas zweite untersuchte Hubbard-Modell beschreibt eine magneto-optische Falle mit einer endlichen Anzahl Gitterplätze, in welcher fermionische Atome platziert sind. Es wird eine z-antiferromagnetische Phase unter Berücksichtigung nicht-lokaler Vielteilchenkorrelationen erhalten, und dabei werden bekannte Ergebnisse einer effektiven Einteilchenbeschreibung verbessert.rnrnDas korrelierte Halbmetall wird im Rahmen einer Mehrbandrechnung im Hinblick auf Korrelationseffekte untersucht. Ausgangspunkt ist eine ab initio-Beschreibung durch die Dichtefunktionaltheorie (DFT), welche dann durch die Hinzunahme lokaler Korrelationen ergänzt wird. Die Vielteilcheneffekte werden an Hand einer einfachen Wechselwirkungsnäherung verdeutlicht, und für ein Wechselwirkungsmodell in sphärischer Symmetrie präzisiert. Es ergibt sich nur eine schwache Quasiteilchenrenormierung. Besonders für röntgenspektroskopische Experimente wird eine gute Übereinstimmung erzielt.rnrnDie numerischen Ergebnisse für das Jz-Modell basieren auf Quanten-Monte-Carlo-Simulationen im Rahmen der dynamischen Molekularfeldtheorie (DMFT). Für alle anderen Systeme wird ein Mehrband-Algorithmus entwickelt und implementiert, welcher explizit nicht-diagonale Mehrbandprozesse berücksichtigt.rnrn

Ab initio quantum mechanical investigation of structural and chemical-physical properties of selected minerals for minero-petrological, structural ceramic and biomaterial applications

The purpose of this thesis is the atomic-scale simulation of the crystal-chemical and physical (phonon, energetic) properties of some strategically important minerals for structural ceramics, biomedical and petrological applications. These properties affect the thermodynamic stability and rule the mineral-environment interface phenomena, with important economical, (bio)technological, petrological and environmental implications. The minerals of interest belong to the family of phyllosilicates (talc, pyrophyllite and muscovite) and apatite (OHAp), chosen for their importance in industrial and biomedical applications (structural ceramics) and petrophysics. In this thesis work we have applicated quantum mechanics methods, formulas and knowledge to the resolution of mineralogical problems ("Quantum Mineralogy”). The chosen theoretical approach is the Density Functional Theory (DFT), along with periodic boundary conditions to limit the portion of the mineral in analysis to the crystallographic cell and the hybrid functional B3LYP. The crystalline orbitals were simulated by linear combination of Gaussian functions (GTO). The dispersive forces, which are important for the structural determination of phyllosilicates and not properly con-sidered in pure DFT method, have been included by means of a semi-empirical correction. The phonon and the mechanical properties were also calculated. The equation of state, both in athermal conditions and in a wide temperature range, has been obtained by means of variations in the volume of the cell and quasi-harmonic approximation. Some thermo-chemical properties of the minerals (isochoric and isobaric thermal capacity) were calculated, because of their considerable applicative importance. For the first time three-dimensional charts related to these properties at different pressures and temperatures were provided. The hydroxylapatite has been studied from the standpoint of structural and phonon properties for its biotechnological role. In fact, biological apatite represents the inorganic phase of vertebrate hard tissues. Numerous carbonated (hydroxyl)apatite structures were modelled by QM to cover the broadest spectrum of possible biological structural variations to fulfil bioceramics applications.