Source author record

Miguel A. L. Marques

Miguel A. L. Marques appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals

The density functional theory (DFT) approximations that are the most accurate for the calculation of band gap of bulk materials are hybrid functionals like HSE06, the MBJ potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGA), like HLE16, or meta-GGAs, like (m)TASK, have proven to be also quite accurate for the band gap. Here, the focus is on 2D materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [Rauch et al., Phys. Rev. B 101, 245163 (2020) and Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the band gap of 2D systems, $G_{0}W_{0}$ results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the band gaps that are the closest to $G_{0}W_{0}$. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.

preprint2021arXiv

Machine-learning correction to density-functional crystal structure optimization

Density functional theory is routinely applied to predict crystal structures. The most common exchange-correlation functionals used to this end are the Perdew-Burke-Ernzerhof (PBE) approximation and its variant PBEsol. We investigate the performance of these functionals for the prediction of lattice parameters and show how to enhance their accuracy using machine learning. Our dataset is constituted by experimental crystal structures of the Inorganic Crystal Structure Database matched with PBE-optmized structures stored in the materials project database. We complement these data with PBEsol calculations. We demonstrate that the accuracy and precision of PBE/PBEsol volume predictions can be noticeably improved a posteriori by employing simple, explainable machine learning models. These models can improve PBE unit cell volumes to match the accuracy of PBEsol calculations, and reduce the error of the latter with respect to experiment by 35%. Further, the error of PBE lattice constants is reduced by a factor of 3--5. A further benefit of our approach is the implicit correction of finite temperature effects without performing phonon calculations.

preprint2021arXiv

Meta-local density functionals: a new rung on Jacob's ladder

The homogeneous electron gas (HEG) is a key ingredient in the construction of most exchange-correlation functionals of density-functional theory. Often, the energy of the HEG is parameterized as a function of its spin density $n$, leading to the local density approximation (LDA) for inhomogeneous systems. However, the connection between the electron density and kinetic energy density of the HEG can be used to generalize the LDA by evaluating it on a weighted geometric average of the local spin density and the spin density of a HEG that has the local kinetic energy density of the inhomogeneous system, with a mixing ratio $x$. This leads to a new family of functionals that we term meta-local density approximations (meta-LDAs), which are still exact for the HEG, which are derived only from properties of the HEG, and which form a new rung of Jacob's ladder of density functionals. The first functional of this ladder, the local $τ$ approximation (LTA) of Ernzerhof and Scuseria that corresponds to $x=1$ is unfortunately not stable enough to be used in self-consistent field calculations, because it leads to divergent potentials as we show in this work. However, a geometric averaging of the LDA and LTA densities with smaller values of $x$ not only leads to numerical stability of the resulting functional, but also yields more accurate exchange energies in atomic calculations than the LDA, the LTA, or the tLDA functional ($x=1/4$) of Eich and Hellgren. We choose $x=0.50$ as it gives the best total energy in self-consistent exchange-only calculations for the argon atom. Atomization energy benchmarks confirm that the choice $x=0.50$ also yields improved energetics in combination with correlation functionals in molecules, almost eliminating the well-known overbinding of the LDA and reducing its error by two thirds.

preprint2020arXiv

Accurate electronic band gaps of two-dimensional materials from the local modified Becke-Johnson potential

The electronic band structures of two-dimensional materials are significantly different from those of their bulk counterparts, due to quantum confinement and strong modifications of electronic screening. An accurate determination of electronic states is a prerequisite to design electronic or optoelectronic applications of two-dimensional materials, however, most of the theoretical methods we have available to compute band gaps are either inaccurate, computationally expensive, or only applicable to bulk systems. Here we show that reliable band structures of nanostructured systems can now be efficiently calculated using density-functional theory with the local modified Becke-Johnson exchange-correlation functional that we recently proposed. After re-optimizing the parameters of this functional specifically for two-dimensional materials, we show, for a test set of almost 300 systems, that the obtained band gaps are of comparable quality as those obtained using the best hybrid functionals, but at a very reduced computational cost. These results open the way for accurate high-throughput studies of band-structures of two-dimensional materials and for the study of van der Waals heterostructures with large unit cells.

preprint2020arXiv

Homogeneous electron gas in arbitrary dimensions

The homogeneous electron gas is one of the most studied model systems in condensed matter physics. It is also at the basis of the large majority of approximations to the functionals of density functional theory. As such, its exchange-correlation energy has been extensively studied, and is well-known for systems of 1, 2, and 3 dimensions. Here, we extend this model and compute the exchange and correlation energy, as a function of the Wigner-Seitz radius $r_s$, for arbitrary dimension $D$. We find a very different behavior for reduced dimensional spaces ($D=1$ and 2), our three dimensional space, and for higher dimensions. In fact, for $D > 3$, the leading term of the correlation energy does not depend on the logarithm of $r_s$ (as for $D=3$), but instead scales polynomialy: $ -c_D /r_s^{γ_D}$, with the exponent $γ_D=(D-3)/(D-1)$. In the large-$D$ limit, the value of $c_D$ is found to depend linearly with the dimension. In this limit, we also find that the concepts of exchange and correlation merge, sharing a common $1/r_s$ dependence.

preprint2020arXiv

Iodine molecule modifications with high pressure

Metallization and dissociation are key transformations in diatomic molecules at high densities particularly significant for modeling giant planets. Using X-ray absorption spectroscopy and atomistic modeling, we demonstrate that in halogens, the formation of a \textit{connected} molecular structure takes place at pressures well below metallization. Here we show that the iodine diatomic molecule first elongates of $\sim$0.007 Å~up to a critical pressure of $P_c$ $\backsim$7~GPa developing bonds between molecules. Then its length continuously decreases with pressure up to 15-20~GPa. Universal trends in halogens are shown and allow to predict for chlorine a pressure of 42$\pm$8~GPa for molecular bond-length reversal. Our findings tackle the molecule invariability paradigm in diatomic molecular phases at high pressures and may be generalized to other abundant diatomic molecules in the universe, including hydrogen.

preprint2020arXiv

Local modified Becke-Johnson exchange-correlation potential for interfaces, surfaces, and two-dimensional materials

The modified Becke-Johnson meta-GGA potential of density functional theory has been shown to be the best exchange-correlation potential to determine band gaps of crystalline solids. However, it cannot be consistently used for the electronic structure of non-periodic or nanostructured systems. We propose an extension of this potential that enables its use to study heterogeneous, finite and low-dimensional systems. This is achieved by using a coordinate-dependent expression for the parameter $c$ that weights the Becke-Russel exchange, in contrast to the original global formulation, where $c$ is just a fitted number. Our potential takes advantage of the excellent description of band gaps provided by the modified Becke-Johnson potential and preserves its modest computational effort. Furthermore, it yields with one single calculation band diagrams and band offsets of heterostructures and surfaces. We exemplify the usefulness and efficiency of our local meta-GGA potential by testing it for a series of interfaces (Si/SiO$_2$, AlAs/GaAs, AlP/GaP, and GaP/Si), a Si surface, and boron nitride monolayer.

preprint2020arXiv

Reduced Density Matrix Functional Theory for Bosons

Based on a generalization of Hohenberg-Kohn's theorem, we propose a ground state theory for bosonic quantum systems. Since it involves the one-particle reduced density matrix $γ$ as a natural variable but still recovers quantum correlations in an exact way it is particularly well-suited for the accurate description of Bose-Einstein condensates. As a proof of principle we study the building block of optical lattices. The solution of the underlying $v$-representability problem is found and its peculiar form identifies the constrained search formalism as the ideal starting point for constructing accurate functional approximations: The exact functionals for this $N$-boson Hubbard dimer and general Bogoliubov-approximated systems are determined. The respective gradient forces are found to diverge in the regime of Bose-Einstein condensation, $\nabla_γ \mathcal{F} \propto 1/\sqrt{1-N_{\mathrm{BEC}}/N}$, providing a natural explanation for the absence of complete BEC in nature.

preprint2020arXiv

Stable ordered phases of cuprous iodide with complexes of copper vacancies

We perform an exhaustive theoretical study of the phase diagram of Cu-I binaries, focusing on Cu-poor compositions, relevant for p-type transparent conduction. We find that the interaction between neighboring Cu vacancies is the determining factor that stabilizes non-stoichiometric zincblende phases. This interaction leads to defect complexes where Cu vacancies align preferentially along the [100] crystallographic direction. It turns out that these defect complexes have an important influence on hole conductivity, as they lead to dispersive conducting $p$-states that extend up to around 0.8 eV above the Fermi level. We furthermore observe a characteristic peak in the density of electronic states, which could provide an experimental signature for this type of defect complexes.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2019arXiv

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials).

