Spectroscopic signatures of nonpolarons: the case of diamond.

Polarons are quasi-particles made from electrons interacting with vibrations in crystal lattices. They derive their name from the strong electron-vibration polar interactions in ionic systems, that induce spectroscopic and optical signatures of such quasi-particles. In this paper, we focus on diamond, a non-polar crystal with inversion symmetry which nevertheless shows interesting signatures stemming from electron-vibration interactions, better denoted "nonpolaron" signatures in this case. The (non)polaronic effects are produced by short-range crystal fields, while long-range quadrupoles only have a small influence. The corresponding many-body spectral function has a characteristic energy dependence, showing a plateau structure that is similar to but distinct from the satellites observed in the polar Fröhlich case. We determine the temperature-dependent spectral function of diamond by two methods: the standard Dyson-Migdal approach, which calculates electron-phonon interactions within the lowest-order expansion of the self-energy, and the cumulant expansion, which includes higher orders of electron-phonon interactions. The latter corrects the nonpolaron energies and broadening, providing a more realistic spectral function, which we examine in detail for both conduction and valence band edges.

Electron-Phonon beyond Fröhlich: Dynamical Quadrupoles in Polar and Covalent Solids.

We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated.

Li diffusion in Si and LiSi: Nuclear quantum effects and anharmonicity.

The diffusion of Li in bulk Si and crystalline LiSi is investigated over a wide range of temperatures employing first-principles calculations based on density functional theory, transition state theory, and the kinetic Monte Carlo method. Nuclear quantum effects are incorporated by computing the vibrational spectrum and its effect on the effective energy barrier. The Li diffusion coefficient in bulk Si calculated with such quantum effects is ∼33% lower than the classical limit near room temperature due to higher effective energy barrier and tends to the classical limit at a high temperature (>1000 K). The presence of anharmonicity, estimated by the quasiharmonic approximation and the cBΩ model, increases the diffusion coefficient by ∼60%. For Li diffusion in LiSi with multiple vacancy jumps, we obtain an effective diffusion barrier of 0.27 eV ± 0.01 eV. In the Li-Si system, the quantum mechanical effects are only marginally significant at room temperature.

ABINIT: Overview and focus on selected capabilities.

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 "temperature-dependent 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.

Computer simulations of cluster impacts: effects of the atomic masses of the projectile and target.

Cluster secondary ion mass spectrometry is now widely used for the characterization of nanostructures. In order to gain a better understanding of the physics of keV cluster bombardment of surfaces and nanoparticles (NPs), the effects of the atomic masses of the projectile and of the target on the energy deposition and induced sputtering have been studied by means of molecular dynamics simulations. 10 keV C60 was used as a model projectile and impacts on both a flat polymer surface and a metal NP were analyzed. In the first case, the mass of the impinging carbon atoms was artificially varied and, in the second case, the mass of the NP atoms was varied. The results can be rationalized on the basis of the different atomic mass ratios of the projectile and target. In general, the emission is at its maximum, when the projectile and target have the same atomic masses. In the case of the supported NP, the emission of the underlying organic material increases as the atomic mass of the NP decreases. However, it is always less than that calculated for the bare organic surface, irrespective of the mass ratio. The results obtained with C60 impacts on the flat polymer are also compared to simulations of C60 and monoatomic Ga impacts on the NP.

Electronic Properties of Functionalized Diamanes for Field-Emission Displays.

Ultrathin diamond films, or diamanes, are promising quasi-2D materials that are characterized by high stiffness, extreme wear resistance, high thermal conductivity, and chemical stability. Surface functionalization of multilayer graphene with different stackings of layers could be an interesting opportunity to induce proper electronic properties into diamanes. Combination of these electronic properties together with extraordinary mechanical ones will lead to their applications as field-emission displays substituting original devices with light-emitting diodes or organic light-emitting diodes. In the present study, we focus on the electronic properties of fluorinated and hydrogenated diamanes with (111), (110), (0001), (101̅0), and (2̅110) crystallographic orientations of surfaces of various thicknesses by using first-principles calculations and Bader analysis of electron density. We see that fluorine induces an occupied surface electronic state, while hydrogen modifies the occupied bulk state and also induces unoccupied surface states. Furthermore, a lower number of layers is necessary for hydrogenated diamanes to achieve the convergence of the work function in comparison with fluorinated diamanes, with the exception of fluorinated (110) and (2̅110) films that achieve rapid convergence and have the same behavior as other hydrogenated surfaces. This induces a modification of the work function with an increase of the number of layers that makes hydrogenated (2̅110) diamanes the most suitable surface for field-emission displays, better than the fluorinated counterparts. In addition, a quasi-quantitative descriptor of surface dipole moment based on the Tantardini-Oganov electronegativity scale is introduced as the average of bond dipole moments between the surface atoms. This new fundamental descriptor is capable of predicting a priori the bond dipole moment and may be considered as a new useful feature for crystal structure prediction based on artificial intelligence.

