Search

Computational models are an essential tool for the design, characterization, and discovery of novel materials. Computationally hard tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their resources for simulation, analysis, and data processing. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of the computational tasks needed for materials science. In order to do that, the quantum technology must interact with conventional high-performance computing in several ways: approximate results validation, identification of hard problems, and synergies in quantum-centric supercomputing. In this paper, we provide a perspective on how quantum-centric supercomputing can help address critical computational problems in materials science, the challenges to face in order to solve representative use cases, and new suggested directions.

In this work, we test a recently developed method to enhance classical auxiliary-field quantum Monte Carlo (AFQMC) calculations with quantum computers against examples from chemistry and material science, representative of classes of industry-relevant systems. As molecular test cases, we calculate the energy curve of H4 and the relative energies of ozone and singlet molecular oxygen with respect to triplet molecular oxygen, which is industrially relevant in organic oxidation reactions. We find that trial wave functions beyond single Slater determinants improve the performance of AFQMC and allow it to generate energies close to chemical accuracy compared to full configuration interaction or experimental results. In the field of material science, we study the electronic structure properties of cuprates through the quasi-1D Fermi–Hubbard model derived from CuBr2, where we find that trial wave functions with both significantly larger fidelities and lower energies over a mean-field solution do not necessarily lead to AFQMC results closer to the exact ground state energy. [ABSTRACT FROM AUTHOR]

Current machine-learning methods to reproduce ab initio potential energy landscapes suffer from an unfavorable computational scaling with respect to the number of chemical species. In this work, we propose a new approach by using optimized symmetry functions to explore similarities of structures in multicomponent systems in order to yield linear complexity. We combine these symmetry functions with the charge equilibration via neural network technique, a reliable artificial neural network potential for ionic materials, and apply this method to study alkali-halide materials MX with 6 chemical species (M = {Li, Na, K} and X = {F, Cl, Br}). Our results show that our approach provides good agreement both with experimental and DFT reference data of many physical and structural properties for any chemical combination. [ABSTRACT FROM AUTHOR]

Measuring similarities/dissimilarities between atomic structures is important for the exploration of potential energy landscapes. However, the cell vectors together with the coordinates of the atoms, which are generally used to describe periodic systems, are quantities not directly suitable as fingerprints to distinguish structures. Based on a characterization of the local environment of all atoms in a cell, we introduce crystal fingerprints that can be calculated easily and define configurational distances between crystalline structures that satisfy the mathematical properties of a metric. This distance between two configurations is a measure of their similarity/dissimilarity and it allows in particular to distinguish structures. The new method can be a useful tool within various energy landscape exploration schemes, such as minima hopping, random search, swarm intelligence algorithms, and high-throughput screenings. [ABSTRACT FROM AUTHOR]

Intermetallic Li-Al compounds are on the one hand key materials for light-weight engineering, and on the other hand, they have been proposed for high-capacity electrodes for Li batteries. We determine from first-principles the phase diagram of Li-Al binary crystals using the minima hopping structural prediction method. Beside reproducing the experimentally reported phases (LiAl, Li3Al2, Li9Al4, LiAl3, and Li2Al), we unveil a structural variety larger than expected by discovering six unreported binary phases likely to be thermodynamically stable. Finally, we discuss the behavior of the elastic constants and of the electric potential profile of all Li-Al stable compounds as a function of their stoichiometry. [ABSTRACT FROM AUTHOR]

The Minima Hopping global optimization method uses physically realizable molecular dynamics moves in combination with an energy feedback that guarantees the escape from any potential energy funnel. For the purpose of finding reaction pathways, we argue that Minima Hopping is particularly suitable as a guide through the potential energy landscape and as a generator for pairs of minima that can be used as input structures for methods capable of finding transition states between two minima. For Lennard-Jones benchmark systems we compared this Minima Hopping guided path search method to a known approach for the exploration of potential energy landscapes that is based on deterministic mode-following. Although we used a stabilized mode-following technique that reliably allows to follow distinct directions when escaping from a local minimum, we observed that Minima Hopping guided path search is far superior in finding lowest-barrier reaction pathways.We, therefore, suggest that Minima Hopping guided path search can be used as a simple and efficient way to identify energetically low-lying chemical reaction pathways. Finally, we applied the Minima Hopping guided path search approach to 75-atom and 102-atom Lennard-Jones systems. For the 75-atom system we found pathways whose highest energies are significantly lower than the highest energy along the previously published lowest-barrier pathway. Furthermore, many of these pathways contain a smaller number of intermediate transition states than the previously publish lowest-barrier pathway. In case of the 102-atom system Minima Hopping guided path search found a previously unknown and energetically low-lying funnel. [ABSTRACT FROM AUTHOR]

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH4)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 ab initio calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with C2221 symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minimahopping 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(BH4)4. [ABSTRACT FROM AUTHOR]

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH4)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 ab initio calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with C2221 symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minimahopping 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(BH4)4. [ABSTRACT FROM AUTHOR]