preprint2016arXiv

Emergence of superconductivity in doped H$_2$O ice at high pressure

We investigate the possibility of achieving high-temperature superconductivity in hydrides under pressure by inducing metallization of otherwise insulating phases through doping, a path previously used to render standard semiconductors superconducting at ambient pressure. Following this idea, we study H$_2$O, one of the most abundant and well-studied substances, we identify nitrogen as the most likely and promising substitution/dopant. We show that for realistic levels of doping of a few percent, the phase X of ice becomes superconducting with a critical temperature of about 60 K at 150GPa. In view of the vast number of hydrides that are strongly covalent bonded, but that remain insulating until rather large pressures, our results open a series of new possibilities in the quest for novel high-temperature superconductors.

preprint2015arXiv

Generalized Pauli constraints in reduced density matrix functional theory

Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble $N$-representability conditions. Recently, the topic of pure-state $N$-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble $N$-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.

preprint2015arXiv

Identification of novel Cu, Ag, and Au ternary oxides from global structural prediction

We use ab initio global structural prediction, and specifically the minima hopping method, to explore the periodic table in search of novel oxide phases. In total, we study 183 different compositions of the form MXO2, where M=(Cu, Ag, Au) and X is an element of the periodic table. This set includes the well-known Cu delafossite compounds that are, up to now, the best p-type transparent conductive oxides known to mankind. Our calculations discover 81 stable compositions, out of which only 36 are included in available databases. Some of these new phases are potentially good candidates for transparent electrodes. These results demonstrate, on one hand, how incomplete is still our knowledge of the phase-space of stable ternary materials. On the other hand, we show that structural prediction combined with high-throughput approaches is a powerful tool to extend that knowledge, paving the way for the experimental discovery of new materials on a large scale.

