No Arabic abstract
Cathodes are critical components of rechargeable batteries. Conventionally, the search for cathode materials relies on experimental trial-and-error and a traversing of existing computational/experimental databases. While these methods have led to the discovery of several commercially-viable cathode materials, the chemical space explored so far is limited and many phases will have been overlooked, in particular those that are metastable. We describe a computational framework for battery cathode exploration, based on ab initio random structure searching (AIRSS), an approach that samples local minima on the potential energy surface to identify new crystal structures. We show that, by delimiting the search space using a number of constraints, including chemically aware minimum interatomic separations, cell volumes, and space group symmetries, AIRSS can efficiently predict both thermodynamically stable and metastable cathode materials. Specifically, we investigate LiCoO2, LiFePO4, and LixCuyFz to demonstrate the efficiency of the method by rediscovering the known crystal structures of these cathode materials. The effect of parameters, such as minimum separations and symmetries, on the efficiency of the sampling is discussed in detail. The adaptation of the minimum interatomic distances, on a species-pair basis, from low-energy optimized structures to efficiently capture the local coordination environment of atoms, is explored. A family of novel cathode materials based, on the transition-metal oxalates,is proposed. They demonstrate superb energy density, oxygen-redox stability, and lithium diffusion properties. This article serves both as an introduction to the computational framework, and as a guide to battery cathode material discovery using AIRSS.
Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, a quantum computer-based AIMD method is presented. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. We compute the energy gradients numerically using the Hellmann-Feynman theorem, finite differences, and a correlated sampling technique. Our method only requires additional classical calculations of electron integrals for each degree of freedom, without any additional computations on a quantum computer beyond the initial VQE run. To achieve comparable accuracy, our gradient calculation method requires three to five orders of magnitude fewer measurements than other brute force methods without correlated sampling. As a proof of concept, AIMD dynamics simulations are demonstrated for the H2 molecule on IBM quantum devices. To the best of our knowledge, it is the first successful attempt to run AIMD on quantum devices for a chemical system. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction (FCI) wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular and material systems.
The concept of machine learning configuration interaction (MLCI) [J. Chem. Theory Comput. 2018, 14, 5739], where an artificial neural network (ANN) learns on the fly to select important configurations, is further developed so that accurate ab initio potential energy curves can be efficiently calculated. This development includes employing the artificial neural network also as a hash function for the efficient deletion of duplicates on the fly so that the singles and doubles space does not need to be stored and this barrier to scalability is removed. In addition configuration state functions are introduced into the approach so that pure spin states are guaranteed, and the transferability of data between geometries is exploited. This improved approach is demonstrated on potential energy curves for the nitrogen molecule, water, and carbon monoxide. The results are compared with full configuration interaction values, when available, and different transfer protocols are investigated. It is shown that, for all of the considered systems, accurate potential energy curves can now be efficiently computed with MLCI. For the potential curves of N$_{2}$ and CO, MLCI can achieve lower errors than stochastically selecting configurations while also using substantially less processor hours.
We examined the reliability of exchange-correlation functionals for molecular encapsulations combined by van der Waals forces, comparing their predictions with those of diffusion Monte Carlo method. We established that functionals with D3 dispersion force correction and including sufficient proportion of exact-exchange in long-ranged interaction can comparatively reliably estimate the binding strength. Our finding agrees with a previous ab initio study on argon dimer. However we found that even such functionals may not be able to distinguish the energy differences among different conformations.
We present four-dimensional ab initio potential energy surfaces for the three spin states of the NH-NH complex. The potentials are partially based on the work of Dhont et al. [J. Chem. Phys. 123, 184302 (2005)]. The surface for the quintet state is obtained at the RCCSD(T)/aug-cc-pVTZ level of theory and the energy diferences with the singlet and triplet states are calculated at the CASPTn/aug-cc-pVTZ (n = 2; 3) level of theory. The ab initio potentials are fitted to coupled spherical harmonics in the angular coordinates, and the long range is further expanded as a power series in 1/R. The RCCSD(T) potential is corrected for a size-consistency error prior to fitting. The long-range coeficients obtained from the fit are found to be in good agreement with perturbation theory calculations.
Reliable quantum chemical methods for the description of molecules with dense-lying frontier orbitals are needed in the context of many chemical compounds and reactions. Here, we review developments that led to our newcomputational toolbo x which implements the quantum chemical density matrix renormalization group in a second-generation algorithm. We present an overview of the different components of this toolbox.