A structure prediction method is presented based on the minima hopping method. To escape local minima, moves on the configurational enthalpy surface are performed by variable cell shape molecular dynamics. To optimize the escape steps the initial atomic and cell velocities are aligned to low curvature directions of the current local minimum. The method is applied to both silicon crystals and well-studied binary Lennard-Jones mixtures. For the latter new putative ground state structures are presented. It is shown that a high success rate is achieved and a reliable prediction of unknown ground state structures is possible. [ABSTRACT FROM AUTHOR]

In this study, we discovered a method to synthesize large single crystals of Cu-BTC MOFs (Cu-s) with the assistance of a surfactant (PF-127) in an ethanol–water mixture solution. Single crystal and powder XRD (X-ray diffraction) measurements revealed that Cu-s, with the appearance of long rods, has the chemical formula Cu(BTC)3·H2O. We proposed the growth mechanism of Cu-s by a series of experiments and demonstrated that the addition of surfactant PF-127 can effectively inhibit the contact between Cu2+ ions and BTC ligands, which favoured the growth of Cu-s along the a-axis and suppressed the fast formation of porous Cu-BTC MOFs (HKUST-1). Considering that the photoactive properties of Cu-s and its relatively large size (an average length of up to ∼3 mm), we constructed a simple photo-detector device, showing excellent photo-sensing properties under illumination. [ABSTRACT FROM AUTHOR]

We report on the discovery of a quasi two-dimensional copper–bismuth nano sheet from ab initio calculations, which we call cubine. According to our predictions, single layers of cubine can be isolated from the recently reported high-pressure CuBi bulk material at an energetic cost of merely ≈20 meV/Ų, comparable to values to separate single layers of graphene from graphite. Our calculations suggest that cubine has remarkable electronic and electrochemical properties: It is a superconductor with a moderate electron–phonon coupling λ = 0.5, leading to Tc ≈ 1 K, and can be readily intercalated with lithium with a high diffusibility, rendering it a promising candidate material to boost the rate capacity of current electrodes in lithium-ion batteries. [ABSTRACT FROM AUTHOR]

Exploration beyond the known phase space of thermodynamically stable compounds into the realm of metastable materials is a frontier of materials chemistry. The application of high pressure in experiment and theory provides a powerful vector by which to explore this uncharted phase space, allowing discovery of complex new structures and bonding in the solid state. We harnessed this approach for the Cu-Bi system, where the realization of new phases offers potential for exotic properties such as superconductivity. This potential is due to the presence of bismuth, which, by virtue of its status as one of the heaviest stable elements, forms a critical component in emergent materials such as superconductors and topological insulators. To fully investigate and understand the Cu-Bi system, we welded theoretical predictions with experiment to probe the Cu-Bi system under high pressures. By employing the powerful approach of in situ X-ray diffraction in a laser-heated diamond anvil cell (LHDAC), we thoroughly explored the high-pressure and high-temperature (high-PT) phase space to gain insight into the formation of intermetallic compounds at these conditions. We employed density functional theory (DFT) calculations to calculate a pressure versus temperature phase diagram, which correctly predicts that CuBi is stabilized at lower pressures than Cu11Bi7, and allows us to uncover the thermodynamic contributions responsible for the stability of each phase. Detailed comparisons between the NiAs structure type and the two high-pressure Cu-Bi phases, Cu11Bi7 and CuBi, reveal the preference for elemental segregation within the Cu-Bi phases, and highlight the unique channels and layers formed by ordered Cu vacancies. The electron localization function from DFT calculations account for the presence of these "voids" as a manifestation of the lone pair orientation on the Bi atoms. Our study demonstrates the power of joint experimental-computational work in exploring the chemistry occurring at high-PT conditions. The existence of multiple high-pressure-stabilized phases in the Cu-Bi binary system, which can be readily identified with in situ techniques, offers promise for other systems in which no ambient pressure phases are known to exist. [ABSTRACT FROM AUTHOR]