preprint2015arXiv

Low-density silicon allotropes for photovoltaic applications

Silicon materials play a key role in many technologically relevant fields, ranging from the electronic to the photovoltaic industry. A systematic search for silicon allotropes was performed by employing a modified ab initio minima hopping crystal structure prediction method. The algorithm was optimized to specifically investigate the hitherto barely explored low-density regime of the silicon phase diagram by imitating the guest-host concept of clathrate compounds. In total 44 metastable phases are presented, of which 11 exhibit direct or quasi-direct band-gaps in the range of $\approx$1.0-1.8 eV, close to the optimal Shockley-Queisser limit of $\approx$1.4 eV, with a stronger overlap of the absorption spectra with the solar spectrum compared to conventional diamond silicon. Due to the structural resemblance to known clathrate compounds it is expected that the predicted phases can be synthesized.

preprint2015arXiv

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

preprint2014arXiv

First-principles predicted low-energy structures of NaSc(BH4)4

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH$_4$)$_4$ crystallizes in the crystallographic space group $Cmcm$ where each sodium (scandium) atom is surrounded by six scandium (sodium) atoms. A careful investigation of this phase based on \textit{ab initio} calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with $C222_1$ symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minima-hopping method we discover a class of new low-energy structures exhibiting a novel structural motif in which each sodium (scandium) atom is surrounded by four scandium (sodium) atoms arranged at the corners of either a rectangle with nearly equal sides or a tetrahedron. These new phases are all predicted to be insulators with band gaps of $7.9-8.2$ eV. Finally, we estimate the influence of these structures on the hydrogen-storage performance of NaSc(BH$_4$)$_4$.