Constrained Density Functional Theory: A Potential-Based Self-Consistency Approach.

Chemical reactions, charge transfer reactions, and magnetic materials are notoriously difficult to describe within Kohn-Sham density functional theory, which is strictly a ground-state technique. However, over the last few decades, an approximate method known as constrained density functional theory (cDFT) has been developed to model low-lying excitations linked to charge transfer or spin fluctuations. Nevertheless, despite becoming very popular due to its versatility, low computational cost, and availability in numerous software applications, none of the previous cDFT implementations is strictly similar to the corresponding ground-state self-consistent density functional theory: the target value of constraints (e.g., local magnetization) is not treated equivalently with atomic positions or lattice parameters. In the present work, by considering a potential-based formulation of the self-consistency problem, the cDFT is recast in the same framework as Kohn-Sham DFT: a new functional of the potential that includes the constraints is proposed, where the constraints, the atomic positions, or the lattice parameters are treated all alike, while all other ingredients of the usual potential-based DFT algorithms are unchanged, thanks to the formulation of the adequate residual. Tests of this approach for the case of spin constraints (collinear and noncollinear) and charge constraints are performed. Expressions for the derivatives with respect to constraints (e.g., the spin torque) for the atomic forces and the stress tensor in cDFT are provided. The latter allows one to study striction effects as a function of the angle between spins. We apply this formalism to body-centered cubic iron and first reproduce the well-known magnetization amplitude as a function of the angle between local magnetizations. We also study stress as a function of such an angle. Then, the local collinear magnetization and the local atomic charge are varied together. Since the atomic spin magnetizations, local atomic charges, atomic positions, and lattice parameters are treated on an equal footing, this formalism is an ideal starting point for the generation of model Hamiltonians and machine-learning potentials, computation of second or third derivatives of the energy as delivered from density-functional perturbation theory, or for second-principles approaches.

Sr-Doped Superionic Hydrogen Glass: Synthesis and Properties of SrH22.

Recently, several research groups announced reaching the point of metallization of hydrogen above 400 GPa. Despite notable progress, detecting superconductivity in compressed hydrogen remains an unsolved problem. Following the mainstream of extensive investigations of compressed metal polyhydrides, here small doping of molecular hydrogen by strontium is demonstrated to lead to a dramatic reduction in the metallization pressure to ≈200 GPa. Studying the high-pressure chemistry of the Sr-H system, the formation of several new phases is observed: C2/m-Sr3 H13 , pseudocubic SrH6 , SrH9 with cubic F 4 ¯ 3 m $F\bar{4}3m$ -Sr sublattice, and pseudo tetragonal superionic P1-SrH22 , the metal hydride with the highest hydrogen content (96 at%) discovered so far. High diffusion coefficients of hydrogen in the latter phase DH  = 0.2-2.1 × 10-9 m2 s-1 indicate an amorphous state of the H-sublattice, whereas the strontium sublattice remains solid. Unlike Ca and Y, strontium forms molecular semiconducting polyhydrides, whereas calcium and yttrium polyhydrides are high-TC superconductors with an atomic H sublattice. The discovered SrH22 , a kind of hydrogen sponge, opens a new class of materials with ultrahigh content of hydrogen.

Computational Modeling of 2D Materials under High Pressure and Their Chemical Bonding: Silicene as Possible Field-Effect Transistor.

To study the possibility for silicene to be employed as a field-effect transistor (FET) pressure sensor, we explore the chemistry of monolayer and multilayered silicene focusing on the change in hybridization under pressure. Ab initio computations show that the effect of pressure depends greatly on the thickness of the silicene film, but also reveals the influence of real experimental conditions, where the pressure is not hydrostatic. For this purpose, we introduce anisotropic strain states. With pure uniaxial stress applied to silicene layers, a path for sp3 silicon to sp3d silicon is found, unlike with pure hydrostatic pressure. Even with mixed-mode stress (in-plane pressure half of the out-of-plane one), we find no such path. In addition to introducing our theoretical approach to study 2D materials, we show how the hybridization change of silicene under pressure makes it a good FET pressure sensor.

