Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPower laws used to extrapolate the coupled cluster correlation energy to the thermodynamic limit
82
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Recent calculations using coupled cluster on solids have raised discussion of using a $N^{-1/3}$ power law to fit the correlation energy when extrapolating to the thermodynamic limit, an approach which differs from the more commonly used $N^{-1}$ power law which is (for example) often used by quantum Monte Carlo methods. In this paper, we present one way to reconcile these viewpoints. Coupled cluster doubles calculations were performed on uniform electron gases reaching system sizes of $922$ electrons for an extremely wide range of densities ($0.1<r_s<100.0$) to study how the correlation energy approaches the thermodynamic limit. The data were corrected for basis set incompleteness error and use a selected twist angle approach to mitigate finite size error from shell filling effects. Analyzing these data, we initially find that a power law of $N^{-1/3}$ appears to fit the data better than a $N^{-1}$ power law in the large system size limit. However, we provide an analysis of the transition structure factor showing that $N^{-1}$ still applies to large system sizes and that the apparent $N^{-1/3}$ power law occurs only at low $N$.
rate research
Read More
Modern electronic structure theories can predict and simulate a wealth of phenomena in surface science and solid-state physics. In order to allow for a direct comparison with experiment, such ab initio predictions have to be made in the thermodynamic limit, substantially increasing the computational cost of many-electron wave-function theories. Here, we present a method that achieves thermodynamic limit results for solids and surfaces using the gold standard coupled cluster ansatz of quantum chemistry with unprecedented efficiency. We study the energy difference between carbon diamond and graphite crystals, adsorption energies of water on h-BN, as well as the cohesive energy of the Ne solid, demonstrating the increased efficiency and accuracy of coupled cluster theory for solids and surfaces.
The Shockley-Queisser (SQ) limit provides a convenient metric for predicting light-to-electricity conversion efficiency of a solar cell based on the band gap of the light-absorbing layer. In reality, few materials approach this radiative limit. We develop a formalism and a computational method to predict the maximum photovoltaic efficiency of imperfect crystals from first principles. Our scheme includes equilibrium populations of native defects, their carrier-capture coefficients, and the associated recombination rates. When applied to kesterite solar cells, we reveal an intrinsic limit of 20% for $mathrm{Cu_2ZnSnSe_4}$, which falls far below the SQ limit of 32%. The effects of atomic substitution and extrinsic doping are studied, leading to pathways for enhanced efficiency of 31%. This approach can be applied to support targeted-materials selection for future solar-energy technologies.
Machine learning models are changing the paradigm of molecular modeling, which is a fundamental tool for material science, chemistry, and computational biology. Of particular interest is the inter-atomic potential energy surface (PES). Here we develop Deep Potential - Smooth Edition (DeepPot-SE), an end-to-end machine learning-based PES model, which is able to efficiently represent the PES for a wide variety of systems with the accuracy of ab initio quantum mechanics models. By construction, DeepPot-SE is extensive and continuously differentiable, scales linearly with system size, and preserves all the natural symmetries of the system. Further, we show that DeepPot-SE describes finite and extended systems including organic molecules, metals, semiconductors, and insulators with high fidelity.
The accuracy of the many-body perturbation theory GW formalism to calculate electron-phonon coupling matrix elements has been recently demonstrated in the case of a few important systems. However, the related computational costs are high and thus represent strong limitations to its widespread application. In the present study, we explore two less demanding alternatives for the calculation of electron-phonon coupling matrix elements on the many-body perturbation theory level. Namely, we test the accuracy of the static Coulomb-hole plus screened-exchange (COHSEX) approximation and further of the constant screening approach, where variations of the screened Coulomb potential W upon small changes of the atomic positions along the vibrational eigenmodes are neglected. We find this latter approximation to be the most reliable, whereas the static COHSEX ansatz leads to substantial errors. Our conclusions are validated in a few paradigmatic cases: diamond, graphene and the C60 fullerene. These findings open the way for combining the present many-body perturbation approach with efficient linear-response theories.
Recent advances in selected CI, including the adaptive sampling configuration interaction (ASCI) algorithm and its heat bath extension, have made the ASCI approach competitive with the most accurate techniques available, and hence an increasingly powerful tool in solving quantum Hamiltonians. In this work, we show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques we present ASCI results applied to a large range of systems and basis sets in order to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-zeta benchmark data for the G1 dataset. The largest of these calculations is Si$_{2}$H$_{6}$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
