Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. ...The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.
We propose the use of preconditioning in FCIQMC which, in combination with perturbative estimators, greatly increases the efficiency of the algorithm. The use of preconditioning allows a time step ...close to unity to be used (without time-step errors), provided that multiple spawning attempts are made per walker. We show that this approach substantially reduces statistical noise on perturbative corrections to initiator error, which improve the accuracy of FCIQMC but which can suffer from significant noise in the original scheme. Therefore, the use of preconditioning and perturbatively corrected estimators in combination leads to a significantly more efficient algorithm. In addition, a simpler approach to sampling variational and perturbative estimators in FCIQMC is presented, which also allows the variance of the energy to be calculated. These developments are investigated and applied to benzene (30e, 108o), an example where accurate treatment is not possible with the original method.
High-quality excitation generators are crucial to the effectiveness of coupled cluster Monte Carlo (CCMC) and full configuration interaction Quantum Monte Carlo (FCIQMC) calculations. The heat bath ...sampling of Holmes et al. Holmes, A. A.; Changlani, H. J.; Umrigar, C. J. J. Chem. Theory Comput. 2016, 12, 1561–1571. dramatically increases the efficiency of the spawn step of such algorithms but requires memory storage scaling quartically with system size which can be prohibitive for large systems. Alternatively, Alavi et al. Smart, S. D.; Booth, G. H.; Alavi, A. Unpublished results. approximated these weights with weights based on Cauchy–Schwarz-like inequalities calculated on-the-fly. While reducing the memory cost, this algorithm scales linearly in system size computationally. We combine both of these ideas with the single-reference nature of many systems studied and introduce a spawn-sampling algorithm that has low memory requirements (quadratic in basis set size) compared to the heat bath algorithm and only scales either independently of system size (CCMC) or linearly in the number of electrons (FCIQMC) that works especially well on localized orbitals. Tests on small water chains with localized orbitals with CCMC and with an initiator point sample in FCIQMC indicate that it can be equally efficient as the other excitation generators. As the system gets larger, calculations with our new algorithm converge faster than the on-the-fly weight algorithm while having a much more favorable memory scaling than the heat bath algorithm.
The Canadian Study of Health and Aging Clinical Frailty Scale (CFS) and the laboratory Frailty Index (FI-lab) are validated tools based on clinical and laboratory data, respectively. Their utility as ...predictors of geriatric trauma outcomes is unknown. Our primary objective was to determine whether pre-admission CFS is associated with adverse discharge destination. Secondary objectives were to evaluate the relationships between CFS and in-hospital complications and between admission FI-lab and discharge destination.
We performed a 4-year (2011 to 2014) retrospective cohort study with patients 65 years and older admitted to a level I trauma center. Admission FI-lab was calculated using 23 variables collected within 48 hours of presentation. The primary outcome was discharge destination, either adverse (death or discharge to a long-term, chronic, or acute care facility) or favorable (home or rehabilitation). The secondary outcome was in-hospital complications. Multivariable logistic regression was used to evaluate the relationship between CFS or FI-lab and outcomes.
There were 266 patients included. Mean age was 76.5 ± 7.8 years and median Injury Severity Score was 17 (interquartile range 13 to 24). There were 260 patients and 221 patients who had sufficient data to determine CFS and FI-lab scores, respectively. Pre-admission frailty as per the CFS (CFS 6 or 7) was independently associated with adverse discharge destination (odds ratio 5.1; 95% CI 2.0 to 13.2; p < 0.001). Severe frailty on admission, as determined by the FI-lab (FI-lab > 0.4), was not associated with adverse outcomes.
Pre-admission clinical frailty independently predicts adverse discharge destination in geriatric trauma patients. The CFS may be used to triage resources to mitigate adverse outcomes in this population. The FI-lab determined on admission for trauma may not be useful.
The convergence of full configuration interaction quantum Monte Carlo (FCIQMC) is accelerated using a quasi-Newton propagation (QN) which can also be applied to coupled cluster Monte Carlo (CCMC). ...The computational scaling of this optimized propagation is O ( 1 ) , keeping the additional computational cost to a bare minimum. Its effects are investigated deterministically and stochastically on a model system, the uniform electron gas, with Hilbert space size up to 1040, and shown to accelerate convergence of the instantaneous projected energy by over an order of magnitude in the FCIQMC test case. Its capabilities are then demonstrated with FCIQMC on an archetypical quantum chemistry problem, the chromium dimer, in an all-electron basis set with Hilbert space size of about 1022 yielding highly accurate FCI energies.