Sampling the reciprocal Coulomb potential in finite anisotropic cells.

We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions: the treatment of the singularity at the Brillouin-zone center and discretization errors, which can cause severe convergence problems in anisotropic cells, necessary for the calculation of low-dimensional systems. We apply our strategy to the Hartree-Fock and coupled cluster (CC) theories and discuss the consequences of different sampling strategies on different theories. We show that sampling the Coulomb potential via the widely used probe-charge Ewald method is unsuitable for CC calculations in anisotropic cells. To demonstrate the applicability of our developed approach, we study two representative, low-dimensional use cases: the infinite carbon chain, for which we report the first periodic CCSD(T) potential energy surface, and a surface slab of lithium hydride, for which we demonstrate the impact of different sampling strategies for calculating surface energies. We find that our Coulomb sampling strategy serves as a vital solution, addressing the critical need for improved accuracy in plane-wave based CC calculations for low-dimensional systems.

Averting the Infrared Catastrophe in the Gold Standard of Quantum Chemistry.

Coupled-cluster theories can be used to compute ab initio electronic correlation energies of real materials with systematically improvable accuracy. However, the widely used coupled cluster singles and doubles plus perturbative triples [CCSD(T)] method is only applicable to insulating materials. For zero-gap materials the truncation of the underlying many-body perturbation expansion leads to an infrared catastrophe. Here, we present a novel perturbative triples formalism denoted as (cT) that yields convergent correlation energies in metallic systems. Furthermore, the computed correlation energies for the three-dimensional uniform electron gas at metallic densities are in good agreement with quantum Monte Carlo results. At the same time the newly proposed method retains all desirable properties of CCSD(T) such as its accuracy for insulating systems as well as its low computational cost compared to a full inclusion of the triples. This paves the way for ab initio calculations of real metals with chemical accuracy.

Clinical, Radiological, and Pathological Diagnosis of Fibro-Osseous Lesions of the Oral and Maxillofacial Region: A Retrospective Study.

BACKGROUND: Fibro-osseous lesions (FOL) of the jaw represent a rare, benign group of lesions that share similar clinical, radiological, and histopathological features and are characterized by progressive, variable replacement of healthy bone tissue by fibrous connective tissue. METHODS: This retrospective study aimed to evaluate the incidence of fibro-osseous lesions and to reassess the efficacy of case-specific treatment management from a clinical, radiological, and histopathological perspective based on 14 years of data. RESULTS: Forty-four patients with a radiological and/or histopathological diagnosis of benign FOLs were identified and re-evaluated. Cemento-osseous dysplasia was the most common group of FOLs present in our patient cohort (45%), followed by ossifying fibroma (39%) and fibrous dysplasia (16%). The diagnostic imaging technique of choice was CBCT (68%), followed by PAN (18%), with most patients (95 %) additionally undergoing biopsy. The mean age of the patients at the time of diagnosis was 40.54 ± 13.7 years, with most lesions being located in the mandible (86%), with females being predominantly affected (73%). CONCLUSION: An interdisciplinary approach that analyzes all case-specific factors, including demographic data, medical history, intraoperative findings, and, most importantly, histopathological and radiological features, is essential for an accurate diagnosis and key to avoiding inappropriate treatment.

Cerium Oxides without U: The Role of Many-Electron Correlation.

Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchmark of ab initio many-electron theory for electron transfer energies and lattice parameters under periodic boundary conditions. We show that the direct random phase approximation clearly outperforms all DFT variations. From this foundation, we, then, systematically improve even further. Periodic second-order Møller-Plesset perturbation theory meanwhile manages to recover standard hybrid functional values. Using these approaches to eliminate parameter bias allows for highly accurate benchmarks of strongly correlated materials, the reliable assessment of various density functionals, and functional fitting via machine-learning.

Local embedding of coupled cluster theory into the random phase approximation using plane waves.

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: The case of water on graphitic carbon nitride.

A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.

Initial morphological symmetry breaking in the foregut and development of the omental bursa in human embryos.

