Search

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.

Polarons are quasi-particles made from electrons interacting with vibrations in crystal lattices. They derive their name from the strong electron-vibration polar interaction in ionic systems, that induces associated spectroscopic and optical signatures of such quasi-particles in these materials. In this paper, we focus on diamond, a non-polar crystal with inversion symmetry which nevertheless shows characteristic signatures of polarons, better denoted "nonpolarons" in this case. The polaronic effects are produced by short-range crystal fields with only a small influence of long-range quadrupoles. The 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\"{o}hlich case. The temperature-dependent spectral function of diamond is determined 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., Comment: 10 pages 16 Figures + Subfigures 2 Tables 82 References Supplementary Information (SI) at the end: SI 10 pages SI 14 Figures + Subfigures SI 1 Table

The electronic structure of semiconductors and insulators is affected by ionic motion through electron-phonon interaction, yielding temperature-dependent band gap energies and zero-point renormalization (ZPR) at absolute zero temperature. For polar materials, the most significant contribution to the band gap ZPR can be understood in terms of the Fr\"ohlich model, which focuses on the nonadiabatic interaction between an electron and the macroscopic electrical polarization created by a long-wavelength optical longitudinal phonon mode. On the other hand, spin-orbit interaction (SOC) modifies the bare electronic structure, which will, in turn, affect the electron-phonon interaction and the ZPR. We present a comparative investigation of the effect of SOC on the band gap ZPR of twenty semiconductors and insulators with cubic symmetry using first-principles calculations. We observe a SOC-induced decrease of the ZPR, up to 30%, driven by the valence band edge, which almost entirely originates from the modification of the bare electronic eigenenergies and the decrease of the hole effective masses near the $\Gamma$ point. We also incorporate SOC into a generalized Fr\"ohlich model, addressing the Dresselhaus splitting which occurs in noncentrosymmetric materials, and confirm that the predominance of nonadiabatic effects on the band gap ZPR of polar materials is unchanged when including SOC. Our generalized Fr\"ohlich model with SOC provides a reliable estimate of the SOC-induced decrease of the polaron formation energy obtained from first principles and brings to light some fundamental subtleties in the numerical evaluation of the effective masses with SOC for noncentrosymmetric materials. We finally warn about a possible breakdown of the parabolic approximation within the physically relevant energy range of the Fr\"ohlich interaction for materials with high phonon frequencies treated with SOC., Comment: 34 pages including Supplemental Material; 10 figures, 1 table, 3 Supplementary figures, 9 Supplementary tables

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.

Starting from recent advances in the first-principles modeling of polarons, variational polaron equations in the strong-coupling adiabatic approximation are formulated in Bloch space. In this framework, polaron formation energy as well as individual electron, phonon and electron-phonon contributions are obtained. We suggest an efficient gradient-based optimization algorithm and apply these equations to the generalized Fr\"ohlich model with anisotropic non-degenerate electronic bands, both in two- and three-dimensional cases. The effect of the divergence of Fr\"ohlich electron-phonon matrix elements at $\Gamma$-point is treated analytically, improving the convergence with respect to the sampling in reciprocal space. The whole methodology is validated by obtaining the known asymptotic solution of the standard Fr\"ohlich model in isotropic scenario and also by comparing our results with the Gaussian ansatz approach, showing the difference between the numerically exact and Gaussian trial wavefunctions. Additionally, decomposition of the energy into individual terms allows one to recover the Pekar's 1:2:3:4 theorem, which is shown to be valid even in the anisotropic case. We expect that the improvements in the formalism and numerical implementation will be applicable beyond the large polaron hypothesis inherent to Fr\"ohlich model., Comment: 14 pages, 5 figures

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-Sr

The zero-point lattice expansion (ZPLE) is a small variation of the lattice parameters induced by the presence of phonons in a material compared to the static lattice picture. It contributes significantly to the zero-point renormalization (ZPR) of the band gap energy, but its consequences have not been investigated as thoroughly as those stemming from electron-phonon interactions. In the usual first-principles approach, one evaluates the ZPLE by minimizing the $T=0$ K Helmholtz free energy. In this work, we show that the formalism based on the Gr\"uneisen parameters, which commonly neglects zero-point effects, can be efficiently used to compute ZPLE for both isotropic and anisotropic materials at much lower computational cost. We systematically test this formalism on 22 cubic and wurtzite materials and obtain excellent agreement with free energy minimization results for both the ZPLE and the resulting band gap ZPR. We use our results to validate an empirical expression estimating the ZPLE-induced ZPR and unveil its sensitivity to the temperature range involved in estimating the ZPLE from experimental data. Our findings finally reveal that the ZPLE contribution to the band gap ZPR can reach 20% to more than 80% of the electron-phonon interaction contribution for heavier or more ionic materials, including materials containing light atoms. Considering both contributions on an equal footing is thus essential should one attempt to compare theoretical ZPR results with experimental data., Comment: 30 pages, 5 figures, 14 tables. Sec. II D expanded for internal consistency, typo corrected in Fig.4., clarified Tables captions, minor corrections to the text; results and conclusions unchanged

