Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs.

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

The impact of intravenous opioid shortages on postoperative pain control in a pediatric cardiac intensive care unit.

BACKGROUND: The national supply of intravenous pain medications has been restricted since 2018 due to shortages. In late 2017, the ongoing shortages combined with the increased demand for alternative formulations resulted in interruptions in the supply of medications. The cascading effect of interruptions in supply caused hospitals to implement restrictive measures to extend the availability of these medications. This study aims to evaluate differences in pain control in postoperative cardiac patients due to intravenous opioid shortage at a standalone pediatric hospital. METHOD: A retrospective chart review of postoperative patients in the cardiac intensive care unit at Children's National Hospital was conducted from January to September 2017 and 2018. The objective of this study is to determine if a difference in postoperative pain scores in the period before and during the national intravenous opioid shortage exists, analyzed by differences in Face, Legs, Activity, Cry, Consolability (FLACC) scores. RESULTS: No differences were found in mean FLACC scores for the 140 patients who were evaluated. Patients in 2017 reported faster decline of pain symptoms compared with those in 2018. No differences in baseline characteristics or secondary outcomes were noted. CONCLUSION: Pain management was better controlled immediately postoperatively in 2018 compared with 2017, but the pain was better controlled for the duration of the postoperative period in 2017.

Hybrid DFT investigation of the energetics of Mg ion diffusion in α-MoO3.

Rechargeable batteries that utilize divalent Mg ions as the charge carrier species can in principle achieve substantially greater volumetric energy densities than conventional Li-ion batteries. One significant impediment to the development of commercially viable Mg-ion batteries is the slow rate of Mg ion diffusion through otherwise promising cathode materials. Accurate prediction of the activation energies associated with this diffusion process using density functional theory (DFT) is especially challenging due to self-interaction errors intrinsic to DFT that lead to over-delocalization of the d-electrons. One effective but highly computationally demanding approach to reducing self-interaction errors is the use of hybrid functionals, which incorporate a fraction of exact Hartree-Fock exchange. In this work, we assess the effects of exact exchange on computed activation energies for ion diffusion in one potential cathode material, α-MoO3. In contrast to previous studies that primarily utilize non-hybrid functionals, we perform nudged elastic band calculations in which the nuclear coordinates are fully converged using both hybrid functionals and k-point sampling. It is found that while non-hybrid functionals indicate the existence of thermodynamically accessible channels for bulk Mg ion diffusion in all three dimensions, hybrid functionals predict that some of these channels are largely inaccessible under typical charge/discharge conditions. Furthermore, it is demonstrated that certain commonly used approximations for incorporating the effects of Hartree-Fock exchange are inadequate for this system, including DFT+U calculations and the use of single-point hybrid calculations using atomic positions obtained using non-hybrid functionals.

Accurate and systematically improvable density functional theory embedding for correlated wavefunctions.

We analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using wavefunction methods, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We test an MP2 correction for this term and demonstrate that the corrected embedding scheme accurately reproduces wavefunction calculations for a series of chemical reactions. Our projector-based embedding method uses localized occupied orbitals to partition the system; as with other local correlation methods, abrupt changes in the character of the localized orbitals along a reaction coordinate can lead to discontinuities in the embedded energy, but we show that these discontinuities are small and can be systematically reduced by increasing the size of the active region. Convergence of reaction energies with respect to the size of the active subsystem is shown to be rapid for all cases where the density functional treatment is able to capture the polarization of the environment, even in conjugated systems, and even when the partition cuts across a double bond.

Accurate basis set truncation for wavefunction embedding.

Density functional theory (DFT) provides a formally exact framework for performing embedded subsystem electronic structure calculations, including DFT-in-DFT and wavefunction theory-in-DFT descriptions. In the interest of efficiency, it is desirable to truncate the atomic orbital basis set in which the subsystem calculation is performed, thus avoiding high-order scaling with respect to the size of the MO virtual space. In this study, we extend a recently introduced projection-based embedding method [F. R. Manby, M. Stella, J. D. Goodpaster, and T. F. Miller III, J. Chem. Theory Comput. 8, 2564 (2012)] to allow for the systematic and accurate truncation of the embedded subsystem basis set. The approach is applied to both covalently and non-covalently bound test cases, including water clusters and polypeptide chains, and it is demonstrated that errors associated with basis set truncation are controllable to well within chemical accuracy. Furthermore, we show that this approach allows for switching between accurate projection-based embedding and DFT embedding with approximate kinetic energy (KE) functionals; in this sense, the approach provides a means of systematically improving upon the use of approximate KE functionals in DFT embedding.

Density functional theory embedding for correlated wavefunctions: improved methods for open-shell systems and transition metal complexes.

Density functional theory (DFT) embedding provides a formally exact framework for interfacing correlated wave-function theory (WFT) methods with lower-level descriptions of electronic structure. Here, we report techniques to improve the accuracy and stability of WFT-in-DFT embedding calculations. In particular, we develop spin-dependent embedding potentials in both restricted and unrestricted orbital formulations to enable WFT-in-DFT embedding for open-shell systems, and develop an orbital-occupation-freezing technique to improve the convergence of optimized effective potential calculations that arise in the evaluation of the embedding potential. The new techniques are demonstrated in applications to the van-der-Waals-bound ethylene-propylene dimer and to the hexa-aquairon(II) transition-metal cation. Calculation of the dissociation curve for the ethylene-propylene dimer reveals that WFT-in-DFT embedding reproduces full CCSD(T) energies to within 0.1 kcal/mol at all distances, eliminating errors in the dispersion interactions due to conventional exchange-correlation (XC) functionals while simultaneously avoiding errors due to subsystem partitioning across covalent bonds. Application of WFT-in-DFT embedding to the calculation of the low-spin/high-spin splitting energy in the hexaaquairon(II) cation reveals that the majority of the dependence on the DFT XC functional can be eliminated by treating only the single transition-metal atom at the WFT level; furthermore, these calculations demonstrate the substantial effects of open-shell contributions to the embedding potential, and they suggest that restricted open-shell WFT-in-DFT embedding provides better accuracy than unrestricted open-shell WFT-in-DFT embedding due to the removal of spin contamination.

Embedded density functional theory for covalently bonded and strongly interacting subsystems.

Embedded density functional theory (e-DFT) is used to describe the electronic structure of strongly interacting molecular subsystems. We present a general implementation of the Exact Embedding (EE) method [J. Chem. Phys. 133, 084103 (2010)] to calculate the large contributions of the nonadditive kinetic potential (NAKP) in such applications. Potential energy curves are computed for the dissociation of Li(+)-Be, CH(3)-CF(3), and hydrogen-bonded water clusters, and e-DFT results obtained using the EE method are compared with those obtained using approximate kinetic energy functionals. In all cases, the EE method preserves excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures. We also demonstrate an accurate pairwise approximation to the NAKP that allows for efficient parallelization of the EE method in large systems; benchmark calculations on molecular crystals reveal ideal, size-independent scaling of wall-clock time with increasing system size.