Search

Lead vanadate PbVO3 is a polar crystal with a P4mm space group at ambient conditions. It is isostructural with the model soft-mode driven ferroelectric PbTiO3, but differs from it by the so-called 'supertetragonal' elongation of its unit cell. In this paper, we report a combined study of the lattice dynamics of PbVO3 by Raman spectroscopy at room temperature and first-principle calculations. All zone-center transverse optical (TO) phonon modes are identified by polarized, angle-dependent Raman spectroscopy and assigned as follows: E modes at 136, 269, 374 and 508 cm-1, A1 modes at 188, 429 and 874 cm-1 and B1 mode at 319 cm-1. The calculations confirm the experimental symmetry assignment and allow to obtain the longitudinal (LO) phonons wavenumbers. Besides, we analyze the mode eigenvectors in detail, in order to identify the atomic displacements associated with each mode and compare them with PbTiO3. In spite of their differences in chemistry and strain, the phonon eigenvectors are found to be remarkably comparable in both compounds. We discuss the position of the ferroelectric soft mode in PbVO3 as compared to PbTiO3. A sizeable splitting of the B1+E modes appears as a characteristic feature of supertetragonal phases. The peculiarity of the vanadyl V-O bond frequency in PbVO3 is also addressed., 11 pages, 8 figures, 6 Tables, Supplemental Information 3 figures and 6 Tables, under review

We present a general relation between the electromechanical couplings of infinitesimal strain and electric field to arbitrary order, measured at fixed voltage and at fixed electric field. We show that the improper response at fixed field can be written as the strain derivative of the $n^{\text{th}}$ order susceptibility tensor, and the proper response at fixed voltage drop can be written as the response at fixed field plus corrections for dilations and 90$^{\circ}$ rotations induced by strain. Our theory correctly reproduces the proper piezoelectric response and we go beyond with the electrostrictive response. We present first-principles calculations of the improper electrostrictive response at fixed field, and illustrate how the correction is used to obtain the proper response at fixed voltage. This distinction is of high importance given the recent interest in giant electrostrictors exhibiting electromechanical responses as large as the piezoelectric ones., Accepted version

Aurivillius oxides have been a research focus due to their ferroelectric properties, but by replacing oxide ions by fluoride, divalent magnetic cations can be introduced, giving Bi

The most important types of non-collinear magnetic orders that are realized in simple perovskite oxides are outlined in relation to multiferroicity. These orders are classified and rationalized in terms of a mimimal spin Hamiltonian, based on which the notion of spin-driven ferroelectricity is illustrated. These concepts find direct application in reference materials such as BiFeO3, GdFeO3 and TbMnO3 whose multiferroic properties are briefly reviewed.

Multiferroics are materials exhibiting the coexistence of ferroelectricity and ideally ferromagnetism. Unfortunately, most known magnetoelectric multiferroics combine ferroelectricity with antiferromagnetism or with weak ferromagnetism. Here, following previous theoretical predictions, we provide clear experimental indications that ferroelectricity can be induced by epitaxial tensile strain in the ferromagnetic simple binary oxide EuO. We investigate the ferroelectric phase transition using infrared reflectance spectroscopy, finding that the frequency of the soft optical phonon reduces with increasing tensile strain and decreasing temperature. We observe such a soft mode anomaly at 100 K in (EuO)2/(BaO)2 superlattices grown epitaxially on (LaAlO3)0.29-(SrAl1/2Ta1/2O3)0.71 substrates, which is a typical signature for a displacive ferroelectric phase transition. The EuO in this superlattice is nominally subjected to 6.4% biaxial tensile strain, i.e., 50% more than believed needed from previously published calculations. We interpret our results with new first-principles density functional calculations using a hybrid functional, which provides a better quantitative agreement with experiment than the previously used local-density approximation and generalized gradient approximation functionals., Communications Materials, 1 (1), ISSN:2662-4443

Funder: Deutsche Forschungsgemeinschaft; Id: http://dx.doi.org/10.13039/501100001659, Outstanding photovoltaic (PV) materials combine a set of advantageous properties including large optical absorption and high charge carrier mobility, facilitated by small effective masses. Halide perovskites (ABX3, where X = I, Br, or Cl) are among the most promising PV materials. Their optoelectronic properties are governed by the B X bond, which is responsible for the pronounced optical absorption and the small effective masses of the charge carriers. These properties are frequently attributed to the ns2 configuration of the B atom, i.e., Pb 6s2 or Sn 5s2 (“lone‐pair”) states. The analysis of the PV properties in conjunction with a quantum‐chemical bond analysis reveals a different scenario. The B X bond differs significantly from ionic, metallic, or conventional 2c 2e covalent bonds. Instead it is better regarded as metavalent, since it shares about one p‐electron between adjacent atoms. The resulting σ‐bond, formally a 2c 1e bond, is half‐filled, causing pronounced optical absorption. Electron transfer between B and X atoms and lattice distortions open a moderate bandgap resulting in charge carriers with small effective masses. Hence, metavalent bonding explains favorable PV properties of halide perovskites, as summarized in a map for different bond types, which provides a blueprint to design PV materials.

Out-of-plane polar domain structures have recently been discovered in strained and twisted bilayers of inversion symmetry broken systems such as hexagonal boron nitride. Here we show that this symmetry breaking also gives rise to an in-plane component of polarization, and the form of the total polarization is determined purely from symmetry considerations. The in-plane component of the polarization makes the polar domains in strained and twisted bilayers topologically non-trivial, forming a network of merons and antimerons (half-skyrmions and half-antiskyrmions). For twisted systems, the merons are of Bloch type whereas for strained systems they are of N\'eel type. We propose that the polar domains in strained or twisted bilayers may serve as a platform for exploring topological physics in layered materials and discuss how control over topological phases and phase transitions may be achieved in such systems., Comment: Accepted version

In this work a new method for the calculation of the electrostrictive properties of materials using density functional theory is presented. The method relies on the thermodynamical equivalence, in a dielectric, of the quadratic mechanical responses (stress or strain) to applied electric stimulus (electric or polarization fields) to the strain or stress dependence of its dielectric susceptibility or stiffness tensors. Comparing with current finite-field methodologies for the calculation of electrostriction, it is demonstrated that this presented methodology offers significant advantages of efficiency, robustness, and ease of use. These advantages render tractable the high throughput theoretical investigation into the largely unknown electrostrictive properties of materials, and the microscopic origins of giant electrostriction.

Ferroelectric Rashba semiconductors (FERSC), in which Rashba spin-splitting can be controlled and reversed by an electric field, have recently emerged as a new class of functional materials useful for spintronic applications. The development of concrete devices based on such materials is, however, still hampered by the lack of robust FERSC compounds. Here, we show that the coexistence of large spontaneous polarization and sizeable spin–orbit coupling is not sufficient to have strong Rashba effects and clarify why simple ferroelectric oxide perovskites with transition metal at the B-site are typically not suitable FERSC candidates. By rationalizing how this limitation can be by-passed through band engineering of the electronic structure in layered perovskites, we identify the Bi$${}_{2}$$ 2 WO$${}_{6}$$ 6 Aurivillius crystal as a robust ferroelectric with large and reversible Rashba spin-splitting, that can even be substantially doped without losing its ferroelectric properties. Importantly, we highlight that a unidirectional spin–orbit field arises in layered Bi$${}_{2}$$ 2 WO$${}_{6}$$ 6 , resulting in a protection against spin-decoherence.

Spin reorientation and magnetization reversal are two important features of the rare-earth orthorhombic perovskites ($RM{\mathrm{O}}_{3}$) that have attracted a lot of attention, though their exact microscopic origin has eluded researchers. Here, using density functional theory and classical atomistic spin dynamics we build a general Heisenberg magnetic model that allows to explore the whole phase diagram of the chromite and ferrite compounds and to scrutinize the microscopic mechanism responsible for spin reorientations and magnetization reversals. We show that the occurrence of a magnetization reversal transition depends on the relative strength and sign of two interactions between rare-earth and transition-metal atoms: superexchange and Dzyaloshinskii-Moriya. We also conclude that the presence of a smooth spin reorientation transition between the so-called ${\mathrm{\ensuremath{\Gamma}}}_{4}$ and the ${\mathrm{\ensuremath{\Gamma}}}_{2}$ phases through a coexisting region, and the temperature range in which it occurs, depends on subtle balance of metal-metal (superexchange and Dzyaloshinskii-Moriya) and metal--rare-earth (Dzyaloshinsky-Moriya) couplings. In particular, we show that the intermediate coexistence region occurs because the spin sublattices rotate at different rates.