Exciton-plasmon States in nanoscale materials: breakdown of the Tamm-Dancoff approximation.

Within the Tamm-Dancoff approximation, ab initio approaches describe excitons as packets of electron-hole pairs propagating only forward in time. However, we show that in nanoscale materials excitons and plasmons hybridize, creating exciton-plasmon states where the electron-hole pairs oscillate back and forth in time. Then, as exemplified by the trans-azobenzene molecule and the carbon nanotubes, the Tamm-Dancoff approximation yields errors larger than the accuracy claimed in ab initio calculations. Instead, we propose a general and efficient approach that avoids the Tamm-Dancoff approximation, correctly describes excitons, plasmons, and exciton-plasmon states, and provides a good agreement with experimental results.

ChemEnv: a fast and robust coordination environment identification tool.

Coordination or local environments have been used to describe, analyze and understand crystal structures for more than a century. Here, a new tool called ChemEnv, which can identify coordination environments in a fast and robust manner, is presented. In contrast to previous tools, the assessment of the coordination environments is not biased by small distortions of the crystal structure. Its robust and fast implementation enables the analysis of large databases of structures. The code is available open source within the pymatgen package and the software can also be used through a web app available on http://crystaltoolkit.org through the Materials Project.

Incipient Metals: Functional Materials with a Unique Bonding Mechanism.

While solid-state materials are commonly classified as covalent, ionic, or metallic, there are cases that defy these iconic bonding mechanisms. Phase-change materials (PCMs) for data storage are a prominent example: they have been claimed to show "resonant bonding," but a clear definition of this mechanism has been lacking. Here, it is shown that these solids are fundamentally different from resonant bonding in the π-orbital systems of benzene and graphene, based on first-principles data for vibrational, optical, and polarizability properties. It is shown that PCMs and related materials exhibit a unique mechanism between covalent and metallic bonding. It is suggested that these materials be called "incipient metals," and their bonding nature "metavalent". Data for a diverse set of 58 materials show that metavalent bonding is not just a superposition of covalent and metallic cases, but instead gives rise to a unique and anomalous set of physical properties. This allows the derivation of a characteristic fingerprint of metavalent bonding, composed of five individual components and firmly rooted in physical properties. These findings are expected to accelerate the discovery and design of functional materials with attractive properties and applications, including nonvolatile memories, thermoelectrics, photonics, and quantum materials.

High-throughput density-functional perturbation theory phonons for inorganic materials.

The knowledge of the vibrational properties of a material is of key importance to understand physical phenomena such as thermal conductivity, superconductivity, and ferroelectricity among others. However, detailed experimental phonon spectra are available only for a limited number of materials, which hinders the large-scale analysis of vibrational properties and their derived quantities. In this work, we perform ab initio calculations of the full phonon dispersion and vibrational density of states for 1521 semiconductor compounds in the harmonic approximation based on density functional perturbation theory. The data is collected along with derived dielectric and thermodynamic properties. We present the procedure used to obtain the results, the details of the provided database and a validation based on the comparison with experimental data.

Kinked silicon nanowires-enabled interweaving electrode configuration for lithium-ion batteries.

A tri-dimensional interweaving kinked silicon nanowires (k-SiNWs) assembly, with a Ni current collector co-integrated, is evaluated as electrode configuration for lithium ion batteries. The large-scale fabrication of k-SiNWs is based on a procedure for continuous metal assisted chemical etching of Si, supported by a chemical peeling step that enables the reuse of the Si substrate. The kinks are triggered by a simple, repetitive etch-quench sequence in a HF and H2O2-based etchant. We find that the inter-locking frameworks of k-SiNWs and multi-walled carbon nanotubes exhibit beneficial mechanical properties with a foam-like behavior amplified by the kinks and a suitable porosity for a minimal electrode deformation upon Li insertion. In addition, ionic liquid electrolyte systems associated with the integrated Ni current collector repress the detrimental effects related to the Si-Li alloying reaction, enabling high cycling stability with 80% capacity retention (1695 mAh/gSi) after 100 cycles. Areal capacities of 2.42 mAh/cm2 (1276 mAh/gelectrode) can be achieved at the maximum evaluated thickness (corresponding to 1.3 mgSi/cm2). This work emphasizes the versatility of the metal assisted chemical etching for the synthesis of advanced Si nanostructures for high performance lithium ion battery electrodes.

