Efficient and accurate treatment of electron correlations with Correlation Matrix Renormalization theory.
Ontology highlight
ABSTRACT:
SUBMITTER: Yao YX
PROVIDER: S-EPMC4551991 | biostudies-other | 2015
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We present a second-order N-electron valence state perturbation theory (NEVPT2) based on a density matrix renormalization group (DMRG) reference wave function that exploits a Cholesky decomposition of the two-electron repulsion integrals (CD-DMRG-NEVPT2). With a parameter-free multireference perturbation theory approach at hand, the latter allows us to efficiently describe static and dynamic correlation in large molecular systems. We demonstrate the applicability of CD-DMRG-NEVPT2 for spin-state energetics of spin-crossover complexes involving calculations with more than 1000 atomic basis functions. We first assess, in a study of a heme model, the accuracy of the strongly and partially contracted variant of CD-DMRG-NEVPT2 before embarking on resolving a controversy about the spin ground state of a cobalt tropocoronand complex.
| S-EPMC5312874 | biostudies-literature
Project description:The density matrix renormalization group (DMRG) is a powerful method to treat static correlation. Here we present an inexpensive way to calculate correlation energy starting from a DMRG wave function using pair-density functional theory (PDFT). We applied this new approach, called DMRG-PDFT, to study singlet-triplet gaps in polyacenes and polyacetylenes that require active spaces larger than the feasibility limit of the conventional complete active-space self-consistent field (CASSCF) method. The results match reasonably well with the most reliable literature values and have only a moderate dependence on the compression of the initial DMRG wave function. Furthermore, DMRG-PDFT is significantly less expensive than other commonly applied ways of adding additional correlation to DMRG, such as DMRG followed by multireference perturbation theory or multireference configuration interaction.
| S-EPMC6368241 | biostudies-literature
Project description:Motivated by multi-hop communication in unreliable wireless networks, we present a percolation theory for time-varying networks. We develop a renormalization group theory for a prototypical network on a regular grid, where individual links switch stochastically between active and inactive states. The question whether a given source node can communicate with a destination node along paths of active links is equivalent to a percolation problem. Our theory maps the temporal existence of multi-hop paths on an effective two-state Markov process. We show analytically how this Markov process converges towards a memoryless Bernoulli process as the hop distance between source and destination node increases. Our work extends classical percolation theory to the dynamic case and elucidates temporal correlations of message losses. Quantification of temporal correlations has implications for the design of wireless communication and control protocols, e.g. in cyber-physical systems such as self-organized swarms of drones or smart traffic networks.
| S-EPMC5964206 | biostudies-literature
Project description:Density-functional theory (DFT) is a rigorous and (in principle) exact framework for the description of the ground state properties of atoms, molecules and solids based on their electron density. While computationally efficient density-functional approximations (DFAs) have become essential tools in computational chemistry, their (semi-)local treatment of electron correlation has a number of well-known pathologies, e.g. related to electron self-interaction. Here, we present a type of machine-learning (ML) based DFA (termed Kernel Density Functional Approximation, KDFA) that is pure, non-local and transferable, and can be efficiently trained with fully quantitative reference methods. The functionals retain the mean-field computational cost of common DFAs and are shown to be applicable to non-covalent, ionic and covalent interactions, as well as across different system sizes. We demonstrate their remarkable possibilities by computing the free energy surface for the protonated water dimer at hitherto unfeasible gold-standard coupled cluster quality on a single commodity workstation.
| S-EPMC7804195 | biostudies-literature
Project description:The development of theories and methods devoted to the accurate calculation of the electronic quasi-particle states and levels of molecules, clusters and solids is of prime importance to interpret the experimental data. These quantum systems are often modelled by using the Born-Oppenheimer approximation where the coupling between the electrons and vibrational modes is not fully taken into account, and the electrons are treated as pure quasi-particles. Here, we show that in small diamond cages, called diamondoids, the electron-vibration coupling leads to the breakdown of the electron quasi-particle picture. More importantly, we demonstrate that the strong electron-vibration coupling is essential to properly describe the overall lineshape of the experimental photoemission spectrum. This cannot be obtained by methods within Born-Oppenheimer approximation. Moreover, we deduce a link between the vibronic states found by our many-body perturbation theory approach and the well-known Jahn-Teller effect.
| S-EPMC4844694 | biostudies-literature
Project description:Face-sharing octahedral dinuclear Cr(III) compounds with d3-d3 electronic configurations represent nontrivial examples of electronic complexity, posing particular challenges for theoretical and computational studies. A tris-hydroxy-bridged Cr(III)-Cr(III) system has proven to be a richly rewarding target for studies of magnetism and electron paramagnetic resonance spectroscopy. It was also reported to be a peculiarly difficult system to treat with density functional theory (DFT). In this work the magnetic coupling problem for this dimer is approached with broken-symmetry (BS)-DFT and multireference calculations that utilize the density matrix renormalization group (DMRG) to handle full-valence active spaces. BS-DFT is shown to recover the correct ordering and energy spacing of Heisenberg spin states if used in conjunction with appropriate spin projection procedures, albeit with pronounced functional sensitivity. The contrasting conclusions of previous studies are traced to incorrect inclusion of electronically excited configurations. Analysis of the direct and differential overlap of corresponding orbital pairs from the BS-DFT solution indicates that metal-metal through-space interaction is the dominant contributor to antiferromagnetic coupling. At the DFT level a procedure that utilizes pseudopotential substitution is demonstrated that allows evaluation of the direct exchange vs superexchange contributions. A complementary description is obtained with DMRG-SCF calculations that enable state-averaged CASSCF calculations with both metal and bridge orbitals in the active space. A localized orbital subspace analysis supports the DFT conclusions that in contrast to doubly bridged isoelectronic analogues, antiferromagnetic coupling in the chromium dimer arises primarily from direct metal-metal interaction but is significantly enhanced by ligand-mediated superexchange.
| S-EPMC6728064 | biostudies-literature
Project description:The cores of globular proteins are densely packed, resulting in complicated networks of structural interactions. These interactions in turn give rise to dynamic structural correlations over a wide range of time scales. Accurate analysis of these complex correlations is crucial for understanding biomolecular mechanisms and for relating structure to function. Here we report a highly accurate technique for inferring the major modes of structural correlation in macromolecules using likelihood-based statistical analysis of sets of structures. This method is generally applicable to any ensemble of related molecules, including families of nuclear magnetic resonance (NMR) models, different crystal forms of a protein, and structural alignments of homologous proteins, as well as molecular dynamics trajectories. Dominant modes of structural correlation are determined using principal components analysis (PCA) of the maximum likelihood estimate of the correlation matrix. The correlations we identify are inherently independent of the statistical uncertainty and dynamic heterogeneity associated with the structural coordinates. We additionally present an easily interpretable method ("PCA plots") for displaying these positional correlations by color-coding them onto a macromolecular structure. Maximum likelihood PCA of structural superpositions, and the structural PCA plots that illustrate the results, will facilitate the accurate determination of dynamic structural correlations analyzed in diverse fields of structural biology.
| S-EPMC2242818 | biostudies-literature
Project description:MOTIVATION:High throughput technologies are widely employed in modern biomedical research. They yield measurements of a large number of biomolecules in a single experiment. The number of experiments usually is much smaller than the number of measurements in each experiment. The simultaneous measurements of biomolecules provide a basis for a comprehensive, systems view for describing relevant biological processes. Often it is necessary to determine correlations between the data matrices under different conditions or pathways. However, the techniques for analyzing the data with a low number of samples for possible correlations within or between conditions are still in development. Earlier developed correlative measures, such as the RV coefficient, use the trace of the product of data matrices as the most relevant characteristic. However, a recent study has shown that the RV coefficient consistently overestimates the correlations in the case of low sample numbers. To correct for this bias, it was suggested to discard the diagonal elements of the outer products of each data matrix. In this work, a principled approach based on the matrix decomposition generates three trace-independent parts for every matrix. These components are unique, and they are used to determine different aspects of correlations between the original datasets. RESULTS:Simulations show that the decomposition results in the removal of high correlation bias and the dependence on the sample number intrinsic to the RV coefficient. We then use the correlations to analyze a real proteomics dataset. AVAILABILITY AND IMPLEMENTATION:The python code can be downloaded from http://dynamic-proteome.utmb.edu/MatrixCorrelations.aspx. SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
| S-EPMC6853692 | biostudies-literature
Project description:Correlation between fluctuations of the square of electron tunneling matrix element TDA2 and the donor-acceptor distance RDA in the electron transfer (ET) reaction from bacteriopheophytin anion to the primary quinone of the reaction center in the photosynthetic bacteria Rhodobacter sphaeroides is investigated by a combined study of molecular dynamics simulations of the protein conformation fluctuation and quantum chemical calculations. We adopted two kinds of RDA ; edge-to-edge distance REE and center-to-center distance RCC . The value of TDA2 distributed over more than 5 orders of magnitude and the fluctuation of the value of RDA distributed over more than 1.8 Å for the 106 instantaneous conformations of 1 ns simulation. We made analysis of the time-averaged correlation step by step as follows. We divide the 106 simulation data into 1000/t parts of small data set to obtain the averaged data points of <TDA2> t and <REE > t or <RCC > t . Plotting the 1000/t sets of log10 <TDA2> t as a function of <REE > t or <RCC > t , we made a principal coordinate analysis for these distributions. The slopes <?E > t and <?C > t of the primary axis are very large at small value of t and they are decreased considerably as t becomes large. The ellipticity for the distribution of <TDA2> tvs <REE > t which can be a measure for the degree of correlation became very small when t is large, while it does not hold for the distribution of <TDA2> tvs <RCC > t . These results indicate that only the correlation between <TDA2> t and <REE > t for large t satisfies the well-known linear relation ("Dutton law"), although the slope is larger than the original value 1.4 Å-1. Based on the present result, we examined the analysis of the dynamic disorder by means of the single-molecule spectroscopy by Xie and co-workers with use of the "Dutton law".
| S-EPMC5036608 | biostudies-literature
Project description:We present a method for computing the resonant inelastic x-ray scattering (RIXS) spectra in one-dimensional systems using the density matrix renormalization group (DMRG) method. By using DMRG to address this problem, we shift the computational bottleneck from the memory requirements associated with exact diagonalization (ED) calculations to the computational time associated with the DMRG algorithm. This approach is then used to obtain RIXS spectra on cluster sizes well beyond state-of-the-art ED techniques. Using this new procedure, we compute the low-energy magnetic excitations observed in Cu L-edge RIXS for the challenging corner shared CuO4 chains, both for large multi-orbital clusters and downfolded t-J chains. We are able to directly compare results obtained from both models defined in clusters with identical momentum resolution. In the strong coupling limit, we find that the downfolded t-J model captures the main features of the magnetic excitations probed by RIXS only after a uniform scaling of the spectra is made.
| S-EPMC6056525 | biostudies-literature