Complex multinary compounds (ternary, quaternary, and higher) offer countless opportunities for discovering new semiconductors for applications such as photovoltaics and thermoelectrics. However, controlling doping has been a major challenge in complex semiconductors as there are many possibilities for charged intrinsic defects (e.g., vacancies, interstitials, antisite defects) whose energy depends on competing impurity phases. Even in compounds with no apparent deviation from a stoichiometric nominal composition, such defects commonly lead to free carrier concentrations in excess of 1020 cm−3. Nevertheless, by slightly altering the nominal composition, these defect concentrations can be tuned with small variation of the chemical potentials (composition) of each element. While the variation of chemical composition is undetectable, it is shown that the changes can be inferred by mapping (in nominal composition space) the boundaries where different competing impurity phases form. In the inexpensive Zintl compound Ca9Zn4+ xSb9, the carrier concentrations can be finely tuned within three different three-phase regions by altering the nominal composition ( x = 0.2-0.8), enabling the doubling of thermoelectric performance (zT). Because of the low thermal conductivity, the zT can reach as high as 1.1 at 875 K, which is one of the highest among the earth abundant p-type thermoelectrics with no ion conducting. [ABSTRACT FROM AUTHOR]

A new intermetallic compound, the first to be structurally identified in the Cu−Bi binary system, is reported. This compound is accessed by high-pressure reaction of the elements. Its detailed characterization, physical property measurements, and ab initio calculations are described. The commensurate crystal structure of Cu11Bi7 is a unique variation of the NiAs structure type. Temperature-dependent electrical resistivity and heat capacity measurements reveal a bulk superconducting transition at Tc=1.36 K. Density functional theory calculations further demonstrate that Cu11Bi7 can be stabilized (relative to decomposition into the elements) at high pressure and temperature. These results highlight the ability of high-pressure syntheses to allow for inroads into heretofore-undiscovered intermetallic systems for which no thermodynamically stable binaries are known. [ABSTRACT FROM AUTHOR]

ZnxSb1-x systems are one of the most promising materials for thermoelectric applications due to their low cost, low toxicity, and high thermoelectric figure of merit. We report the discovery and characterization of novel ZnSb polymorphs with improved thermoelectric properties from a systematic ab initio structural search. The putative ground state at low temperature is a structure with Pbca symmetry. Corresponding electron and phonon transport properties are computed, where the lattice thermal conductivity was accurately obtained through third-order interatomic force constants. At 500 K, both the reduced figure of merit S2/ L as well as the lattice thermal conductivity are superior to any previously reported numbers on ZnSb polymorphs, rendering the Pbca phase a promising candidate for an exceptionally good thermoelectric material. [ABSTRACT FROM AUTHOR]

Hydrogen-rich compounds have been extensively studied both theoretically and experimentally in the quest for novel high-temperature superconductors. Reports on sulfur hydride attaining metallicity under pressure and exhibiting superconductivity at temperatures as high as 200 K have spurred an intense search for room-temperature superconductors in hydride materials. Recently, compressed phosphine was reported to metallize at pressures above 45 GPa, reaching a superconducting transition temperature (TC) of 100 K at 200 GPa. However, neither the exact composition nor the crystal structure of the superconducting phase have been conclusively determined. In this work, the phase diagram of PHn (n = 1,2,3,4,5,6) was extensively explored by means of ab initio crystal structure predictions using the minima hopping method (MHM). The results do not support the existence of thermodynamically stable PHn compounds, which exhibit a tendency for elemental decomposition at high pressure even when vibrational contributions to the free energies are taken into account. Although the lowest energy phases of PH1,2,3 display TC's comparable to experiments, it remains uncertain if the measured values of TC can be fully attributed to a phase-pure compound of PHn. [ABSTRACT FROM AUTHOR]

Intermetallic compounds made of alkali metals and gold have intriguing electronic and structural properties that have not been extensively explored. We perform a systematic study of the phase diagram of one binary system belonging to this family, namely NaxAu1−x, using the ab initio minima hopping structural prediction method. We discover that the most stable composition is NaAu2, in agreement with available experimental data. We also confirm the crystal structures of NaAu2 and Na2Au, that were fully characterized in experiments, and identify a candidate ground-state structure for the experimental stoichiometry NaAu. Moreover, we obtain three other stoichiometries, namely Na3Au2, Na3Au and Na5Au, that could be thermodynamically stable. We do not find any evidence for the existence of the experimentally proposed composition NaAu5. Finally, we perform phonon calculations to check the dynamical stability of all reported phases and we simulate x-ray diffraction spectra for comparison with future experimental data. [ABSTRACT FROM AUTHOR]

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, monolayer 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 nanostructures. [ABSTRACT FROM AUTHOR]

