Quantum-centric high performance computing for quantum chemistry.

High performance computing (HPC) is renowned for its capacity to tackle complex problems. Meanwhile, quantum computing (QC) provides a potential way to accurately and efficiently solve quantum chemistry problems. The emerging field of quantum-centric high performance computing (QCHPC), which merges these two powerful technologies, is anticipated to enhance computational capabilities for solving challenging problems in quantum chemistry. The implementation of QCHPC for quantum chemistry requires interdisciplinary research and collaboration across multiple fields, including quantum chemistry, quantum physics, computer science and so on. This perspective provides an introduction to the quantum algorithms that are suitable for deployment in QCHPC, focusing on conceptual insights rather than technical details. Parallel strategies to implement these algorithms on quantum-centric supercomputers are discussed. We also summarize high performance quantum emulating simulators, which are considered a viable tool to explore QCHPC. We conclude with challenges and outlooks in this field.

Transient Photoinduced Pb2+ Disproportionation for Exciton Self-Trapping and Broadband Emission in Low-Dimensional Lead Halide Perovskites.

Low-dimensional lead halide perovskites with broadband emission hold great promise for single-component white-light-emitting (WLE) devices. The origin of their broadband emission has been commonly attributed to self-trapped excitons (STEs) composed of localized electronic polarization with a distorted lattice. Unfortunately, the exact electronic and structural nature of the STE species in these WLE materials remains elusive, hindering the rational design of high-efficiency WLE materials. In this study, by combining ultrafast transient absorption spectroscopy and ab initio calculations, we uncover surprisingly similar STE features in two prototypical low dimensional WLE perovskite single crystals: 1D (DMEDA)PbBr4 and 2D (EDBE)PbBr4, despite of their different dimensionalities. Photoexcited excitons rapidly localize to intrinsic STEs within ∼250 fs, contributing to the white light emission. Crucially, STEs in both systems exhibit characteristic absorption features akin to those of Pb+ and Pb3+. Further atomic level theoretical simulations confirm photoexcited electrons and holes are localized on the Pb2+ site to form Pb+- and Pb3+-like species, resembling transient photoinduced Pb2+ disproportionation. This study provides conclusive evidence on the key excited state species for exciton self-trapping and broadband emission in low dimensional lead halide WLE perovskites and paves the way for the rational design of high-efficiency WLE materials.

Efficient implementation of analytical gradients for periodic hybrid functional calculations within fitted numerical atomic orbitals from NAO2GTO.

The NAO2GTO scheme provides an efficient way to evaluate the electron repulsion integrals (ERIs) over numerical atomic orbitals (NAOs) with auxiliary Gaussian-type orbitals (GTOs). However, the NAO2GTO fitting will significantly impact the accuracy and convergence of hybrid functional calculations. To address this issue, here we propose to use the fitted orbitals as a new numerical basis to properly handle the mismatch between NAOs and fitted GTOs. We present an efficient and linear-scaling implementation of analytical gradients of Hartree-Fock exchange (HFX) energy for periodic HSE06 calculations with fitted NAOs in the HONPAS package. In our implementation, the ERIs and their derivatives for HFX matrix and forces are evaluated analytically with the auxiliary GTOs, while other terms are calculated using numerically discretized GTOs. Several integral screening techniques are employed to reduce the number of required ERI derivatives. We benchmark the accuracy and efficiency of our implementation and demonstrate that our results of lattice constants, bulk moduli, and band gaps of several typical semiconductors are in good agreement with the experimental values. We also show that the calculation of HFX forces based on a master-worker dynamic parallel scheme has a very high efficiency and scales linearly with respect to system size. Finally, we study the geometry optimization and polaron formation due to an excess electron in rutile TiO2 by means of HSE06 calculations to further validate the applicability of our implementation.

Corrigendum: OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources.

[This corrects the article DOI: 10.3389/fchem.2023.1156891.].

OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources.

We have proposed, for the first time, an OpenCL implementation for the all-electron density-functional perturbation theory (DFPT) calculations in FHI-aims, which can effectively compute all its time-consuming simulation stages, i.e., the real-space integration of the response density, the Poisson solver for the calculation of the electrostatic potential, and the response Hamiltonian matrix, by utilizing various heterogeneous accelerators. Furthermore, to fully exploit the massively parallel computing capabilities, we have performed a series of general-purpose graphics processing unit (GPGPU)-targeted optimizations that significantly improved the execution efficiency by reducing register requirements, branch divergence, and memory transactions. Evaluations on the Sugon supercomputer have shown that notable speedups can be achieved across various materials.

