Browse
Submit Data
Databases
API
Help

Dataset Information

18 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package.

ABSTRACT: Constrained density functional theory (CDFT) is a powerful tool for the prediction of electron transfer parameters in condensed phase simulations at a reasonable computational cost. In this work we present an extension to CDFT in the popular mixed Gaussian/plane wave electronic structure package CP2K, implementing the additional force terms arising from a constraint based on Hirshfeld charge partitioning. This improves upon the existing Becke partitioning scheme, which is prone to give unphysical atomic charges. We verify this implementation for a variety of systems: electron transfer in (H₂O)₂⁺ in a vacuum, electron tunnelling between oxygen vacancy centers in solid MgO, and electron self-exchange in aqueous Ru²⁺-Ru³⁺. We find good agreement with previous plane-wave CDFT results for the same systems, but at a significantly lower computational cost, and we discuss the general reliability of condensed phase CDFT calculations.

SUBMITTER: Ahart CS

PROVIDER: S-EPMC9281399 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Photochemical Ring Opening of Oxirane Modeled by Constrained Density Functional Theory.

Project description:A constrained density functional theory/classical trajectory surface hopping study of the photochemical dissociation of oxirane (CH2)2O is presented. The calculations confirm the Gomer-Noyes mechanism for the initial reaction and agree largely with experimental photolysis data including reaction yields. The calculated yields, however, depend both on temperature and its modeling. The timescales of the various reaction steps are well below 100 fs, similar to previous time-dependent density functional calculations. At variance with those, however, the present calculations obey Kasha's rule, i.e., the photoreaction is initiated in the energetically lowest excited state.

| S-EPMC7513336 | biostudies-literature

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

Project description:Orthogonality constrained density functional theory (OCDFT) [F. A. Evangelista, P. Shushkov and J. C. Tully, J. Phys. Chem. A, 2013, 117, 7378] is a variational time-independent approach for the computation of electronic excited states. In this work we extend OCDFT to compute core-excited states and generalize the original formalism to determine multiple excited states. Benchmark computations on a set of 13 small molecules and 40 excited states show that unshifted OCDFT/B3LYP excitation energies have a mean absolute error of 1.0 eV. Contrary to time-dependent DFT, OCDFT excitation energies for first- and second-row elements are computed with near-uniform accuracy. OCDFT core excitation energies are insensitive to the choice of the functional and the amount of Hartree-Fock exchange. We show that OCDFT is a powerful tool for the assignment of X-ray absorption spectra of large molecules by simulating the gas-phase near-edge spectrum of adenine and thymine.

| S-EPMC4449318 | biostudies-literature

Accurate and Transferable Reactive Molecular Dynamics Models from Constrained Density Functional Theory.

Project description:Chemical reactions constitute the central feature of many liquid, material, and biomolecular processes. Conventional molecular dynamics (MD) is inadequate for simulating chemical reactions given the fixed bonding topology of most force fields, while modeling chemical reactions using ab initio molecular dynamics is limited to shorter time and length scales given its high computational cost. As such, the multiscale reactive molecular dynamics method provides one promising alternative for simulating complex chemical systems at atomistic detail on a reactive potential energy surface. However, the parametrization of such models is a key barrier to their applicability and success. In this work, we present reactive MD models derived from constrained density functional theory that are both accurate and transferable. We illustrate the features of these models for proton dissociation reactions of amino acids in both aqueous and protein environments. Specifically, we present models for ionizable glutamate and lysine that predict accurate absolute pKa values in water as well as their significantly shifted pKa in staphylococcal nuclease (SNase) without any modification of the models. As one outcome of the new methodology, the simulations show that the deprotonation of ionizable residues in SNase can be closely coupled with side chain rotations, which is a concept likely generalizable to many other proteins. Furthermore, the present approach is not limited to only pKa prediction but can enable the fully atomistic simulation of many other reactive systems along with a determination of the key aspects of the reaction mechanisms.

| S-EPMC8480781 | biostudies-literature

Quantitative theory of entropic forces acting on constrained nucleotide sequences applied to viruses.

Project description:We outline a theory to quantify the interplay of entropic and selective forces on nucleotide organization and apply it to the genomes of single-stranded RNA viruses. We quantify these forces as intensive variables that can easily be compared between sequences, outline a computationally efficient transfer-matrix method for their calculation, and apply this method to influenza and HIV viruses. We find viruses altering their dinucleotide motif use under selective forces, with these forces on CpG dinucleotides growing stronger in influenza the longer it replicates in humans. For a subset of genes in the human genome, many involved in antiviral innate immunity, the forces acting on CpG dinucleotides are even greater than the forces observed in viruses, suggesting that both effects are in response to similar selective forces involving the innate immune system. We further find that the dynamics of entropic forces balancing selective forces can be used to predict how long it will take a virus to adapt to a new host, and that it would take H1N1 several centuries to adapt to humans from birds, typically contributing many of its synonymous substitutions to the forcible removal of CpG dinucleotides. By examining the probability landscape of dinucleotide motifs, we predict where motifs are likely to appear using only a single-force parameter and uncover the localization of UpU motifs in HIV. Essentially, we extend the natural language and concepts of statistical physics, such as entropy and conjugated forces, to understanding viral sequences and, more generally, constrained genome evolution.

| S-EPMC3977288 | biostudies-literature