Combining ferroelectricity with other properties such as visible light absorption or long-range magnetic order requires the discovery of new families of ferroelectric materials. Here, through the analysis of a high-throughput database of phonon band structures, we identify a new structural family of anti-Ruddlesden-Popper phases A$_4$X$_2$O (A=Ca, Sr, Ba, Eu, X=Sb, P, As, Bi) showing ferroelectric and anti-ferroelectric behaviors. The discovered ferroelectrics belong to the new class of hyperferroelectrics which polarize even under open-circuit boundary conditions. The polar distortion involves the movement of O anions against apical A cations and is driven by geometric effects resulting from internal chemical strains. Within this new structural family, we show that Eu$_4$Sb$_2$O combines coupled ferromagnetic and ferroelectric order at the same atomic site, a very rare occurrence in materials physics., 21 pages, 4 figures with supplemental information attached

Raman spectroscopy is an advantageous method for studying the local structure of materials, but the interpretation of measured spectra is complicated by the presence of oblique phonons in polycrystals of polar materials. Whilst group theory considerations and standard ab initio calculations are helpful, they are often valid only for single crystals. In this paper, we introduce a method for computing Raman spectra of polycrystalline materials from first principles. We start from the standard approach based on the (Placzek) rotation invariants of the Raman tensors and extend it to include the effect of the coupling between the lattice vibrations and the induced electric field, and the electro-optic contribution, relevant for polar materials like ferroelectrics. As exemplified by applying the method to rhombohedral BaTiO3, AlN, and LiNbO3, such an extension brings the simulated Raman spectrum to a much better correspondence with the experimental one. Additional advantages of the method are that it is general, permits automation, and thus can be used in high-throughput fashion., 23 pages, 4 figures

Multiferroics are a unique class of materials where magnetic and ferroelectric orders coexist. The research on multiferroics contributes significantly to the fundamental understanding of the strong correlations between different material degrees of freedom and provides an energy-efficient route toward the electrical control of magnetism. While multiple ABO3 oxide perovskites have been identified as being multiferroic, their magnetoelectric coupling strength is often weak, necessitating the material search in different compounds. Here, we report the observation of room-temperature multiferroic orders in multi-anion SrNbO3-xNx thin films. In these samples, the multi-anion state enables the room-temperature ferromagnetic ordering of the Nb d-electrons. Simultaneously, we find ferroelectric responses that originate from the structural symmetry breaking associated with both the off-center displacements of Nb and the geometric displacements of Sr, depending on the relative O-N arrangements within the Nb-centered octahedra. Our findings not only diversify the available multiferroic material pool but also demonstrate a new multiferroism design strategy via multi-anion engineering.

We report a theoretical study of the non-linear magnetoelectric response of GdFeO$_3$ through an analytical approach combined with a Heisenberg model which is fitted against first-principles calculations. Our theory reproduces the non-linear change of polarization under applied magnetic field reported experimentally such that it allows to analyze the origin of the large responses in the different directions. We show that the non-linear character of the response in these materials originates from the fact that the antiferromagnetic order of Gd atoms changes non-linearly with respect to the applied magnetic field. Our model can be generalized to other materials in which the antiferromagnetic ordering breaks inversion symmetry.

The structural and electronic properties of ${\mathrm{KTaO}}_{3}$/${\mathrm{KZnF}}_{3}$ and ${\mathrm{KTaO}}_{3}$/${\mathrm{KNiF}}_{3}$ oxyfluoride superlattices are studied from first-principles density functional theory calculations. We highlight, that beyond a critical layer thickness, these systems exhibit an insulator to metal transition that gives rise to the appearance of two-dimensional electron and hole gas, confined both, due the band alignment, within the oxide layer. The origin of the insulator to metal transition is related to the polar discontinuity at the interfaces. The behavior is discussed in terms of a simple electrostatic model and compared to that of the prototypical ${\mathrm{LaAlO}}_{3}$/${\mathrm{SrTiO}}_{3}$ oxide system. The magnetic properties ${\mathrm{KTaO}}_{3}$/${\mathrm{KNiF}}_{3}$ superlattices are further discussed, revealing a sizable Rashba-type spin splitting at these interfaces, much larger than in similar oxide/oxide systems.