Bilaterally symmetrical primordia of visceral organs undergo asymmetrical morphogenesis leading to typical arrangement of visceral organs in the adult. Asymmetrical morphogenesis within the upper abdomen leads, among others, to the formation of the omental bursa dorsally to the rotated stomach. A widespread view of this process assumes kinking of thin mesenteries as a main mechanism. This view is based on a theory proposed already by Johannes Müller in 1830 and was repeatedly criticized, but some of the most plausible alternative views (initially proposed by Swaen in 1897 and Broman in 1904) still remain to be proven. Here, we analyzed serial histological sections of human embryos between stages 12 and 15 at high light microscopical resolution to reveal the succession of events giving rise to the development of the omental bursa and its relation to the emerging stomach asymmetry. Our analysis indicates that morphological symmetry breaking in the upper abdomen occurs within a wide mesenchymal plate called here mesenteric septum and is based on differential behavior of the coelomic epithelium which causes asymmetric paragastric recess formation and, importantly, precedes initial rotation of stomach. Our results thus provide the first histological evidence of breaking the symmetry of the early foregut anlage in the human embryo and pave the way for experimental studies of left-right symmetry breaking in the upper abdomen in experimental model organisms.

A shortcut to the thermodynamic limit for quantum many-body calculations of metals.

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.

Human C-terminal CUBN variants associate with chronic proteinuria and normal renal function.

BACKGROUNDProteinuria is considered an unfavorable clinical condition that accelerates renal and cardiovascular disease. However, it is not clear whether all forms of proteinuria are damaging. Mutations in CUBN cause Imerslund-Gräsbeck syndrome (IGS), which is characterized by intestinal malabsorption of vitamin B12 and in some cases proteinuria. CUBN encodes for cubilin, an intestinal and proximal tubular uptake receptor containing 27 CUB domains for ligand binding.METHODSWe used next-generation sequencing for renal disease genes to genotype cohorts of patients with suspected hereditary renal disease and chronic proteinuria. CUBN variants were analyzed using bioinformatics, structural modeling, and epidemiological methods.RESULTSWe identified 39 patients, in whom biallelic pathogenic variants in the CUBN gene were associated with chronic isolated proteinuria and early childhood onset. Since the proteinuria in these patients had a high proportion of albuminuria, glomerular diseases such as steroid-resistant nephrotic syndrome or Alport syndrome were often the primary clinical diagnosis, motivating renal biopsies and the use of proteinuria-lowering treatments. However, renal function was normal in all cases. By contrast, we did not found any biallelic CUBN variants in proteinuric patients with reduced renal function or focal segmental glomerulosclerosis. Unlike the more N-terminal IGS mutations, 37 of the 41 proteinuria-associated CUBN variants led to modifications or truncations after the vitamin B12-binding domain. Finally, we show that 4 C-terminal CUBN variants are associated with albuminuria and slightly increased GFR in meta-analyses of large population-based cohorts.CONCLUSIONCollectively, our data suggest an important role for the C-terminal half of cubilin in renal albumin reabsorption. Albuminuria due to reduced cubilin function could be an unexpectedly common benign condition in humans that may not require any proteinuria-lowering treatment or renal biopsy.FUNDINGATIP-Avenir program, Fondation Bettencourt-Schueller (Liliane Bettencourt Chair of Developmental Biology), Agence Nationale de la Recherche (ANR) Investissements d'avenir program (ANR-10-IAHU-01) and NEPHROFLY (ANR-14-ACHN-0013, to MS), Steno Collaborative Grant 2018 (NNF18OC0052457, to TSA and MS), Heisenberg Professorship of the German Research Foundation (KO 3598/5-1, to AK), Deutsche Forschungsgemeinschaft (DFG) Collaborative Research Centre (SFB) KIDGEM 1140 (project 246781735, to CB), and Federal Ministry of Education and Research (BMB) (01GM1515C, to CB).

RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods.

We present a method to approximate post-Hartree-Fock correlation energies by using approximate natural orbitals obtained by the random phase approximation (RPA). We demonstrate the method by applying it to the helium atom, the hydrogen and fluorine molecule, and to diamond as an example of a periodic system. For these benchmark systems, we show that RPA natural orbitals converge the MP2 correlation energy rapidly. Additionally, we calculated full configuration interaction energies for He and H2, which are in excellent agreement with the literature and experimental values. We conclude that the proposed method may serve as a compromise to reach good approximations to correlation energies at moderate computational cost, and we expect the method to be especially useful for theoretical studies on surface chemistry by providing an efficient basis to correlated wave function based methods.

Novel biologics targeting the P2X7 ion channel.

Targeting the P2X7 ion channel, a danger sensor for extracellular nucleotides, improves outcomes in models of inflammation, cancer, and brain-diseases. Antibodies and nanobodies have been developed that antagonize or potentiate gating of P2X7. Their potential advantages over small-molecule drugs include high specificity, lower off-target effects, and tunable in vivo half-life. Genetic fusion of P2X7-specific biologics to binding modules may enable targeting of specific cell subsets. Besides directly modulating P2X7 function, antibodies can also initiate specific depletion of P2X7-expressing cells. Adeno-associated viral vectors (AAV) can be used to express P2X7-specific antibodies in vivo to achieve long-lasting biological effects. Furthermore, if successfully targeted to P2X7-expressing cells, AAVs may enable modulation of the function of P2X7-expressing immune cells via encoded transgenic RNA or proteins.

