Catalysis has entered everyday life through a number of technological processes relying on different catalytic systems. The increasing demand for such systems requires rationalization of the use of their expensive components, like noble metal catalysts. As such, a catalyst with low noble metal concentration, in which each one of the noble atoms is active, would reach the lowest price possible. Nevertheless, there are no reactivity descriptors outlined for this type of low coordinated supported atoms. Using DFT calculations, we consider three diverse systems as models of single atom catalysts. We investigate monomers and bimetallic dimers of Ru, Rh, Pd, Ir and Pt on MgO(001), Cu adatom on thin Mo(001)-supported films (NaF, MgO and ScN) and single Pt adatoms on oxidized graphene surfaces. Reactivity of these metal atoms was probed by CO. In each case we see the interaction through the donation-backdonation mechanism. In some cases the CO adsorption energies can be linked to the position of the d-band center and the charge of the adatom. Higher positioned d-band center and less charged supported single atoms bind CO weaker. Also, in some cases metal atoms less strongly bonded to the substrate bind CO more strongly. The results suggest that the identification of common activity descriptor(s) for single metal atoms on foreign supports is a difficult task with no unique solution. However, it is also suggested that the stability of adatoms and strong anchoring to the support are prerequisites for the application of descriptor-based search for novel single atom catalysts.

For the Cu(100), Cu(110), and Cu(111) surfaces varying asymmetric line shapes are found for the atomic 3d84s2 multiplet two-hole final state binding energies reached in MVV Auger photoelectron coincidence spectroscopy. Higher asymmetry for Cu(111) and Cu(110) in comparison to Cu(100) is caused by reduced dynamic screening for Cu(111) and Cu(110) in contrast to free electron like Cu(100). This is a consequence of the surface projected band gaps in Cu(111) and Cu(110) not present in Cu(100). We describe the distinct tailing in the experimental line shapes of the three Cu surfaces with first principles calculations of layer-dependent two-hole binding energy shifts, depth-dependent intensity distribution and Doniach-Sunjic asymmetry parametrization. These fundamental insights into the surface-specific electronic structure can advance the understanding of structure-reactivity relationships in Copper-based surfaces and catalysts.

We show results from first-principles calculations for cerium at very high compressions. These reveal a most remarkable behavior in a material; depending on atomic volume, cerium adopts three distinct face-centered cubic (fcc) phases driven by different physical mechanisms. The two well-known a and phases are vigorously debated in the literature, but we focus on the a phase as a metal with delocalized character of the 4f electron. The ultimate high compression fcc phase, here named ω, is driven partly by electrostatics. Our density-functional theory (DFT) study excellently reproduces the experimentally known compression behavior of cerium up to a few Mbar but goes beyond those pressures with structural transitions to tetragonal, hexagonal, and cubic (fcc) phases occurring before 100 Mbar (10000 GPa or 10 TPa). The 4f-electron contribution to the chemical bonding is shown to rule phase transitions and compressibility. The change of 4f occupation nicely explains the pressure dependence of the structural axial ratio in the tetragonal phase. At very high pressure, structures known at low pressures return because of band broadening, electrostatic ion repulsion, and an increase in hybridization between states that under normal conditions can be considered core (atomic like) states and the valence-band states.

The aim of this paper is to demonstrate how economic growth stimulates business-service providers to develop new service varieties, which, in turn, enhance the productivity of business-service buyers. This creates a coevolutionary process where service suppliers and customers interact, leading to an increase in the number of differentiated service offerings. We introduce a framework for local economies, wherein business-service sectors evolve in response to local demand potential, while non-business-service sectors grow based on each economy’s supply potential. Business service growth is more rapid in local economies with higher demand potential, while non-business-service sectors expand faster in areas where the business-service supply potential is greater. A key assumption is that business service firms operate in a monopolistic competition environment, where an increase in business-service capacity leads to an expansion in the variety of services offered. This, in turn, enhances the diversity of service offerings in municipalities with strong demand potential. Additionally, service providers not only deliver innovation-related information to client firms but also unintentionally disseminate knowledge within the region, fostering knowledge spillovers among firms.