In the simplest picture, a "self-trapped" polaron forms when an excess electron or hole deforms a crystal lattice, creating a potential well with bound states. Properties of self-trapped polarons in methylammonium lead iodide perovskite (MAPbI3), which is widely used as solar cell absorber, are of great interest, and are a subject of ongoing investigations and debates concerning the existence of large polarons with the co-presence of metastable self-trapping. Herein, we employ a self-interaction-free density functional theory method to investigate the stability of small polarons in tetragonal MAPbI3 phase. The electron small polaron is found to be unstable, while the hole small polaron is found to be metastable at realistic operation temperatures of solar cells. Further, the hole polaron is found to have a hole band close to the conduction band, which in conjunction with its metastability suggests that small polarons will have an appreciable effect on charge-carrier recombinations in MAPbI3. Further, we posit that the existence of the metastable polarons in addition to the large polarons may explain the experimentally observed non-monotonic temperature dependence of bimolecular charge-carrier recombination rate in tetragonal MAPbI3 phase.

Polarons, that is, charge carriers correlated with lattice deformations, are ubiquitous quasiparticles in semiconductors, and play an important role in electrical conductivity. To date most theoretical studies of so-called large polarons, in which the lattice can be considered as a continuum, have focused on the original Fr\"ohlich model: a simple (non-degenerate) parabolic isotropic electronic band coupled to one dispersionless longitudinal optical phonon branch. The Fr\"ohlich model allows one to understand characteristics such as polaron formation energy, radius, effective mass and mobility. Real cubic materials, instead, have electronic band extrema that are often degenerate or anisotropic and present several phonon modes. In the present work, we address such issues. We keep the continuum hypothesis inherent to the large polaron Fr\"ohlich model, but waive the isotropic and non-degeneracy hypotheses, and also include multiple phonon branches. For polaron effective masses, working at the lowest order of perturbation theory, we provide analytical results for the case of anisotropic electronic energy dispersion, with two distinct effective masses (uniaxial) and numerical simulations for the degenerate 3-band case, typical of III-V and II-VI semiconductor valence bands. We also deal with the strong-coupling limit, using a variational treatment: we propose trial wavefunctions for the above-mentioned cases, providing polaron radii and energies. Then, we evaluate the polaron formation energies, effective masses and localisation lengths using parameters representative of a dozen II-VI, III-V and oxide semiconductors, for both electron and hole polarons...In the non-degenerate case, we compare the perturbative approach with the Feynman path integral approach in characterisizing polarons in the weak coupling limit..., Comment: 18 pages, 10 figures in main text; 11 pages, 1 figure in supplementary material

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.

Electronic and optical properties of materials are affected by atomic motion through the electron-phonon interaction: not only band gaps change with temperature, but even at absolute zero temperature, zero-point motion causes band-gap renormalization. We present a large-scale first-principles evaluation of the zero-point renormalization of band edges beyond the adiabatic approximation. For materials with light elements, the band gap renormalization is often larger than 0.3 eV, and up to 0.7 eV. This effect cannot be ignored if accurate band gaps are sought. For infrared-active materials, global agreement with available experimental data is obtained only when non-adiabatic effects are taken into account. They even dominate zero-point renormalization for many materials, as shown by a generalized Fr\"ohlich model that includes multiple phonon branches, anisotropic and degenerate electronic extrema, whose range of validity is established by comparison with first-principles results., Comment: Note that the supplementary material, included in the main file, is structured differently than in the published version, with a different numbering scheme

A new tool called ChemEnv, which can identify coordination environments in a fast and robust manner, is presented., 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.

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., G. B., G.-M. R., and G. H. acknowledge financial support from F.R.S.-FNRS. H. P. C. M. acknowledges financial support from F.R.S.-FNRS through the PDR Grants HTBaSE (T.1071.15). M. G., M. J. V., and X. G. acknowledge financial support from F.R.S.-FNRS through the PDR Grants AIXPHO (T.0238.13) and ALPS (T.0103.19). M. J. V. thanks FNRS and ULiege for a sabbatical grant in ICN2. M. S. and M. R. acknowledge the support of Ministerio de Economia, Industria y Competitividad (MINECO-Spain) through Grants No. MAT2016-77100-C2-2-P and No. SEV-2015-0496; of Generalitat de Catalunya (Grant No. 2017 SGR1506); and of the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (Grant Agreement No. 724529). The present research benefited from computational resources made available on the Tier-1 supercomputer of the Federation WallonieBruxelles, infrastructure funded by the Walloon Region under Grant Agreement No. 1117545.

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.

We describe a new approach to compute the electron-phonon self-energy and carrier mobilities in semiconductors. Our implementation does not require a localized basis set to interpolate the electron-phonon matrix elements, with the advantage that computations can be easily automated. Scattering potentials are interpolated on dense q meshes using Fourier transforms and ab initio models to describe the long-range potentials generated by dipoles and quadrupoles. To reduce significantly the computational cost, we take advantage of crystal symmetries and employ the linear tetrahedron method and double-grid integration schemes, in conjunction with filtering techniques in the Brillouin zone. We report results for the electron mobility in Si, GaAs, and GaP obtained with this new methodology., G.B., G.-M.R., and G.H. acknowledge financial support from F.R.S.-FNRS. H.P.C.M. acknowledges financial support from F.R.S.-FNRS through the PDR Grants HTBaSE (T.1071.15). M.G., M.J.V., and X.G. acknowledge financial support from F.R.S.-FNRS through the PDR Grants AIXPHO (T.0238.13) and ALPS (T.0103.19). M.J.V. thanks FNRS and ULiege for a sabbatical grant in ICN2. M.S. and M.R. acknowledge the support of Ministerio de Economia, Industria y Competitividad (MINECO-Spain) through Grants No. MAT2016-77100-C2-2-P and No. SEV-2015-0496; of Generalitat de Catalunya (Grant No. 2017 SGR1506); and of the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (Grant Agreement No. 724529).

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.

Ga ( 1 − x ) In x N alloys, widely employed to produce light-emitting diodes, exhibit a bowing of the band gap as a function of concentration x , and a luminescence spectral width which differs from the expected value of 1.8 kT. Through first-principles calculations, based on many-body perturbation theory and density-functional theory with a meta-GGA exchange–correlation functional, we explore jointly these effects, in an exhaustive set of Ga ( 1 − x ) In x N supercells with 16 atoms. We disentangle the bowing due to the average volume change with the one due local atomic configuration and local relaxation. The first one account for about 40% of the bowing, despite that fact that the change of volume with respect to concentration is nearly linear (Vegard’s law). The computed bowing parameter is 1.39 eV. The experimental broadening between 3 kT and 8 kT, not examined theoretically until now, is well accounted by local atomic configuration changes and lifting of the degeneracy of the top of the valence band.

We study from first principles the stability of neutral and charged oxygen vacancies in lutetium oxyorthosilicate, Lu2SiO5, as well as its possible modification due to the presence of Ce 3 + , as present in commercial scintillators. We show that the neutral oxygen vacancy with the lowest formation energy forms at the oxygen sites within the [SiO4] tetrahedra instead of the interstitial oxygen site bonded exclusively to lutetium atoms. The discrepancy with a previous study is attributed to the quality of the pseudopotential. Support for these results is found by performing a bonding analysis of the oxygen sites, as well as oxygen vacancy calculations in the iso-structural Y2SiO5 compound. In addition, we find that the incorporation of Ce 3 + ion does not affect the stability of oxygen vacancies in Lu2SiO5.

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.

The diffusion of large databases collecting different kind of material properties from high-throughput density functional theory calculations has opened new paths in the study of materials science thanks to data mining and machine learning techniques. Phonon calculations have already been employed successfully to predict materials properties and interpret experimental data, e.g. phase stability, ferroelectricity and Raman spectra, so their availability for a large set of materials will further increase the analytical and predictive power at hand. Moving to a larger scale with density functional perturbation calculations, however, requires the presence of a robust framework to handle this challenging task. In light of this, we automatized the phonon calculation and applied the result to the analysis of the convergence trends for several materials. This allowed to identify and tackle some common problems emerging in this kind of simulations and to lay out the basis to obtain reliable phonon band structures from high-throughput calculations, as well as optimizing the approach to standard phonon simulations.

The recently discovered Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ red phosphor, candidate for the next generation of eco-efficient white light-emitting diodes, exhibits excellent emission spectral position and exceptionally small linewidth. It belongs to the UCr$_4$C$_4$-structure family of phosphors, containing many potential candidates commercial phosphors, whose small linewidth, tentatively ascribed to the high-symmetry cuboid environment of the doping site, has drawn the attention of researchers in the last five years. We use density functional theory, $\Delta$SCF method and configuration coordinate models (CCM) to provide a complete characterization of this material. Using a multi-dimensional CCM, an accurate description of the coupling of the vibronic structure with the electronic 5d$\rightarrow$4f transition is obtained, including the partial Huang-Rhys factors and frequency of the dominant modes. We show that, in addition to the first-coordination shell cuboid deformation mode, low-frequency phonon modes involving chains of strontium atoms along the tetragonal axis shape the emission linewidth in Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$. This finding sheds new light on the emission properties of UCr$_4$C$_4$-structure phosphors, possessing similar Ca/Sr/Ba channel. Our approach provides a robust theoretical framework to systematically study the emission spectra of such Eu-doped phosphors, and predict candidates with expected similar or even sharper linewidth.

