scispace - formally typeset
Search or ask a question

Showing papers on "Ab initio quantum chemistry methods published in 2015"


Journal ArticleDOI
TL;DR: The high pressure structures, metallization, and superconductivity of recently synthesized H2-containing compounds (H2S)2H2 are elucidated by ab initio calculations and the ordered crystal structure with P1 symmetry is determined.
Abstract: The high pressure structures, metallization and superconductivity of recently synthesized H2-containing compounds (H2S)2H2 are elucidated by ab initio calculations. The ordered crystal structure with P1 symmetry is determined, supported by the good agreement between theoretical and experimental X-ray diffraction data, equation of states and Raman spectra. The Cccm structure is favorable with partial hydrogen bond symmetrization above 37 GPa. Upon further compression, H2 molecules disappear and two intriguing metallic structures with R3m and Im-3m symmetries are reconstructive above 111 and 180 GPa, respectively. The predicted metallization pressure is 111 GPa, which is approximately one-third of the currently suggested metallization pressure of bulk molecular hydrogen. Application of the Allen-Dynes-modified McMillan equation for the Im-3m structure yields high Tc values of 191 K to 204 K at 200 GPa, which is among the highest values reported for H2-rich van der Waals compounds and MH3 type hydride thus far.

703 citations


Journal ArticleDOI
TL;DR: The transferability of the approach is demonstrated, using semiempirical quantum chemistry and machine learning models trained on 1 and 10% of 134k organic molecules, to reproduce enthalpies of all remaining molecules at density functional theory level of accuracy.
Abstract: Chemically accurate and comprehensive studies of the virtual space of all possible molecules are severely limited by the computational cost of quantum chemistry. We introduce a composite strategy that adds machine learning corrections to computationally inexpensive approximate legacy quantum methods. After training, highly accurate predictions of enthalpies, free energies, entropies, and electron correlation energies are possible, for significantly larger molecular sets than used for training. For thermochemical properties of up to 16k isomers of C7H10O2 we present numerical evidence that chemical accuracy can be reached. We also predict electron correlation energy in post Hartree-Fock methods, at the computational cost of Hartree-Fock, and we establish a qualitative relationship between molecular entropy and electron correlation. The transferability of our approach is demonstrated, using semiempirical quantum chemistry and machine learning models trained on 1 and 10% of 134k organic molecules, to reproduce enthalpies of all remaining molecules at density functional theory level of accuracy.

595 citations


Journal ArticleDOI
TL;DR: The present formalism generalizes the Fröhlich vertex to the case of anisotropic materials and multiple phonon branches, and enables ab initio calculations of carrier mobilities, lifetimes, mass enhancement, and pairing in polar materials.
Abstract: We develop a method for calculating the electron-phonon vertex in polar semiconductors and insulators from first principles. The present formalism generalizes the Frohlich vertex to the case of anisotropic materials and multiple phonon branches, and can be used either as a postprocessing correction to standard electron-phonon calculations, or in conjunction with ab initio interpolation based on maximally localized Wannier functions. We demonstrate this formalism by investigating the electron-phonon interactions in anatase TiO(2), and show that the polar vertex significantly reduces the electron lifetimes and enhances the anisotropy of the coupling. The present work enables ab initio calculations of carrier mobilities, lifetimes, mass enhancement, and pairing in polar materials.

260 citations


Journal ArticleDOI
TL;DR: It is shown that an electric field perpendicular to the layers can induce a semiconductor to metal transition in this family of compounds.
Abstract: We use first-principle calculations to investigate the electronic structure of InSe and In2Se3. The interlayer binding energy is found to be in the same range as for other 2D systems, and the monolayers are found to be dynamically stable, which suggest the possibility to obtain them as isolated layers. The GW approximation including spin-orbit is used to obtain the bandgaps, which are in the range relevant for application in electronics. Also, it is shown that an electric field perpendicular to the layers can induce a semiconductor to metal transition in this family of compounds.

184 citations


