Self-Consistent Calculation of the Localized Orbital Scaling Correction for Correct Electron Densities and Energy-Level Alignments in Density Functional Theory.
Ontology highlight
ABSTRACT:
SUBMITTER: Mei Y
PROVIDER: S-EPMC8341159 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.
| S-EPMC5135379 | biostudies-literature
Project description:This work describes the extension of a previously reported empirical localized orbital correction model for density functional theory (DFT-LOC) for atomization energies, ionization potentials, electron affinities, and reaction enthalpies to the correction of barrier heights for chemical reactions of various types including cycloadditions, cycloreversions, dipolar cycloadditions, S(N)2's, carbon radical reactions, hydrogen radical reactions, sigmatropic shifts, and electrocyclizations. The B3LYP localized orbital correction version of the model (B3LYP-LOC) reduces the number of outliers and overall mean unsigned error (MUE) vs. experiment or ab initio values from 3.2 to 1.3 kcal/mole for barrier heights and from 5.1 to 1.1 kcal/mole for reaction enthalpies versus B3LYP. Furthermore, the new model has essentially zero additional computational cost beyond standard DFT calculations. Although the model is heuristic and is based on multiple linear regression to experimental or ab initio data, each of the parameters is justified on chemical grounds and provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation.
| S-EPMC2808043 | biostudies-other
Project description:Our previous works have demonstrated the ability of our localized orbital correction (LOC) methodology to greatly improve the accuracy of various thermochemical properties at the stationary points of the Density Functional Theory (DFT) reaction coordinate (RC). Herein we extend this methodology from stationary points to the entire RC connecting any stationary points by developing continuous localized orbital corrections (CLOCs). We show that the resultant method, DFT-CLOC, is capable of producing RCs with far greater accuracy than uncorrected DFT and yet requires negligible computational cost beyond the uncorrected DFT calculations. Various post-Hartree-Fock (post-HF) reaction coordinate profiles were used, including a sigmatropic shift, Diels-Alder reaction, electrocyclization, carbon radical and three hydrogen radical reactions to show that this method is robust across multiple reaction types of general interest.
| S-EPMC3030200 | biostudies-literature
Project description:A new self-consistent procedure for calculating the total energy with an orbital-dependent density functional approximation (DFA), the generalized optimized effective potential (GOEP), is developed in the present work. The GOEP is a nonlocal Hermitian potential that delivers the sets of occupied and virtual orbitals and minimizes the total energy. The GOEP optimization leads to the same minimum as does the orbital optimization. The GOEP method is promising as an effective optimization approach for orbital-dependent functionals, as demonstrated for the self-consistent calculations of the random phase approximation (RPA) to the correlation functionals in the particle-hole (ph) and particle-particle (pp) channels. The results show that the accuracy in describing the weakly interacting van der Waals systems is significantly improved in the self-consistent calculations. In particular, the important single excitations contribution in non-self-consistent RPA calculations can be captured self-consistently through the GOEP optimization, leading to orbital renormalization, without using the single excitations in the energy functional.
| S-EPMC6209318 | biostudies-literature
Project description:The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier's wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine-cytosine (GC) base pairs and among 3 adenine-thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π-π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.
| S-EPMC7456463 | biostudies-literature
Project description:Sequence similarity searches performed with BLAST, SSEARCH and FASTA achieve high sensitivity by using scoring matrices (e.g. BLOSUM62) that target low identity (<33%) alignments. Although such scoring matrices can effectively identify distant homologs, they can also produce local alignments that extend beyond the homologous regions.We measured local alignment start/stop boundary accuracy using a set of queries where the correct alignment boundaries were known, and found that 7% of BLASTP and 8% of SSEARCH alignment boundaries were overextended. Overextended alignments include non-homologous sequences; they occur most frequently between sequences that are more closely related (>33% identity). Adjusting the scoring matrix to reflect the identity of the homologous sequence can correct higher identity overextended alignment boundaries. In addition, the scoring matrix that produced a correct alignment could be reliably predicted based on the sequence identity seen in the original BLOSUM62 alignment. Realigning with the predicted scoring matrix corrected 37% of all overextended alignments, resulting in more correct alignments than using BLOSUM62 alone.
| S-EPMC3834790 | biostudies-literature
Project description:Networks distribute energy, materials and information to the components of a variety of natural and human-engineered systems, including organisms, brains, the Internet and microprocessors. Distribution networks enable the integrated and coordinated functioning of these systems, and they also constrain their design. The similar hierarchical branching networks observed in organisms and microprocessors are striking, given that the structure of organisms has evolved via natural selection, while microprocessors are designed by engineers. Metabolic scaling theory (MST) shows that the rate at which networks deliver energy to an organism is proportional to its mass raised to the 3/4 power. We show that computational systems are also characterized by nonlinear network scaling and use MST principles to characterize how information networks scale, focusing on how MST predicts properties of clock distribution networks in microprocessors. The MST equations are modified to account for variation in the size and density of transistors and terminal wires in microprocessors. Based on the scaling of the clock distribution network, we predict a set of trade-offs and performance properties that scale with chip size and the number of transistors. However, there are systematic deviations between power requirements on microprocessors and predictions derived directly from MST. These deviations are addressed by augmenting the model to account for decentralized flow in some microprocessor networks (e.g. in logic networks). More generally, we hypothesize a set of constraints between the size, power and performance of networked information systems including transistors on chips, hosts on the Internet and neurons in the brain.
| S-EPMC2607348 | biostudies-literature
Project description:Metabolic rate, heart rate, lifespan, and many other physiological properties vary with body mass in systematic and interrelated ways. Present empirical data suggest that these scaling relationships take the form of power laws with exponents that are simple multiples of one quarter. A compelling explanation of this observation was put forward a decade ago by West, Brown, and Enquist (WBE). Their framework elucidates the link between metabolic rate and body mass by focusing on the dynamics and structure of resource distribution networks-the cardiovascular system in the case of mammals. Within this framework the WBE model is based on eight assumptions from which it derives the well-known observed scaling exponent of 3/4. In this paper we clarify that this result only holds in the limit of infinite network size (body mass) and that the actual exponent predicted by the model depends on the sizes of the organisms being studied. Failure to clarify and to explore the nature of this approximation has led to debates about the WBE model that were at cross purposes. We compute analytical expressions for the finite-size corrections to the 3/4 exponent, resulting in a spectrum of scaling exponents as a function of absolute network size. When accounting for these corrections over a size range spanning the eight orders of magnitude observed in mammals, the WBE model predicts a scaling exponent of 0.81, seemingly at odds with data. We then proceed to study the sensitivity of the scaling exponent with respect to variations in several assumptions that underlie the WBE model, always in the context of finite-size corrections. Here too, the trends we derive from the model seem at odds with trends detectable in empirical data. Our work illustrates the utility of the WBE framework in reasoning about allometric scaling, while at the same time suggesting that the current canonical model may need amendments to bring its predictions fully in line with available datasets.
| S-EPMC2518954 | biostudies-literature
Project description:Water plays a unique role in all living organisms. Not only is it nature's ubiquitous solvent, but it also actively takes part in many cellular processes. In particular, the structure and properties of interfacial water near biomolecules such as proteins are often related to the function of the respective molecule. It can therefore be highly instructive to study the local water density around solutes in cellular systems, particularly when solvent-mediated forces such as the hydrophobic effect are relevant. Computational methods such as molecular dynamics (MD) simulations seem well suited to study these systems at the atomic level. However, due to sampling requirements, it is not clear that MD simulations are, indeed, the method of choice to obtain converged densities at a given level of precision. We here compare the calculation of local water densities with two different methods: MD simulations and the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). In particular, we investigate the convergence of the local water density to assess the required simulation times for different levels of resolution. Moreover, we provide a quantitative comparison of the densities calculated with MD and with 3D-RISM-KH and investigate the effect of the choice of the water model for both methods. Our results show that 3D-RISM-KH yields density distributions that are very similar to those from MD up to a 0.5 Å resolution, but for significantly reduced computational cost. The combined use of MD and 3D-RISM-KH emerges as an auspicious perspective for efficient solvent sampling in dynamical systems.
| S-EPMC3407544 | biostudies-literature
Project description:BackgroundAccurate measurement and reconstruction of orbital soft tissue is important to diagnosis and treatment of orbital diseases. This study applied an interactive graph cut method to orbital soft tissue precise segmentation and calculation in computerized tomography (CT) images, and to estimate its application in orbital reconstruction.MethodsThe interactive graph cut method was introduced to segment extraocular muscle and intraorbital fat in CT images. Intra- and inter-observer variability of tissue volume measured by graph cut segmentation was validated. Accuracy and reliability of the method was accessed by comparing with manual delineation and commercial medical image software. Intraorbital structure of 10 patients after enucleation surgery was reconstructed based on graph cut segmentation and soft tissue volume were compared within two different surgical techniques.ResultsBoth muscle and fat tissue segmentation results of graph cut method showed good consistency with ground truth in phantom data. There were no significant differences in muscle calculations between observers or segmental methods (p > 0.05). Graph cut results of fat tissue had coincidental variable trend with ground truth which could identify 0.1cm3 variation. The mean performance time of graph cut segmentation was significantly shorter than manual delineation and commercial software (p < 0.001). Jaccard similarity and Dice coefficient of graph cut method were 0.767 ± 0.045 and 0.836 ± 0.032 for human normal extraocular muscle segmentation. The measurements of fat tissue were significantly better in graph cut than those in commercial software (p < 0.05). Orbital soft tissue volume was decreased in post-enucleation orbit than that in normal orbit (p < 0.05).ConclusionThe graph cut method was validated to have good accuracy, reliability and efficiency in orbit soft tissue segmentation. It could discern minor volume changes of soft tissue. The interactive segmenting technique would be a valuable tool for dynamic analysis and prediction of therapeutic effect and orbital reconstruction.
| S-EPMC6916112 | biostudies-literature