Towards practical and massively parallel quantum computing emulation for quantum chemistry.

Quantum computing is moving beyond its early stage and seeking for commercial applications in chemical and biomedical sciences. In the current noisy intermediate-scale quantum computing era, the quantum resource is too scarce to support these explorations. Therefore, it is valuable to emulate quantum computing on classical computers for developing quantum algorithms and validating quantum hardware. However, existing simulators mostly suffer from the memory bottleneck so developing the approaches for large-scale quantum chemistry calculations remains challenging. Here we demonstrate a high-performance and massively parallel variational quantum eigensolver (VQE) simulator based on matrix product states, combined with embedding theory for solving large-scale quantum computing emulation for quantum chemistry on HPC platforms. We apply this method to study the torsional barrier of ethane and the quantification of the protein-ligand interactions. Our largest simulation reaches 1000 qubits, and a performance of 216.9 PFLOP/s is achieved on a new Sunway supercomputer, which sets the state-of-the-art for quantum computing emulation for quantum chemistry.

The electron-phonon renormalization in the electronic structure calculation: Fundamentals, current status, and challenges.

Electron-phonon (e-ph) interaction plays a crucial role in determining many physical properties of the materials, such as the superconducting transition temperature, the relaxation time and mean free path of hot carriers, the temperature dependence of the electronic structure, and the formation of the vibrational polaritons. In the past two decades, the calculations of e-ph properties from first-principles has become possible. In particular, the renormalization of electronic structures due to e-ph interaction can be evaluated, providing greater insight into the quantum zero-point motion effect and the temperature dependence behavior. In this perspective, we briefly overview the basic theory, outline the computational challenges, and describe the recent progress in this field, as well as future directions and opportunities of the e-ph coupling calculations.

Multiscale quantum algorithms for quantum chemistry.

Exploring the potential applications of quantum computers in material design and drug discovery is attracting more and more attention after quantum advantage has been demonstrated using Gaussian boson sampling. However, quantum resource requirements in material and (bio)molecular simulations are far beyond the capacity of near-term quantum devices. In this work, multiscale quantum computing is proposed for quantum simulations of complex systems by integrating multiple computational methods at different scales of resolution. In this framework, most computational methods can be implemented in an efficient way on classical computers, leaving the critical portion of the computation to quantum computers. The simulation scale of quantum computing strongly depends on available quantum resources. As a near-term scheme, we integrate adaptive variational quantum eigensolver algorithms, second-order Møller-Plesset perturbation theory and Hartree-Fock theory within the framework of the many-body expansion fragmentation approach. This new algorithm is applied to model systems consisting of hundreds of orbitals with decent accuracy on the classical simulator. This work should encourage further studies on quantum computing for solving practical material and biochemistry problems.

Ultrafast Spontaneous Localization of a Jahn-Teller Exciton Polaron in Two-Dimensional Semiconducting CrI3 by Symmetry Breaking.

The excited state species and properties in low-dimensional semiconductors can be completely redefined by electron-lattice coupling or a polaronic effect. Here, by combining ultrafast broadband pump-probe spectroscopy and first-principles GW and Bethe-Salpeter equation calculations, we show semiconducting CrI3 as a prototypical 2D polaronic system with characteristic Jahn-Teller exciton polaron induced by symmetry breaking. A photogenerated electron and hole in CrI3 localize spontaneously in ∼0.9 ps and pair geminately to a Jahn-Teller exciton polaron with elongated Cr-I octahedra, large binding energy, and an unprecedentedly small exciton-exciton annihilation rate constant (∼10-20 cm3 s-1). Coherent phonon dynamics indicates the localization is mainly triggered by the coherent nuclear vibration of the I-Cr-I out-of-plane stretch mode at 128.5 ± 0.1 cm-1. The excited state Jahn-Teller exciton polaron in CrI3 broadens the realm of 2D polaron systems and reveals the decisive role of coupled electron-lattice motion on excited state properties and exciton physics in 2D semiconductors.

