shadedancer1

CO2 adsorption and activation on a catalyst are key elementary steps for CO2 conversion to various valuable products. In the present computational study, we screened different Cu-ZrO2 interface structures and analyzed the influence of the interface structure on CO2 binding strength using density functional theory calculations. Our results demonstrate that a Cu nanorod favors one position on both tetragonal and monoclinic ZrO2 surfaces, where the bottom Cu atoms are placed close to the lattice oxygens. In agreement with previous calculations, we find that CO2 prefers a bent bidentate configuration at the Cu-ZrO2 interface and the molecule is clearly activated being negatively charged. Straining of the Cu nanorod influences CO2 adsorption energy but does not change the preferred nanorod position on zirconia. Altogether, our results highlight that CO2 adsorption and activation depend sensitively on the chemical composition and atomic structure of the interface used in the calculations. This structure sensitivity may potentially impact further catalytic steps and the overall computed reactivity profile.In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation.Quantum dephasing of excitonic transitions in CsPbBr3 nanocrystals has been studied using two-dimensional electronic spectroscopy at cryogenic temperatures. The exciton-phonon interactions for acoustic and optical modes exhibit different effects on the coherent dynamics of excitonic transitions. The homogeneous linewidth shows a proportional dependence on the temperature, suggesting the primary dephasing channel of the elastic scattering between exciton and acoustic modes. The exciton-optical mode interaction is manifested as the beatings of off-diagonal signals in the population time domain at the frequencies of 29 and 51 cm-1, indicating phonon replicas of excitonic transitions arising from coherent exciton-phonon interaction. The insight information of exciton homogeneous broadening in perovskite nanocrystals is essential for the potential application of quantum light sources.Diamine-appended metal-organic frameworks (MOFs) of the form Mg2(dobpdc)(diamine)2 adsorb CO2 in a cooperative fashion, exhibiting an abrupt change in CO2 occupancy with pressure or temperature. This change is accompanied by hysteresis. While hysteresis is suggestive of a first-order phase transition, we show that hysteretic temperature-occupancy curves associated with this material are qualitatively unlike the curves seen in the presence of a phase transition; they are instead consistent with CO2 chain polymerization, within one-dimensional channels in the MOF, in the absence of a phase transition. Our simulations of a microscopic model reproduce this dynamics, providing a physical understanding of cooperative adsorption in this industrially important class of materials.Transition metal-catalyzed reactions invariably include steps where ligands associate or dissociate. In order to obtain reliable energies for such reactions, sufficiently large basis sets need to be employed. In this paper, we have used high-precision multiwavelet calculations to compute the metal-ligand association energies for 27 transition metal complexes with common ligands, such as H2, CO, olefins, and solvent molecules. By comparing our multiwavelet results to a variety of frequently used Gaussian-type basis sets, we show that counterpoise corrections, which are widely employed to correct for basis set superposition errors, often lead to underbinding. Additionally, counterpoise corrections are difficult to employ when the association step also involves a chemical transformation. Multiwavelets, which can be conveniently applied to all types of reactions, provide a promising alternative for computing electronic interaction energies free from any basis set errors.Correlation-driven phenomena in molecular periodic systems are challenging to predict computationally not only because such systems are periodically infinite but also because they are typically strongly correlated. Here, we generalize the variational two-electron reduced density matrix (2-RDM) theory to compute the energies and properties of strongly correlated periodic systems. The 2-RDM of the unit cell is directly computed subject to necessary N-representability conditions such that the unit-cell 2-RDM represents at least one N-electron density matrix. Two canonical but non-trivial systems, periodic metallic hydrogen chains and periodic acenes, are treated to demonstrate the methodology. We show that while single-reference correlation theories do not capture the strong (static) correlation effects in either of these molecular systems, the periodic variational 2-RDM theory predicts the Mott metal-to-insulator transition in the hydrogen chains and the length-dependent polyradical formation in acenes. For both hydrogen chains and acenes, the periodic calculations are compared with previous non-periodic calculations with the results showing a significant change in energies and increase in the electron correlation from the periodic boundary conditions. The 2-RDM theory, which allows for much larger active spaces than are traditionally possible, is applicable to studying correlation-driven phenomena in general periodic molecular solids and materials.Computational software workflows are emerging as all-in-one solutions to speed up the discovery of new materials. Many computational approaches require the generation of realistic structural models for property prediction and candidate screening. However, molecular and supramolecular materials represent classes of materials with many potential applications for which there is no go-to database of existing structures or general protocol for generating structures. Here, we report a new version of the supramolecular toolkit, stk, an open-source, extendable, and modular Python framework for general structure generation of (supra)molecular structures. Our construction approach works on arbitrary building blocks and topologies and minimizes the input required from the user, making stk user-friendly and applicable to many material classes. This version of stk includes metal-containing structures and rotaxanes as well as general implementation and interface improvements. Additionally, this version includes built-in tools for exploring chemical space with an evolutionary algorithm and tools for database generation and visualization. The latest version of stk is freely available at github.com/lukasturcani/stk.The surface of ice in contact with water contains sites that undergo deprotonation and protonation and can act as adsorption sites for aqueous ions. Therefore, an electrical double layer should form at this interface and existing models for describing the electrical double layer at metal oxide-water interfaces should be able to be modified to describe the surface charge, surface potential, and ionic occupancy at the ice-water interface. I used a surface complexation model along with literature measurements of the zeta potential of ice in brines of various strength and pH to constrain equilibrium constants. I then made predictions of ion site occupancy, surface charge density, and partitioning of counterions between the Stern and diffuse layers. The equilibrium constant for cation adsorption is more than 5 orders of magnitude larger than the other constants, indicating that this reaction dominates even at low salinity. Deprotonated OH sites are predicted to be slightly more abundant than dangling O sites, consistent with previous work. Surface charge densities are on the order of ±0.001 C/m2 and are always negative at the moderate pH values of interest to atmospheric and geophysical applications (6-9). In this pH range, over 99% of the counterions are contained in the Stern layer. This suggests that diffuse layer polarization will not occur because the ionic concentrations in the diffuse layer are nearly identical to those in the bulk electrolyte and that electrical conduction and polarization in the Stern layer will be negligible due to reduced ion mobility.We present a version of the T-moves approach for treating nonlocal pseudopotentials in diffusion Monte Carlo, which has much smaller time-step errors than the existing T-moves approaches, while at the same time preserving desirable features such as the upper-bound property for the energy. In addition, we modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error. The latter is applicable not only to pseudopotential calculations but also to all-electron calculations.Photoisomerization in the retinal leads to a channel opening in rhodopsins that triggers translocation or pumping of ions/protons. Crystal structures of rhodopsins contain several structurally conserved water molecules. It has been suggested that water plays an active role in facilitating the ion pumping/translocation process by acting as a lubricant in these systems. In this paper, we systematically investigate the localization, structure, dynamics, and energetics of the water molecules along the channel for the resting/dark state of KR2 rhodopsin. By employing several microseconds long atomistic molecular dynamics simulation of this trans-membrane protein system, we demonstrate the presence of five distinct water containing pockets/cavities separated by gateways controlled by protein side-chains. There exists a strong hydrogen bonded network involving these buried water molecules and functionally important key residues. We present evidence of significant structural and dynamical heterogeneity in the water molecules present in these cavities, with very rare exchange between them. The exchange time scale of such buried water with the bulk has an extremely wide range, from tens of nanoseconds to >1.5 ?s. The translational and rotational dynamics of buried water are found to be strongly dependent on the protein cavity size and local interactions with a classic signature of trapped diffusion and rotational anisotropy.