T Cell Phenotype and T Cell Receptor Repertoire in Patients with Major Depressive Disorder.

While a link between inflammation and the development of neuropsychiatric disorders, including major depressive disorder (MDD) is supported by a growing body of evidence, little is known about the contribution of aberrant adaptive immunity in this context. Here, we conducted in-depth characterization of T cell phenotype and T cell receptor (TCR) repertoire in MDD. For this cross-sectional case-control study, we recruited antidepressant-free patients with MDD without any somatic or psychiatric comorbidities (n = 20), who were individually matched for sex, age, body mass index, and smoking status to a non-depressed control subject (n = 20). T cell phenotype and repertoire were interrogated using a combination of flow cytometry, gene expression analysis, and next generation sequencing. T cells from MDD patients showed significantly lower surface expression of the chemokine receptors CXCR3 and CCR6, which are known to be central to T cell differentiation and trafficking. In addition, we observed a shift within the CD4+ T cell compartment characterized by a higher frequency of CD4+CD25highCD127low/- cells and higher FOXP3 mRNA expression in purified CD4+ T cells obtained from patients with MDD. Finally, flow cytometry-based TCR Vß repertoire analysis indicated a less diverse CD4+ T cell repertoire in MDD, which was corroborated by next generation sequencing of the TCR ß chain CDR3 region. Overall, these results suggest that T cell phenotype and TCR utilization are skewed on several levels in patients with MDD. Our study identifies putative cellular and molecular signatures of dysregulated adaptive immunity and reinforces the notion that T cells are a pathophysiologically relevant cell population in this disorder.

Laplace transformed MP2 for three dimensional periodic materials using stochastic orbitals in the plane wave basis and correlated sampling.

We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three-dimensional periodic systems in second-order Møller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are required, as it is necessary for calculations of adsorption, binding, or surface energies. The algorithm is based on the Laplace transformed MP2 (LTMP2) formulation in the plane wave basis. The time-dependent Hartree-Fock orbitals, appearing in the LTMP2 formulation, are stochastically rotated in the occupied and unoccupied Hilbert space. This avoids a full summation over all combinations of occupied and unoccupied orbitals, as inspired by the work of Neuhauser, Rabani, and Baer [J. Chem. Theory Comput. 9, 24 (2013)]. Additionally, correlated sampling is introduced, accelerating the statistical convergence significantly.

Assessing Density Functionals Using Many Body Theory for Hybrid Perovskites.

Which density functional is the "best" for structure simulations of a particular material? A concise, first principles, approach to answer this question is presented. The random phase approximation (RPA)-an accurate many body theory-is used to evaluate various density functionals. To demonstrate and verify the method, we apply it to the hybrid perovskite MAPbI_{3}, a promising new solar cell material. The evaluation is done by first creating finite temperature ensembles for small supercells using RPA molecular dynamics, and then evaluating the variance between the RPA and various approximate density functionals for these ensembles. We find that, contrary to recent suggestions, van der Waals functionals do not improve the description of the material, whereas hybrid functionals and the strongly constrained appropriately normed (SCAN) density functional yield very good agreement with the RPA. Finally, our study shows that in the room temperature tetragonal phase of MAPbI_{3}, the molecules are preferentially parallel to the shorter lattice vectors but reorientation on ps time scales is still possible.

Quartic scaling MP2 for solids: A highly parallelized algorithm in the plane wave basis.

We present a low-complexity algorithm to calculate the correlation energy of periodic systems in second-order Møller-Plesset (MP2) perturbation theory. In contrast to previous approximation-free MP2 codes, our implementation possesses a quartic scaling, O(N4), with respect to the system size N and offers an almost ideal parallelization efficiency. The general issue that the correlation energy converges slowly with the number of basis functions is eased by an internal basis set extrapolation. The key concept to reduce the scaling is to eliminate all summations over virtual orbitals which can be elegantly achieved in the Laplace transformed MP2 formulation using plane wave basis sets and fast Fourier transforms. Analogously, this approach could allow us to calculate second order screened exchange as well as particle-hole ladder diagrams with a similar low complexity. Hence, the presented method can be considered as a step towards systematically improved correlation energies.

