Efficient Generation of Permutationally Invariant Potential Energy Surfaces for Large Molecules.
Ontology highlight
ABSTRACT:
SUBMITTER: Conte R
PROVIDER: S-EPMC7997398 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The electronic energy in the Hartree-Fock (HF) theory is the trace of the product of the charge density matrix (CDM) with the one-electron and two-electron matrices represented in an atomic orbital basis, where the two-electron matrix is also a function of the same CDM. In this work, we examine a formalism of analytic representation of a generic molecular potential energy surface (PES) as a sum of a linearly parameterized HF and a correction term, the latter formally representing the electron correlation energy, also linearly parameterized, by expressing the elements of CDM using permutationally invariant polynomials (PIPs). We show on a variety of numerical examples, ranging from exemplary two-electron systems HeH+ and H3+ to the more challenging cases of methanium (CH5+) fragmentation and high-energy tautomerization of formamide to formimidic acid that such a formulation requires significantly fewer, 10-20% of PIPs, to accomplish the same accuracy of the fit as the conventional representation at practically the same computational cost.
| S-EPMC10501011 | biostudies-literature
Project description:Accurate conformational energetics of molecules are of great significance to understand maby chemical properties. They are also fundamental for high-quality parameterization of force fields. Traditionally, accurate conformational profiles are obtained with density functional theory (DFT) methods. However, obtaining a reliable energy profile can be time-consuming when the molecular sizes are relatively large or when there are many molecules of interest. Furthermore, incorporation of data-driven deep learning methods into force field development has great requirements for high-quality geometry and energy data. To this end, we compared several possible alternatives to the traditional DFT methods for conformational scans, including the semi-empirical method GFN2-xTB and the neural network potential ANI-2x. It was found that a sequential protocol of geometry optimization with the semi-empirical method and single-point energy calculation with high-level DFT methods can provide satisfactory conformational energy profiles hundreds of times faster in terms of optimization.
| S-EPMC9738817 | biostudies-literature
Project description:Despite the continuous development of theoretical methodologies for describing nonadiabatic dynamics of molecular systems, there is a lack of approaches for processes where the norm of the wave function is not conserved, i.e., when an imaginary potential accounts for some irreversible decaying mechanism. Current approaches rely on building potential energy surfaces of reduced dimensionality, which is not optimal for more involving and realistic multidimensional problems. Here, we present a novel methodology for describing the dynamics of complex-valued molecular Hamiltonians, which is a generalisation of the trajectory surface hopping method. As a first application, the complex surface fewest switches surface hopping (CS-FSSH) method was employed to survey the relaxation mechanisms of the shape resonant anions of iodoethene. We have provided the first detailed and dynamical picture of the π*/σ* mechanism of dissociative electron attachment in halogenated unsaturated compounds, which is believed to underlie electron-induced reactions of several molecules of interest. Electron capture into the π* orbital promotes C[double bond, length as m-dash]C stretching and out-of-plane vibrations, followed by charge transfer from the double bond into the σ* orbital at the C-I bond, and, finally, release of the iodine ion, all within only 15 fs. On-the-fly dynamics simulations of a vast class of processes can be envisioned with the CS-FSSH methodology, including autoionisation from transient anions, core-ionised and superexcited states, Auger and interatomic coulombic decay, and time-dependent luminescence.
| S-EPMC8162122 | biostudies-literature
Project description:Energy migrations at metal nanomaterial surfaces are fundamentally important to heterogeneous reactions. Here we report two distinctive energy migration pathways of monolayer adsorbate molecules on differently sized metal nanoparticle surfaces investigated with ultrafast vibrational spectroscopy. On a 5?nm platinum particle, within a few picoseconds the vibrational energy of a carbon monoxide adsorbate rapidly dissipates into the particle through electron/hole pair excitations, generating heat that quickly migrates on surface. In contrast, the lack of vibration-electron coupling on approximately 1?nm particles results in vibrational energy migration among adsorbates that occurs on a twenty times slower timescale. Further investigations reveal that the rapid carbon monoxide energy relaxation is also affected by the adsorption sites and the nature of the metal but to a lesser extent. These findings reflect the dependence of electron/vibration coupling on the metallic nature, size and surface site of nanoparticles and its significance in mediating energy relaxations and migrations on nanoparticle surfaces.
| S-EPMC4757789 | biostudies-literature
Project description:Adversarial generative models are becoming an essential tool in molecular design and discovery due to their efficiency in exploring the desired chemical space with the assistance of deep learning. In this article, we introduce an integrated framework by combining the modules of algorithmic synthesis, deep prediction, adversarial generation, and fine screening for the purpose of effective design of the thermally activated delayed fluorescence (TADF) molecules that can be used in the organic light-emitting diode devices. The retrosynthetic rules are employed to algorithmically synthesize the D-A complex based on the empirically defined donor and acceptor moieties, which is followed by the high-throughput labeling and prediction with the deep neural network. The new D-A molecules are subsequently generated via the adversarial autoencoder, with the excited-state property distributions perfectly matching those of the original samples. Fine screening of the generated molecules, including the spin-orbital coupling calculation and the excited-state optimization, is eventually implemented to select the qualified TADF candidates within the novel chemical space. Further investigation shows that the created structures fully mimic the original D-A samples by maintaining a significant charge transfer characteristic, a minimal adiabatic singlet-triplet gap, and a moderate spin-orbital coupling that are desirable for the delayed fluorescence.
| S-EPMC9161419 | biostudies-literature
Project description:In macromolecular design, conformational energies are sensitive to small changes in atom coordinates; thus, modeling the small, continuous motions of atoms around low-energy wells confers a substantial advantage in structural accuracy. However, modeling these motions comes at the cost of a very large number of energy function calls, which form the bottleneck in the design calculations. In this work, we remove this bottleneck by consolidating all conformational energy evaluations into the pre-computation of a local polynomial expansion of the energy about the "ideal" conformation for each low-energy, "rotameric" state of each residue pair. This expansion is called "energy as polynomials in internal coordinates" (EPIC), where the internal coordinates can be side-chain dihedrals, backrub angles, and/or any other continuous degrees of freedom of a macromolecule, and any energy function can be used without adding any asymptotic complexity to the design. We demonstrate that EPIC efficiently represents the energy surface for both molecular-mechanics and quantum-mechanical energy functions, and apply it specifically to protein design for modeling both side chain and backbone degrees of freedom.
| S-EPMC4469035 | biostudies-literature
Project description:The structure and dynamics of a molecular system is governed by its potential energy surface (PES), representing the total energy as a function of the nuclear coordinates. Obtaining accurate potential energy surfaces is limited by the exponential scaling of Hilbert space, restricting quantitative predictions of experimental observables from first principles to small molecules with just a few electrons. Here, we present an explicitly physics-informed approach for improving and assessing the quality of families of PESs by modifying them through linear coordinate transformations based on experimental data. We demonstrate this "morphing" of the PES for the He - H2+ complex using recent comprehensive Feshbach resonance (FR) measurements for reference PESs at three different levels of quantum chemistry. In all cases, the positions and intensities of peaks in the energy distributions are improved. We find these observables to be mainly sensitive to the long-range part of the PES.
| S-EPMC10906917 | biostudies-literature
Project description:Photophysical and photochemical processes are ruled by the interplay between transient vibrational and electronic degrees of freedom, which are ultimately determined by the multidimensional potential energy surfaces (PESs). Differences between ground and excited PESs are encoded in the relative intensities of resonant Raman bands, but they are experimentally challenging to access, requiring measurements at multiple wavelengths under identical conditions. Herein, we perform a two-color impulsive vibrational scattering experiment to launch nuclear wavepacket motions by an impulsive pump and record their coupling with a targeted excited-state potential by resonant Raman processes with a delayed probe, generating in a single measurement background-free vibrational spectra across the entire sample absorption. Building on the interference between the multiple pathways resonant with the excited-state manifold that generate the Raman signal, we show how to experimentally tune their relative phase by varying the probe chirp, decoding nuclear displacements along different normal modes and revealing the multidimensional PESs. Our results are validated against time-dependent density functional theory.
| S-EPMC8488957 | biostudies-literature
Project description:We present Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF), an automated methodology for the generation of global potential energy surfaces (PES), molecular properties surfaces, e.g., dipole, polarizabilities, etc. using a single random geometry as input. The SRP-MGPF workflow integrates: (i) a fully automated procedure for the global topographical characterization of a (intermolecular) PES based on the Transition State Search Using Chemical Dynamical Simulations (TSSCDS) family of methods;i (ii) the global optimization of the parameters of a semiempirical Hamiltonian in order to reproduce a given level of electronic structure theory; and (iii) a tensor decomposition algorithm which turns the resulting SRP-PES into sum of products (Tucker) form with the Multigrid POTFIT algorithm. The latter is necessary for quantum dynamical studies within the Multiconfiguration Time-Dependent Hartree (MCTDH) quantum dynamics method. To demonstrate our approach, we have applied our methodology to the cis-trans isomerization reaction in HONO in full dimensionality (6D). The resulting SRP-PES has been validated through the computation of classical on-the-fly dynamical calculations as well as calculations of the lowest vibrational eigenstates of HONO as well as high-energy wavepacket propagations.
| S-EPMC6702682 | biostudies-literature
Project description:We develop a rapid and efficient method for the comparison of protein local surface similarities using geometric invariants (fingerprints). By combining fast fingerprint comparison with explicit alignment, we successfully screen the entire Protein Data Bank for proteins that possess local surface similarities. Our method is independent of sequence and fold similarities, and has potential application to protein structure annotation and protein-protein interface design.
| S-EPMC2757799 | biostudies-literature