Massively parallel GPU enabled third-order cluster perturbation excitation energies for cost-effective large scale excitation energy calculations.

We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

Multilayer Divide-Expand-Consolidate Coupled-Cluster Method: Demonstrative Calculations of the Adsorption Energy of Carbon Dioxide in the Mg-MOF-74 Metal-Organic Framework.

The implementation and evaluation of a multilayer extension of the divide-expand-consolidate (DEC) scheme within the LSDalton program is presented. The DEC scheme is a linear-scaling, fragmentation-based local coupled-cluster (CC) method that provides a means of overcoming the scaling wall associated with canonical CC electronic structure calculations on large molecular systems. Taking advantage of the local nature of correlation effects, the correlation energy for the full molecule is calculated from a set of independent fragments using localized molecular orbitals. However, when only a small subsystem of a larger system is of interest, for example, adsorption sites or catalytically active sites, the majority of the computational time may be spent evaluating the correlation energy of fragments which have little effect on the properties in the area of interest (AOI). The multilayer DEC (ML-DEC) scheme addresses this by taking advantage of the independent nature of the fragments in order to evaluate the correlation energy of various regions of the system at different levels of theory. Regions far from the AOI are evaluated at lower (cheaper) levels of theory such as Hartree-Fock (HF) or Møller-Plesset second-order perturbation theory (MP2), while the area immediately surrounding the AOI is treated with a higher level CC model. Through the ML-DEC scheme, the computational cost of CC calculations on these types of systems can be significantly reduced while maintaining the accuracy of higher-level calculations. Results from HF/RI-MP2 and RI-MP2/CCSD ML-DEC calculations of the binding energy of a fatty acid dimer are presented. We find that the ML-DEC scheme is capable of reproducing DEC energy differences at a target level of theory, provided that the region treated at the target level of theory is chosen to be sufficiently large. Time-to-solution is found to be significantly reduced, particularly in the RI-MP2/CCSD calculations. Finally, the ML-DEC scheme is applied to the calculation of CO2 adsorption in a Mg-MOF-74 channel.

Three-body interactions and the elastic constants of hcp solid 4He.

The effect of three-body interactions on the elastic properties of hexagonal close packed solid 4He is investigated using variational path integral (VPI) Monte Carlo simulations. The solid's nonzero elastic constants are calculated, at T = 0 K and for a range of molar volumes from 7.88 cm3/mol to 20.78 cm3/mol, from the bulk modulus and the three pure shear constants C0, C66, and C44. Three-body interactions are accounted for using our recently reported perturbative treatment based on the nonadditive three-body potential of Cencek et al. Previous studies have attempted to account for the effect of three-body interactions on the elastic properties of solid 4He; however, these calculations have treated zero point motions using either the Einstein or Debye approximations, which are insufficient in the molar volume range where solid 4He is characterized as a quantum solid. Our VPI calculations allow for a more accurate treatment of the zero point motions which include atomic correlation. From these calculations, we find that agreement with the experimental bulk modulus is significantly improved when three-body interactions are considered. In addition, three-body interactions result in non-negligible differences in the calculated pure shear constants and nonzero elastic constants, particularly at higher densities, where differences of up to 26.5% are observed when three-body interactions are included. We compare to the available experimental data and find that our results are generally in as good or better agreement with experiment as previous theoretical investigations.

Search for anisotropy in the Debye-Waller factor of HCP solid (4)He.

The properties of hexagonal close packed (hcp) solid (4)He are dominated by large atomic zero point motions. An accurate description of these motions is therefore necessary in order to accurately calculate the properties of the system, such as the Debye-Waller (DW) factors. A recent neutron scattering experiment reported significant anisotropy in the in-plane and out-of-plane DW factors for hcp solid (4)He at low temperatures, where thermal effects are negligible and only zero-point motions are expected to contribute. By contrast, no such anisotropy was observed either in earlier experiments or in path integral Monte Carlo (PIMC) simulations of solid hcp (4)He. However, the earlier experiments and the PIMC simulations were both carried out at higher temperatures where thermal effects could be substantial. We seek to understand the cause of this discrepancy through variational quantum Monte Carlo simulations utilizing an accurate pair potential and a modified trial wavefunction which allows for anisotropy. Near the melting density, we find no anisotropy in an ideal hcp (4)He crystal. A theoretical equation of state is derived from the calculated energies of the ideal crystal over a range of molar volumes from 7.88 to 21.3 cm(3), and is found to be in good qualitative agreement with experimental data.