This paper analyzes city system dynamics, based on a theoretical framework relating interaction potentials to agglomeration economies and density externalities. It employs new historical time series data on population size of cities in Sweden over two centuries (1810–2010) and introduces two schematic growth factors: (i) the intra-city potential and (ii) the extra-city potential located in in rings encircling each city. The first factor is measured by each city’s population size, while the second is a vector of distance-discounted population size for each of a city’s urban rings. In this way, we can explain a city’s growth as a function of its interaction potential inside the city, as well as inside the first, second hand third ring. A robust finding is that cities with large ring potentials follow different development paths than those with small ring potentials. We also find clear evidence of structural change between the two centuries 1810–1910 and 1910–2010. In the first period, city growth is positively impacted by the size of the intra-city potential, whereas the same potential dampens or reduces the growth in the second period. Moreover, the ring potentials outside the city tend to switch from having negative growth stimulation in the first period to having positive stimulation in the second period.

The development of new electrochromic materials and devices, like smart windows, has an enormous impact on the energy efficiency of modern society. One of the crucial materials in this technology is nickel oxide. Ni-deficient NiO shows anodic electrochromism, whose mechanism is still under debate. We use DFT+U calculations to show that Ni vacancy generation results in the formation of hole polarons localized at the two oxygens next to the vacancy. In the case of NiO bulk, upon Li insertion or injection of an extra electron into Ni-deficient NiO, one hole gets filled, and the hole bipolaron is converted into a hole polaron well-localized at one O atom, resulting from the transition between oxidized (colored) to reduced (bleached) state. In the case of the Ni-deficient NiO(001) surface, the qualitatively same picture is obtained upon embedding Li, Na, and K into the Ni surface vacancy, reinforcing the conclusion that the electron injection, resulting in the filling of the hole states, is responsible for the modulation of the optical properties of NiO. Hence, our results suggest a new mechanism of Ni-deficient NiO electrochromism not related to the change of the Ni oxidation states, i.e., the Ni2+/Ni3+ transition, but based on the formation and annihilation of hole polarons in oxygen p-states.

High temperature properties of ceria surfaces are important for many applications. Here we report the temperature dependences of surface energy for the (111) and (110) CeO 2 obtained in the framework of the extended two-stage up-sampled thermodynamic integration using Langevin dynamics (TU-TILD). The method was used together with machine-learning potentials called moment tensor potentials (MTPs), which were fitted to the results of the ab initio MD calculations for (111) and (110) CeO 2 at different temperatures. The parameters of MTPs training and fitting were tested and the optimal algorithm for the ceria systems was proposed. We found that the temperature increases from 0 K to 2100 K led to the decrease of the Helmholtz free energy of (111) CeO 2 from 0.78 J/m 2 to 0.64 J/m 2 . The energy of (110) CeO 2 dropped from 1.19 J/m 2 at 0 K to 0.92 J/m 2 at 1800 K. We show that it is important to take anharmonicity into account as simple consideration of volume expansion gives wrong temperature dependences of the surface energies.

justification Igor A. Pašti1,2*, Ana S. Dobrota1, Dmitri Migas,3,4,5 Börje Johansson,2,6 Natalia V. Skorodumova2* 1 University of Belgrade – Faculty of Physical Chemistry, Belgrade, Serbia 2 Department of Materials Science and Engineering, School of Industrial Engineering and Management, KTH – Royal Institute of Technology, Stockholm, Sweden 3 Belarusian State University of Informatics and Radioelectronics, Minsk, Belarus 4 National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Moscow, Russia 5 Institute for Nuclear Problems of Belarusian State University, Minsk, Belarus 6 Department of Physics and Astronomy, Uppsala University, Uppsala, Sweden