Localized Resolution of Identity Approach to the Analytical Gradients of Random-Phase Approximation Ground-State Energy: Algorithm and Benchmarks.

We develop and implement a formalism which enables calculating the analytical gradients of particle-hole random-phase approximation (RPA) ground-state energy with respect to the atomic positions within the atomic orbital basis set framework. Our approach is based on a localized resolution of identity (LRI) approximation for evaluating the two-electron Coulomb integrals and their derivatives, and the density functional perturbation theory for computing the first-order derivatives of the Kohn-Sham (KS) orbitals and orbital energies. Our implementation allows one to relax molecular structures at the RPA level using both Gaussian-type orbitals (GTOs) and numerical atomic orbitals (NAOs). Benchmark calculations against previous implementations show that our approach delivers adequate numerical precision, highlighting the usefulness of LRI in the context of RPA gradient evaluations. A careful assessment of the quality of RPA geometries for small molecules reveals that post-KS RPA systematically overestimates the bond lengths. We furthermore optimized the geometries of the four low-lying water hexamers-cage, prism, cyclic, and book isomers, and determined the energy hierarchy of these four isomers using RPA. The obtained RPA energy ordering is in good agreement with that yielded by the coupled cluster method with single, double and perturbative triple excitations, despite that the dissociation energies themselves are appreciably underestimated. The underestimation of the dissociation energies by RPA is well corrected by the renormalized single excitation correction.

Coupled Electronic and Anharmonic Structural Dynamics for Carrier Self-Trapping in Photovoltaic Antimony Chalcogenides.

V-VI antimony chalcogenide semiconductors have shown exciting potentials for thin film photovoltaic applications. However, their solar cell efficiencies are strongly hampered by anomalously large voltage loss (>0.6 V), whose origin remains controversial so far. Herein, by combining ultrafast pump-probe spectroscopy and density functional theory (DFT) calculation, the coupled electronic and structural dynamics leading to excited state self-trapping in antimony chalcogenides with atomic level characterizations is reported. The electronic dynamics in Sb2 Se3 indicates a ≈20 ps barrierless intrinsic self-trapping, with electron localization and accompanied lattice distortion given by DFT calculations. Furthermore, impulsive vibrational coherences unveil key SbSe vibrational modes and their real-time interplay that drive initial excited state relaxation and energy dissipation toward stabilized small polaron through electron-phonon and subsequent phonon-phonon coupling. This study's findings provide conclusive evidence of carrier self-trapping arising from intrinsic lattice anharmonicity and polaronic effect in antimony chalcogenides and a new understanding on the coupled electronic and structural dynamics for redefining excited state properties in soft semiconductor materials.

Growth, Raman Scattering Investigation and Photodetector Properties of 2D SnP.

As an important metal phosphides material, 2D tin phosphides (SnPx 0 <  x ≤ 3) have been theoretically predicted to have intriguing physicochemical properties and potential applications in electronics, optoelectronics, and energy fields. However, the synthesis of high-quality 2D SnP single crystal has not been reported due to the lack of efficiency and reliable growth method. Here, a facile atmospheric pressure chemical vapor deposition (APCVD) method is developed to realize the growth of high-quality 2D SnP nanosheets, by employing tin (Sn) foil as both liquid metal substrates and reaction precursor. Temperature-dependent and angle-resolved polarization Raman spectra observed Raman peaks located at 142.6, 303.3, and 444.2 cm-1 are concluded to belong to A1g mode, which are consistent with the theoretical calculation results. Moreover, the field-effect transistor (FET) devices based on SnP nanosheets show a typical n-type characteristic with an on/off ratio of 103 at 200 K. SnP nanosheets also demonstrate excellent photoresponse performance under the illumination of 473, 532, and 639 nm lasers, which can be tuned by Vgs , Vds , and light power density. It is believed that these findings can provide the first-hand experimental information for the future study of 2D SnP nanosheets.

Controlled synthesis and Raman study of a 2D antiferromagnetic P-type semiconductor: α-MnSe.

