Researcher profile

Miguel A. L. Marques

Miguel A. L. Marques contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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).