Synthesis and biochemical evaluation of benzoylbenzophenone thiosemicarbazone analogues as potent and selective inhibitors of cathepsin L.

Upregulation of cathepsin L in a variety of tumors and its ability to promote cancer cell invasion and migration through degradation of the extracellular matrix suggest that cathepsin L is a promising biological target for the development of anti-metastatic agents. Based on encouraging results from studies on benzophenone thiosemicarbazone cathepsin inhibitors, a series of fourteen benzoylbenzophenone thiosemicarbazone analogues were designed, synthesized, and evaluated for their inhibitory activity against cathepsins L and B. Thiosemicarbazone inhibitors 3-benzoylbenzophenone thiosemicarbazone 1, 1,3-bis(4-fluorobenzoyl)benzene thiosemicarbazone 8, and 1,3-bis(2-fluorobenzoyl)-5-bromobenzene thiosemicarbazone 32 displayed the greatest potency against cathepsin L with low IC50 values of 9.9 nM, 14.4 nM, and 8.1 nM, respectively. The benzoylbenzophenone thiosemicarbazone analogues evaluated were selective in their inhibition of cathepsin L compared to cathepsin B. Thiosemicarbazone analogue 32 inhibited invasion through Matrigel of MDA-MB-231 breast cancer cells by 70% at 10 µM. Thiosemicarbazone analogue 8 significantly inhibited the invasive potential of PC-3ML prostate cancer cells by 92% at 5 µM. The most active cathepsin L inhibitors from this benzoylbenzophenone thiosemicarbazone series (1, 8, and 32) displayed low cytotoxicity toward normal primary cells [in this case human umbilical vein endothelial cells (HUVECs)]. In an initial in vivo study, 3-benzoylbenzophenone thiosemicarbazone (1) was well-tolerated in a CDF1 mouse model bearing an implanted C3H mammary carcinoma, and showed efficacy in tumor growth delay. Low cytotoxicity, inhibition of cell invasion, and in vivo tolerability are desirable characteristics for anti-metastatic agents functioning through an inhibition of cathepsin L. Active members of this structurally diverse group of benzoylbenzophenone thiosemicarbazone cathepsin L inhibitors show promise as potential anti-metastatic, pre-clinical drug candidates.

Testing novel strategies for patients hospitalised with HIV-associated disseminated tuberculosis (NewStrat-TB): protocol for a randomised controlled trial.

BACKGROUND: HIV-associated tuberculosis (TB) contributes disproportionately to global tuberculosis mortality. Patients hospitalised at the time of the diagnosis of HIV-associated disseminated TB are typically severely ill and have a high mortality risk despite initiation of tuberculosis treatment. The objective of the study is to assess the safety and efficacy of both intensified TB treatment (high dose rifampicin plus levofloxacin) and immunomodulation with corticosteroids as interventions to reduce early mortality in hospitalised patients with HIV-associated disseminated TB. METHODS: This is a phase III randomised controlled superiority trial, evaluating two interventions in a 2 × 2 factorial design: (1) high dose rifampicin (35 mg/kg/day) plus levofloxacin added to standard TB treatment for the first 14 days versus standard tuberculosis treatment and (2) adjunctive corticosteroids (prednisone 1.5 mg/kg/day) versus identical placebo for the first 14 days of TB treatment. The study population is HIV-positive patients diagnosed with disseminated TB (defined as being positive by at least one of the following assays: urine Alere LAM, urine Xpert MTB/RIF Ultra or blood Xpert MTB/RIF Ultra) during a hospital admission. The primary endpoint is all-cause mortality at 12 weeks comparing, first, patients receiving intensified TB treatment to standard of care and, second, patients receiving corticosteroids to those receiving placebo. Analysis of the primary endpoint will be by intention to treat. Secondary endpoints include all-cause mortality at 2 and 24 weeks. Safety and tolerability endpoints include hepatoxicity evaluations and corticosteroid-related adverse events. DISCUSSION: Disseminated TB is characterised by a high mycobacterial load and patients are often critically ill at presentation, with features of sepsis, which carries a high mortality risk. Interventions that reduce this high mycobacterial load or modulate associated immune activation could potentially reduce mortality. If found to be safe and effective, the interventions being evaluated in this trial could be easily implemented in clinical practice. TRIAL REGISTRATION: ClinicalTrials.gov NCT04951986. Registered on 7 July 2021 https://clinicaltrials.gov/study/NCT04951986.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers.

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Coupled cluster theory on modern heterogeneous supercomputers.

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.