We assess the predictive accuracy of perturbation theory based estimates of changes in covalent bonding due to linear alchemical interpolations among molecules. We have investigated $\sigma$ bonding to hydrogen, as well as $\sigma$ and $\pi$ bonding between main-group elements, occurring in small sets of iso-valence-electronic molecular species with elements drawn from second to fourth rows in the $p$-block of the periodic table. Numerical evidence suggests that first order estimates of covalent bonding potentials can achieve chemical accuracy if (i) the alchemical interpolation is vertical (fixed geometry), (ii) involves molecules containing elements in the third and fourth row of the periodic table, and (iii) a reference geometry is optimized. In this case, changes in the bonding potential become near-linear in coupling parameter, resulting in analytical predictions with very high accuracy ($\sim$1 kcal/mol). Second order estimates deteriorate the prediction. If initial and final molecules differ not only in composition but also in geometry, all estimates become substantially worse, with second order being slightly more accurate than first order. The independent particle approximation to the second order perturbation performs poorly when compared to the coupled perturbed or finite difference approach. Taylor series expansions up to fourth order of the potential energy curve of highly symmetric systems indicate a finite radius of convergence, as illustrated for the alchemical stretching of H$_2^+$. Numerical results are presented for covalent bonds to hydrogen in 12 molecules with 8 valence electrons; (ii) main-group single bonds in 9 molecules with 14 valence electrons; (iii) main-group double bonds in 9 molecules with 12 valence electrons; (iv) main-group triple bonds in 9 molecules with 10 valence electrons; (v) H$_2^+$ single bond with 1 electron

In silico design of new molecules and materials with desirable quantum properties by high-throughput screening is a major challenge due to the high dimensionality of chemical space. To facilitate its navigation, we present a unification of coordinate and composition space in terms of alchemical normal modes (ANMs) which result from second order perturbation theory. ANMs assume a predominantly smooth nature of chemical space and form a basis in which new compounds can be expanded and identified. We showcase the use of ANMs for the energetics of the isoelectronic series of diatomics with 14 electrons, BN doped benzene derivatives (C 6−2x (BN) x H 6 with x = 0, 1, 2, 3), predictions for over 1.8 million BN doped coronene derivatives, and genetic energy optimizations in the entire BN doped coronene space. Using Ge lattice scans as reference, the applicability ANMs across the periodic table is demonstrated for III-V and IV-IV-semiconductors Si, Sn, SiGe, SnGe, SiSn, as well as AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs, and InSb. Analysis of our results indicates simple qualitative structure property rules for estimating energetic rankings among isomers. Useful quantitative estimates can also be obtained when few atoms are changed to neighboring or lower lying elements in the periodic table. The quality of the predictions often increases with the symmetry of system chosen as reference due to cancellation of odd order terms. Rooted in perturbation theory the ANM approach promises to generally enable unbiased compound exploration campaigns at reduced computational cost. arXiv:1809.03302v2 [physics.chem-ph]

We use alchemical first order derivatives for the rapid yet robust prediction of band structures. The power of the approach is demonstrated for the design challenge of finding AlxGa1−xAs semiconductor alloys with large direct band gap using computational alchemy within a genetic algorithm. Dozens of crystal polymorphs are identified for x > 2 3 with direct band gaps larger than 2 eV according to HSE approximated density functional theory. Based on a single generalized gradient approximated density functional theory band structure calculation of pure GaAs we observe convergence after visiting only ∼800 crystal candidates. The general applicability of alchemical gradients is demonstrated for band structure estimates in III-V and IV-IV crystals as well as for H2 uptake in Sr and Ca-alanate crystals.

'Alchemical' interpolation paths, i.e. coupling systems along fictitious paths without realistic correspondence, are frequently used within materials and molecular modeling simulation protocols for the estimation of changes in state functions such as free energies. We discuss alchemical changes in the context of quantum chemistry, and present illustrative numerical results for the changes of HOMO eigenvalue of the He atom due to alchemical teleportation - the simultaneous annihilation and creation of nuclear charges at different locations. To demonstrate the predictive power of alchemical first order derivatives (Hellmann-Feynman) the covalent bond potential of hydrogen fluoride and hydrogen chloride is investigated, as well as the hydrogen bond in the water-water and water-hydrogen fluoride dimer, respectively. Based on converged electron densities for one configuration, the versatility of alchemical derivatives is exemplified for the screening of entire binding potentials with reasonable accuracy. Finally, we discuss new constraints for the identification of non-linear coupling potentials for which the energy's Hellmann-Feynman derivative will yield accurate predictions.

No abstract