Double zero tillage and foliar phosphorus fertilization coupled with microbial inoculants enhance maize productivity and quality in a maize-wheat rotation.

Maize is an important industrial crop where yield and quality enhancement both assume greater importance. Clean production technologies like conservation agriculture and integrated nutrient management hold the key to enhance productivity and quality besides improving soil health and environment. Hence, maize productivity and quality were assessed under a maize-wheat cropping system (MWCS) using four crop-establishment and tillage management practices [FBCT-FBCT (Flat bed-conventional tillage both in maize and wheat); RBCT-RBZT (Raised bed-CT in maize and raised bed-zero tillage in wheat); FBZT-FBZT (FBZT both in maize and wheat); PRBZT-PRBZT (Permanent raised bed-ZT both in maize and wheat], and five P-fertilization practices [P100 (100% soil applied-P); P50 + 2FSP (50% soil applied-P + 2 foliar-sprays of P through 2% DAP both in maize and wheat); P50 + PSB + AM-fungi; P50 + PSB + AMF + 2FSP; and P0 (100% NK with no-P)] in split-plot design replicated-thrice. Double zero-tilled PRBZT-PRBZT system significantly enhanced the maize grain, starch, protein and oil yield by 13.1-19% over conventional FBCT-FBCT. P50 + PSB + AMF + 2FSP, integrating soil applied-P, microbial-inoculants and foliar-P, had significantly higher grain, starch, protein and oil yield by 12.5-17.2% over P100 besides saving 34.7% fertilizer-P both in maize and on cropping-system basis. P50 + PSB + AMF + 2FSP again had significantly higher starch, lysine and tryptophan content by 4.6-10.4% over P100 due to sustained and synchronized P-bioavailability. Higher amylose content (24.1%) was observed in grains under P50 + PSB + AMF + 2FSP, a beneficial trait due to its lower glycemic-index highly required for diabetic patients, where current COVID-19 pandemic further necessitated the use of such dietary ingredients. Double zero-tilled PRBZT-PRBZT reported greater MUFA (oleic acid, 37.1%), MUFA: PUFA ratio and P/S index with 6.9% higher P/S index in corn-oil (an oil quality parameter highly required for heart-health) over RBCT-RBCT. MUFA, MUFA: PUFA ratio and P/S index were also higher under P50 + PSB + AMF + 2FSP; avowing the obvious role of foliar-P and microbial-inoculants in influencing maize fatty acid composition. Overall, double zero-tilled PRBZT-PRBZT with crop residue retention at 6 t/ha per year along with P50 + PSB + AMF + 2FSP while saving 34.7% fertilizer-P in MWCS, may prove beneficial in enhancing maize productivity and quality so as to reinforce the food and nutritional security besides boosting food, corn-oil and starch industry in south-Asia and collateral arid agro-ecologies across the globe.

Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments.

The method of increments and frozen natural orbital (MI-FNO) framework is introduced to help expedite the application of noisy, intermediate-scale quantum (NISQ) devices for quantum chemistry simulations. The MI-FNO framework provides a systematic reduction of the occupied and virtual orbital spaces for quantum chemistry simulations. The correlation energies of the resulting increments from the MI-FNO reduction can then be solved by various algorithms, including quantum algorithms such as the phase estimation algorithm and the variational quantum eigensolver (VQE). The unitary coupled-cluster singles and doubles VQE framework is used to obtain correlation energies for the case of small molecules (i.e., BeH2, CH4, NH3, H2O, and HF) using the cc-pVDZ basis set. The quantum resource requirements are estimated for a constrained geometry complex catalyst that is utilized in industrial settings for the polymerization of α-olefins. We show that the MI-FNO approach provides a significant reduction in the quantum bit (qubit) requirements relative to the full system simulations. We propose that the MI-FNO framework can create scalable examples of quantum chemistry problems that are appropriate for assessing the progress of NISQ devices.

Modeling protic dicationic ionic liquids based on quaternary ammonium, imidazolium or pyrrolidinium cations and bis(trifluoromethanesulfonyl)imide anion: Structure and spectral characteristics.

Protic dicationic ionic liquids (PDILs) have attracted growing attention owing to their applications in domains of electrochemistry, proton conducting materials and other diverse areas. In the present work protic dicationic ionic liquids (PDILs) comprising of quaternary ammonium-, imidazolium- or pyrrolidinium-dications and bis(trifluoromethanesulfonyl)imide (Tf2N‾) anion have been modelled as the dication-(Tf2N)2 complexes. Electronic structure, vibrational and 1H NMR spectra of these complexes have been derived employing the M06-2x density functional theory. Theoretical calculations have shown that the strength of cation-anion binding follows the order: methylpyrrolidinum > quaternary ammonium > butylpyrrolidinium > imidazolium, which can be attributed to number and strength of N-H---O and C-H---O interactions. The dication-(Tf2N)2 complexes emerge with signature in frequency up-shift of the characteristic N-H stretching in their infrared spectra. Underlying molecular interactions are unveiled through natural bond orbital analyses, Quantum theory of atoms in molecules (QTAIM) and noncovalent interaction reduced density gradient method. The calculations have shown that cation-anion binding energies increase linearly with kinetic energy density component G(r) in QTAIM analysis and proton affinities in the PDILs. A correlation between change in free energies accompanying the dication-(Tf2N)2 complexes and proton affinities has also been established.

Unveiling Noncovalent Interactions in Imidazolium, Pyrrolidinium, or Quaternary Ammonium Cation and Acetate Anion Based Protic Ionic Liquids: Structure and Spectral Characteristics.

In the present work protic ionic liquids (PILs) composed of imidazolium-, quaternary ammonium-, or pyrrolidinium-dications and acetate (OAc-) anion have been modeled as the dication-anion complexes through the M06-2x based density functional theory. It has been shown that cation-anion interaction energies are larger for the PILs containing the quaternary ammonium cation, which can be attributed to strong hydrogen bonding from the terminal ammonium protons. Underlying N-H···O and C-H···O hydrogen bonding, electrostatic, and van der Waals interactions are unraveled using the natural bond orbital analyses in conjunction with the quantum theory of atoms in molecules (QTAIM) and noncovalent interaction index reduced density gradient methods. The ramifications of noncovalent binding to 1H NMR and vibrational spectra are presented. The calculations further demonstrate a linear correlation of the kinetic energy density parameter G( r) in QTAIM analysis with the characteristic frequency shift of -NH3+ stretching in the dication-anion complexes. Moreover, the chemical shifts (δH) in 1H NMR spectra from theory reveal larger deshielding; the corresponding δH value correlates well with proton affinities and cation-anion binding energies as well. Effect of solvent (DMSO) on structure, binding energies, and 1H NMR are presented. The shifts of the characteristic carbonyl and the terminal ammonium stretching vibrations accompanying the dication-anion complexes from gas phase calculations are in consonance with the self-consistent reaction field theory.

Sudden cardiac arrest as a rare presentation of myxedema coma: case report.

Myxedema coma is a decompensated hypothyroidism which occurs due to long-standing, undiagnosed, or untreated hypothyroidism. Untreated hypothyroidism is known to affect almost all organs including the heart. It is associated with a decrease in cardiac output, stroke volume due to decreased myocardial contractility, and an increase in systemic vascular resistance. It can cause cardiac arrhythmias and the most commonly seen conduction abnormalities are sinus bradycardia, heart block, ventricular tachycardia, and torsade de pointes. The authors report a case of an elderly man who presented with sudden cardiac arrest and myxedema coma and who was successfully revived.

Primary Pulmonary Lymphoma Presenting with Superior Vena Cava Syndrome in a Young Female.

Primary Pulmonary Diffuse Large B Cell Lymphoma (PPDLBCL) is an extremely rare entity, which exhibits an aggressive behavior by compressing local blood vessels. It represents only 0.04% of all lymphoma cases and is extremely rare in young age. We present a case of a primary pulmonary lymphoma with superior vena cava syndrome (SVCS) in a young female. 27-year-old African American female presented with fever, cough, and facial puffiness for 2 weeks and unintentional weight loss. Chest examination showed decreased breath sounds and dullness on percussion on right side. Labs were normal except for mild leukocytosis, high lactate, and lactate dehydrogenase. Chest X-ray showed a large right side infiltrate with pleural effusion but chest CT showed 10 × 14 × 16 cm mass in the right lung without hilar and mediastinal lymphadenopathy. CT guided biopsy of the right lung mass was done and large B cell lymphoma was diagnosed. She received "involved field radiation" because of the bulky tumor size and superior vena cava involvement prior to R-CHOP to which she responded well. PPDLBCL should be considered as one of the differentials in a young patient with a large lung mass, which needs timely diagnosis and management.

A low-cost approach to electronic excitation energies based on the driven similarity renormalization group.

We propose an economical state-specific approach to evaluate electronic excitation energies based on the driven similarity renormalization group truncated to second order (DSRG-PT2). Starting from a closed-shell Hartree-Fock wave function, a model space is constructed that includes all single or single and double excitations within a given set of active orbitals. The resulting VCIS-DSRG-PT2 and VCISD-DSRG-PT2 methods are introduced and benchmarked on a set of 28 organic molecules [M. Schreiber et al., J. Chem. Phys. 128, 134110 (2008)]. Taking CC3 results as reference values, mean absolute deviations of 0.32 and 0.22 eV are observed for VCIS-DSRG-PT2 and VCISD-DSRG-PT2 excitation energies, respectively. Overall, VCIS-DSRG-PT2 yields results with accuracy comparable to those from time-dependent density functional theory using the B3LYP functional, while VCISD-DSRG-PT2 gives excitation energies comparable to those from equation-of-motion coupled cluster with singles and doubles.

Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability.

Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.

Benchmark coupled-cluster g-tensor calculations with full inclusion of the two-particle spin-orbit contributions.

We present a parallel implementation to compute electron spin resonance g-tensors at the coupled-cluster singles and doubles (CCSD) level which employs the ACES III domain-specific software tools for scalable parallel programming, i.e., the super instruction architecture language and processor (SIAL and SIP), respectively. A unique feature of the present implementation is the exact (not approximated) inclusion of the five one- and two-particle contributions to the g-tensor [i.e., the mass correction, one- and two-particle paramagnetic spin-orbit, and one- and two-particle diamagnetic spin-orbit terms]. Like in a previous implementation with effective one-electron operators [J. Gauss et al., J. Phys. Chem. A 113, 11541-11549 (2009)], our implementation utilizes analytic CC second derivatives and, therefore, classifies as a true CC linear-response treatment. Therefore, our implementation can unambiguously appraise the accuracy of less costly effective one-particle schemes and provide a rationale for their widespread use. We have considered a large selection of radicals used previously for benchmarking purposes including those studied in earlier work and conclude that at the CCSD level, the effective one-particle scheme satisfactorily captures the two-particle effects less costly than the rigorous two-particle scheme. With respect to the performance of density functional theory (DFT), we note that results obtained with the B3LYP functional exhibit the best agreement with our CCSD results. However, in general, the CCSD results agree better with the experimental data than the best DFT/B3LYP results, although in most cases within the rather large experimental error bars.

Ionization potential optimized double-hybrid density functional approximations.

Double-hybrid density functional approximations (DH-DFAs) provide an accurate description of the electronic structure of molecules by semiempirically mixing density functional and wavefunction theory. In this paper, we investigate the properties of the potential used in such approximations. By using the optimized effective potential approach, the consistent Kohn-Sham (KS) potential for a double-hybrid functional (including the second-order perturbational contribution) can be generated. This potential is shown to provide an improved description of orbital energies as vertical ionization potentials (IPs), relative to the perturbation-free KS potential typically used. Based on this observation, we suggest that DH-DFAs should be constructed in such a way that the potential provides accurate orbital energies. As a proof of principle, the B2-PLYP functional is reparameterized to obtain the IP-optimized B2IP-PLYP functional, using a small set of vertical IPs and atomization energies as reference data. This functional is shown to outperform B2-PLYP in a wide range of benchmarks and is en par with the related B2GP-PLYP. In particular, it is shown to be the most reliable choice in electronically difficult and multireference cases.

Probing Molecular Interactions in Functionalized Asymmetric Quaternary Ammonium-Based Dicationic Ionic Liquids.

Electronic structure, binding energies, and spectral characteristics of functionalized asymmetric dicationic ionic liquids (DILs) composed of quaternary ammonium cations substituted with the ethoxyethyl and allyl/3-phenylpropyl/methoxyethoxyethyl/pentyl functionalities on two different nitrogen centers of the dication and the bis(trifluoromethanesulfonyl)imide (Tf2N-) anion were derived employing the dispersion-corrected density functional theory. DILs based on methoxyethoxyethyl-substituted cation reveal stronger binding toward the Tf2N- anion. The measured glass transition temperatures are found to be strongly dependent on the cation-anion binding facilitated through noncovalent interactions with dominant contributions from the electrostatics and hydrogen bonding. The manifestations of these interactions to vibrational spectra, in particular, to SO2 and CF3 stretchings in the complexes are presented. It has been demonstrated that the frequency down (red)-shift of the SO2 stretching in these DILs with varying substituent follows the order: methoxyethoxyethyl (35 cm-1) > allyl (23 cm-1) > pentyl (20 cm-1) > 3-phenylpropyl (5 cm-1), which is consistent with the strength of cation-anion binding. The CF3 stretching of the anion exhibits the frequency shift in the opposite direction with its hierarchy being reversed to that of SO2 stretchings; the largest upshift (blue shift) of 60 cm-1 was predicted for the DILs composed of 3-phenlpropyl substituted dications. The direction of such frequency shift has been rationalized through the difference molecular electron density maps in conjunction with the electron density at the bond critical point in the quantum theory of atoms in molecules. The underlying cation-anion binding has been analyzed through charge distribution analysis characterized in terms of molecular electrostatic potential topography. Furthermore, the observed decomposition temperatures of DILs are shown to correlate well with the maximum surface electrostatic potential parameters in quantum theory of atoms in molecules.

Increasing the applicability of density functional theory. V. X-ray absorption spectra with ionization potential corrected exchange and correlation potentials.

Core excitation energies are computed with time-dependent density functional theory (TD-DFT) using the ionization energy corrected exchange and correlation potential QTP(0,0). QTP(0,0) provides C, N, and O K-edge spectra to about an electron volt. A mean absolute error (MAE) of 0.77 and a maximum error of 2.6 eV is observed for QTP(0,0) for many small molecules. TD-DFT based on QTP (0,0) is then used to describe the core-excitation spectra of the 22 amino acids. TD-DFT with conventional functionals greatly underestimates core excitation energies, largely due to the significant error in the Kohn-Sham occupied eigenvalues. To the contrary, the ionization energy corrected potential, QTP(0,0), provides excellent approximations (MAE of 0.53 eV) for core ionization energies as eigenvalues of the Kohn-Sham equations. As a consequence, core excitation energies are accurately described with QTP(0,0), as are the core ionization energies important in X-ray photoionization spectra or electron spectroscopy for chemical analysis.

Predicting Near Edge X-ray Absorption Spectra with the Spin-Free Exact-Two-Component Hamiltonian and Orthogonality Constrained Density Functional Theory.

Orthogonality constrained density functional theory (OCDFT) provides near-edge X-ray absorption (NEXAS) spectra of first-row elements within one electronvolt from experimental values. However, with increasing atomic number, scalar relativistic effects become the dominant source of error in a nonrelativistic OCDFT treatment of core-valence excitations. In this work we report a novel implementation of the spin-free exact-two-component (X2C) one-electron treatment of scalar relativistic effects and its combination with a recently developed OCDFT approach to compute a manifold of core-valence excited states. The inclusion of scalar relativistic effects in OCDFT reduces the mean absolute error of second-row elements core-valence excitations from 10.3 to 2.3 eV. For all the excitations considered, the results from X2C calculations are also found to be in excellent agreement with those from low-order spin-free Douglas-Kroll-Hess relativistic Hamiltonians. The X2C-OCDFT NEXAS spectra of three organotitanium complexes (TiCl4, TiCpCl3, TiCp2Cl2) are in very good agreement with unshifted experimental results and show a maximum absolute error of 5-6 eV. In addition, a decomposition of the total transition dipole moment into partial atomic contributions is proposed and applied to analyze the nature of the Ti pre-edge transitions in the three organotitanium complexes.

Increasing the applicability of density functional theory. IV. Consequences of ionization-potential improved exchange-correlation potentials.

This paper's objective is to create a "consistent" mean-field based Kohn-Sham (KS) density functional theory (DFT) meaning the functional should not only provide good total energy properties, but also the corresponding KS eigenvalues should be accurate approximations to the vertical ionization potentials (VIPs) of the molecule, as the latter condition attests to the viability of the exchange-correlation potential (VXC). None of the prominently used DFT approaches show these properties: the optimized effective potential VXC based ab initio dft does. A local, range-separated hybrid potential cam-QTP-00 is introduced as the basis for a "consistent" KS DFT approach. The computed VIPs as the negative of KS eigenvalue have a mean absolute error of 0.8 eV for an extensive set of molecule's electron ionizations, including the core. Barrier heights, equilibrium geometries, and magnetic properties obtained from the potential are in good agreement with experiment. A similar accuracy with less computational efforts can be achieved by using a non-variational global hybrid variant of the QTP-00 approach.

Massively parallel implementations of coupled-cluster methods for electron spin resonance spectra. I. Isotropic hyperfine coupling tensors in large radicals.

Coupled cluster (CC) methods provide highly accurate predictions of molecular properties, but their high computational cost has precluded their routine application to large systems. Fortunately, recent computational developments in the ACES III program by the Bartlett group [the OED∕ERD atomic integral package, the super instruction processor, and the super instruction architecture language] permit overcoming that limitation by providing a framework for massively parallel CC implementations. In that scheme, we are further extending those parallel CC efforts to systematically predict the three main electron spin resonance (ESR) tensors (A-, g-, and D-tensors) to be reported in a series of papers. In this paper inaugurating that series, we report our new ACES III parallel capabilities that calculate isotropic hyperfine coupling constants in 38 neutral, cationic, and anionic radicals that include the (11)B, (17)O, (9)Be, (19)F, (1)H, (13)C, (35)Cl, (33)S,(14)N, (31)P, and (67)Zn nuclei. Present parallel calculations are conducted at the Hartree-Fock (HF), second-order many-body perturbation theory [MBPT(2)], CC singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] levels using Roos augmented double- and triple-zeta atomic natural orbitals basis sets. HF results consistently overestimate isotropic hyperfine coupling constants. However, inclusion of electron correlation effects in the simplest way via MBPT(2) provides significant improvements in the predictions, but not without occasional failures. In contrast, CCSD results are consistently in very good agreement with experimental results. Inclusion of perturbative triples to CCSD via CCSD(T) leads to small improvements in the predictions, which might not compensate for the extra computational effort at a non-iterative N(7)-scaling in CCSD(T). The importance of these accurate computations of isotropic hyperfine coupling constants to elucidate experimental ESR spectra, to interpret spin-density distributions, and to characterize and identify radical species is illustrated with our results from large organic radicals. Those include species relevant for organic chemistry, petroleum industry, and biochemistry, such as the cyclo-hexyl, 1-adamatyl, and Zn-porphycene anion radicals, inter alia.