The Dion-Jacobson (DJ) family of perovskite-related materials have recently attracted interest due to their polar structures and properties, resulting from hybrid-improper mechanisms for ferroelectricity in n=2 systems, and from proper mechanisms in n=3 CsBi\(_2\)Ti\(_2\)NbO\(_{10}\). We report here a combined experimental and computational study on analogous n=3 CsLn\(_2\)Ti\(_2\)NbO\(_{10}\) (Ln = La, Nd) materials. Density functional theory calculations reveal the shallow energy landscape in these systems and give an understanding of the competing structural models suggested by neutron and electron diffraction studies. The structural disorder resulting from the shallow energy landscape breaks inversion symmetry at a local level, consistent with the observed second-harmonic generation. This study reveals the potential to tune between proper and hybrid-improper mechanisms by composition in the DJ family. The disorder and shallow energy landscape have implications for designing functional materials with properties reliant on competing low-energy phases such as relaxors and antiferroelectrics.

We present the use and implementation of the firefly algorithm to help in scanning the multiple metastable minima of orbital occupations in density functional theory (DFT) plus Hubbard U correction and to identify the ground state occupations in strongly correlated materials. We show the application of this implementation with the Abinit code on KCoF3 and UO2 crystals, which are typical d and f electron systems with numerous occupation minima. We demonstrate the validity and performance of the method by comparing with previous methodologies. The method is general and can be applied to any code using constrained occupation matrices.

Ferroelectric tungsten-based Aurivillius oxides are naturally stable superlattice structures, in which A-site deficient perovskite blocks [WnO3n+1]−2 (n=1,2,3,⋯) interleave with fluorite-like bismuth oxide layers [Bi2O2]+2 along the c-axis. In the n=2Bi2W2O9 phase, an in-plane antipolar distortion dominates but there has been controversy as to the ground-state symmetry. Here we show, using a combination of first-principles density functional theory calculations and experiments, that the ground state is a nonpolar phase of Pnab symmetry. We explore the energetics of metastable phases and the potential for antiferroelectricity in this n=2 Aurivillius phase.

Escherichia coli is one of the major pathogens in humans and animals causing localized and systemic infections, which often lead to acute inflammation, watery diarrhea, and hemorrhagic colitis. Bacterial lipopolysaccharide (LPS) and Shiga exotoxins (Stx) are mostly responsible for such clinical signs. Therefore, highly effective treatment of E. coli infections should include both eradication of bacteria and neutralization of their toxins. Here, for the first time, we compared the in vitro ability of common antibiotics to decrease LPS- and Stx-mediated cytotoxicity: colistin, amoxicillin (used separately or combined), enrofloxacin, and its metabolite ciprofloxacin. Three experimental scenarios were realized as follows: (a) the direct effect of antibiotics on endotoxin, (b) the effect of antibiotic treatment on LPS-mediated cytotoxicity in an experiment mimicking "natural infection," (c) the effect of antibiotics to decrease Stx2e-mediated cytotoxicity. Two cell lines, A549 and Vero cells, were used to perform cytotoxic assays with the methyl tetrazolium (MTT) and lactate dehydrogenase leakage (LDH) methods, respectively. Colistin and amoxicillin, especially used in combination, were able to attenuate LPS toxic effect, which was reflected by increase in A549 cell viability. In comparison with other antibiotics, the combination of colistin and amoxicillin exhibited the highest boster or additive effect in protecting cells against LPS- and Stx2e-induced toxicity. In summary, in comparison with fluoroquinolones, the combination of colistin and amoxicillin at concentrations similar to those achieved in plasma of treated animals exhibited the highest ability to attenuate LPS- and Stx2e-mediated cytotoxicity.