The development of new electrochromic materials and devices, like smart windows, has an enormous impact on the energy efficiency of modern society. One of the crucial materials in this technology is nickel-oxide. Ni-deficient NiO shows anodic electrochromism whose mechanism is still under debate. Using DFT+U calculations, we show that Ni vacancy generation results in the formation of hole polarons localised at the two oxygens next to the vacancy. Upon Li insertion or injection of an extra electron into Ni-deficient NiO, one hole gets filled, and the hole bipolaron is converted into a hole polaron well-localized at one O atom. Furthermore, the calculated absorption coefficients demonstrate that Li insertion/extraction or rather the addition/removal of an extra electron into Ni-deficient NiO can lead to switching between the oxidized (colored) and the reduced (bleached) states. Hence, our results suggest a new mechanism of Ni-deficient NiO electrochromism not related to the Ni2+/Ni3+ transition but based on the formation and annihilation of hole polarons in oxygen p-states.

For carbon-based materials, in contrast to metal surfaces, a general relationship between strain and reactivity is not yet established, even though there are literature reports on strained graphene. Knowledge of such relationships would be extremely beneficial for understanding the reactivity of graphene-based surfaces and finding optimisation strategies which would make these materials more suitable for targeted applications. Here we investigate the effects of compressive and tensile strain (up to +/- 5%) on the structure, electronic properties and reactivity of pure, N-doped and P-doped graphene, using DFT calculations. We demonstrate the possibility of tuning the topology of the graphene surface by strain, as well as by the choice of the dopant atom. The reactivity of (doped) strained graphene is probed using H and Na as simple adsorbates of great practical importance. Strain can both enhance and weaken H and Na adsorption on (doped) graphene. In case of Na adsorption, a linear relationship is observed between the Na adsorption energy on P-doped graphene and the phosphorus charge. A linear relationship between the Na adsorption energy on flat graphene surfaces and strain is found. Based on the adsorption energies and electrical conductivity, potentially good candidates for hydrogen storage and sodiumion battery electrodes are discussed.

Understanding the reactions in M-O2 cells (M = Li or Na) is of great importance for further advancement of this promising technology. Computational modelling can be helpful along this way, but an adequate approach is needed to model such complex systems. We propose a new scheme for modelling processes in M-O2 cells, where reference energies are obtained from high-level theory, CCSD(T), while the interactions of reaction intermediates with catalyst surfaces are extracted from computationally less expensive DFT. The approach is demonstrated for the case of graphene-based surfaces as model catalysts in Li-O2 and Na-O2 cells using the minimum viable mechanism. B-doped graphene was identified as the best catalyst among considered surfaces, while pristine graphene performs poorly. Moreover, we show that the inclusion of dispersion corrections for DFT has a significant impact on calculated discharge and charge potentials and suggests that long-range dispersion interactions should always be considered when graphene-based materials are modelled as electrocatalysts. Finally, we offer general guidelines for designing new ORR catalysts for M-O2 cells in terms of the optimization of the interactions of catalyst surface with reaction intermediates.

Due to its unique physical and chemical properties, graphene is being considered as a promising material for energy conversion and storage applications. Introduction of functional groups and dopants on/in graphene is a useful strategy for tuning its properties. In order to fully exploit its potential, atomic-level understanding of its interaction with species of importance for such applications is required. We present a DFT study of the interaction of sodium atoms with epoxy-graphene and analyze how this interaction is affected upon doping with boron and nitrogen. We demonstrate how the dopants, combined with oxygen-containing groups alter the reactivity of graphene towards Na. Dopants act as attractors of epoxy groups, enhancing the sodium adsorption on doped oxygen-functionalized graphene when compared to the case of non-doped epoxy-graphene. Furthermore, by considering thermodynamics of the Na interaction with doped epoxy-graphene it has been concluded that such materials are good candidates for Na storage applications. Therefore, we suggest that controlled oxidation of doped carbon materials could lead to the development of advanced anode materials for rechargeable Na-ion batteries.

