Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient heat-bath sampling in Fock space
130
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce an algorithm for sampling many-body quantum states in Fock space. The algorithm efficiently samples states with probability approximately proportional to an arbitrary function of the second-quantized Hamiltonian matrix element connecting the sampled state to the current state. We apply the new sampling algorithm to the recently-developed Semistochastic Full Configuration Interaction Quantum Monte Carlo method (S-FCIQMC), a semistochastic implementation of the power method for projecting out the ground state energy in a basis of Slater determinants. The heat-bath sampling requires modest additional computational time and memory compared to uniform sampling but results in newly-spawned weights that are approximately of the same magnitude, thereby greatly improving the efficiency of projection. A comparison in efficiency between uniform and approximate heat-bath sampling is performed on the all-electron nitrogen dimer at equilibrium in Dunnings cc-pVXZ basis sets with X in {D, T, Q, 5}, demonstrating a large gain in efficiency that increases with basis set size. In addition, a comparison in efficiency is performed on three all-electron first-row dimers, B_2, N_2, and F_2, in a cc-pVQZ basis, demonstrating that the gain in efficiency compared to uniform sampling also increases dramatically with the number of electrons.
rate research
Read More
We extend our recently-developed heat-bath configuration interaction (HCI) algorithm, and our semistochastic algorithm for performing multireference perturbation theory, to the calculation of excited-state wavefunctions and energies. We employ time-reversal symmetry, which reduces the memory requirements by more than a factor of two. An extrapolation technique is introduced to reliably extrapolate HCI energies to the Full CI limit. The resulting algorithm is used to compute the twelve lowest-lying potential energy surfaces of the carbon dimer using the cc-pV5Z basis set, with an estimated error in energy of 30-50 {mu}Ha compared to Full CI. The excitation energies obtained using our algorithm have a mean absolute deviation of 0.02 eV compared to experimental values. We also calculate the complete active-space (CAS) energies of the S0, S1, and T0 states of tetracene, which are of relevance to singlet fission, by fully correlating active spaces as large as 18 electrons in 36 orbitals.
We introduce a new selected configuration interaction plus perturbation theory algorithm that is based on a deterministic analog of our recent efficient heat-bath sampling algorithm. This Heat-bath Configuration Interaction (HCI) algorithm makes use of two parameters that control the tradeoff between speed and accuracy, one which controls the selection of determinants to add to a variational wavefunction, and one which controls the the selection of determinants used to compute the perturbative correction to the variational energy. We show that HCI provides an accurate treatment of both static and dynamic correlation by computing the potential energy curve of the multireference carbon dimer in the cc-pVDZ basis. We then demonstrate the speed and accuracy of HCI by recovering the full configuration interaction energy of both the carbon dimer in the cc-pVTZ basis and the strongly-correlated chromium dimer in the Ahlrichs VDZ basis, correlating all electrons, to an accuracy of better than 1 mHa, in just a few minutes on a single core. These systems have full variational spaces of 3x10^14 and 2x10^22 determinants respectively.
In molecular devices electronic degrees of freedom are coupled to vibrational modes of the molecule, offering an opportunity to study fundamental aspects of this coupling between at the nanoscale. To this end we consider the nonequilibrium heat exchange between a conduction band and a bosonic bath mediated by a single molecule. For molecules large enough so that on-site interactions can be dropped we carry out an asymptotically exact calculation of the heat current, governed by the smallness of the electron-phonon coupling, and obtain the steady state heat current driven by a finite temperature drop. At low temperatures the heat current is found to have a power-law behavior with respect to the temperature difference with the power depending on the nature of the bosonic bath. At high temperatures, on the other hand, the current is linear in the temperature difference for all types of bosonic baths. The crossover between these behaviors is described. Some of the results are given a physical explanation by comparing to a perturbative Master equation calculation (whose limitation we examine).
By using Poissons summation formula, we calculate periodic integrals over Gaussian basis functions by partitioning the lattice summations between the real and reciprocal space, where both sums converge exponentially fast with a large exponent. We demonstrate that the summation can be performed efficiently to calculate 2-center Gaussian integrals over various kernels including overlap, kinetic, and Coulomb. The summation in real space is performed using an efficient flavor of the McMurchie-Davidson Recurrence Relation (MDRR). The expressions for performing summation in the reciprocal space are also derived and implemented. The algorithm for reciprocal space summation allows us to reuse several terms and leads to significant improvement in efficiency when highly contracted basis functions with large exponents are used. We find that the resulting algorithm is only between a factor of 5 to 15 slower than that for molecular integrals, indicating the very small number of terms needed in both the real and reciprocal space summations. An outline of the algorithm for calculating 3-center Coulomb integrals is also provided.
We introduce vibrational heat-bath configuration interaction (VHCI) as an accurate and efficient method for calculating vibrational eigenstates of anharmonic systems. Inspired by its origin in electronic structure theory, VHCI is a selected CI approach that uses a simple criterion to identify important basis states with a pre-sorted list of anharmonic force constants. Screened second-order perturbation theory and simple extrapolation techniques provide significant improvements to variational energy estimates. We benchmark VHCI on four molecules with 12 to 48 degrees of freedom and use anharmonic potential energy surfaces truncated at fourth and sixth order. For all molecules studied, VHCI produces vibrational spectra of tens or hundreds of states with sub-wavenumber accuracy at low computational cost.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