Third-generation photovoltaic (PV) materials combine many advantageous properties, including a high optical absorption together with a large charge carrier mobility, facilitated by small effective masses. Halide perovskites (ABX3, where X = I, Br or Cl) appear to be the most promising third-generation PV materials at present. Their opto-electronic properties are governed by the B-X bond. A quantum-chemical bond analysis reveals that this bond differs significantly from ionic, metallic or covalent bonds. Instead, it is better regarded as metavalent, since it shares approximately one p-electron between adjacent atoms. The resulting s–bond is half-filled, which causes pronounced optical absorption. Electron transfer and lattice distortions open a moderate band gap, resulting in charge carriers with small effective masses. Hence metavalent bonding explains the favorable PV properties of halide perovskites. This is summarized in a map for different bond types, which provides a blueprint to design third-generation PV materials.

Paper published as part of the special topic on Electronic Structure Software, ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperaturedependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library LIBPAW. ABINIT has strong links with many other software projects that are briefly mentioned., This work (A.H.R.) was supported by the DMREF-NSF Grant No. 1434897, National Science Foundation OAC-1740111, and U.S. Department of Energy DE-SC0016176 and DE-SC0019491 projects. N.A.P. and M.J.V. gratefully acknowledge funding from the Belgian Fonds National de la Recherche Scientifique (FNRS) under Grant No. PDR T.1077.15-1/7. M.J.V. also acknowledges a sabbatical “OUT” grant at ICN2 Barcelona as well as ULiège and the Communauté Française de Belgique (Grant No. ARC AIMED G.A. 15/19-09). X.G. and M.J.V. acknowledge funding from the FNRS under Grant No. T.0103.19-ALPS. X.G. and G.-M. R. acknowledge support from the Communauté française de Belgique through the SURFASCOPE Project (No. ARC 19/24-057). X.G. acknowledges the hospitality of the CEA DAM-DIF during the year 2017. G.H. acknowledges support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division under Contract No. DE-AC02-05-CH11231 (Materials Project Program No. KC23MP). The Belgian authors acknowledge computational resources from supercomputing facilities of the University of Liège, the Consortium des Equipements de Calcul Intensif (Grant No. FRS-FNRS G.A. 2.5020.11), and Zenobe/CENAERO funded by the Walloon Region under Grant No. G.A. 1117545. M.C. and O.G. acknowledge support from the Fonds de Recherche du Québec Nature et Technologie (FRQ-NT), Canada, and the Natural Sciences and Engineering Research Council of Canada (NSERC) under Grant No. RGPIN-2016-06666. The implementation of the libpaw library (M.T., T.R., and D.C.) was supported by the ANR NEWCASTLE project (Grant No. ANR-2010-COSI-005-01) of the French National Research Agency. M.R. and M.S. acknowledge funding from Ministerio de Economia, Industria y Competitividad (MINECO-Spain) (Grants Nos. MAT2016-77100-C2-2-P and SEV-2015-0496) and Generalitat de Catalunya (Grant No. 2017 SGR1506). This work has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation program (Grant Agreement No. 724529). P.G. acknowledges support from FNRS Belgium through PDR (Grant No. HiT4FiT), ULiège and the Communauté française de Belgique through the ARC project AIMED, the EU and FNRS through M.ERA.NET project SIOX, and the European Funds for Regional Developments (FEDER) and the Walloon Region in the framework of the operational program “Wallonie-2020.EU” through the project Multifunctional thin films/LoCoTED. The Flatiron Institute is a division of the Simons Foundation. A large part of the data presented in this paper is available directly from the Abinit Web page www.abinit.org. Any other data not appearing in this web page can be provided by the corresponding author upon reasonable request.

