Orbital Correction Method Research Articles

Bond length alternation of π-conjugated polymers predicted by the Fermi-Löwdin orbital self-interaction correction method.

π-conjugated polymers have been used in a wide range of practical applications, partly due to their unique properties that originate in the delocalization of electrons through the polymer backbone. The level of delocalization can be characterized by the induced bond length alternation (BLA), with shorter BLA connected with strong delocalization and vice versa. The accurate description of this structural parameter can be considered a benchmark for testing the capability of different electronic structure methods for self-interaction error (SIE) removal and electron correlation inclusion. Density functional theory (DFT), in its local or semi-local flavors, suffers from SIE and, thus, underestimates the BLA compared to self-interaction-free methods. In this work, we utilize the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for one-electron self-interaction removal to characterize the BLA of five oligomers with increasing length extrapolated to the polymeric limit. We compare the self-interaction-free BLA to several DFT approximations, Møller-Plesset second-order perturbation theory (MP2), and the BLA obtained with the domain based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] approximation. Our findings show that FLOSIC corrects for the small BLA given by (semi-)local DFT approximations, but it tends to overcorrect with respect to CAM-B3LYP, MP2, and DLPNO-CCSD(T).

Study the closed orbit correction of a compact hadron driver for cancer therapy

A design of a compact hadron driver for future cancer therapies based on the induction synchrotron concept is given. In order to realize a slow extraction technique in a fast cycling synchrotron, which allows the energy sweep beam scanning, the zero momentum-dispersion D(s) region and high flat D(s) region are necessary. The present design meets both requirements. The lattice has the two-fold symmetry with a circumference of 52.8 m, 2 m-long dispersion-free straight section, and 3 m-long large flat dispersion straight section. Assuming a 1.5 T bending magnet, the ring can deliver heavy ions of 200 MeV/au at 10 Hz. Details of the lattice parameters and orbit correction approach and method are discussed.

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi-Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi-Löwdin orbitals which are used to remove the self-interaction error in the Perdew-Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi-Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn-Sham with one additional term in the Kohn-Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn-Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. We analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville-von Neumann equation to calculate excitation energies.

How Good Is the Density-Corrected SCAN Functional for Neutral and Ionic Aqueous Systems, and What Is So Right about the Hartree-Fock Density?

Density functional theory (DFT) is the most widely used electronic structure method, due to its simplicity and cost effectiveness. The accuracy of a DFT calculation depends not only on the choice of the density functional approximation (DFA) adopted but also on the electron density produced by the DFA. SCAN is a modern functional that satisfies all known constraints for meta-GGA functionals. The density-driven errors, defined as energy errors arising from errors of the self-consistent DFA electron density, can hinder SCAN from achieving chemical accuracy in some systems, including water. Density-corrected DFT (DC-DFT) can alleviate this shortcoming by adopting a more accurate electron density which, in most applications, is the electron density obtained at the Hartree-Fock level of theory due to its relatively low computational cost. In this work, we present extensive calculations aimed at determining the accuracy of the DC-SCAN functional for various aqueous systems. DC-SCAN (SCAN@HF) shows remarkable consistency in reproducing reference data obtained at the coupled cluster level of theory, with minimal loss of accuracy. Density-driven errors in the description of ionic aqueous clusters are thoroughly investigated. By comparison with the orbital-optimized CCD density in the water dimer, we find that the self-consistent SCAN density transfers a spurious fraction of an electron across the hydrogen bond to the hydrogen atom (H*, covalently bound to the donor oxygen atom) from the acceptor (OA) and donor (OD) oxygen atoms, while HF makes a much smaller spurious transfer in the opposite direction, consistent with DC-SCAN (SCAN@HF) reduction of SCAN overbinding due to delocalization error. While LDA seems to be the conventional extreme of density delocalization error, and HF the conventional extreme of (usually much smaller) density localization error, these two densities do not quite yield the conventional range of density-driven error in energy differences. Finally, comparisons of the DC-SCAN results with those obtained with the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method show that DC-SCAN represents a more accurate approach to reducing density-driven errors in SCAN calculations of ionic aqueous clusters. While the HF density is superior to that of SCAN for noncompact water clusters, the opposite is true for the compact water molecule with exactly 10 electrons.

Describing Chemical Reactivity with Frontier Molecular Orbitalets.