We study the thermodynamic stability at low temperatures of a series of alkali-metal-zinc double-cation borohydrides, including LiZn(BH4)3, LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5, KZn(BH4)3, and KZn2(BH4)5. While LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5, and KZn(BH4)3 were recently synthesized, LiZn(BH4)3 and KZn2(BH4)5 are hypothetical compounds. Using the minima-hopping method, we discover two new lowest-energy structures for NaZn(BH4)3 and KZn2(BH4)5 which belong to the C2/c and P2 space groups, respectively. These structures are predicted to be both thermodynamically stable and dynamically stable, implying that their existence may be possible. On the other hand, we predict that the lowest-energy P1 structure of LiZn(BH4)3 is unstable, suggesting a possible reason elucidating why this compound has not been experimentally identified. In exploring the low-energy structures of these compounds, we find that their energetic ordering is sensitive to the inclusion of the van der Waals interactions. We also find that a proper treatment of these interactions, e.g., as given by a nonlocal density functional such as vdW-DF2, is necessary to address the stability of the low-energy structures of these compounds. [ABSTRACT FROM AUTHOR]

We performed a systematic structural search of high-pressure carbon allotropes for unit cells containing from 6 to 24 atoms using the minima hopping method. We discovered a series of new structures that are consistently lower in enthalpy than the ones previously reported. Most of these include (5+7)- or (4+8)-membered rings and can therefore be placed in the families proposed by H. Niu et al. [ Phys. Rev. Lett. 108 135501 (2012)]. However, we also found three more families with competitive enthalpies that contain (5+5+8)-membered rings, sp2 motives, or buckled hexagons. These structures are likely to play an important role in dislocation planes and structural defects of diamond and hexagonal diamond. [ABSTRACT FROM AUTHOR]

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

We present a systematic study of the low-energy structures of zinc borohydride, a crystalline material proposed for the purpose of hydrogen storage. In addition to previously proposed structures, many new low-energy structures of zinc borohydride are found by utilizing the minima-hopping method. We identify a new dynamically stable structure which belongs to the I4122 space group as the lowest-energy phase of zinc borohydride at low temperatures. A low transition barrier between I4122 and P1, the two lowest-lying phases of zinc borohydride, is predicted, implying that a coexistence of low-energy phases of zinc borohydride is possible at ambient conditions. An analysis based on the simulated x-ray-diffraction pattern reveals that the I4122 structure exhibits the same major features as the experimentally synthesized zinc borohydnde samples. [ABSTRACT FROM AUTHOR]

We identify a class of novel low energy phases of the hydrogen storage material LiAlH4 by using the ab initio minima hopping crystal structure prediction method. These phases are, unlike previous predictions and known structures of similar materials, characterized by polymeric networks consisting of Al atoms interlinked with H atoms. The most stable structure is a layered ionic crystal with P21/c symmetry, and it has lower free energy than the previously reported structure over a wide range of temperatures. Furthermore, we carry out x-ray diffraction, phonon, and GW band-structure analysis in order to characterize this phase. Its experimental synthesis would have profound implications for the study of dehydrogenation and rehydrogenation processes and the stability problem of LiAlH4 for hydrogen storage applications. [ABSTRACT FROM AUTHOR]

A systematic ab initio search for low-enthalpy phases of disilane (Si2H6) at high pressures was performed based on the minima hopping method. We found a novel metallic phase of disilane with Cmcm symmetry, which is enthalpically more favorable than the recently proposed structures of disilane up to 280 GPa, but revealing compositional instability below 190 GPa. The Cmcm phase has a moderate electron-phonon coupling yielding a superconducting transition temperature Tc of around 20 K at 100 GPa, decreasing to 13 K at 220 GPa. These values are significantly smaller than previously predicted Tc's for disilane at equivalent pressure. This shows that similar but different crystalline structures of a material can result in dramatically different Tc's and stresses the need for a systematic search for a crystalline ground state. [ABSTRACT FROM AUTHOR]

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 it provides an explanation to several 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. [ABSTRACT FROM AUTHOR]