preprint2013arXiv

Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations

We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $α$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $α$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant).

preprint2013arXiv

Conducting boron sheets formed by the reconstruction of the α-boron (111) surface

Systematic ab initio structure prediction was applied for the first time to predict low energy surface reconstructions by employing the minima hopping method on the α-boron (111) surface. Novel reconstruction geometries were identified and carefully characterized in terms of structural and electronic properties. Our calculations predict the formation of a planar, mono-layer sheet at the surface, which is responsible for conductive surface states. Furthermore, the isolated boron sheet is shown to be the ground state 2D-structure in vacuum at a hole density of η=1/5 and is therefore a potential candidate as a precursor for boron nano-structures.

preprint2013arXiv

Low-Energy Polymeric Phases of Alanates

Low-energy structures of alanates are currently known to be described by patterns of isolated, nearly ideal tetrahedral [AlH$_4$] anions and metal cations. We discover that the novel polymeric motif recently proposed for LiAlH$_4$ plays a dominant role in a series of alanates, including LiAlH$_4$, NaAlH$_4$, KAlH$_4$, Mg(AlH$_4$)$_2$, Ca(AlH$_4$)$_2$ and Sr(AlH$_4$)$_2$. In particular, most of the low-energy structures discovered for the whole series are characterized by networks of corner-sharing [AlH$_6$] octahedra, forming wires and/or planes throughout the materials. Finally, for Mg(AlH$_4$)$_2$ and Sr(AlH$_4$)$_2$, we identify two polymeric phases to be lowest in energy at low temperatures.

preprint2013arXiv

Strong renormalization of the electronic band gap due to lattice polarization in the GW formalism

The self-consistent GW band gaps are known to be significantly overestimated. We show that this overestimation is, to a large extent, due to the neglect of the contribution of the lattice polarization to the screening of the electron-electron interaction. To solve this problem, we derive within the GW formalism a generalized plasmon-pole model that accounts for lattice polarization. The resulting GW self-energy is used to calculate the band structures of a set of binary semiconductors and insulators. The lattice contribution always decreases the band gap. The shrinkage increases with the size of the longitudinal-transverse optical splitting and it can represent more than 15% of the band gap in highly polar compounds, reducing the band-gap percentage error by a factor of three.

preprint2012arXiv

Libxc: a library of exchange and correlation functionals for density functional theory

The central quantity of density functional theory is the so-called exchange-correlation functional. This quantity encompasses all non-trivial many-body effects of the ground-state and has to be approximated in any practical application of the theory. For the past 50 years, hundreds of such approximations have appeared, with many successfully persisting in the electronic structure community and literature. Here, we present a library that contains routines to evaluate many of these functionals (around 180) and their derivatives.

preprint2012arXiv

Low-energy silicon allotropes with strong absorption in the visible for photovoltaic applications

We present state-of-the-art first-principle calculations of the electronic and optical properties of silicon allotropes with interesting characteristics for applications in thin-film solar cells. These new phases consist of distorted sp$^3$ silicon networks and have a lower formation energy than other experimentally produced silicon phases. Some of these structures turned out to have quasi-direct and dipole-allowed band gaps in the range 0.8--1.5\,eV, and to display absorption coefficients comparable with those of chalcopyrites used in thin-film record solar cells.

preprint2012arXiv

Prediction of a novel monoclinic carbon allotrope

A novel allotrope of carbon with $P2/m$ symmetry was identified during an \emph{ab-initio} minima-hopping structural search which we call $M10$-carbon. This structure is predicted to be more stable than graphite at pressures above 14.4 GPa and consists purely of $sp^3$ bonds. It has a high bulk modulus and is almost as hard as diamond. A comparison of the simulated X-ray diffraction pattern shows a good agreement with experimental results from cold compressed graphite.

preprint2012arXiv