Journal ArticleDOI
TL;DR: The method is used to illustrate the 2D-3D transition of the dielectric function of multilayer MoS2 crystals, the hybridization of quantum plasmons in thick graphene/hBN heterostructures, and to demonstrate the intricate effect of substrate screening on the non-Rydberg exciton series in supported WS2.
Abstract: Vertical stacking of two-dimensional (2D) crystals, such as graphene and hexagonal boron nitride, has recently lead to a new class of materials known as van der Waals heterostructures (vdWHs) with unique and highly tunable electronic properties. Ab initio calculations should in principle provide a powerful tool for modeling and guiding the design of vdWHs, but in their traditional form such calculations are only feasible for commensurable structures with a few layers. Here we show that the dielectric properties of realistic, incommensurable vdWHs comprising hundreds of layers can be efficiently calculated using a multiscale approach where the dielectric functions of the individual layers (the dielectric building blocks) are computed ab initio and coupled together via the long-range Coulomb interaction. We use the method to illustrate the 2D–3D transition of the dielectric function of multilayer MoS2 crystals, the hybridization of quantum plasmons in thick graphene/hBN heterostructures, and to demonstrate ...

182 citations


Journal ArticleDOI
TL;DR: It is shown that doping of boron clusters with a transition metal atom induces an earlier two-dimensional to three-dimensional structural transition and the CoB16− cluster is tested as a building block in a triple-decker sandwich, suggesting a promising route for its realization in the solid state.
Abstract: The electron deficiency and strong bonding capacity of boron have led to a vast variety of molecular structures in chemistry and materials science. Here we report the observation of highly symmetric cobalt-centered boron drum-like structures of CoB16(-), characterized by photoelectron spectroscopy and ab initio calculations. The photoelectron spectra display a relatively simple spectral pattern, suggesting a high symmetry structure. Two nearly degenerate isomers with D8d (I) and C4v (II) symmetries are found computationally to compete for the global minimum. These drum-like structures consist of two B8 rings sandwiching a cobalt atom, which has the highest coordination number known heretofore in chemistry. We show that doping of boron clusters with a transition metal atom induces an earlier two-dimensional to three-dimensional structural transition. The CoB16(-) cluster is tested as a building block in a triple-decker sandwich, suggesting a promising route for its realization in the solid state.

178 citations


Journal ArticleDOI
TL;DR: The recent advent of the density matrix renormalization group (DMRG) theory has delivered a new capability to compute multireference (MR) wave function with large configuration space, which far exceeds the limitation of conventional approaches as discussed by the authors.
Abstract: The recent advent of the density matrix renormalization group (DMRG) theory has delivered a new capability to compute multireference (MR) wave function with large configuration space, which far exceeds the limitation of conventional approaches. Here, we provide an overview of our recent work on the developments of ab initio DMRG methods in the context of the active space approaches and their applications to MR chemical systems. © 2014 Wiley Periodicals, Inc.

174 citations


Journal ArticleDOI
TL;DR: In this article, the authors proposed an ab initio model for disordered carbon and investigated the intercalation mechanism of Na into the layered domains, which revealed that a larger interlayer distance and the presence of defects can effectively overcome the van der Waals interaction between graphene sheets and help Na inter calcation to form NaC8.
Abstract: Graphite, a predominantly chosen anode material for commercial lithium ion batteries (LIBs), has been reported to have negligible intercalation capacity as an anode for sodium ion batteries (NIBs). Disordered carbon exhibits high Na intercalation capacity and emerges as a leading candidate for NIB applications. However, the mechanism of Na+ ion insertion into disordered carbon is still controversial. Here, we propose an ab initio model for disordered carbon and investigate the intercalation mechanism of Na into the layered domains. Our ab initio calculations reveal that a larger interlayer distance and the presence of defects can effectively overcome the van der Waals interaction between graphene sheets and help Na intercalation to form NaC8. The calculation results clarify the mechanism of the Na intercalation and account for the presence of sloping and flat regions of charge–discharge curves in disordered carbon reported in numerous experiments. This reveals new prospects for helping Na intercalation into graphite.

166 citations


Journal Article
TL;DR: Applications to dimers outside the database suggests that the same scaling procedure significantly improves the FodFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.

161 citations


Journal ArticleDOI
TL;DR: In this paper, a new high-κ dielectric material was identified by automated ab initio calculations on ~1800 oxides and the cubic BeO was found to possess an unprecedented material property of 10
Abstract: Novel high-κ dielectric materials are identified by automated ab initio calculations on ~1800 oxides. The cubic BeO is found to possess an unprecedented material property of 10 eV for band gap and 275 for dielectric constant. Candidate high-κ oxides are suggested for microelectronic devices such as CPU, DRAM and flash memory.