Locality in physical space is critical in understanding chemical reactivity in the analysis of various phenomena and processes in chemistry, biology, and materials science, as exemplified in the concepts of reactive functional groups and active sites. Frontier molecular orbitals (FMOs) pinpoint the locality of chemical bonds that are chemically reactive because of the associated orbital energies and thus have achieved great success in describing chemical reactivity, mainly for small systems. For large systems, however, the delocalization nature of canonical molecular orbitals makes it difficult for FMOs to highlight the locality of the chemical reactivity. To obtain localized molecular orbitals that also reflect the frontier nature of the chemical processes, we develop the concept of frontier molecular orbitalets (FMOLs) for describing the reactivity of large systems. The concept of orbitalets was developed recently in the localized orbital scaling correction method, which aims for eliminating the delocalization error in common density functional approximations. Orbitalets are localized in both physical and energy spaces and thus contain both orbital locality and energy information. The FMOLs are thus the orbitalets with energies highest among occupied orbitalets and lowest among unoccupied ones. The applications of FMOLs to hexadeca-1,3,5,7,9,11,13,15-octaene in its equilibrium geometry, inter- and intra-molecular charge-transfer systems, and two transition states of a bifurcating reaction demonstrate that FMOLs can connect quantum mechanical treatments of chemical systems and chemical reactivities by locating the reactive region of large chemical systems. Therefore, FMOLs extend the role of FMOs for small systems and describe the chemical reactivity of large systems with energy and locality insight, with potentially broad applications.

Surface Displacement and Source Model Separation of the Two Strongest Earthquakes During the 2019 Ridgecrest Sequence: Insights From InSAR, GPS, and Optical Data

AbstractSince the occurrence of the 2019 Mw 6.4 and Mw 7.1 Ridgecrest earthquake sequence in the Eastern California Shear Zone, coseismic deformation following the two earthquakes has been intensively studied, and source modeling with interferometric synthetic aperture radar (InSAR), Global Positioning System (GPS) and seismological datasets has been favored. However, we recently found that the coseismic modeling of the two earthquakes constrained by the dense near‐field Planet‐Lab optical measurements can be more detailed than the previously estimated, which requires an accurate Planet‐Lab displacement. In this study, we improve the long wavelength orbital ramps correction method for obtaining a more accurate Planet‐Lab horizontal displacement field. The corrected dense near‐field Planet‐Lab data are firstly used together with the intermediate‐field GPS data to invert and distinguish the fault slip distribution for the two earthquakes. The same scheme is used to simultaneously invert the InSAR, SAR, GPS, and optical datasets, to improve the constraints on seismic source parameters. Our inversion results show that the joint‐event slip model is rougher than the two single‐event slip models, but it has a more concentrated slip pattern and larger slip amplitude in some zones. We show that adding a near‐field constraint of the Planet‐Lab data in the combined‐data inversion can reduce the slip parameter uncertainty and enhance the model resolution. The Coulomb failure stress changes on the southeastern Blackwater, the southern Owens Valley, and the central Panamint Valley faults are enhanced by about 0.4–0.8 bar by the 2019 Ridgecrest earthquake sequence.

LibSC: Library for Scaling Correction Methods in Density Functional Theory.

In recent years, a series of scaling correction (SC) methods have been developed in the Yang laboratory to reduce and eliminate the delocalization error, which is an intrinsic and systematic error existing in conventional density functional approximations (DFAs) within density functional theory (DFT). On the basis of extensive numerical results, the SC methods have been demonstrated to be capable of reducing the delocalization error effectively and producing accurate descriptions for many critical and challenging problems, including the fundamental gap, photoemission spectroscopy, charge transfer excitations, and polarizability. In the development of SC methods, the SC methods were mainly implemented in the QM4D package that was developed in the Yang laboratory for research development. The heavy dependency on the QM4D package hinders the SC methods from access by researchers for broad applications. In this work, we developed a reliable and efficient implementation, LibSC, for the global scaling correction (GSC) method and the localized orbital scaling correction (LOSC) method. LibSC will serve as a lightweight and open-source library that can be easily accessed by the quantum chemistry community. The implementation of LibSC is carefully modularized to provide the essential functionalities for conducting calculations of the SC methods. In addition, LibSC provides simple and consistent interfaces to support multiple popular programing languages, including C, C++, and Python. In addition to the development of the library, we also integrated LibSC with two popular and open-source quantum chemistry packages, the Psi4 package and the PySCF package, which provides immediate access for general users to perform calculations with SC methods.