Multilevel Density Functional Theory.

Project description:Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn-Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.

| S-EPMC7880574 | biostudies-literature

Using Constrained Density Functional Theory to Track Proton Transfers and to Sample Their Associated Free Energy Surface.

Project description:Ab initio molecular dynamics (AIMD) and quantum mechanics/molecular mechanics (QM/MM) methods are powerful tools for studying proton solvation, transfer, and transport processes in various environments. However, due to the high computational cost of such methods, achieving sufficient sampling of rare events involving excess proton motion-especially when Grotthuss proton shuttling is involved-usually requires enhanced free energy sampling methods to obtain informative results. Moreover, an appropriate collective variable (CV) that describes the effective position of the net positive charge defect associated with an excess proton is essential both for tracking the trajectory of the defect and for the free energy sampling of the processes associated with the resulting proton transfer and transport. In this work, such a CV is derived from first principles using constrained density functional theory (CDFT). This CV is applicable to a broad array of proton transport and transfer processes as studied via AIMD and QM/MM simulations.

| S-EPMC8444337 | biostudies-literature

Validation of experimental molecular crystal structures with dispersion-corrected density functional theory calculations.

Project description:This paper describes the validation of a dispersion-corrected density functional theory (d-DFT) method for the purpose of assessing the correctness of experimental organic crystal structures and enhancing the information content of purely experimental data. 241 experimental organic crystal structures from the August 2008 issue of Acta Cryst. Section E were energy-minimized in full, including unit-cell parameters. The differences between the experimental and the minimized crystal structures were subjected to statistical analysis. The r.m.s. Cartesian displacement excluding H atoms upon energy minimization with flexible unit-cell parameters is selected as a pertinent indicator of the correctness of a crystal structure. All 241 experimental crystal structures are reproduced very well: the average r.m.s. Cartesian displacement for the 241 crystal structures, including 16 disordered structures, is only 0.095 Å (0.084 Å for the 225 ordered structures). R.m.s. Cartesian displacements above 0.25?A either indicate incorrect experimental crystal structures or reveal interesting structural features such as exceptionally large temperature effects, incorrectly modelled disorder or symmetry breaking H atoms. After validation, the method is applied to nine examples that are known to be ambiguous or subtly incorrect.

| S-EPMC2940256 | biostudies-literature

Density-functional fluctuation theory of crowds.

Project description:A primary goal of collective population behavior studies is to determine the rules governing crowd distributions in order to predict future behaviors in new environments. Current top-down modeling approaches describe, instead of predict, specific emergent behaviors, whereas bottom-up approaches must postulate, instead of directly determine, rules for individual behaviors. Here, we employ classical density functional theory (DFT) to quantify, directly from observations of local crowd density, the rules that predict mass behaviors under new circumstances. To demonstrate our theory-based, data-driven approach, we use a model crowd consisting of walking fruit flies and extract two functions that separately describe spatial and social preferences. The resulting theory accurately predicts experimental fly distributions in new environments and provides quantification of the crowd "mood". Should this approach generalize beyond milling crowds, it may find powerful applications in fields ranging from spatial ecology and active matter to demography and economics.

| S-EPMC6117271 | biostudies-literature

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations With Numerical Atomic Orbitals in HONPAS.

Project description:Linear-scaling density functional theory (DFT) is an efficient method to describe the electronic structures of molecules, semiconductors, and insulators to avoid the high cubic-scaling cost in conventional DFT calculations. Here, we present a parallel implementation of linear-scaling density matrix trace correcting (TC) purification algorithm to solve the Kohn-Sham (KS) equations with the numerical atomic orbitals in the HONPAS package. Such a linear-scaling density matrix purification algorithm is based on the Kohn's nearsightedness principle, resulting in a sparse Hamiltonian matrix with localized basis sets in the DFT calculations. Therefore, sparse matrix multiplication is the most time-consuming step in the density matrix purification algorithm for linear-scaling DFT calculations. We propose to use the MPI_Allgather function for parallel programming to deal with the sparse matrix multiplication within the compressed sparse row (CSR) format, which can scale up to hundreds of processing cores on modern heterogeneous supercomputers. We demonstrate the computational accuracy and efficiency of this parallel density matrix purification algorithm by performing large-scale DFT calculations on boron nitrogen nanotubes containing tens of thousands of atoms.

| S-EPMC7726133 | biostudies-literature

Surface Induced Phenytoin Polymorph. 2. Structure Validation by Comparing Experimental and Density Functional Theory Raman Spectra.

Project description:A method for structure solution in thin films that combines grazing incidence X-ray diffraction data analysis and crystal structure prediction was presented in a recent work (Braun et al. Cryst. Growth Des.2019, DOI: 10.1021/acs.cgd.9b00857). Applied to phenytoin form II, which is only detected in films, the approach gave a very reasonable, but not fully confirmed, candidate structure with Z = 4 and Z' = 2. In the present work, we demonstrate how, by calculating and measuring the crystal Raman spectrum in the low wavenumber energy region with the aim of validating the candidate structure, this can be further refined. In fact, we find it to correspond to a saddle point of the energy landscape of the system, from which a minimum of lower symmetry may be reached. With the new structure, with Z = 4 and Z' = 2, we finally obtain an excellent agreement between experimental and calculated Raman spectra.

| S-EPMC8016182 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data