Search

Accurate charge densities are essential for reliable electronic structure calculations because they significantly impact predictions of various chemical properties and in particular, according to the Hellmann-Feynman theorem, atomic forces. This study examines the accuracy of charge densities obtained from different DFT exchange-correlation functionals in comparison with coupled cluster calculations with single and double excitations. We find that modern DFT functionals can provide highly accurate charge densities, particularly in case of meta-GGA and hybrid functionals. In connection with Gaussian basis sets, it is necessary to use the largest basis sets available to obtain densitites that are nearly basis set error free. These findings highlight the importance of selecting appropriate computational methods for generating high-precision charge densities, which are for instance needed to generate reference data for training modern machine learned potentials.

We introduce a method for computing quantum mechanical forces through surface integrals over the stress tensor within the framework of density functional theory. This approach avoids the inaccuracies of traditional force calculations using the Hellmann-Feynman theorem when applied to multiresolution wavelet representations of orbitals. By integrating the quantum mechanical stress tensor over surfaces that enclose individual nuclei, we achieve highly accurate forces that exhibit superior consistency with the potential energy surface. Extensive benchmarks show that surface integrals over the stress tensor offer a robust and reliable alternative to the direct use of the Hellmann-Feynman theorem for force computations in DFT with discontinuous basis sets, particularly in cases where wavelet-based methods are employed. In addition, we integrate this approach with machine learning techniques, demonstrating that the forces obtained through surface integrals are sufficiently accurate to be used as training data for machine-learned potentials. This stands in contrast to forces calculated using the Hellmann-Feynman theorem, which do not offer this level of accuracy., Comment: 8 pages, 4 figures, 1 table

Machine learning potentials (MLP) have revolutionized the field of atomistic simulations by describing the atomic interactions with the accuracy of electronic structure methods at a small fraction of the costs. Most current MLPs construct the energy of a system as a sum of atomic energies, which depend on information about the atomic environments provided in form of predefined or learnable feature vectors. If, in addition, non-local phenomena like long-range charge transfer are important, fourth-generation MLPs need to be used, which include a charge equilibration (Qeq) step to take the global structure of the system into account. This Qeq can significantly increase the computational cost and thus can become the computational bottleneck for large systems. In this paper we present a highly efficient formulation of Qeq that does not require the explicit computation of the Coulomb matrix elements resulting in a quasi-linearly scaling method. Moreover, our approach also allows for the efficient calculation of energy derivatives, which explicitly consider the global structure-dependence of the atomic charges as obtained from Qeq. Due to its generality, the method is not restricted to MLPs but can also be applied within a variety of other force fields.

In the dynamic field of materials science, the quest to find optimal structures with low potential energy is of great significance. Over the past two decades, the minima hopping algorithm has emerged as a successful tool in this pursuit. We present a robust, user friendly and efficient implementation of the minima hopping algorithm as a Python library, enhancing in this way the global structure optimization simulations significantly. Our implementation significantly accelerates the exploration the potential energy surfaces, leveraging an MPI parallelization scheme that allows for multi level parallelization. In this scheme, multiple minima hopping processes are running simultaneously communicating their findings to a single database and, therefore, sharing information with each other about which parts of the potential energy surface have already been explored. Also multiple features from several existing implementations such as variable cell shape molecular dynamics and combined atomic position and cell geometry optimization for bulk systems, enhanced temperature feedback and fragmentation fixing for clusters are included in this implementation. Finally, this implementation takes advantage of the Atomic Simulation Environment (ASE) Python library allowing for high flexibility regarding the underlying energy and force evaluation.

This paper presents the results of an extensive structural search of ternary solids containing lutetium, nitrogen and hydrogen. Based on thousands of thermodynamically stable structures, available online, the convex hull of the formation enthalpies is constructed. To obtain the correct energetic ordering, the highly accurate RSCAN DFT functional is used in high quality all-electron calculations. In this way possible pseudopotential errors are eliminated. A novel lutetium hydride structure (HLu$_2$) that is on the convex hull is found in our search. An electron phonon analysis however shows that it is not a candidate structure for near ambient superconductivity. Besides this structure, which appears to have been missed in previous searches, possibly due to different DFT methodologies, our results agree closely with the results of previously published structure search efforts. This shows, that the field of crystal structure prediction has matured to a state where independent methodologies produce consistent and reproducible results, underlining the trustworthiness of modern crystal structure predictions. Hence it is quite unlikely that a structure, that would give rise within standard BCS theory to the superconducting properties, claimed to have been observed by Dasenbrock-Gammon et al. 10.1038/s41586-023-05742-0 , exists. This solidifies the evidence that no structure with conventional superconducting properties exists that could explain the experimental observation made by Dasenbrock-Gammon et al. 10.1038/s41586-023-05742-0

Ground state structures found in nature are in many cases of high symmetry. But structure prediction methods typically render only a small fraction of high symmetry structures. Especially for large crystalline unit cells there are many low energy defect structures. For this reason methods have been developed where either preferentially high symmetry structures are used as input or where the whole structural search is done within a certain symmetry group. In both cases it is necessary to specify the correct symmetry group beforehand. However it can in general not be predicted which symmetry group is the correct one leading to the ground state. For this reason we introduce a potential energy biasing scheme that favors symmetry and where it is not necessary to specify any symmetry group beforehand. On this biased potential energy surface, high symmetry structures will be found much faster than on an unbiased surface and independently of the symmetry group to which they belong. For our two test cases, a $C_{60}$ fullerene and bulk silicon carbide, we get a speedups of 25 and 63. In our data we also find a clear correlation between the similarity of the atomic environments and the energy. In low energy structures all the atoms of a species tend to have similar environments., Comment: 8 pages, 4 figures, 1 table

A fast and reliable geometry optimization algorithm is presented that optimizes atomic positions and lattice vectors simultaneously. Using a series of benchmarks, it is shown that the method presented in this paper outperforms in most cases the standard optimization methods implemented in popular codes such as QUANTUM ESPRESSO and VASP. To motivate the variable cell shape optimization method presented in here, the eigenvalues of the lattice Hessian matrix are investigated thoroughly. It is shown that they change depending on the shape of the cell and the number of particles inside the cell. For certain cell shapes the resulting condition number of the lattice matrix can grow quadratically with respect to the number of particles. By a coordinate transformation which can be applied to all variable cell shape optimization methods, the undesirable conditioning of the lattice Hessian matrix is eliminated.

To elucidate the geometric structure of the putative room temperature superconductor, carbonaceous sulfur hydride, at high pressure, we present the results of an extensive computational structure search of bulk C-S-H at 250 gigapascals. Using the minima hopping structure prediction method coupled to the GPU accelerated Sirius library, more than 17,000 local minima with different stochiometries in large simulation cells were investigated. Only 24 stochiometries are favourable against elemental decomposition, all of them are carbon doped H$_3$S crystals. The absence of van Hove singularities or similar peaks in the electronic density of states of more than 3.000 candidate phases rules out conventional superconductivity in C-S-H at room-temperature.