Application of machine learning in orbital correction of storage ring

Synchrotron light source is one of the most powerful tools in modern science and technology. Shanghai Synchrotron Radiation Facility (SSRF), located in Shanghai, China, is an advanced 3.5 GeV 3rd-generation medium energy light source. The 3rd-generation synchrotron radiation light source will provide high brilliance and high stability synchrotron radiation to fulfill the advanced experimental conditions in frontier researches. To achieve highly stable radiation, it is important to have highly stable beam orbit. Thus we adopted machine learning method to control and feedback the orbit. Using this neural network-based orbit correction method, which doesn’t rely on the response matrix, we can establish a nonlinear mapping relationship between correctors and the orbit distortions and perform continuous online retraining. This new method can significantly improve the orbit stability of SSRF.

Self-interaction correction in water-ion clusters.

We study the importance of self-interaction errors in density functional approximations for various water-ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with the local spin-density approximation, Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriately normed (SCAN) meta-GGA to describe binding energies of hydrogen-bonded water-ion clusters, i.e., water-hydronium, water-hydroxide, water-halide, and non-hydrogen-bonded water-alkali clusters. In the hydrogen-bonded water-ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water-ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher-level calculations. The non-hydrogen-bonded water-alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤50 meV/H2O) for hydrogen-bonded water-ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

NOAA-AVHRR Orbital Drift Correction: Validating Methods Using MSG-SEVIRI Data as a Benchmark Dataset

National Oceanic and Atmospheric Administration–Advanced Very High Resolution Radiometer (NOAA-AVHRR) data provides the possibility to build the longest Land Surface Temperature (LST) dataset to date, starting in 1981 up to the present. However, due to the orbital drift of the NOAA platforms, no LST dataset is available before 2000 and the arrival of newer platforms. Although numerous methods have been developed to correct this orbital drift effect on the LST, a lack of validation has prevented their application. This is the gap we bridge here by using the 15 min temporal resolution of Meteosat Second Generation–Spinning Enhanced Visible and Infra-Red Imager (MSG-SEVIRI) data to simulate drifted and reference LST time series. We then use these time series to validate an orbital drift correction method based on solar zenith angle (SZA) anomalies that we presented in a previous work (C1), as well as two variations of this approach (C0 and C2). Our results show that the C0 method performs better than the two others, although its overall bias absolute value ranges up to 1 K, while standard deviation values remain around 3 K. This is verified for most land covers, for all NOAA platforms, and these statistics remain mostly stable with noise on SZA time series (from 0° to ±10°). With this study, we show that orbital drift correction methods can be thoroughly validated and that such validation should aim toward bias absolute values below 0.1 K and standard deviation values around 1.4 K at coarse spatial resolution. To validate other orbital drift correction approaches, the drifted and reference time series used in this work are freely available for download from the first author’s webpage. This will be the first step toward the building of an orbital-drift-corrected long-term LST dataset.

Self-Consistent Calculation of the Localized Orbital Scaling Correction for Correct Electron Densities and Energy-Level Alignments in Density Functional Theory.

The recently developed localized orbital scaling correction (LOSC) method shows the ability to systematically and size-consistently reduce the delocalization error existing in conventional density functional approximations (DFAs). However, the application of LOSC to DFAs was mainly through a post self-consistent field (SCF) manner, and few results from applying LOSC to DFAs in an SCF manner have been reported. The reason is that the originally proposed SCF approach to SCF-LOSC calculation uses an approximate Hamiltonian and encounters convergence problems easily in practice. In this work, we develop a new SCF approach with a correct Hamiltonian and achieve reliable SCF-LOSC calculations. We demonstrate the capability of the new SCF approach for SCF-LOSC to correctly describe the electron densities, total energies, and energy-level alignment for the molecular dissociation process, while conventional DFAs or LOSC-DFAs with post-SCF calculations show large errors. This work demonstrates that the new SCF approach for SCF-LOSC would be a promising method for studying problems for correct electron densities and energy-level alignments in large systems.

Fermi-Löwdin orbital self-interaction correction using the strongly constrained and appropriately normed meta-GGA functional.

Despite the success of density functional approximations (DFAs) in describing the electronic properties of many-electron systems, the most widely used approximations suffer from self-interaction errors (SIEs) that limit their predictive power. Here, we describe the effects of removing SIE from the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation using the Fermi-Löwdin Orbital Self-Interaction Correction (FLOSIC) method. FLOSIC is a size-extensive implementation of the Perdew-Zunger self-interaction correction (PZ-SIC) formalism. We find that FLOSIC-SCAN calculations require careful treatment of numerical details and an integration grid that yields reliable accuracy with this approach. We investigate the performance of FLOSIC-SCAN for predicting a wide array of properties and find that it provides better results than FLOSIC-LDA and FLOSIC-PBE in nearly all cases. It also gives better predictions than SCAN for orbital energies and dissociation energies where self-interaction effects are known to be important, but total energies and atomization energies are made worse. For these properties, we also investigate the use of the self-consistent FLOSIC-SCAN density in the SCAN functional and find that this DFA@FLOSIC-DFA approach yields improved results compared to pure, self-consistent SCAN calculations. Thus, FLOSIC-SCAN provides improved results over the parent SCAN functional in cases where SIEs are dominant, and even when they are not, if the SCAN@FLOSIC-SCAN method is used.

Interpretation and Automatic Generation of Fermi-Orbital Descriptors.

We present an interpretation of Fermi-orbital descriptors (FODs) and argue that these descriptors carry chemical bonding information. We show that a bond order derived from these FODs agrees well with reference values, and highlight that optimized FOD positions used within the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) method correspond to expectations from Linnett's double-quartet theory, which is an extension of Lewis theory. This observation is independent of the underlying exchange-correlation functional, which is shown using the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA), and the strongly constrained and appropriately normed meta-GGA. To make FOD positions generally accessible, we propose and discuss four independent methods for the generation of Fermi-orbital descriptors, their implementation as well as their advantages and drawbacks. In particular, we introduce a re-implementation of the electron force field, an approach based on the centers of mass of orbital densities, a Monte Carlo-based algorithm, and a method based on Lewis-like bonding information. All results are summarized with respect to future developments of FLO-SIC and related methods. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Self-interaction-free electric dipole polarizabilities for atoms and their ions using the Fermi-Löwdin self-interaction correction

The static electric dipole polarizability of a system is a measure of the binding of its electrons. In density functional theory calculations, this binding is weakened by the presence of unphysical self-interaction in the density functional approximation (DFA), leading to overestimates of polarizabilities. To investigate this systematically we compare polarizabilities for the atoms from H to Ar and their anions and cations calculated in several DFAs and the corresponding self-interaction-corrected (SIC) DFAs with experiment and with high-level quantum chemistry reference values. The SIC results are obtained using the Fermi-L\"owdin orbital self-interaction correction (FLO-SIC) method. Removing self-interaction generally leads to smaller polarizabilities that agree significantly better with reference values. We find that FLO-SIC improves the performance of the local spin density approximation and the generalized gradient approximation (GGA) for polarizabilities to a quality that is comparable to so-called rung 4 functionals, but slightly degrades the performance of the strongly constrained and appropriately normed meta-GGA functional.

New Orbit Correction Method Based on SVDC Algorithm for Ring Based Light Sources

New Orbit Correction Method Based on SVDC Algorithm for Ring Based Light Sources

Analytic atomic gradients in the fermi-löwdin orbital self-interaction correction.

We derived, implemented, and thoroughly tested the complete analytic expression for atomic forces, consisting of the Hellmann-Feynman term and the Pulay correction, for the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) method. Analytic forces are shown to be numerically accurate through an extensive comparison to forces obtained from finite differences. Using the analytic forces, equilibrium structures for a small set of molecules were obtained. This work opens the possibility of routine self-interaction free geometrical relaxations of molecules using the FLO-SIC method. © 2018 Wiley Periodicals, Inc.

Shrinking Self-Interaction Errors with the Fermi–Löwdin Orbital Self-Interaction-Corrected Density Functional Approximation

The self-interaction error (SIE) is one of the major drawbacks of practical exchange-correlation functionals for Kohn-Sham density functional theory. Despite this, the use of methods that explicitly remove SIE from approximate density functionals is scarce in the literature due to their relatively high computational cost and lack of consistent improvement over standard modern functionals. In this article we assess the performance of a novel approach recently proposed by Pederson, Ruzsinszky, and Perdew [ J. Chem. Phys. 2014, 140, 121103] for performing self-interaction free calculations in density functional theory based on Fermi orbitals. To this end, we employ test sets consisting of reaction energies that are considered particularly sensitive to SIE. We found that the parameter-free Fermi-Löwdin orbital self-interaction correction method combined with the standard local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof (PBE) functionals gives a much better estimate of reaction energies compared to their parent LSDA and PBE functionals for most of the reactions in these two sets. They also perform on par with the global PBE0 and range-separated LC-ωPBE hybrids, which partially eliminate the SIE by including Hartree-Fock exchange. This shows the potential of the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for practical density functional calculations without SIE.

Fermi-Löwdin orbital self-interaction correction to magnetic exchange couplings

We analyze the effect of removing self-interaction error on magnetic exchange couplings using the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in the framework of density functional theory (DFT). We compare magnetic exchange couplings obtained from self-interaction-free FLOSIC calculations with the local spin density approximation (LSDA) with several widely used DFT realizations and wave function based methods. To this end, we employ the linear H-He-H model system, six organic radical molecules, and [Cu2Cl6]2- as representatives of different types of magnetic interactions. We show that the simple self-interaction-free version of LSDA improves calculated couplings with respect to LSDA in all cases, even though the nature of the exchange interaction varies across the test set, and in most cases, it yields results comparable to modern hybrids and range-separated approximate functionals.

Impacts of Earth rotation parameters on GNSS ultra-rapid orbit prediction: Derivation and real-time correction

Analysis centers (ACs) for global navigation satellite systems (GNSSs) cannot accurately obtain real-time Earth rotation parameters (ERPs). Thus, the prediction of ultra-rapid orbits in the international terrestrial reference system (ITRS) has to utilize the predicted ERPs issued by the International Earth Rotation and Reference Systems Service (IERS) or the International GNSS Service (IGS). In this study, the accuracy of ERPs predicted by IERS and IGS is analyzed. The error of the ERPs predicted for one day can reach 0.15mas and 0.053ms in polar motion and UT1-UTC direction, respectively. Then, the impact of ERP errors on ultra-rapid orbit prediction by GNSS is studied. The methods for orbit integration and frame transformation in orbit prediction with introduced ERP errors dominate the accuracy of the predicted orbit. Experimental results show that the transformation from the geocentric celestial references system (GCRS) to ITRS exerts the strongest effect on the accuracy of the predicted ultra-rapid orbit. To obtain the most accurate predicted ultra-rapid orbit, a corresponding real-time orbit correction method is developed. First, orbits without ERP-related errors are predicted on the basis of ITRS observed part of ultra-rapid orbit for use as reference. Then, the corresponding predicted orbit is transformed from GCRS to ITRS to adjust for the predicted ERPs. Finally, the corrected ERPs with error slopes are re-introduced to correct the predicted orbit in ITRS. To validate the proposed method, three experimental schemes are designed: function extrapolation, simulation experiments, and experiments with predicted ultra-rapid orbits and international GNSS Monitoring and Assessment System (iGMAS) products. Experimental results show that using the proposed correction method with IERS products considerably improved the accuracy of ultra-rapid orbit prediction (except the geosynchronous BeiDou orbits). The accuracy of orbit prediction is enhanced by at least 50% (error related to ERP) when a highly accurate observed orbit is used with the correction method. For iGMAS-predicted orbits, the accuracy improvement ranges from 8.5% for the inclined BeiDou orbits to 17.99% for the GPS orbits. This demonstrates that the correction method proposed by this study can optimize the ultra-rapid orbit prediction.

Analysis of RDSS positioning accuracy based on RNSS wide area differential technique

The BeiDou Navigation Satellite System (BDS) provides Radio Navigation Service System (RNSS) as well as Radio Determination Service System (RDSS). RDSS users can obtain positioning by responding the Master Control Center (MCC) inquiries to signal transmitted via GEO satellite transponder. The positioning result can be calculated with elevation constraint by MCC. The primary error sources affecting the RDSS positioning accuracy are the RDSS signal transceiver delay, atmospheric trans-mission delay and GEO satellite position error. During GEO orbit maneuver, poor orbit forecast accuracy significantly impacts RDSS services. A real-time 3-D orbital correction method based on wide-area differential technique is raised to correct the orbital error. Results from the observation shows that the method can successfully improve positioning precision during orbital maneuver, independent from the RDSS reference station. This improvement can reach 50% in maximum. Accurate calibration of the RDSS signal transceiver delay precision and digital elevation map may have a critical role in high precise RDSS positioning services.