Coordination or local environments (e.g., tetrahedra and octahedra) are powerful descriptors of the crystalline structure of materials. These structural descriptors are essential to the understanding of crystal chemistry and the design of new materials. However, extensive statistics on the occurence of local environment are not available even on common chemistries such as oxides. Here, we present the first large-scale statistical analysis of the coordination environments of cations in oxides using a large set of experimentally observed compounds (about 8,000). Using a newly developed method, we provide the distribution of local environment for each cation in oxides. We discuss our results highlighting previously know trends, unexpected coordination environments as well as compounds presenting very rare coordinations. Our work complements the know-how of the solid state chemist with a statistically sound analysis and pave the way for further data mining efforts linking for instance coordination environment...

We describe recent advances in condensed matter physics that make it possible to study non-linear optical properties of periodic solids from first-principles density functional theory. These methods are either based on perturbation theory in the framework of the 2n + 1 theorem, finite electric fields or effective hamiltonians. We apply them to the study of non-linear optical properties of various ferroelectric ABO3 compounds such as BaTiO3, PbTiO3 or LiNbO3.

The one-dimensional configuration coordinate model (1D-CCM) is widely used for the analysis of photoluminescence in molecules and doped solids, and relies on a linear combination of the equilibrium nuclear configurations of ground and excited states. It delivers an estimation of the energy barrier at which ground and excited state curves cross, semi-classically linked to non-radiative transition rate and thermal quenching. To assess its predictive power for the latter properties, we propose a new \textit{optimized configuration path (OCP) method in which} the ground-state and excited-state forces are mixed instead of their configurations. We also define another one-parameter model thanks a double energy parabola hypothesis (DEPH). We compare the OCP method and the DEPH reference with the 1D-CCM for three paradigmatic 4f-5d phosphors Y$_3$Al$_5$O$_{12}$:Ce, Lu$_2$SiO$_5$:Ce, and YAlO$_3$:Ce. We find that the OCP and DEPH methods yield similar results with geometries that have significantly lower ground-state energies than the 1D-CCM for the same 4f-5d energy difference. However, the OCP method suffers from the appearance of multiple local minima, rendering the clear determination of the optimal geometry very difficult in practice. Still the OCP method allows one to quantify the deviations from the 1D-CCM, therefore increasing confidence in the lower bound obtained from the DEPH for the 4f-5d crossing barrier, and its comparison with the energy of the auto-ionization thermal quenching mechanism. We expect the OCP approach to be applicable to other luminescent materials or molecules., Comment: 12 pages, 6 figure (Accepted by Phys. Rev. B)

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 equilibri...

Understanding the physical mechanisms behind thermal effects in phosphors is crucial for white light-emitting device (WLEDs) applications, as thermal quenching of their photoluminescence might render them useless. We analyze from first-principles, before and after absorption/emission of light, two chemically close Eu-doped Ba3Si6O12N2 and Ba3Si6O9N4 crystals for WLEDs. The first one has an almost constant emission intensity with increasing temperature whereas the other one does not. Our results, in which the Eu-5d levels are obtained inside the band gap thanks to the removal of an electron from the 4f7 shell, and the atomic neighborhood properly relaxed in the excited state, attributes the above-mentioned experimental difference to an autoionization model of the thermal quenching, based on the energy difference between Eu5d and the conduction band minimum. Our depleted-shifted 4f method can identify luminescent centers and therefore allows for effective crystal site engineering of luminescent centers in p...

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.

International audience; In a series of publications, Hardy Gross and co-workers have highlighted the interest of an “exact factorization” approach to the interacting electron-nuclei problem, be it time-independent or time-dependent. In this approach, an effective potential governs the dynamics of the nuclei such that the resulting N-body nuclear density is in principle exact. This contrasts with the more usual adiabatic approach, where the effective potential leads to an approximate nuclear density. Inspired by discussions with Hardy, we explore the factorization idea for arbitrary many-body Hamiltonians, generalizing the electron-nuclei case, with a focus on the static case. While the exact equations do not lead to any practical advantage, they are illuminating, and may therefore constitute a suitable starting point for approximations. In particular, we find that unitary transformations that diagonalize the coupling term for one of the sub-systems make exact factorization appealing. The algorithms by which the equations for the separate subsystems can be solved in the time-independent case are also explored. We illustrate our discussions using the two-site Holstein model and the quantum Rabi model. Two factorization schemes are possible: one where the boson field feels a potential determined by the electrons, and the reverse exact factorization, where the electrons feel a potential determined by the bosons; both are explored in this work. A comparison with a self-energy approach is also presented.