A new ab initio potential energy surface for the NeCO complex with the vibrational coordinate dependence.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang Z
PROVIDER: S-EPMC3036545 | biostudies-other | 2011 Jan
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We built a full-dimensional analytical potential energy surface of the ground electronic state of Li₂H from ca. 20,000 ab initio multi-reference configuration interaction calculations, including core⁻valence correlation effects. The surface is flexible enough to accurately describe the three dissociation channels: Li (2s ²S) + LiH (¹Σ⁺), Li₂ (¹Σg⁺) + H (1s ²S) and 2Li (2s ²S) + H (1s ²S). Using a local fit of this surface, we calculated pure (J = 0) vibrational states of Li₂H up to the barrier to linearity (ca. 3400 cm-1 above the global minimum) using a vibrational self-consistent field/virtual state configuration interaction method. We found 18 vibrational states below this barrier, with a maximum of 6 quanta in the bending mode, which indicates that Li₂H could be spectroscopically observable. Moreover, we show that some of these vibrational states are highly correlated already ca. 1000 cm-1 below the height of the barrier. We hope these calculations can help the assignment of experimental spectra. In addition, the first low-lying excited states of each B₁, B₂ and A₂ symmetry of Li₂H were characterized.
| S-EPMC6337435 | biostudies-literature
Project description:The rates of intramolecular proton transfer are calculated on a full-dimensional reactive electronic potential energy surface that incorporates high level ab initio calculations along the reaction path and by using classical Transition State theory, Path-Integral Quantum Transition State Theory, and the Quantum Instanton approach. The specific example problem studied is malonaldehyde. Estimates of the kinetic isotope effect using the latter two methods are found to be in reasonable agreement with each other. Improvements and extensions of this practical, yet chemically accurate framework for the calculations of quantized, reactive dynamics are also discussed.
| S-EPMC2992356 | biostudies-literature
Project description:We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50?000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.
| S-EPMC5472367 | biostudies-literature
Project description:The routine prediction of three-dimensional protein structure from sequence remains a challenge in computational biochemistry. It has been intuited that calculated energies from physics-based scoring functions are able to distinguish native from nonnative folds based on previous performance with small proteins and that conformational sampling is the fundamental bottleneck to successful folding. We demonstrate that as protein size increases, errors in the computed energies become a significant problem. We show, by using error probability density functions, that physics-based scores contain significant systematic and random errors relative to accurate reference energies. These errors propagate throughout an entire protein and distort its energy landscape to such an extent that modern scoring functions should have little chance of success in finding the free energy minima of large proteins. Nonetheless, by understanding errors in physics-based score functions, they can be reduced in a post-hoc manner, improving accuracy in energy computation and fold discrimination.
| S-EPMC3081830 | biostudies-literature
Project description:Accurate prediction of the placement and comformations of protein side chains given only the backbone trace has a wide range of uses in protein design, structure prediction, and functional analysis. Prediction has most often relied on discrete rotamer libraries so that rapid fitness of side-chain rotamers can be assessed against some scoring function. Scoring functions are generally based on experimental parameters from small-molecule studies or empirical parameters based on determined protein structures. Here, we describe the NCN algorithm for predicting the placement of side chains. A predominantly first-principles approach was taken to develop the potential energy function incorporating van der Waals and electrostatics based on the OPLS parameters, and a hydrogen bonding term. The only empirical knowledge used is the frequency of rotameric states from the PDB. The rotamer library includes nearly 50,000 rotamers, and is the most extensive discrete library used to date. Although the computational time tends to be longer than most other algorithms, the overall accuracy exceeds all algorithms in the literature when placing rotamers on an accurate backbone trace. Considering only the most buried residues, 80% of the total residues tested, the placement accuracy reaches 92% for chi(1), and 83% for chi(1 + 2), and an overall RMS deviation of 1 A. Additionally, we show that if information is available to restrict chi(1) to one rotamer well, then this algorithm can generate structures with an average RMS deviation of 1.0 A for all heavy side-chains atoms and a corresponding overall chi(1 + 2) accuracy of 85.0%.
| S-EPMC2286725 | biostudies-literature
Project description:In "Lattice dynamics and structure of the new langasites Ln3CrGe3Be2O14 (Ln = La, Pr, Nd): vibrational spectra and ab initio calculations" [1], experimental and calculated results on lattice dynamics of the recently discovered new compounds La3CrGe3Be2O14, Pr3CrGe3Be2O14, and Nd3CrGe3Be2O14 are reported. These compounds belong to the langasite series and constitute a new class of low-dimensional antiferromagnets. The data presented in this article includes IR diffuse transmission spectra of powder samples of Ln3CrGe3Be2O14 (Ln = La, Pr, Nd) registered at room temperature with a Bruker 125HR Fourier spectrometer, Raman spectra taken in the backscattering geometry (also at room temperature) with a triple monochromator using the line 514, 5 nm of an argon laser as an excitation, results of the DFT calculations with the B3LYP and PBE0 hybrid functionals on the optimized crystal structures, eigenfrequencies and eigenvectors of the normal vibrational modes. These data can be used to analyse electron-phonon interaction and multiferroic properties of the new langasites and to compare the lattice dynamics of different langasites. The dataset is available on mendeley data public repository at https://doi.org/10.17632/32grbb4p82.1.
| S-EPMC6939064 | biostudies-literature
Project description:Ab initio intermolecular potential energy surfaces (PES) of N2-NO have been constructed at the level of theory CCSD(T) with the augmented correlation-consistent basis sets aug-cc-pVmZ (with m = 2, 3, 4). The nitrogen in the closed-shell electronic configuration X1?+ and nitric oxide in the open-shell electronic configuration A2?+ were employed to calculate ab initio intermolecular interaction energies. The two new ab initio 5-site intermolecular pair potentials at the theoretical level CCSD(T)/aug-cc-pVmZ (with m = 4, 24) were developed appropriately and are suitable for N2-NO dimer by using the nonlinear least-squares fitting method combining MIO and Levenberg-Marquardt algorithms. The correlation quality of these two potentials was found to be very good with R 2 values in the range of 0.98372 to 0.99775. The cross second virial coefficients B 12(T) of the N2-NO dimer were calculated in the temperature range of 100 to 470 K using the two ab initio 5-site potentials. The discrepancies between the calculated results and the experimental data can be acceptable.
| S-EPMC7271377 | biostudies-literature
Project description:This work reports on the theoretical equilibrium crystal shapes of GaAs and InAs as a function of temperature and pressure, taking into account the contribution of the surface vibration, using ab-initio thermodynamic calculations. For this purpose, new (111)B reconstructions, which are energetically stable at a high temperature, are suggested. It was found that there was a feasible correspondence between the calculated equilibrium shapes and the experimental shapes, which implied that the previous experimental growth was performed under conditions that were close to equilibrium. In this study, GaAs and InAs were selected as prototype compound semiconductors, but the developed calculation methodology can also be applied to other III-V compound semiconductor materials.
| S-EPMC6361998 | biostudies-literature
Project description:Although the thermoelectric figure of merit zT above 300 K has seen significant improvement recently, the progress at lower temperatures has been slow, mainly limited by the relatively low Seebeck coefficient and high thermal conductivity. Here we report, for the first time to our knowledge, success in first-principles computation of the phonon drag effect--a coupling phenomenon between electrons and nonequilibrium phonons--in heavily doped region and its optimization to enhance the Seebeck coefficient while reducing the phonon thermal conductivity by nanostructuring. Our simulation quantitatively identifies the major phonons contributing to the phonon drag, which are spectrally distinct from those carrying heat, and further reveals that although the phonon drag is reduced in heavily doped samples, a significant contribution to Seebeck coefficient still exists. An ideal phonon filter is proposed to enhance zT of silicon at room temperature by a factor of 20 to ? 0.25, and the enhancement can reach 70 times at 100 K. This work opens up a new venue toward better thermoelectrics by harnessing nonequilibrium phonons.
| S-EPMC4672777 | biostudies-literature
Project description:We have developed a new combined approach for ab initio protein structure prediction. The protein conformation is described as a lattice chain connecting C(alpha) atoms, with attached C(beta) atoms and side-chain centers of mass. The model force field includes various short-range and long-range knowledge-based potentials derived from a statistical analysis of the regularities of protein structures. The combination of these energy terms is optimized through the maximization of correlation for 30 x 60,000 decoys between the root mean square deviation (RMSD) to native and energies, as well as the energy gap between native and the decoy ensemble. To accelerate the conformational search, a newly developed parallel hyperbolic sampling algorithm with a composite movement set is used in the Monte Carlo simulation processes. We exploit this strategy to successfully fold 41/100 small proteins (36 approximately 120 residues) with predicted structures having a RMSD from native below 6.5 A in the top five cluster centroids. To fold larger-size proteins as well as to improve the folding yield of small proteins, we incorporate into the basic force field side-chain contact predictions from our threading program PROSPECTOR where homologous proteins were excluded from the data base. With these threading-based restraints, the program can fold 83/125 test proteins (36 approximately 174 residues) with structures having a RMSD to native below 6.5 A in the top five cluster centroids. This shows the significant improvement of folding by using predicted tertiary restraints, especially when the accuracy of side-chain contact prediction is >20%. For native fold selection, we introduce quantities dependent on the cluster density and the combination of energy and free energy, which show a higher discriminative power to select the native structure than the previously used cluster energy or cluster size, and which can be used in native structure identification in blind simulations. These procedures are readily automated and are being implemented on a genomic scale.
| S-EPMC1303233 | biostudies-literature