Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExplaining and Fixing DFT Failures for Torsional Barriers
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Most torsional barriers are predicted to high accuracy (about 1kJ/mol) by standard semilocal functionals, but a small subset has been found to have much larger errors. We create a database of almost 300 carbon-carbon torsional barriers, including 12 poorly behaved barriers, all stemming from Y=C-X group, where X is O or S, and Y is a halide. Functionals with enhanced exchange mixing (about 50%) work well for all barriers. We find that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.
rate research
Read More
We present a simple torsional potential for graphene to accurately describe its out-of-plane deformations. The parameters of the potential are derived through appropriate fitting with suitable DFT calculations regarding the deformation energy of graphene sheets folded around two different folding axes, along an armchair or along a zig-zag direction. Removing the energetic contribution of bending angles, using a previously introduced angle bending potential, we isolate the purely torsional deformation energy, which is then fitted to simple torsional force fields. The presented out-of-plane torsional potential can accurately fit the deformation energy for relatively large torsional angles up to 0.5 rad. To test our proposed potential, we apply it to the problem of the vertical displacement of a single carbon atom out of the graphene plane and compare the obtained deformation energy with corresponding DFT calculations. The dependence of the deformation energy on the vertical displacement of the pulled carbon atom is indistinguishable in these two cases, for displacements up to about 0.5 $AA$. The presented potential is applicable to other sp$^2$ carbon structures.
We analyze theoretically the interplay between the torsional and the rotational motion of an aligned biphenyl-like molecule. To do so, we consider a transition between two electronic states with different internal torsional potentials, induced by means of a resonant laser pulse. The change in the internal torsional potential provokes the motion of the torsional wavepacket in the excited electronic state, modifying the structure of the molecule, and hence, its inertia tensor. We find that this process has a strong impact on the rotational wave function, displaying different behavior depending on the electronic states involved and their associated torsional potentials. We describe the dynamics of the system by considering the degree of alignment and the expectations values of the angular momentum operators for the overall rotation of the molecule.
Projection-based embedding provides a simple, robust, and accurate approach for describing a small part of a chemical system at the level of a correlated wavefunction method while the remainder of the system is described at the level of density functional theory. Here, we present the derivation, implementation, and numerical demonstration of analytical nuclear gradients for projection-based wavefunction-in-density functional theory (WF-in-DFT) embedding. The gradients are formulated in the Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. An important aspect of the gradient theory is that WF contributions to the total WF-in-DFT gradient can be simply evaluated using existing WF gradient implementations without modification. Another simplifying aspect is that Kohn-Sham (KS) DFT contributions to the projection-based embedding gradient do not require knowledge of the WF calculation beyond the relaxed WF density. Projection-based WF-in-DFT embedding gradients are thus easily generalized to any combination of WF and KS-DFT methods. We provide numerical demonstration of the method for several applications, including calculation of a minimum energy pathway for a hydride transfer in a cobalt-based molecular catalyst using the nudged-elastic-band method at the CCSD-in-DFT level of theory, which reveals large differences from the transition state geometry predicted using DFT.
Using a set of oscillator strengths and excited-state dipole moments of near full configuration interaction (FCI) quality determined for small compounds, we benchmark the performances of several single-reference wave function methods (CC2, CCSD, CC3, CCSDT, ADC(2), and ADC(3/2)) and time-dependent density-functional theory (TD-DFT) with various functionals (B3LYP, PBE0, M06-2X, CAM-B3LYP, and $omega$B97X-D). We consider the impact of various gauges (length, velocity, and mixed) and formalisms: equation of motion (EOM) emph{vs} linear response (LR), relaxed emph{vs} unrelaxed orbitals, etc. Beyond the expected accuracy improvements and a neat decrease of formalism sensitivy when using higher-order wave function methods, the present contribution shows that, for both ADC(2) and CC2, the choice of gauge impacts more significantly the magnitude of the oscillator strengths than the choice of formalism, and that CCSD yields a notable improvement on this transition property as compared to CC2. For the excited-state dipole moments, switching on orbital relaxation appreciably improves the accuracy of both ADC(2) and CC2, but has a rather small effect at the CCSD level. Going from ground to excited states, the typical errors on dipole moments for a given method tend to roughly triple. Interestingly, the ADC(3/2) oscillator strengths and dipoles are significantly more accurate than their ADC(2) counterparts, whereas the two models do deliver rather similar absolute errors for transition energies. Concerning TD-DFT, one finds: i) a rather negligible impact of the gauge on oscillator strengths for all tested functionals (except for M06-2X); ii) deviations of ca.~0.10 D on ground-state dipoles for all functionals; iii) the better overall performance of CAM-B3LYP for the two considered excited-state properties.
Bottom-up coarse-grained molecular dynamics models are parameterized using complex effective Hamiltonians. These models are typically optimized to approximate high dimensional data from atomistic simulations. In contrast, human validation of these models is often limited to low dimensional statistics that do not necessarily differentiate between the CG model and said atomistic simulations. We propose that explainable machine learning can directly convey high-dimensional error to scientists and use Shapley additive explanations do so in two coarse-grained protein models.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