151 citations


Journal ArticleDOI
TL;DR: Schreiber et al. as mentioned in this paper studied the performance of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, for the optical gaps and spectra of solid-state systems.
Abstract: The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green’s function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel’s widely used time-dependent density functional theory benchmark set [Schreiber et al., J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

Journal ArticleDOI
TL;DR: In this paper, the authors conduct an extensive investigation into the Hf-Ta-C system, which includes the compounds that have the highest melting points known to date, and identify three major chemical factors that contribute to the high melting temperatures.
Abstract: Using electronic structure calculations, we conduct an extensive investigation into the Hf-Ta-C system, which includes the compounds that have the highest melting points known to date. We identify three major chemical factors that contribute to the high melting temperatures. Based on these factors, we propose a class of materials that may possess even higher melting temperatures and explore it via efficient ab initio molecular dynamics calculations in order to identify the composition maximizing the melting point. This study demonstrates the feasibility of automated and high-throughput materials screening and discovery via ab initio calculations for the optimization of “higher-level” properties, such as melting points, whose determination requires extensive sampling of atomic configuration space.

Journal ArticleDOI
TL;DR: In this paper, a hybrid high-level (MP2/CBS): low-level method is used to determine adsorption structures and energies, and vibrational entropies and thermal enthalpy contributions are obtained from vibrational partition functions for the DFT+dispersion potential energy surface.
Abstract: Heats of adsorption of methane, ethane, and propane in H-chabazite (Si/Al = 14.4) have been measured and entropies have been derived from adsorption isotherms. For these systems quantum chemical ab initio calculations of Gibbs free energies have been performed. The deviations from the experimental values for methane, ethane, and propane are below 3 kJ/mol for the enthalpy, and the Gibbs free energy. A hybrid high-level (MP2/CBS): low-level (DFT+dispersion) method is used to determine adsorption structures and energies. Vibrational entropies and thermal enthalpy contributions are obtained from vibrational partition functions for the DFT+dispersion potential energy surface. Anharmonic corrections have been evaluated for each normal mode separately. One-dimensional Schrodinger equations are solved for potentials obtained by (curvilinear) distortions of the normal modes using a representation in internal coordinates.

Journal ArticleDOI
TL;DR: Theoretical energy levels of the n = 1 and n = 2 states of hydrogen-like atoms with the nuclear charge numbers 1 ≤ Z ≤ 110 are tabulated in this article.
Abstract: Theoretical energy levels of the n = 1 and n = 2 states of hydrogen-like atoms with the nuclear charge numbers 1 ≤ Z ≤ 110 are tabulated. The tabulation is based on ab initio quantum electrodynamics calculations performed to all orders in the nuclear binding strength parameter Zα, where α is the fine structure constant. Theoretical errors due to various effects are critically examined and estimated.

Journal ArticleDOI
TL;DR: In this article, a direct molecular simulation (DMS) approach is used to predict the internal energy relaxation and dissociation dynamics of high-temperature nitrogen molecules by providing forces between the four atoms.
Abstract: The direct molecular simulation (DMS) approach is used to predict the internal energy relaxation and dissociation dynamics of high-temperature nitrogen. An ab initio potential energy surface (PES) is used to calculate the dynamics of two interacting nitrogen molecules by providing forces between the four atoms. In the near-equilibrium limit, it is shown that DMS reproduces the results obtained from well-established quasiclassical trajectory (QCT) analysis, verifying the validity of the approach. DMS is used to predict the vibrational relaxation time constant for N2–N2 collisions and its temperature dependence, which are in close agreement with existing experiments and theory. Using both QCT and DMS with the same PES, we find that dissociation significantly depletes the upper vibrational energy levels. As a result, across a wide temperature range, the dissociation rate is found to be approximately 4–5 times lower compared to the rates computed using QCT with Boltzmann energy distributions. DMS calculations predict a quasi-steady-state distribution of rotational and vibrational energies in which the rate of depletion of high-energy states due to dissociation is balanced by their rate of repopulation due to collisional processes. The DMS approach simulates the evolution of internal energy distributions and their coupling to dissociation without the need to precompute rates or cross sections for all possible energy transitions. These benchmark results could be used to develop new computational fluid dynamics models for high-enthalpy flow applications.

Journal ArticleDOI
TL;DR: This work presents ab initio calculations of hot electrons in gallium arsenide (GaAs) using density functional theory and many-body perturbation theory and finds a dominant role of acoustic vibrations, challenging the common notion that optical vibrations dominate excited electron energy loss.
Abstract: Hot carrier dynamics critically impacts the performance of electronic, optoelectronic, photovoltaic, and plasmonic devices. Hot carriers lose energy over nanometer lengths and picosecond timescales and thus are challenging to study experimentally, whereas calculations of hot carrier dynamics are cumbersome and dominated by empirical approaches. In this work, we present ab initio calculations of hot electrons in gallium arsenide (GaAs) using density functional theory and many-body perturbation theory. Our computed electron–phonon relaxation times at the onset of the Γ, L, and X valleys are in excellent agreement with ultrafast optical experiments and show that the ultrafast (tens of femtoseconds) hot electron decay times observed experimentally arise from electron–phonon scattering. This result is an important advance to resolve a controversy on hot electron cooling in GaAs. We further find that, contrary to common notions, all optical and acoustic modes contribute substantially to electron–phonon scattering, with a dominant contribution from transverse acoustic modes. This work provides definitive microscopic insight into hot electrons in GaAs and enables accurate ab initio computation of hot carriers in advanced materials.

Journal ArticleDOI
TL;DR: In this article, the authors review a rather extended field of ab initio calculations of mechanical properties of materials and make an overview of advanced approaches to crystal stability, focusing on results of studies on interfaces and grain boundaries.

Journal ArticleDOI
TL;DR: The original force field for clay materials (ClayFF) developed by Cygan et al. (J. Phys. Chem. B 2004, 108, 1255) is modified to describe negative charging of the (101) quartz surface above its point of zero charge (pH ≈ 2.0-4.5) as mentioned in this paper.
Abstract: The original force field for clay materials (ClayFF) developed by Cygan et al. (J. Phys. Chem. B 2004, 108, 1255) is modified to describe negative charging of the (101) quartz surface above its point of zero charge (pH ≈ 2.0–4.5). The modified force field adopts the scaled natural bond orbital charges derived by the quantum mechanical calculations which are used to obtain the desired surface charge density and to determine the delocalization of the charge after deprotonation of surface silanol groups. Classical molecular dynamics simulations (CMD) of the (101) surface of α-quartz with different surface charge densities (0, −0.03, −0.06, and −0.12 C m–2) are performed to evaluate the influence of the negative surface charge on interfacial water and adsorption of Na+, Rb+, and Sr2+ ions. The CMD results are compared with ab initio calculations, X-ray experiment, and the triple-layer model. The modified force field can be easily implemented in common molecular dynamics packages and used for simulations of in...

Journal ArticleDOI
TL;DR: It is demonstrated that the O-H stretch SFG spectra at the water-air interface calculated by using the ssVVCF formalism well reproduce those calculated by use of the dipole moment-polarizability time correlation function.
Abstract: Interfacial water structures have been studied intensively by probing the O-H stretch mode of water molecules using sum-frequency generation (SFG) spectroscopy. This surface-specific technique is finding increasingly widespread use, and accordingly, computational approaches to calculate SFG spectra using molecular dynamics (MD) trajectories of interfacial water molecules have been developed and employed to correlate specific spectral signatures with distinct interfacial water structures. Such simulations typically require relatively long (several nanoseconds) MD trajectories to allow reliable calculation of the SFG response functions through the dipole moment-polarizability time correlation function. These long trajectories limit the use of computationally expensive MD techniques such as ab initio MD and centroid MD simulations. Here, we present an efficient algorithm determining the SFG response from the surface-specific velocity-velocity correlation function (ssVVCF). This ssVVCF formalism allows us to calculate SFG spectra using a MD trajectory of only ∼100 ps, resulting in the substantial reduction of the computational costs, by almost an order of magnitude. We demonstrate that the O-H stretch SFG spectra at the water-air interface calculated by using the ssVVCF formalism well reproduce those calculated by using the dipole moment-polarizability time correlation function. Furthermore, we applied this ssVVCF technique for computing the SFG spectra from the ab initio MD trajectories with various density functionals. We report that the SFG responses computed from both ab initio MD simulations and MD simulations with an ab initio based force field model do not show a positive feature in its imaginary component at 3100 cm(-1).

Journal ArticleDOI
TL;DR: In this article, the authors report on the observation of multiphoton interband absorption processes for quantum gases in shaken light crystals and investigate the occurrence of such excitation features with respect to the potential depth and the driving amplitude.
Abstract: We report on the observation of multiphoton interband absorption processes for quantum gases in shaken light crystals. Periodic inertial forcing, induced by a spatial motion of the lattice potential, drives multiphoton interband excitations of up to the ninth order. The occurrence of such excitation features is systematically investigated with respect to the potential depth and the driving amplitude. Ab initio calculations of resonance positions as well as numerical evaluation of their strengths exhibit good agreement with experimental data. In addition our findings could make it possible to reach novel phases of quantum matter by tailoring appropriate driving schemes.

Journal ArticleDOI
TL;DR: Support is found for the π-hole interaction involving the xenon atom from the solid state architecture of several X-ray structures retrieved from the crystal structural depot.
Abstract: In this manuscript we combine high level ab initio calculations (RI-MP2/aug-cc-pVTZ) and the analysis of several crystal structures to demonstrate the existence of π-hole aerogen bonding interactions in Xe(IV) compounds. The ability of XeF4 and Xe(OMe)4 to interact with electron rich molecules is rationalized using several computational tools, including molecular electrostatic potential surfaces, energetic and geometric features of the complexes and “atoms in molecules” (AIM) and Natural Bond Orbital (NBO) analyses. We have found support for the π-hole interaction involving the xenon atom from the solid state architecture of several X-ray structures retrieved from the crystal structural depot. Particularly, π-hole aerogen bonding interactions are quite common in the solid state of Xe(IV) compounds.

Journal ArticleDOI
TL;DR: An extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems is reported and it is shown that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space.
Abstract: The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green’s function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.

Journal ArticleDOI
TL;DR: The agreement between experiments and calculations is quite good; however, the band gaps are underestimated by calculations, with the exception of the bang gap of SrCrO4, for which theory and calculations agree.
Abstract: The crystal structures, lattice vibrations, and electronic band structures of PbCrO4, PbSeO4, SrCrO4, and SrSeO4 were studied by ab initio calculations, Raman spectroscopy, X-ray diffraction, and optical-absorption measurements. Calculations properly describe the crystal structures of the four compounds, which are isomorphic to the monazite structure and were confirmed by X-ray diffraction. Information is also obtained on the Raman- and IR-active phonons, with all of the vibrational modes assigned. In addition, the band structures and electronic densities of states of the four compounds were determined. All are indirect-gap semiconductors. In particular, chromates are found to have band gaps smaller than 2.5 eV and selenates higher than 4.3 eV. In the chromates (selenates), the upper part of the valence band is dominated by O 2p states and the lower part of the conduction band is composed primarily of electronic states associated with the Cr 3d and O 2p (Se 4s and O 2p) states. Calculations also show that the band gap of PbCrO4 (PbSeO4) is smaller than the band gap of SrCrO4 (SrSeO4). This phenomenon is caused by Pb states, which, to some extent, also contribute to the top of the valence band and the bottom of the conduction band. The agreement between experiments and calculations is quite good; however, the band gaps are underestimated by calculations, with the exception of the bang gap of SrCrO4, for which theory and calculations agree. Calculations also provide predictions of the bulk modulus of the studied compounds.

Journal ArticleDOI
TL;DR: As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including two- and three-dimensional solids.
Abstract: We propose a multireference linearized coupled cluster theory using matrix product states (MPSs-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles, for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual Density Matrix Renormalization Group (DMRG) algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here, MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including two- and three-dimensional solids.

Journal ArticleDOI
TL;DR: In this article, the Hartree-Fock method was combined with generalized gradient approximation (GGA) to tune the electronic and optically active energy gap of partially reduced graphene oxide.
Abstract: Controlled reduction of graphene oxide is an alternative and promising method to tune the electronic and optically active energy gap of this two-dimensional material in the energy range of the visible light spectrum. By means of ab initio calculations, based on hybrid density functional theory, that combine the Hartree–Fock method with the generalized gradient approximation (GGA), we investigated the electronic, optical, and radiative recombination properties of partially reduced graphene oxide, modelled as small islands of pristine graphene formed in an infinite sheet of graphene oxide. We predict that tuning of optically active gaps, in the wide range from ∼6.5 eV to ∼0.25 eV, followed by the electron radiative transition times in the range from ns to μs, can be effected by controlling the level of oxidization.

Journal ArticleDOI
TL;DR: In this article, aberrationcorrected scanning transmission electron microscopy and ab initio calculations are used to study the crystal structure of few-layer black phosphorus, and the effects of oxidation on both the atomic and electronic structure of black phosphorus are analyzed to explain observed device degradation.
Abstract: Black phosphorus, a layered two-dimensional crystal with tunable electronic properties and high hole mobility, is quickly emerging as a promising candidate for future electronic and photonic devices. Although theoretical studies using ab initio calculations have tried to predict its atomic and electronic structure, uncertainty in its fundamental properties due to a lack of clear experimental evidence continues to stymie our full understanding and application of this novel material. In this work, aberration-corrected scanning transmission electron microscopy and ab initio calculations are used to study the crystal structure of few-layer black phosphorus. Directly interpretable annular dark-field images provide a three-dimensional atomic-resolution view of this layered material in which its stacking order and all three lattice parameters can be unambiguously identified. In addition, electron energy-loss spectroscopy (EELS) is used to measure the conduction band density of states of black phosphorus, which agrees well with the results of density functional theory calculations performed for the experimentally determined crystal. Furthermore, experimental EELS measurements of interband transitions and surface plasmon excitations are also consistent with simulated results. Finally, the effects of oxidation on both the atomic and electronic structure of black phosphorus are analyzed to explain observed device degradation. The transformation of black phosphorus into amorphous PO3 or H3PO3 during oxidation may ultimately be responsible for the degradation of devices exposed to atmosphere over time.

Journal ArticleDOI
TL;DR: In this article, the authors investigated the behavior of oxygen vacancies in three different metal-oxide semiconductors (rutile and anatase TiO2, monoclinic WO3, and tetragonal ZrO2) using a hybrid density-functional method in which the fraction of exact exchange is material-dependent but obtained ab initio in a self-consistent scheme.
Abstract: We investigate the behavior of oxygen vacancies in three different metal-oxide semiconductors (rutile and anatase TiO2, monoclinic WO3, and tetragonal ZrO2) using a recently proposed hybrid density-functional method in which the fraction of exact exchange is material-dependent but obtained ab initio in a self-consistent scheme. In particular, we calculate charge-transition levels relative to the oxygen-vacancy defect and compare computed optical and thermal excitation/emission energies with the available experimental results, shedding light on the underlying excitation mechanisms and related materials properties. We find that this novel approach is able to reproduce not only ground-state properties and band structures of perfect bulk oxide materials but also provides results consistent with the optical and electrical behavior observed in the corresponding substoichiometric defective systems.

Journal ArticleDOI
TL;DR: The results obtained showed that the new method allows easy evaluation of lattice energies of inorganic ionic crystals without the need for ab initio calculations and complex calculations.
Abstract: This paper presents a new technique for estimation of lattice energies of inorganic ionic compounds using a simple formula. This new method demonstrates the relationship between chemical hardness and lattice energies of ionic compounds. Here chemical hardness values of ionic compounds are calculated via our molecular hardness equation. The results obtained using the present method and comparisons made by considering experimental data and the results from other theoretical methods in the literature showed that the new method allows easy evaluation of lattice energies of inorganic ionic crystals without the need for ab initio calculations and complex calculations.

Journal ArticleDOI
TL;DR: In this article, a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods is presented.
Abstract: Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

Journal ArticleDOI
TL;DR: In this article, the authors formulate and apply Bogoliubov coupled cluster theory to the description of open-shell nuclei up to mass ε ≈ 1.5.
Abstract: Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass $A\ensuremath{\approx}130$ on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking.Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory.Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in $m$ scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation.Results: Proof-of-principle calculations in an ${N}_{\text{max}}=6$ spherical harmonic oscillator basis for $^{16,18}\mathrm{O}$ and $^{18}\mathrm{Ne}$ in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while $^{20}\mathrm{O}$ and $^{20}\mathrm{Mg}$ display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations.Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.