Ab Initio Approach to Second-order Resonant Raman Scattering Including Exciton-Phonon Interaction.

Raman spectra obtained by the inelastic scattering of light by crystalline solids contain contributions from first-order vibrational processes (e.g. the emission or absorption of one phonon, a quantum of vibration) as well as higher-order processes with at least two phonons being involved. At second order, coupling with the entire phonon spectrum induces a response that may strongly depend on the excitation energy, and reflects complex processes more difficult to interpret. In particular, excitons (i.e. bound electron-hole pairs) may enhance the absorption and emission of light, and couple strongly with phonons in resonance conditions. We design and implement a first-principles methodology to compute second-order Raman scattering, incorporating dielectric responses and phonon eigenstates obtained from density-functional theory and many-body theory. We demonstrate our approach for the case of silicon, relating frequency-dependent relative Raman intensities, that are in excellent agreement with experiment, to different vibrations and regions of the Brillouin zone. We show that exciton-phonon coupling, computed from first principles, indeed strongly affects the spectrum in resonance conditions. The ability to analyze second-order Raman spectra thus provides direct insight into this interaction.

When Density Functional Approximations Meet Iron Oxides.

Three density functional approximations (DFAs), PBE, PBE+U, and Heyd-Scuseria-Ernzerhof screened hybrid functional (HSE), were employed to investigate the geometric, electronic, magnetic, and thermodynamic properties of four iron oxides, namely, α-FeOOH, α-Fe2O3, Fe3O4, and FeO. Comparing our calculated results with available experimental data, we found that HSE (a = 0.15) (containing 15% "screened" Hartree-Fock exchange) can provide reliable values of lattice constants, Fe magnetic moments, band gaps, and formation energies of all four iron oxides, while standard HSE (a = 0.25) seriously overestimates the band gaps and formation energies. For PBE+U, a suitable U value can give quite good results for the electronic properties of each iron oxide, but it is challenging to accurately get other properties of the four iron oxides using the same U value. Subsequently, we calculated the Gibbs free energies of transformation reactions among iron oxides using the HSE (a = 0.15) functional and plotted the equilibrium phase diagrams of the iron oxide system under various conditions, which provide reliable theoretical insight into the phase transformations of iron oxides.

Reproducibility in density functional theory calculations of solids.

The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Identification and design principles of low hole effective mass p-type transparent conducting oxides.

The development of high-performance transparent conducting oxides is critical to many technologies from transparent electronics to solar cells. Whereas n-type transparent conducting oxides are present in many devices, their p-type counterparts are not largely commercialized, as they exhibit much lower carrier mobilities due to the large hole effective masses of most oxides. Here we conduct a high-throughput computational search on thousands of binary and ternary oxides and identify several highly promising compounds displaying exceptionally low hole effective masses (up to an order of magnitude lower than state-of-the-art p-type transparent conducting oxides), as well as wide band gaps. In addition to the discovery of specific compounds, the chemical rationalization of our findings opens new directions, beyond current Cu-based chemistries, for the design and development of future p-type transparent conducting oxides.

First-principles study of the electro-optic effect in ferroelectric oxides.

We first present a method to compute the electro-optic tensor from first principles, explicitly taking into account the electronic, ionic and piezoelectric contributions. We then study the nonlinear optic behavior of three paradigmatic ferroelectric oxides. Our calculations reveal the dominant contribution of the soft mode to the electro-optic coefficients in LiNbO3 and BaTiO3 and its minor role in PbTiO3. We identify the coupling between the electric field and the polar atomic displacements along the B-O chains as the origin of the large electro-optic response in perovskite ABO3 compounds.

Concentration of small ring structures in vitreous silica from a first-principles analysis of the Raman spectrum.

Using a first-principles approach, we calculate Raman spectra for a model structure of vitreous silica. We develop a perturbational method for calculating the dielectric tensor in an ultrasoft pseudopotential scheme and obtain Raman coupling tensors by finite differences with respect to atomic displacements. For frequencies below 1000 cm(-1), the parallel-polarized Raman spectrum of vitreous silica is dominated by oxygen bending motions, showing a strong sensitivity to the intermediate range structure. By modeling the Raman coupling, we derive estimates for the concentrations of three- and four-membered rings from the experimental intensities of the Raman defect lines.