Abinit is a material- and nanostructure-oriented package that implements density-functional theory (DFT) and many-body perturbation theory (MBPT) to find, from first principles, numerous properties including total energy, electronic structure, vibrational and thermodynamic properties, different dielectric and non-linear optical properties, and related spectra. In the special issue to celebrate the 40th anniversary of CPC, published in 2009, a detailed account of Abinit was included [Gonze et al. (2009)], and has been amply cited. The present article comes as a follow-up to this 2009 publication. It includes an analysis of the impact that Abinit has had, through for example the bibliometric indicators of the 2009 publication. Links with several other computational materials science projects are described. This article also covers the new capabilities of Abinit that have been implemented during the last three years, complementing a recent update of the 2009 article published in 2016. Physical and technical developments inside the abinit application are covered, as well as developments provided with the Abinit package, such as the multibinit and a-tdep projects, and related Abinit organization developments such as AbiPy . The new developments are described with relevant references, input variables, tests, and tutorials. Program summary Program Title: Abinit Program Files doi: http://dx.doi.org/10.17632/csvdrr4d68.1 Licensing provisions: GPLv3 Programming language: Fortran2003, Python Journal reference of previous version: X .Gonze et al, Comput. Phys. Commun. 205 (2016) 106–131 Does the new version supersede the previous version?: Yes. The present 8.10.3 version is now the up-to-date stable version of abinit , and supercedes the 7.10.5 version. Reasons for the new version: New developments Summary of revisions: • Many new capabilities of the main abinit application, related to density-functional theory, density-functional perturbation theory, GW, the Bethe-Salpeter equation, dynamical mean-field theory, etc. • New applications in the package: multibinit (second-principles calculations)and tdep (temperature-dependent properties) Nature of problem: Computing accurately material and nanostructure properties: electronic structure, bond lengths, bond angles, primitive cell, cohesive energy, dielectric properties, vibrational properties, elastic properties, optical properties, magnetic properties, non-linear couplings, electronic and vibrational lifetimes, etc. For large-scale systems, second-principles calculations, building upon the first-principles results, are also possible. Solution method: Software application based on density-functional theory and many-body perturbation theory, pseudopotentials, with plane waves or wavelets as basis functions. Different real-time algorithms are implemented for second-principles calculations.

In parallel with the revival of interest for magneto-electric multiferroic materials in the beginning of the century, first-principles simulations have grown incredibly in efficiency during the last two decades. Density functional theory calculations, in particular, have so become a must-have tool for physicists and chemists in the multiferroic community. While these calculations were originally used to support and explain experimental behaviour, their interest has progressively moved to the design of novel magneto-electric multiferroic materials. In this article, we mainly focus on oxide perovskites, an important class of multifunctional material, and review some significant advances to which contributed first-principles calculations. We also briefly introduce the various theoretical developments that were at the core of all these advances.

Despite of the tremendous efforts spent on the oxygen vacancy migration in determining the property optimization of oxygen-vacancy enrichment transition metal oxides, few has focused on their dynamic behaviors non-equilibrium states. In this work, we performed multi-timescale ultrafast X-ray diffraction measurements by using picosecond synchrotron X-ray pulses and femtosecond table-top X-ray pulses to monitor the structural dynamics in the oxygen-vacancy ordered SrCoO2.5 thin films. A giant photoinduced strain ({\Delta}c/c > 1%) was observed, whose distinct correlation with the pump photon energy indicates a non-thermal origin of the photoinduced strain. The sub-picosecond resolution X-ray diffraction reveals the formation and propagation of the coherent acoustic phonons inside the film. We also simulate the effect of photoexcited electron-hole pairs and the resulting lattice changes using the Density Function Theory method to obtain further insight on the microscopic mechanism of the measured photostriction effect. Comparable photostrictive responses and the strong dependence on excitation wavelength are predicted, revealing a bonding to anti-bonding charge transfer or high spin to intermediate spin crossover induced lattice expansion in the oxygen-vacancy films., Comment: 12 pages, 4 figures, support material

Entropy is a fundamental thermodynamic quantity that is a measure of the accessible microstates available to a system, with the stability of a system determined by the magnitude of the total entropy of the system. This is valid across truly mind boggling length scales - from nanoparticles to galaxies. However, quantitative measurements of entropy change using calorimetry are predominantly macroscopic, with direct atomic scale measurements being exceedingly rare. Here for the first time, we experimentally quantify the polar configurational entropy (in meV/K) using sub-\r{a}ngstr\"{o}m resolution aberration corrected scanning transmission electron microscopy. This is performed in a single crystal of the prototypical ferroelectric $\mathsf{LiNbO_3}$ through the quantification of the niobium and oxygen atom column deviations from their paraelectric positions. Significant excursions of the niobium - oxygen polar displacement away from its symmetry constrained direction is seen in single domain regions which increases in the proximity of domain walls. Combined with first principles theory plus mean field effective Hamiltonian methods, we demonstrate the variability in the polar order parameter, which is stabilized by an increase in the magnitude of the configurational entropy. This study presents a powerful tool to quantify entropy from atomic displacements and demonstrates its dominant role in local symmetry breaking at finite temperatures in classic, nominally Ising ferroelectrics., Comment: 23 pages, 21 figures (5 main, 16 supplemental)

