Selection of terms for a CI wavefunction to preserve potential surface features

A cumulative selection procedure for choosing configuration functions for inclusion in CI calculations is described. The objective of the method is to obtain equal energy loss, relative to unselected calculations, for different states and different regions of the potential surface. Results obtained from calculations on the BH molecule indicate an overall advantage in comparison to the threshold selection procedure, particularly with regard to molecular geometry changes.     

A low-energy passage for C+ + H2 → CH+(1Σ+) + H

Ab initio (SCF and CI) calculations have been performed for the linear approach (C_∝v) of C⁺ to H₂. For the ² Σ⁺ surface the saddle point and barrier height were determined. The interaction of the ²Σ⁺ and ²Π surfaces was investigated, leading to the conclusion that in near-C_∝v symmetry a low-energy path exists by which CH⁺¹ Σ⁺ can be formed, with negligible barrier in excess of endothermicity.     

Disappearing Polymorphs in Metal–Organic Framework Chemistry: Unexpected Stabilization of a Layered Polymorph over an Interpenetrated Three-Dimensional Structure in Mercury Imidazolate

The “disappearing polymorph” phenomenon is well established in organic solids, and has had a profound effect in pharmaceutical materials science. The first example of this effect in metal-containing systems in general, and in coordination-network solids in particular, is here reported. Specifically, attempts to mechanochemically synthesize a known interpenetrated diamondoid (dia) mercury(II) imidazolate metal–organic framework (MOF) yielded a novel, more stable polymorph based on square-grid (sql) layers. Simultaneously, the dia-form was found to be highly elusive, observed only as a short-lived intermediate in monitoring solvent-free synthesis and not at all from solution. The destabilization of a dense dia-framework relative to a lower dimensionality one is in contrast to the behavior of other imidazolate MOFs, with periodic density functional theory (DFT) calculations showing that it arises from weak interactions, including structure-stabilizing agostic C−H⋅⋅⋅Hg contacts. While providing a new link between MOFs and crystal engineering of organic solids, these findings highlight a possible role for agostic interactions in directing topology and stability of MOF polymorphs.     

Order-Unity 13C Nuclear Polarization of [1-13C]Pyruvate in Seconds and the Interplay of Water and SABRE Enhancement

Signal Amplification By Reversible Exchange in SHield Enabled Alignment Transfer (SABRE-SHEATH) is investigated to achieve rapid hyperpolarization of ¹³C₁ spins of [1-¹³C]pyruvate, using parahydrogen as the source of nuclear spin order. Pyruvate exchange with an iridium polarization transfer complex can be modulated via a sensitive interplay between temperature and co-ligation of DMSO and H₂O. Order-unity ¹³C (>50 %) polarization of catalyst-bound [1-¹³C]pyruvate is achieved in less than 30 s by restricting the chemical exchange of [1-¹³C]pyruvate at lower temperatures. On the catalyst bound pyruvate, 39 % polarization is measured using a 1.4 T NMR spectrometer, and extrapolated to >50 % at the end of build-up in situ. The highest measured polarization of a 30-mM pyruvate sample, including free and bound pyruvate is 13 % when using 20 mM DMSO and 0.5 M water in CD₃OD. Efficient ¹³C polarization is also enabled by favorable relaxation dynamics in sub-microtesla magnetic fields, as indicated by fast polarization buildup rates compared to the T₁ spin-relaxation rates (e. g., ∼0.2 s⁻¹ versus ∼0.1 s⁻¹, respectively, for a 6 mM catalyst-[1-¹³C]pyruvate sample). Finally, the catalyst-bound hyperpolarized [1-¹³C]pyruvate can be released rapidly by cycling the temperature and/or by optimizing the amount of water, paving the way to future biomedical applications of hyperpolarized [1-¹³C]pyruvate produced via comparatively fast and simple SABRE-SHEATH-based approaches.     

Experimental measurement of the strength of a C alpha-H...O bond in a lipid bilayer

Due to the apolar nature of the lipid bilayer, the weak Calpha-H...O H-bond is thought to contribute significantly toward the stability of transmembrane helical bundles such as glycophorin A (GPA). Here for the first time we measured the strength of such a bond, using vibrational frequency shifts of a dimeric and nondimeric variants of GPA containing a Gly CD2 label. Although the resulting estimated bond strength of 0.88 kcal/mol is relatively weak, several such bonds could contribute significantly toward bundle stabilization.     