Speeding up the solution of the Bethe-Salpeter equation by a double-grid method and Wannier interpolation

The Bethe-Salpeter equation is a widely used approach to describe optical excitations in bulk semiconductors. It leads to spectra that are in very good agreement with experiment, but the price to pay for such accuracy is a very high computational burden. One of the main bottlenecks is the large number of k-points required to obtain converged spectra. In order to circumvent this problem we propose a strategy to solve the Bethe-Salpeter equation based on a double-grid technique coupled to a Wannier interpolation of the Kohn-Sham band structure. This strategy is then benchmarked for a particularly difficult case, the calculation of the absorption spectrum of GaAs, and for the well studied case of Si. The considerable gains observed in these cases fully validate our approach, and open the way for the application of the Bethe-Salpeter equation to large and complex systems.

preprint2011arXiv

Optimal Control of the Electronic Current Density - An application to one- and two-dimensional one-electron systems

Quantum optimal control theory is a powerful tool for engineering quantum systems subject to external fields such as the ones created by intense lasers. The formulation relies on a suitable definition for a target functional, that translates the intended physical objective to a mathematical form. We propose the use of target functionals defined in terms of the one-particle density and its current. A strong motivation for this is the possibility of using time-dependent density-functional theory for the description of the system dynamics. We exemplify this idea by defining an objective functional that on one hand attempts a large overlap with a target density and on the other hand minimizes the current. The latter requirement leads to optimized states with increased stability, which we prove with a few examples of one- and two-dimensional one-electron systems.

preprint2011arXiv

The crystal structure of cold compressed graphite

Through a systematic structural search we found an allotrope of carbon with Cmmm symmetry which we predict to be more stable than graphite for pressures above 10 GPa. This material, which we refer to as Z-carbon, is formed by pure sp3 bonds and is the only carbon allotrope which provides an excellent match to unexplained features in experimental X-ray diffraction and Raman spectra of graphite under pressure. The transition from graphite to Z-carbon can occur through simple sliding and buckling of graphene sheets. Our calculations predict that Z-carbon is a transparent wide band gap semiconductor with a hardness comparable to diamond.

preprint2010arXiv

Bioheterojunction Effect on Fluorescence Origin and Efficiency Improvement of Firefly Chromophores

We propose the heterojunction effect in the analysis of the fluorescence mechanism of the firefly chromophore. Following this analysis, and with respect to the HOMO-LUMO gap alignment between the chromophore's functional fragments, three main heterojunction types (I, II, and I*) are identified. Time-dependent density-functional theory optical absorption calculations for the firefly chromophore show that the strongest excitation appears in the deprotonated anion state of the keto form. This can be explained by its high HOMO-LUMO overlap due to strong bio-heterojunction confinement. It is also found that the nitrogen atom in the thiazolyl rings, due to its larger electronegativity, plays a key role in the emission process, its importance growing when HOMO and LUMO overlap at its location. This principle is applied to enhance the chromophore's fluorescence efficiency and to guide the functionalization of molecular optoelectronic devices.

preprint2010arXiv

Density-based mixing parameter for hybrid functionals

A very popular ab-initio scheme to calculate electronic properties in solids is the use of hybrid functionals in density functional theory (DFT) that mixes a portion of Fock exchange with DFT functionals. In spite of their success, a major problem still remains, related to the use of one single mixing parameter for all materials. Guided by physical arguments that connect the mixing parameter to the dielectric properties of the solid, and ultimately to its band gap, we propose a method to calculate this parameter from the electronic density alone. This method is able to cut significantly the error of traditional hybrid functionals for large and small gap materials, while retaining a good description of structural properties. Moreover, its implementation is simple and leads to a negligible increase of the computational time.

preprint2010arXiv

Effects of polarization on the band-structure of delafossite transparent conductive oxides