Relativistic Density Functional Calculations of Hyperfine Coupling with Variational versus Perturbational Treatment of Spin-Orbit Coupling.

Different approaches are compared for relativistic density functional theory (DFT) and Hartree-Fock (HF) calculations of electron-nucleus hyperfine coupling (HFC) in molecules with light atoms, in transition metal complexes, and in selected actinide halide complexes with a formal metal 5f(1) configuration. The comparison includes hybrid density functionals with range-separated exchange. Within the variationally stable zeroth-order regular approximation (ZORA) relativistic framework, the HFC is obtained (i) with a linear response (LR) method where spin-orbit (SO) coupling is treated as a linear perturbation, (ii) with a spin-polarized approach closely related to a DFT method for calculating magnetic anisotropy (MA) previously devised by van Wüllen et al. where SO coupling is included variationally, (iii) with a quasi-restricted variational SO method previously devised by van Lenthe, van der Avoird, and Wormer (LWA). The MA and LWA approaches for HFC calculations were implemented in the open-source NWChem quantum chemistry package as part of this study. The methodology extends recent implementations for calculations of electronic g-factors (J. Chem. Theor. Comput.2013, 9, 1052). The impact of electron correlation (DFT vs HF) and DFT delocalization errors, the effects of spin-polarization, the importance of treating spin-orbit coupling beyond first-order, and the magnitude of finite-nucleus effects, are investigated. Similar to calculations of g-factors, the MA approach in conjunction with hybrid functionals performs reasonably well for theoretical predictions of HFC in a wide range of scenarios.

