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.

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.

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.

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.

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.

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.

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.

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.

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.

The nephronophthisis gene product NPHP2/Inversin interacts with Aurora A and interferes with HDAC6-mediated cilia disassembly.

BACKGROUND: Nephronophthisis (NPH) is a rare recessive disease caused by several different gene mutations. Most gene products localize to the cilium, and thus, the various NPH manifestations including kidney cysts and situs inversus have been linked to ciliary defects. RESULTS: Here, we describe that targeted knockdown of NPHP2 significantly reduced the number of cilia on polarized MDCK cells. As one of the underlying molecular mechanisms, we identified a direct interaction between NPHP2 and Aurora A, a cell cycle kinase that promotes ciliary disassembly after activation by Hef1. NPHP2 inhibited the phosphorylation and activation of Aurora A, and reduced its kinase activity in vitro. Aurora A and histone deacetylase inhibitors ameliorated the ciliogenesis defect in NPHP2-deficient MDCK cells, supporting our hypothesis that NPHP2 is involved in the control of ciliary disassembly. Furthermore, we observed that nephrocystin (NPHP1), an interaction partner of NPHP2, also binds Aurora A, exerting very similar inhibitory effects on Hef1-mediated Aurora A activation. CONCLUSIONS: Taken together, these findings suggest that NPHP gene products can interfere with ciliary disassembly through interaction with the Hef1/Aurora A module, thereby modulating cell cycle control and cell proliferation.

Inversin relays Frizzled-8 signals to promote proximal pronephros development.

Mutations of inversin cause type II nephronophthisis, an infantile autosomal recessive disease characterized by cystic kidney disease and developmental defects. Inversin regulates Wnt signaling and is required for convergent extension movements during early embryogenesis. We now show that Inversin is essential for Xenopus pronephros formation, involving two distinct and opposing forms of cell movements. Knockdown of Inversin abrogated both proximal pronephros extension and distal tubule differentiation, phenotypes similar to that of Xenopus deficient in Frizzled-8. Exogenous Inversin rescued the pronephric defects caused by lack of Frizzled-8, indicating that Inversin acts downstream of Frizzled-8 in pronephros morphogenesis. Depletion of Inversin prevents the recruitment of Dishevelled in response to Frizzled-8 and impeded the accumulation of Dishevelled at the apical membrane of tubular epithelial cells in vivo. Thus, defective tubule morphogenesis seems to contribute to the renal pathology observed in patients with nephronophthisis type II.

Regulation of ciliary polarity by the APC/C.

Planar cell polarity signaling controls a variety of polarized cell behaviors. In multiciliated Xenopus epidermal cells, recruitment of Dishevelled (Dvl) to the basal body and its localization to the center of the ciliary rootlet are required to correctly position the motile cilia. We now report that the anaphase-promoting complex (APC/C) recognizes a D-box motif of Dvl and ubiquitylates Dvl on a highly conserved lysine residue. Inhibition of APC/C function by knockdown of the ANAPC2 subunit disrupts the polarity of motile cilia and alters the directionality of the fluid movement along the epidermis of the Xenopus embryo. Our results suggest that the APC/C activity enables cilia to correctly polarize in Xenopus epidermal cells.

Scribble participates in Hippo signaling and is required for normal zebrafish pronephros development.

Spatial organization of cells and their appendages is controlled by the planar cell polarity pathway, a signaling cascade initiated by the protocadherin Fat in Drosophila. Vertebrates express 4 Fat molecules, Fat1-4. We found that depletion of Fat1 caused cyst formation in the zebrafish pronephros. Knockdown of the PDZ domain containing the adaptor protein Scribble intensified the cyst-promoting phenotype of Fat1 depletion, suggesting that Fat1 and Scribble act in overlapping signaling cascades during zebrafish pronephros development. Supporting the genetic interaction with Fat1, Scribble recognized the PDZ-binding site of Fat1. Depletion of Yes-associated protein 1 (YAP1), a transcriptional co-activator inhibited by Hippo signaling, ameliorated the cyst formation in Fat1-deficient zebrafish, whereas Scribble inhibited the YAP1-induced cyst formation. Thus, reduced Hippo signaling and subsequent YAP1 disinhibition seem to play a role in the development of pronephric cysts after depletion of Fat1 or Scribble. We hypothesize that Hippo signaling is required for normal pronephros development in zebrafish and that Scribble is a candidate link between Fat and the Hippo signaling cascade in vertebrates.

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.

Loss of nephrocystin-3 function can cause embryonic lethality, Meckel-Gruber-like syndrome, situs inversus, and renal-hepatic-pancreatic dysplasia.

Many genetic diseases have been linked to the dysfunction of primary cilia, which occur nearly ubiquitously in the body and act as solitary cellular mechanosensory organelles. The list of clinical manifestations and affected tissues in cilia-related disorders (ciliopathies) such as nephronophthisis is broad and has been attributed to the wide expression pattern of ciliary proteins. However, little is known about the molecular mechanisms leading to this dramatic diversity of phenotypes. We recently reported hypomorphic NPHP3 mutations in children and young adults with isolated nephronophthisis and associated hepatic fibrosis or tapetoretinal degeneration. Here, we chose a combinatorial approach in mice and humans to define the phenotypic spectrum of NPHP3/Nphp3 mutations and the role of the nephrocystin-3 protein. We demonstrate that the pcy mutation generates a hypomorphic Nphp3 allele that is responsible for the cystic kidney disease phenotype, whereas complete loss of Nphp3 function results in situs inversus, congenital heart defects, and embryonic lethality in mice. In humans, we show that NPHP3 mutations can cause a broad clinical spectrum of early embryonic patterning defects comprising situs inversus, polydactyly, central nervous system malformations, structural heart defects, preauricular fistulas, and a wide range of congenital anomalies of the kidney and urinary tract (CAKUT). On the functional level, we show that nephrocystin-3 directly interacts with inversin and can inhibit like inversin canonical Wnt signaling, whereas nephrocystin-3 deficiency leads in Xenopus laevis to typical planar cell polarity defects, suggesting a role in the control of canonical and noncanonical (planar cell polarity) Wnt signaling.

A clinical trial alert tool to recruit large patient samples and assess selection bias in general practice research.

BACKGROUND: Many research projects in general practice face problems when recruiting patients, often resulting in low recruitment rates and an unknown selection bias, thus limiting their value for health services research. The objective of the study is to evaluate the recruitment performance of the practice staff in 25 participating general practices when using a clinical trial alert (CTA) tool. METHODS: The CTA tool was developed for an osteoporosis survey of patients at risk for osteoporosis and fractures. The tool used data from electronic patient records (EPRs) to automatically identify the population at risk (net sample), to apply eligibility criteria, to contact eligible patients, to enroll and survey at least 200 patients per practice. The effects of the CTA intervention were evaluated on the basis of recruitment efficiency and selection bias. RESULTS: The CTA tool identified a net sample of 16,067 patients (range 162 to 1,316 per practice), of which the practice staff reviewed 5,161 (32%) cases for eligibility. They excluded 3,248 patients and contacted 1,913 patients. Of these, 1,526 patients (range 4 to 202 per practice) were successfully enrolled and surveyed. This made up 9% of the net sample and 80% of the patients contacted. Men and older patients were underrepresented in the study population. CONCLUSION: Although the recruitment target was unreachable for most practices, the practice staff in the participating practices used the CTA tool successfully to identify, document and survey a large patient sample. The tool also helped the research team to precisely determine a slight selection bias.

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.