A buffer layer formed by depositing a ½ monolayer of Sr on Si(0 0 1) is known to passivate the Si surface, while its surface structure constitutes a suitable template for the integration of various functional oxides with the existing Si platform. We used Pulsed Laser Deposition (PLD) to prepare a Sr/Si(0 0 1)(1 × 2) surface and analysed it using in-situ Reflection High-Energy Diffraction (RHEED) in combination with low-temperature Scanning Tunneling Microscopy (STM). The STM images reveal an atomically ordered surface with terraces composed of one-dimensional (1D) chains running along perpendicular directions on neighbouring terraces. The 1D chains are separated by 0.78 nm and exhibit a low-amplitude corrugation with a period of 0.39 nm. The measured values agree well with the size of the (1 × 2) unit cell observed for similar MBE-grown surfaces, while the density of the surface defects is somewhat higher in the presented case. According to simulated STM images based on DFT calculations, two types of surface defects were identified and explored: arrays of Sr vacancies and Sr adatoms. These results show that PLD can offer precise control for the preparation of high-quality Sr-buffered Si(0 0 1) surfaces.

This work explores the possibility of tuning the structural chemistry of A-site deficient perovskite materials Ln1/3NbO3, and the balance between proper and hybrid-improper mechanisms for polar behaviour in these systems.

We extend the density functional perturbation theory formalism to the case of non-collinear magnetism. The main problem comes with the exchange-correlation (XC) potential derivatives, which are the only ones that are affected by the non-collinearity of the system. Most of the present XC functionals are constructed at the collinear level, such that the off-diagonal (containing magnetization densities along $x$ and $y$ directions) derivatives cannot be calculated simply in the non-collinear framework. To solve this problem, we consider here possibilities to transform the non-collinear XC derivatives to a local collinear basis, where the $z$ axis is aligned with the local magnetization at each point. The two methods we explore are i) expanding the spin rotation matrix as a Taylor series, ii) evaluating explicitly the XC for the local density approximation through an analytical expression of the expansion terms. We compare the two methods and describe their practical implementation. We show their application for atomic displacement and electric field perturbations at the second order, within the norm-conserving pseudopotential methods.

The PyProcar Python package plots the band structure and the Fermi surface as a function of site and/or s,p,d,f - projected wavefunctions obtained for each $k$-point in the Brillouin zone and band in an electronic structure calculation. This can be performed on top of any electronic structure code, as long as the band and projection information is written in the PROCAR format, as done by the VASP and ABINIT codes. PyProcar can be easily modified to read other formats as well. This package is particularly suitable for understanding atomic effects into the band structure, Fermi surface, spin texture, etc. PyProcar can be conveniently used in a command line mode, where each one of the parameters define a plot property. In the case of Fermi-surfaces, the package is able to plot the surface with colors depending on other properties such as the electron velocity or spin projection. The mesh used to calculate the property does not need to be the same as the one used to obtain the Fermi surface. A file with a specific property evaluated for each $k$-point in a $k-$mesh and for each band can be used to project other properties such as electron-phonon mean path, Fermi velocity, electron effective mass, etc. Another existing feature refers to the band unfolding of supercell calculations into predefined unit cells., 32 pages, 14 figures, Submitted to "Computer Physics Communications", project url: "https://github.com/romerogroup/pyprocar"

In most noncollinear crystal magnets, the number of metastable states is quite large and any calculation that tries to predict the ground state can fall into one of the possible metastable phases. In this work, we generalize the population based meta-heuristic firefly algorithm to the problem of the noncollinear magnetic phase ground state prediction within density functional theory (DFT). We extend the different steps in the firefly algorithm to this specific problem by using polarized constrained DFT calculations, whereby using Lagrange multipliers the directions of the atom magnetic moments remain fixed. By locking the directions of the magnetic moments at each search iteration, the method allows one to explore the entire Born–Oppenheimer energy surface of existing and physically plausible noncollinear configurations present in a crystal. We demonstrate that the number of minima can be large, which restrains the use of exhaustive searches.