ترغب بنشر مسار تعليمي؟ اضغط هنا

All-electron quantum Monte Carlo with Jastrow single determinant Ansatz: application to the sodium dimer

29   0   0.0 ( 0 )
 نشر من قبل Kousuke Nakano
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

In this work, we report potential energy surfaces (PESs) of the sodium dimer calculated by variational (VMC) and lattice regularized diffusion Monte Carlo (LRDMC). The VMC calculation is accurate for determining the equilibrium distance and the qualitative shape of the experimental PES. Remarkably, after the application of the LRDMC projection to this single determinant ansatz, namely the Jastrow Antisymmetrized Geminal Product (JAGP), chemical accuracy (~ 1kcal/mol) is reached, and the obtained dissociation energy, equilibrium internuclear distance, and harmonic vibrational frequency are in very good agreement with the experimental ones. This outcome crucially depends on the quality of the optimization used to determine the best possible trial function within the chosen ansatz. The strategy adopted in this work is to minimize the variational energy by initializing the trial function with the DFT single determinant ansatz expanded exactly in the same atomic basis used for the corresponding VMC and LRDMC calculations. This atomic basis is ad-hoc reshaped for QMC calculations. Indeed, we multiply the standard Gaussian type atomic orbitals by a one-body Jastrow factor, satisfying in this way the electron-ion cusp conditions. This allows us to use a very small basis almost converged in the complete basis set limit, by reducing the computational effort as well as the statistical fluctuations on the total energy. In order to achieve these important advantages, we have defined a very efficient DFT algorithm in the mentioned basis, by estimating the corresponding matrix elements on a mesh, and by using a much finer mesh grid in the vicinity of nuclei.

قيم البحث

اقرأ أيضاً

66 - Sam Azadi , R. E. Cohen 2015
We report an accurate study of interactions between Benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory (DFT) using different van der W aals (vdW) functionals. In our QMC calculations, we use accurate correlated trial wave functions including three-body Jastrow factors, and backflow transformations. We consider two benzene molecules in the parallel displaced (PD) geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the CCSD(T)/CBS limit is -2.65(2) kcal/mol [E. Miliordos et al, J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
We report results of both Diffusion Quantum Monte Carlo(DMC) method and Reptation Quantum Monte Carlo(RMC) method on the potential energy curve of the helium dimer. We show that it is possible to obtain a highly accurate description of the helium dim er. An improved stochastic reconfiguration technique is employed to optimize the many-body wave function, which is the starting point for highly accurate simulations based on the Diffusion Quantum Monte Carlo(DMC) and Reptation Quantum Monte Carlo (RMC) methods. We find that the results of these methods are in excellent agreement with the best theoretical results at short range, especially recently developed Reptation Quantum Monte Carlo(RMC) method, yield practically accurate results with reduced statistical error, which gives very excellent agreement across the whole potential. For the equilibrium internuclear distance of 5.6 bohr, the calculated electronic energy with Reptation Quantum Monte Carlo(RMC) method is 5.807483599$pm$0.000000015 hartrees and the corresponding well depth is -11.003$pm$0.005 K.
67 - Carsten Bauer 2020
We assess numerical stabilization methods employed in fermion many-body quantum Monte Carlo simulations. In particular, we empirically compare various matrix decomposition and inversion schemes to gain control over numerical instabilities arising in the computation of equal-time and time-displaced Greens functions within the determinant quantum Monte Carlo (DQMC) framework. Based on this comparison, we identify a procedure based on pivoted QR decompositions which is both efficient and accurate to machine precision. The Julia programming language is used for the assessment and implementations of all discussed algorithms are provided in the open-source software library StableDQMC.jl [http://github.com/crstnbr/StableDQMC.jl].
We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. These quantities are derived from the changes in electr onic polarization as a function of applied finite electric field - an approach we recently introduced and made possible by the use of a Berry-phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities are found to be in excellent agreement with the best estimates available from quantum chemistry - usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is found that while exchange effects dominate the proper description of the susceptibilities, second hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different approximations to the nodal surface of the many-body wavefunction affect the accuracy of the calculated susceptibilities.
122 - Y. F. Kung , C.-C. Chen , Yao Wang 2016
We characterize the three-orbital Hubbard model using state-of-the-art determinant quantum Monte Carlo (DQMC) simulations with parameters relevant to the cuprate high-temperature superconductors. The simulations find that doped holes preferentially r eside on oxygen orbitals and that the ({pi},{pi}) antiferromagnetic ordering vector dominates in the vicinity of the undoped system, as known from experiments. The orbitally-resolved spectral functions agree well with photoemission spectroscopy studies and enable identification of orbital content in the bands. A comparison of DQMC results with exact diagonalization and cluster perturbation theory studies elucidates how these different numerical techniques complement one another to produce a more complete understanding of the model and the cuprates. Interestingly, our DQMC simulations predict a charge-transfer gap that is significantly smaller than the direct (optical) gap measured in experiment. Most likely, it corresponds to the indirect gap that has recently been suggested to be on the order of 0.8 eV, and demonstrates the subtlety in identifying charge gaps.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا