Removal regarding large Dzyaloshinskii-Moriya conversation from abdominal initio data: First-order s
Author : Mcmahon Merritt | Published On : 15 Nov 2024
We present a relativistic correction scheme to improve the accuracy of 1s core-level binding energies calculated from Green's function theory in the GW approximation, which does not add computational overhead. An element-specific corrective term is derived as the difference between the 1s eigenvalues obtained from the self-consistent solutions to the non- or scalar-relativistic Kohn-Sham equations and the four-component Dirac-Kohn-Sham equations for a free neutral atom. We examine the dependence of this corrective term on the molecular environment and the amount of exact exchange in hybrid exchange-correlation functionals. This corrective term is then added as a perturbation to the quasiparticle energies from partially self-consistent and single-shot GW calculations. We show that this element-specific relativistic correction, when applied to a previously reported benchmark set of 65 core-state excitations [D. Golze et al., J. Phys. Chem. Lett. 11, 1840-1847 (2020)], reduces the mean absolute error (MAE) with respect to the experiment from 0.55 eV to 0.30 eV and eliminates the species dependence of the MAE, which otherwise increases with the atomic number. The relativistic corrections also reduce the species dependence for the optimal amount of exact exchange in the hybrid functional used as a starting point for the single-shot G0W0 calculations. Our correction scheme can be transferred to other methods, which we demonstrate for the delta self-consistent field (ΔSCF) approach based on density functional theory.Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems-a particle solvated in a Lennard-Jones fluid, a Diels-Alder reaction, and crystallization of liquid sodium-to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.A computational expression for the Faraday A term of magnetic circular dichroism (MCD) is derived within coupled cluster response theory and alternative computational expressions for the B term are discussed. Moreover, an approach to compute the (temperature-independent) MCD ellipticity in the context of coupled cluster damped response is presented, and its equivalence with the stick-spectrum approach in the limit of infinite lifetimes is demonstrated. The damped response approach has advantages for molecular systems or spectral ranges with a high density of states. click here Illustrative results are reported at the coupled cluster singles and doubles level and compared to time-dependent density functional theory results.The radiative association (RA) rate constant is computed for the formation of the diatomic sodium chloride (NaCl) molecule in the temperature interval 1 K-30 K. At these temperatures, RA of NaCl through non-adiabatic dynamics is important. A scattering program has been implemented to carry out calculations of RA cross sections, accounting for coupled dynamics on the lowest ionic and the lowest neutral diabatic 1Σ+ states. The study shows that the non-adiabatic treatment gives a cross section that exceeds that of conventional adiabatic dynamics by one to four orders of magnitude. The contribution to the RA rate constant from Na and Cl approaching each other in the A1Π state has also been computed using an established quantum mechanical method. Ab initio data from the literature have been used for the potential energy curves, the diabatic coupling, and the electric dipole moments of NaCl.In this paper, we present an overview of crystal imperfections in ice Ih. Due to its molecular nature, the fundamental asymmetry of the hydrogen bond, and proton disorder, crystal defects in this condensed form of water reveal a complexity not usually seen in atomic crystalline solids. The discussion is organized in terms of the spatial extent of the defects. We start with zero-dimensional imperfections such as the molecular vacancy and interstitial, Bjerrum, and ionic defects, as well as possible defect complexes that can be formed from them. Subsequently, we turn to the properties of dislocations, which are the one-dimensional disturbances that carry plastic deformation in crystalline solids. Finally, we discuss two-dimensional defects such as stacking faults and grain boundaries and discuss to what extent the latter are similar to other interfaces in ice Ih such as the free surface. We conclude with an outlook at the road ahead, discussing future challenges toward understanding the role of crystal defects in the macroscopic behavior of ice Ih.