Spin-phonon coupling effects in transition-metal perovskites:a DFT+UU and hybrid-functional study

Spin-phonon coupling effects, as reflected in phonon frequency shifts between ferromagnetic (FM) and G-type antiferromagnetic (AFM) configurations in cubic CaMnO$_3$, SrMnO$_3$, BaMnO$_3$, LaCrO$_3$, LaFeO$_3$ and La$_2$(CrFe)O$_6$, are investigated using density-functional methods. The calculations are carried out both with a hybrid-functional (HSE) approach and with a DFT+$U$ approach using a $U$ that has been fitted to HSE calculations. The phonon frequency shifts obtained in going from the FM to the AFM spin configuration agree well with those computed directly from the more accurate HSE approach, but are obtained with much less computational effort. We find that in the $A$MnO$_3$ materials class with $A$=Ca, Sr, and Ba, this frequency shift decreases as the A cation radius increases for the $\Gamma$ phonons, while it increases for R-point phonons. In La$M$O$_3$ with $M$=Cr, Fe, and Cr/Fe, the phonon frequencies at $\Gamma$ decrease as the spin order changes from AFM to FM for LaCrO$_3$ and LaFeO$_3$, but they increase for the double perovskite La$_2$(CrFe)O$_6$. We discuss these results and the prospects for bulk and superlattice forms of these materials to be useful as multiferroics.Comment: 13 pages, 7 figures, 9 table

The multiferroic phase of DyFeO3_{3}:an ab--initio study

By performing accurate ab-initio density functional theory calculations, we study the role of $4f$ electrons in stabilizing the magnetic-field-induced ferroelectric state of DyFeO$_{3}$. We confirm that the ferroelectric polarization is driven by an exchange-strictive mechanism, working between adjacent spin-polarized Fe and Dy layers, as suggested by Y. Tokunaga [Phys. Rev. Lett, \textbf{101}, 097205 (2008)]. A careful electronic structure analysis suggests that coupling between Dy and Fe spin sublattices is mediated by Dy-$d$ and O-$2p$ hybridization. Our results are robust with respect to the different computational schemes used for $d$ and $f$ localized states, such as the DFT+$U$ method, the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional and the GW approach. Our findings indicate that the interaction between the $f$ and $d$ sublattice might be used to tailor ferroelectric and magnetic properties of multiferroic compounds.Comment: 6 pages, 4 figures-Revised versio

ZnSe/GaAs(001) heterostructures with defected interfaces: structural, thermodynamic and electronic properties

We have performed accurate \emph{ab--initio} pseudopotential calculations for the structural and electronic properties of ZnSe/GaAs(001) heterostructures with interface configurations accounting for charge neutrality prescriptions. Beside the simplest configurations with atomic interdiffusion we consider also some configurations characterized by As depletion and cation vacancies, motivated by the recent successfull growth of ZnSe/GaAs pseudomorphic structures with minimum stacking fault density characterized by the presence of a defected (Zn,Ga)Se alloy in the interface region. We find that--under particular thermodynamic conditions--some defected configurations are favoured with respect to undefected ones with simple anion or cation mixing, and that the calculated band offsets for some defected structures are compatible with those measured. Although it is not possible to extract indications about the precise interface composition and vacancy concentration, our results support the experimental indication of (Zn,Ga)Se defected compounds in high-quality ZnSe/GaAs(001) heterojunctions with low native stacking fault density. The range of measured band offset suggests that different atoms at interfaces rearrange, with possible presence of vacancies, in such a way that not only local charges but also ionic dipoles are vanishing.Comment: 26 pages. 5 figures, revised version, in press (Physical Review B

The shortcomings of semi-local and hybrid functionals: what we can learn from surface science studies

A study of the adsorption of CO on late 4d and $5d$ transition metal (111) surfaces (Ru, Rh, Pd, Ag, Os, Ir, and Pt) considering atop and hollow site adsorption is presented. The applied functionals include the gradient corrected PBE and BLYP functional, and the corresponding hybrid Hartree-Fock density functionals HSE and B3LYP. We find that PBE based hybrid functionals (specifically HSE) yield, with the exception of Pt, the correct site order on all considered metals, but they also considerably overestimate the adsorption energies compared to experiment. On the other hand, the semi-local BLYP functional and the corresponding hybrid functional B3LYP yield very satisfactory adsorption energies and the correct adsorption site for all surfaces. We are thus faced with a Procrustean problem: the B3LYP and BLYP functionals seem to be the overall best choice for describing adsorption on metal surfaces, but they simultaneously fail to account well for the properties of the metal, vastly overestimating the equilibrium volume and underestimating the atomization energies. Setting out from these observations, general conclusions are drawn on the relative merits and drawbacks of various semi-local and hybrid functionals. The discussion includes a revised version of the PBE functional specifically optimized for bulk properties and surface energies (PBEsol), a revised version of the PBE functional specifically optimized to predict accurate adsorption energies (rPBE), as well as the aforementioned BLYP functional. We conclude that no semi-local functional is capable to describe all aspects properly, and including non-local exchange also only improves some, but worsens other properties.Comment: 12 pages, 6 figures; to be published in New Journal of Physic

Large ferroelectric polarization in the new double perovskite NaLaMnWO6_{6} induced by non-polar instabilities

Based on density functional theory calculations and group theoretical analysis, we have studied NaLaMnWO$_{6}$ compound which has been recently synthesized [Phys. Rev. B 79, 224428 (2009)] and belongs to the $AA'BB'{\rm O}_{6}$ family of double perovskites. At low temperature, the structure has monoclinic $P2_{1}$ symmetry, with layered ordering of the Na and La ions and rocksalt ordering of Mn and W ions. The Mn atoms show an antiferromagnetic (AFM) collinear spin ordering, and the compound has been reported as a potential multiferroic. By comparing the low symmetry structure with a parent phase of $P4/nmm$ symmetry, two distortion modes are found dominant. They correspond to MnO$_{6}$ and WO$_{6}$ octahedron \textit{tilt} modes, often found in many simple perovskites. While in the latter these common tilting instabilities yield non-polar phases, in NaLaMnWO$_{6}$ the additional presence of the $A$-$A^{'}$ cation ordering is sufficient to make these rigid unit modes as a source of the ferroelectricity. Through a trilinear coupling with the two unstable tilting modes, a significant polar distortion is induced, although the system has no intrinsic polar instability. The calculated electric polarization resulting from this polar distortion is as large as $\sim$ 16 ${\mu}{\rm C/cm^{2}}$. Despite its secondary character, this polarization is coupled with the dominant tilting modes and its switching is bound to produce the switching of one of two tilts, enhancing in this way a possible interaction with the magnetic ordering. The transformation of common non-polar purely steric instabilities into sources of ferroelectricity through a controlled modification of the parent structure, as done here by the cation ordering, is a phenomenon to be further explored.Comment: Physical Chemistry Chemical physics (in press