Variational versus Perturbational Treatment of Spin-Orbit Coupling in Relativistic Density Functional Calculations of Electronic g Factors: Effects from Spin-Polarization and Exact Exchange.

Different approaches are compared for relativistic calculations of electronic g factors of molecules with light atoms, transition metal complexes, and selected complexes with actinides, using density functional theory (DFT) and Hartree-Fock (HF) theory. The comparison includes functionals with range-separated exchange. Within the variationally stable zeroth-order regular approximation (ZORA) relativistic framework, g factors are obtained with a linear response (LR) method where spin-orbit (SO) coupling is treated as a linear perturbation, a spin-polarized approach based on magnetic anisotropy (MA) that includes SO coupling variationally, and a quasi-restricted variational SO method previously devised by van Lenthe, van der Avoird, and Wormer (LWA). The MA and LWA approaches were implemented in the open-source NWChem quantum chemistry package. We address the importance of electron correlation (DFT vs HF), the importance of including spin polarization in the g tensor methodology, the question of whether the use of nonrelativistic spin density functionals is adequate for such calculations, and the importance of treating spin-orbit coupling beyond first-order. For selected systems, the extent of the DFT delocalization error is explicitly investigated via calculations of the energy as a function of fractional electron numbers. For a test set of small molecules with light main group atoms, all levels of calculation perform adequately as long as there is no energetic near-degeneracy among occupied and unoccupied orbitals. The interplay between different factors determining the accuracy of calculated g factors becomes more complex for systems with heavy elements such as third row transition metals and actinides. The MA approach is shown to perform acceptably well for a wide range of scenarios.

Increasing the applicability of density functional theory. III. Do consistent Kohn-Sham density functional methods exist?

The concept of a "consistent," Kohn-Sham (KS) density functional theory (DFT) is discussed, where the functional is able to provide good total energies and its self-consistent potential is such that the KS eigenvalues correspond to accurate approximations to the principal ionization potentials for the molecule. Today, none of the vast number of DFT approximations show this property. The one exception is the ab initio dft method built upon the optimized effective potential strategy for exchange and correlation. This qualifies as a DFT method because it represents the correlated density as a single determinant and by imposing that condition, generates local exchange and correlation operators which are used in self-consistent solutions of the orbitals and eigenvalues. Such a "consistent" DFT shares many of the properties of the Dyson equation, but without its frequency dependence and associated complications. The relationship between ab initio dft based on MBPT2 functional and GW method is discussed. Ab initio dft provides a self-consistent, frequency independent, effective independent particle alternative with a local correlation potential.

Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond.

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.