The authors discuss research on a metastable cubic silicon phase that was theoretically predicted by the modified particle swarm optimization method. They reference the study "Towards Direct-Gap Silicon Phases by the Inverse Band Structure Design Approach" by H. J. Xiang and colleagues published in the journal in 2013. They believe that optimizing a single physical property alone when designing novel functional materials is not sufficient to predict promising structures.

Semiconducting half and, to a lesser extent, full Heusler compounds are promising thermoelectric materials due to their compelling electronic properties with large power factors. However, intrinsically high thermal conductivity resulting in a limited thermoelectric efficiency has so far impeded their widespread use in practical applications. Here, we report the computational discovery of a class of hitherto unknown stable semiconducting full Heusler compounds with ten valence electrons (X2YZ, X=Ca, Sr, and Ba; Y=Au and Hg; Z=Sn, Pb, As, Sb, and Bi) through high-throughput ab initio screening. These new compounds exhibit ultralow lattice thermal conductivity κL close to the theoretical minimum due to strong anharmonic rattling of the heavy noble metals, while preserving high power factors, thus resulting in excellent phonon-glass electron-crystal materials. [ABSTRACT FROM AUTHOR]

Novel ZnSb polymorphs with improved thermoelectric (TE) properties are discovered by a systematic ab initio structural search based on DFT calculations revealing the glass-like energy landscape and polymorphism of the system. [ABSTRACT FROM AUTHOR]

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 quasidirect band gaps in the range of 1.0-1.8 eV, close to the optimal Shockley-Queisser limit of 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. [ABSTRACT FROM AUTHOR]

We present a systematic study on hydrogenated silicon m vacancies above saturation. For each system a global geometry optimization search for low-lying local minima is performed using a newly developed SiH tight binding model. Subsequently a large number of low-energy structures are examined by density functional calculations using a minimal basis set. Finally the energetically favorable structures are reexamined using a systematically extendable basis set with local, semilocal, and hybrid exchange-correlation functionals. Particular attention is paid to the divacancy to which the Raman peak at 3822 cm-1 associated with the H2 molecule had previously been assigned. Both the energetics and vibrational analysis of divacancy-related stable configurations suggest a revision of the above conclusion. [ABSTRACT FROM AUTHOR]

We present unexpected structural motifs for boron-carbon nanocages of the stoichiometries B12C48 and B12C50, based on first-principles calculations. These configurations are distinct from those proposed so far because the boron atoms are not isolated and distributed over the entire surface of the cages, but rather aggregate at one location to form a patch. Our putative ground state of B12C48 is 1.8 eV lower in energy than the previously proposed ground state and violates all the suggested empirical rules for constructing low-energy fullerenes. The B12C50 configuration is energetically even more favorable than B12C48, showing that structures derived from the C60 buckminsterfullerene are not necessarily magic sizes for heterofullerenes. [ABSTRACT FROM AUTHOR]

We present state-of-the-art first-principles 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 sp3 silicon networks and have a lower formation energy than other experimentally produced silicon phases. Some of these structures turned out to have quasidirect 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. [ABSTRACT FROM AUTHOR]

Raman spectroscopy is a powerful tool to study the intrinsic vibrational characteristics of crystals, and, therefore, it is an adequate technique to explore phase transitions of carbon under pressure. However, the diamond-anvil cell, which is used in experiments to apply pressure, appears as a broad intense feature in the spectra. This feature lies, unfortunately, in the same range as the principal modes of recently proposed sp3 carbon structures. As these modes are hard to distinguish from the diamond cell background, we analyze all Raman-active modes present in the sp3 carbon structures in order to find detectable fingerprint features for an experimental identification. [ABSTRACT FROM AUTHOR]

We report the ab initio discovery of a novel putative ground state for quasi two-dimensional TiO2 through a structural search using the minima hopping method with an artificial neural network potential. The structure is based on a honeycomb lattice and is energetically lower than the experimentally reported lepidocrocite sheet by 7 meV/atom and merely 13 meV/atom higher in energy than the rutile bulk structure. According to our calculations, the hexagonal sheet is stable against mechanical stress, chemically inert, and can be deposited on various substrates without disrupting the structure. Its properties differ significantly from all known TiO2 bulk phases with a large gap of 5.05 eV that can be tuned through strain and defect engineering. [ABSTRACT FROM AUTHOR]