Photoluminescence spectra of MOCVD-grown P-doped GaAS/Alx Ga1-x As MQW

PhotoLuminescence (PL) measurement techniques have been used to investigate on MOCVD grown P-doped GaAs/Al_xGa_1mx As (x=0.3) Multiple Quantum Wells (MQW). The spectra reveal extrinsic luminescence characteristics of e-A⁰ transitions for interface and centre of well acceptors in addition to both bound and free exciton emissions.     

A conserved asparagine in a ubiquitin-conjugating enzyme positions the substrate for nucleophilic attack

The mechanism used by the ubiquitin-conjugating enzyme, Ubc13, to catalyze ubiquitination is probed with three computational techniques: Born–Oppenheimer molecular dynamics, single point quantum mechanics/molecular mechanics energies, and classical molecular dynamics. These simulations support a long-held hypothesis and show that Ubc13-catalyzed ubiquitination uses a stepwise, nucleophilic attack mechanism. Furthermore, they show that the first step—the formation of a tetrahedral, zwitterionic intermediate—is rate limiting. However, these simulations contradict another popular hypothesis that supposes that the negative charge on the intermediate is stabilized by a highly conserved asparagine (Asn79 in Ubc13). Instead, calculated reaction profiles of the N79A mutant illustrate how charge stabilization actually increases the barrier to product formation. Finally, an alternate role for Asn79 is suggested by simulations of wild-type, N79A, N79D, and H77A Ubc13: it stabilizes the motion of the electrophile prior to the reaction, positioning it for nucleophilic attack. © 2019 Wiley Periodicals, Inc.     

Comparison of high-order many-body perturbation theory and configuration interaction for H2O

Diagrammatic many-body perturbation theory, coupled with a recursive computational procedure, is employed to obtain the correlation energy of H₂O within a 39-STO basis set by evaluating all double-excitation diagrams through twelfth order without any approximations. This provides, for the first time, the complete double-excitation diagrams contributions to the correlation energy, which is ?0.28826 hartree, compared with a correlation energy of ?0.27402 hartree obtained from a configuration interaction calculation which includes all double excitations. The difference of 0.0142 hartree includes the “size consistency” correction to the all-double-excitations CI energy, due to the “pathological” unliked-diagram terms remaining in that result, but also involves certain fourth- and higher-order rearrangement diagrams. Contrary to conventional belief, the unshifted, or Møller-Plesset partitioning of the hamiltonian provides a much more rapid convergence of the perturbation series that does the shifted, or Epstein-Nesbet partitioning. In both cases. Padé approximants enhance the convergence of the series considerably. A simple variation-perturbation scheme based on the first-order MBPT wavefunction is sufficient to provide 97.5% of the all-doubles CI correlation energy.     

Analysis of Hydrogen Tunneling in an Enzyme Active Site using von Neumann Measurements

We build on our earlier quantum wavepacket study of hydrogen transfer in the biological enzyme, soybean lipoxygenase-1, by using von Neumann quantum measurement theory to gain qualitative insights into the transfer event. We treat the enzyme active site as a measurement device which acts on the tunneling hydrogen nucleus via the potential it exerts at each configuration. A series of changing active site geometries during the tunneling process effects a sequential projection of the initial, reactant state onto the final, product state. We study this process using several different kinds of von Neumann measurements and show how a discrete sequence of such measurements not only progressively increases the projection of the hydrogen nuclear wavepacket onto the product side but also favors proton over deuteron transfer. Several qualitative features of the hydrogen tunneling problem found in wavepacket dynamics studies are also recovered here. These include the shift in the "transition state" towards the reactant as a result of nuclear quantization, greater participation of excited states in the case of deuterium, and presence of critical points along the reaction coordinate that facilitate hydrogen and deuterium transfer and coincide with surface crossings. To further "tailor" the dynamics, we construct a perturbation to the sequence of measurements, that is a perturbation to the dynamical sequence of active site geometry evolution, which leads us to insight on the existence of sensitive regions of the reaction profile where subtle changes to the dynamics of the active site can have an effect on the hydrogen and deuterium transfer process.