Two-dimensional (2D) non-van der Waals magnetic materials have attracted considerable attention due to their high-temperature ferromagnetism, active surface/interface properties originating from dangling bonds, and good stability under ambient conditions. Here, we demonstrate the controlled synthesis and systematic Raman investigation of ultrathin non-van der Waals antiferromagnetic α-MnSe single crystals. Square and triangular nanosheets with different growth orientations can be achieved by introducing different precursors via the atmospheric chemical vapor deposition (APCVD) method. The temperature-dependent resonant enhancement in the Raman intensity of two peaks at 233.8 cm-1 and 459.9 cm-1 gives obvious evidence that the antiferromagnetic spin-ordering is below TN∼ 160 K. Besides, a new peak located at 254.2 cm-1, gradually appearing as the temperature decreased from 180 K to 100 K, may also be a signature of phase transition from paramagnetic to antiferromagnetic. The phonon dispersion spectra of α-MnSe simulated by density functional perturbation theory (DFPT) match well with the observed Raman signals. Moreover, a fabricated α-MnSe phototransistor exhibits p-type conducting behavior and high photodetection performance. We believe that these findings will be beneficial for the applications of 2D α-MnSe in magnetic and semiconducting fields.

Capturing the Electron-Phonon Renormalization in Molecules from First-Principles.

Since the interaction between electrons and atomic nuclei can affect the electronic structure, in recent years, first-principles-based electron-phonon renormalization methods have been applied in the condensed matter physics community to account for the influence of the electron-phonon coupling in solid systems. However, little is yet known about the behavior and trends of the electron-phonon renormalization in the molecules. In this work, the method for the electron-phonon renormalization in molecules has been derived, using which, we exhaustively investigate the zero-point renormalization in 32 molecules with three different density functions. We find that the renormalization of the highest occupied molecular orbital-lowest unoccupied molecular orbital gap due to electron-vibration coupling does not relate to the atomic masses but quite relates to the electronic structure properties of the molecules.

Implementation of Laplace Transformed MP2 for Periodic Systems With Numerical Atomic Orbitals.

We present an implementation of the canonical and Laplace-transformed formulation of the second-order Møller-Plesset perturbation theory under periodic boundary conditions using numerical atomic orbitals. To validate our approach, we show that our results of the Laplace-transformed MP2 correlation correction for the total energy and the band gap are in excellent agreement with the results of the canonical MP2 formulation. We have calculated the binding energy curve for the stacked trans-polyacetylene at the Hartree-Fock + MP2 level as a preliminary application.

The Moving-Grid Effect in the Harmonic Vibrational Frequency Calculations with Numeric Atom-Centered Orbitals.

When using atom-centered integration grids, the portion of the grid that belongs to a certain atom also moves when this atom is displaced. In the paper, we investigate the moving-grid effect in the calculation of the harmonic vibrational frequencies when using all-electron full-potential numeric atomic-centered orbitals as the basis set. We find that, unlike the first-order derivative (i.e., forces), the moving-grid effect plays an essential role for the second-order derivatives (i.e., vibrational frequencies). Further analysis reveals that the predominantly diagonal force constant terms are affected, which can be bypassed efficiently by invoking translational symmetry. Our approaches have been demonstrated in both finite (molecules) and extended (periodic) systems.

Evidence for photogenerated intermediate hole polarons in ZnO.

Despite their pronounced importance for oxide-based photochemistry, optoelectronics and photovoltaics, only fairly little is known about the polaron lifetimes and binding energies. Polarons represent a crucial intermediate step populated immediately after dissociation of the excitons formed in the primary photoabsorption process. Here we present a novel approach to studying photoexcited polarons in an important photoactive oxide, ZnO, using infrared (IR) reflection-absorption spectroscopy (IRRAS) with a time resolution of 100 ms. For well-defined (10-10) oriented ZnO single-crystal substrates, we observe intense IR absorption bands at around 200 meV exhibiting a pronounced temperature dependence. On the basis of first-principles-based electronic structure calculations, we assign these features to hole polarons of intermediate coupling strength.

Implementation of screened hybrid density functional for periodic systems with numerical atomic orbitals: basis function fitting and integral screening.

We present an efficient O(N) implementation of screened hybrid density functional for periodic systems with numerical atomic orbitals (NAOs). NAOs of valence electrons are fitted with gaussian-type orbitals, which is convenient for the calculation of electron repulsion integrals and the construction of Hartree-Fock exchange matrix elements. All other parts of Hamiltonian matrix elements are constructed directly with NAOs. The strict locality of NAOs is adopted as an efficient two-electron integral screening technique to speed up calculations.