Analytic Interatomic Forces in the Random Phase Approximation.

We discuss that in the random phase approximation (RPA) the first derivative of the energy with respect to the Green's function is the self-energy in the GW approximation. This relationship allows us to derive compact equations for the RPA interatomic forces. We also show that position dependent overlap operators are elegantly incorporated in the present framework. The RPA force equations have been implemented in the projector augmented wave formalism, and we present illustrative applications, including ab initio molecular dynamics simulations, the calculation of phonon dispersion relations for diamond and graphite, as well as structural relaxations for water on boron nitride. The present derivation establishes a concise framework for forces within perturbative approaches and is also applicable to more involved approximations for the correlation energy.

Dithiocarbamate Self-Assembled Monolayers as Efficient Surface Modifiers for Low Work Function Noble Metals.

Tuning the work function of the electrode is one of the crucial steps to improve charge extraction in organic electronic devices. Here, we show that N,N-dialkyl dithiocarbamates (DTC) can be effectively employed to produce low work function noble metal electrodes. Work functions between 3.1 and 3.5 eV are observed for all metals investigated (Cu, Ag, and Au). Ultraviolet photoemission spectroscopy (UPS) reveals a maximum decrease in work function by 2.1 eV as compared to the bare metal surface. Electronic structure calculations elucidate how the complex interplay between intrinsic dipoles and dipoles induced by bond formation generates such large work function shifts. Subsequently, we quantify the improvement in contact resistance of organic thin film transistor devices with DTC coated source and drain electrodes. These findings demonstrate that DTC molecules can be employed as universal surface modifiers to produce stable electrodes for electron injection in high performance hybrid organic optoelectronics.

SUMOylation Blocks the Ubiquitin-Mediated Degradation of the Nephronophthisis Gene Product Glis2/NPHP7.

Glis2/NPHP7 is a transcriptional regulator mutated in type 7 nephronophthisis, an autosomal recessive ciliopathy associated with cystic and fibrotic kidney disease as well as characteristic extrarenal manifestations. While most ciliopathy-associated molecules are found in the cilium, Glis2/NPHP7 presumably localizes to the nucleus. However, the detection of endogenous Glis2/NPHP7 has remained unsuccessful, potentially due to its ubiquitylation-dependent rapid degradation. We report now that Glis2/NPHP7 is also SUMOylated, preferentially by PIAS4, which conjugates Glis2/NPHP7 to SUMO3. SUMOylation interferes with ubiquitylation and degradation of Glis2/NPHP7, suggesting that Glis2/NPHP7 protein levels are regulated by competing ubiquitylation/ SUMOylation. SUMOylation also alters the transcriptional activity of Glis2/NPHP7. While Glis2/NPHP7 activates the mouse insulin-2-promotor (mIns2), SUMOylated Glis2/NPHP7 lacks this property, and seems to act as a repressor. Taken together, our data reveal that Glis2/NPHP7 is extensively modified by post-translational modifications, suggesting that a tight control of this transcriptional regulator is required for normal development and tissue homeostasis.

Interaction with the Bardet-Biedl gene product TRIM32/BBS11 modifies the half-life and localization of Glis2/NPHP7.

Although the two ciliopathies Bardet-Biedl syndrome and nephronophthisis share multiple clinical manifestations, the molecular basis for this overlap remains largely unknown. Both BBS11 and NPHP7 are unusual members of their respective gene families. Although BBS11/TRIM32 represents a RING finger E3 ubiquitin ligase also involved in hereditary forms of muscular dystrophy, NPHP7/Glis2 is a Gli-like transcriptional repressor that localizes to the nucleus, deviating from the ciliary localization of most other ciliopathy-associated gene products. We found that BBS11/TRIM32 and NPHP7/Glis2 can physically interact with each other, suggesting that both proteins form a functionally relevant protein complex in vivo. This hypothesis was further supported by the genetic interaction and synergist cyst formation in the zebrafish pronephros model. However, contrary to our expectation, the E3 ubiquitin ligase BBS11/TRIM32 was not responsible for the short half-life of NPHP7/Glis2 but instead promoted the accumulation of mixed Lys(48)/Lys(63)-polyubiquitylated NPHP7/Glis2 species. This modification not only prolonged the half-life of NPHP7/Glis2, but also altered the subnuclear localization and the transcriptional activity of NPHP7/Glis2. Thus, physical and functional interactions between NPHP and Bardet-Biedl syndrome gene products, demonstrated for Glis2 and TRIM32, may help to explain the phenotypic similarities between these two syndromes.