We use hybrid functionals and restricted self-consistent GW, state-of-the-art theoretical approaches for quasiparticle band structures, to study the electronic states of delafossite Cu(Al,In)O$_2$, the first p-type and bipolar transparent conductive oxides. We show that self-consistent GW gives remarkably wider band gaps than all the other approaches used so far. Accounting for polaronic effects in the GW scheme we recover a very nice agreement with experiments. Furthermore, the modifications with respect to the Kohn-Sham bands are strongly k-dependent, which makes questionable the common practice of using a scissor operator. Finally, our results support the view that the low energy structures found in optical experiments, and initially attributed to an indirect transition, are due to intrinsic defects in the samples.

preprint2010arXiv

Enhancing the superconducting transition temperature of BaSi2 by structural tuning

We present a joint experimental and theoretical study of the superconducting phase of the layered binary silicide BaSi2. Compared with the layered AlB2 structure of graphite or diboride-like superconductors, in the hexagonal structure of binary silicides the sp3 arrangement of silicon atoms leads to corrugated sheets. Through a high-pressure synthesis procedure we are able to modify the buckling of these sheets, obtaining the enhancement of the superconducting transition temperature from 4 K to 8.7 K when the silicon planes flatten out. By performing ab-initio calculations based on density functional theory we explain how the electronic and phononic properties of the system are strongly affected by changes in the buckling. This mechanism is likely present in other intercalated layered superconductors, opening the way to the tuning of superconductivity through the control of internal structural parameters.

preprint2010arXiv

Toward an all-round semi-local potential for the electronic exchange

We test local and semi-local density functionals for the electronic exchange for a variety of systems including atoms, molecules, and atomic chains. In particular, we focus on a recent universal extension of the Becke-Johnson exchange potential [E. Räsänen, S. Pittalis, and C. R. Proetto, J. Chem. Phys. 132, 044112 (2010)]. It is shown that when this potential is used together with the Becke-Roussel approximation to the Slater potential [A. D. Becke and M. R. Roussel, Phys. Rev. A 39, 3761 (1989)], a good overall agreement is obtained with experimental and numerically exact results for several systems, and with a moderate computational cost. Thus, this approximation is a very promising candidate in the quest for a simple and all-round semi-local potential.

preprint2009arXiv

Alloying effects on the optical properties of Ge$_{1-x}$Si$_x$ nanocrystals from TDDFT and comparison with effective-medium theory

We present the optical spectra of Ge$_{1-x}$Si$_x$ alloy nanocrystals calculated with time-dependent density-functional theory in the adiabatic local-density ap proximation (TDLDA). The spectra change smoothly as a function of the compositio n $x$. On the Ge side of the composition range, the lowest excitations at the ab sorption edge are almost pure Kohn-Sham independent-particle HOMO-LUMO transitio ns, while for higher Si contents strong mixing of transitions is found. Within T DLDA the first peak is slightly higher in energy than in earlier independent-par ticle calculations. However, the absorption onset and in particular its composit ion dependence is similar to independent-particle results. Moreover, classical depolarization effects are responsible for a very strong suppression of the abs orption intensity. We show that they can be taken into account in a simpler way using Maxwell-Garnett classical effective-medium theory. Emission spectra are in vestigated by calculating the absorption of excited nanocrystals at their relaxe d geometry. The structural contribution to the Stokes shift is about 0.5 eV. Th e decomposition of the emission spectra in terms of independent-particle transit ions is similar to what is found for absorption. For the emission, very weak tra nsitions are found in Ge-rich clusters well below the strong absorption onset.

preprint2007arXiv

Identification of fullerene-like CdSe nanoparticles from optical spectroscopy calculations

Semiconducting nanoparticles are the building blocks of optical nanodevices as their electronic states, and therefore light absorption and emission, can be controlled by modifying their size and shape. CdSe is perhaps the most studied of these nanoparticles, due to the efficiency of its synthesis, the high quality of the resulting samples, and the fact that the optical gap is in the visible range. In this article, we study light absorption of CdSe nanostructures with sizes up to 1.5 nm within density functional theory. We study both bulk fragments with wurtzite symmetry and novel fullerene-like core-cage structures. The comparison with recent experimental optical spectra allows us to confirm the synthesis of these fullerene-like CdSe clusters.