A density-functional theory investigation of the (100) and (110) surfaces of the body-centered cubic (bcc) Fe1-xbCrxb binary alloys, x(b) <= 15 at.%, is reported. The energies and segregation energies of these surfaces were calculated for chemically homogeneous concentration profiles and for Cr surface contents deviating from the nominal one of the bulk. The implications of these results for the surface alloy phase diagram are discussed. The surface chemistry of Fe-Cr(100) is characterized by a transition from Cr depletion to Cr enrichment in a critical bulk Cr composition window of 6 < x(b) < 9 at.%. In contrast, such threshold behavior of the surface Cr content is absent for Fe-Cr(110) and a nearly homogeneous Cr concentration profile is energetically favorable. The strongly suppressed surface-layer relaxation at both surfaces is shown to be of magnetic origin. The compressive, magnetic contribution to the surface relaxation stress is found to correlate well with the surface magnetic moment squared at both surface terminations. The stability of the Cr surface magnetic moments against bulk Cr content is clarified based on the surface electronic structure.

Interaction across geographical space is made possible and facilitated by different categories of interdependent networks, in urban regions and between regions. Changes in the pertinent spatial interdependencies will have consequences for spatial equilibrium patterns of urban regions and systems of regions, for the distribution of economic activities and population, for productivity change and growth, the size of urban regions and a series of social conditions. The described interdependence between fast and slow adjustments concern the development of new networks for capital markets, stock exchange organizations, mobile phone and Internet systems. In a superficial discussion one would definitely argue that globalization of markets is an issue that refers to a region's external markets. Short distance to external markets is of course not completely independent of geographical distance. The chapter also presents an overview of the key concepts discussed in this book.

This chapter examines a model structure that can project multiregional job and populations dynamics. Such dynamics include the attraction of firms and households to each functional region. The method to identify a region's internal market potential is in a fundamental way related to the delineation of the geographic extension of a functional region. A functional region is distinguished by its concentration of activities and of its infrastructure which facilitates particularly high factor mobility within its interaction borders. The resource-based theory has been challenged in recent decades by scale-based models, which explain location and trade patterns in a context of internal and external economies of scale. The chapter presents model for long term projection of employment and population growth in Swedish municipalities. In Sweden a functional region normally consists of several interacting municipalities. Such a group of municipalities is considered to be integrated into a region when the interaction between the municipalities is intensive.

In 2017, we discovered quaternary i-MAX phases atomically layered solids, where M is an early transition metal, A is an A group element, and X is C-with a ((M2/3M1/32)-M-1)(2)AC chemistry, where the M-1 and M-2 atoms are in-plane ordered. Herein, we report the discovery of a class of magnetic i-MAX phases in which bilayers of a quasi-2D magnetic frustrated triangular lattice overlay a Mo honeycomb arrangement and an Al Kagome lattice. The chemistry of this family is (Mo2/3RE1/3)(2)AlC, and the rare-earth, RE, elements are Ce, Pr, Nd, Sm, Gd, Tb, Dy, Ho, Er, Tm, and Lu. The magnetic properties were characterized and found to display a plethora of ground states, resulting from an interplay of competing magnetic interactions in the presence of magnetocrystalline anisotropy.

Using the Exact Muffin-Tin Orbitals method within the Perdew-Burke-Ernzerhof exchange-correlation approximation for solids and solid surfaces (PBEso1), we study the single crystal elastic constants of 4d transition metals (atomic number Z between 39 and 47) and their binary alloys in the body centered cubic (bcc) and face centered cubic (fcc) structures. Alloys between the first neighbors Z(Z + 1) and between the second neighbors Z(Z + 2) are considered. The lattice constants, bulk moduli and elastic constants are found in good agreement with the available experimental and theoretical data. It is shown that the correlation between the relative tetragonal shear elastic constant C-fcc'-2C(bcc)' and the structural energy difference between the fcc and bcc lattices Delta E is superior to the previously considered models. For a given crystal structure, the equiatomic Z(Z + 2) alloys turn out to have similar structural and elastic properties as the pure elements with atomic number (Z + 1). Furthermore, alloys with composition Z(1-x)(Z + 2)(x) possess similar properties as Z(1-2x)(Z + 1)(2x). The present theoretical data on the structural and the elastic properties of 4d transition metal alloys provides consistent input for coarse scale modeling of material properties.