RESUMEN
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.
RESUMEN
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.
RESUMEN
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.
RESUMEN
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.