Predicting the redox properties of uranyl complexes using electronic structure calculations

A plethora of chemical reactions is redox driven processes. The conversion of toxic and highly soluble U(VI) complexes to nontoxic and insoluble U(IV) form are carried out through proton coupled electron transfer by iron containing cytochromes and mineral surfaces such as machinawite. This redox process takes place through the formation of U(V) species which is unstable and immediately undergo the disproportionation reaction. Thus, theoretical methods are extremely useful to understand the reduction process of U(VI) to U(V) species. We here have carried out the structures and reduction properties of several U(VI) to U(V) complexes using a variety of electronic structure methods. Due to the lack of experimental ionization energies for uranyl (UO₂(V)‐UO₂(VI)) couple, we have benchmarked the current and popularly used density functionals and cost effective ab initio methods against the experimental electron detachment energies of [UO₂F₄]^1‐/2‐ and [UO₂Cl₄]^1‐/2‐. We find that electron detachment energy of U(VI) predicted by RI‐MP2 level on the BP86 geometries correlate nicely with the experimental and CCSD(T) data. Based on our benchmark studies, we have predicted the structures and electron detachment energies of U(V) to U(VI) species for a series of uranium complexes at the RI‐MP2//BP86 level which are experimentally inaccessible till date. We find that the redox active molecular orbital is ligand centered for the oxidation of U(VI) species, where it is metal centered (primarily f‐orbital) for the oxidation of U(V) species. Finally, we have also calculated the detachment energies of a known uranyl [UO₂]¹⁺ complex whose X‐ray crystal structures of both oxidation states are available. The large bulky nature of the ligand stabilizing the uncommon U(V) species which cannot be routinely studied by present day CCSD(T) methods as the system size are more than 20–30 atoms. The success of our efficient computational strategy can be experimentally verified in the near future for the complex as the structures are stable in gas phase which can undergo oxidation.     

Electronic structure,magnetic properties,and mixed valence character of Ce2Ni3Si5 from first principles calculations

Cerium intermetallic compounds exhibit anomalous physical properties such as heavy fermion and Kondo behaviors. Here, an ab initio study of the electronic structure, magnetic properties, and mixed valence character of Ce₂Ni₃Si₅ using density functional theory (DFT) is presented. Two theoretical methods, including pure Perdew–Burke–Ernzerhof (PBE) and PBE + U , are used. In this study, Ce³⁺ and Ce⁴⁺ are considered as two different constituents in the unit cell. The formation energy calculations on the DFT level propose that Ce is in a stable mixed valence of 3.379 at 0 K. The calculated electronic structure shows that Ce₂Ni₃Si₅ is a metallic compound with a contribution at the Fermi level from Ce 4f and Ni 3d states. With the inclusion of the effective Hubbard parameter (U _eff), the five valence electrons of 5 Ce³⁺ ions are distributed only on Ce³⁺ 4f orbitals. Therefore, the occupied Ce³⁺ 4f band is located in the valence band (VB) while Ce⁴⁺ 4f orbitals are empty and Located at the Fermi level. The calculated magnetic moment in Ce₂Ni₃Si₅ is only due to cerium (Ce³⁺) in good agreement with the experimental results. The U _eff value of 5.4 eV provides a reasonable magnetic moment of 0.981 for the unpaired electron per Ce³⁺ ion. These results may serve as a guide for studying present mixed valence cerium‐based compounds. © 2017 Wiley Periodicals, Inc.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

Origin of the Magnetic Anisotropy in Heptacoordinate NiII and CoII Complexes

The nature and magnitude of the magnetic anisotropy of heptacoordinate mononuclear Ni^II and Co^II complexes were investigated by a combination of experiment and ab initio calculations. The zero‐field splitting (ZFS) parameters D of [Ni(H₂DAPBH)(H₂O)₂](NO₃)₂ ? 2 H₂O ( 1 ) and [Co(H₂DAPBH)(H₂O)(NO₃)](NO₃) [ 2 ; H₂DAPBH=2,6‐diacetylpyridine bis‐ (benzoyl hydrazone)] were determined by means of magnetization measurements and high‐field high‐frequency EPR spectroscopy. The negative D value, and hence an easy axis of magnetization, found for the Ni^II complex indicates stabilization of the highest M_S value of the S=1 ground spin state, while a large and positive D value, and hence an easy plane of magnetization, found for Co^II indicates stabilization of the M_S=±1/2 sublevels of the S=3/2 spin state. Ab initio calculations were performed to rationalize the magnitude and the sign of D, by elucidating the chemical parameters that govern the magnitude of the anisotropy in these complexes. The negative D value for the Ni^II complex is due largely to a first excited triplet state that is close in energy to the ground state. This relatively small energy gap between the ground and the first excited state is the result of a small energy difference between the d_xy and ${{\rm{d}}_{x^2 - y^2 } }$ orbitals owing to the pseudo‐pentagonal‐bipyramidal symmetry of the complex. For Co^II, all of the excited states contribute to a positive D value, which accounts for the large magnitude of the anisotropy for this complex.     

Magnetic Properties and Electronic Structure of Neptunyl(VI) Complexes: Wavefunctions,Orbitals, and Crystal‐Field Models

The electronic structure and magnetic properties of neptunyl(VI), NpO₂²⁺, and two neptunyl complexes, [NpO₂(NO₃)₃]^? and [NpO₂Cl₄]^2?, were studied with a combination of theoretical methods: ab initio relativistic wavefunction methods and density functional theory (DFT), as well as crystal‐field (CF) models with parameters extracted from the ab initio calculations. Natural orbitals for electron density and spin magnetization from wavefunctions including spin–orbit coupling were employed to analyze the connection between the electronic structure and magnetic properties, and to link the results from CF models to the ab initio data. Free complex ions and systems embedded in a crystal environment were studied. Of prime interest were the electron paramagnetic resonance g‐factors and their relation to the complex geometry, ligand coordination, and nature of the nonbonding 5f orbitals. The g‐factors were calculated for the ground and excited states. For [NpO₂Cl₄]^2?, a strong influence of the environment of the complex on its magnetic behavior was demonstrated. Kohn–Sham DFT with standard functionals can produce reasonable g‐factors as long as the calculation converges to a solution resembling the electronic state of interest. However, this is not always straightforward.     

Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine.

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.     

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f‐shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self‐interaction correction and LDA + U, the f‐(non‐f)‐mixing interaction is decreased by spectral weights of the many‐body states of the f‐ion. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ferromagnetic interaction in an asymmetric end-to-end azido double-bridged copper(II) dinuclear complex: a combined structure, magnetic, polarized neutron diffraction and theoretical study

A new end-to-end azido double-bridged copper(II) complex [Cu(2)L(2)(N(3))2] (1) was synthesized and characterized (L=1,1,1-trifluoro-7-(dimethylamino)-4-methyl-5-aza-3-hepten-2-onato). Despite the rather long Cu-Cu distance (5.105(1) A), the magnetic interaction is ferromagnetic with J= +16 cm(-1) (H=-JS(1)S(2)), a value that has been confirmed by DFT and high-level correlated ab initio calculations. The spin distribution was studied by using the results from polarized neutron diffraction. This is the first such study on an end-to-end system. The experimental spin density was found to be localized mainly on the copper(II) ions, with a small degree of delocalization on the ligand (L) and terminal azido nitrogens. There was zero delocalization on the central nitrogen, in agreement with DFT calculations. Such a picture corresponds to an important contribution of the d(x2-y2) orbital and a small population of the d(z2) orbital, in agreement with our calculations. Based on a correlated wavefunction analysis, the ferromagnetic behavior results from a dominant double spin polarization contribution and vanishingly small ionic forms.     

Ab initio study of the complexes of halogen-containing molecules RX (X=Cl, Br, and I) and NH3: towards understanding the nature of halogen bonding and the electron-accepting propensities of covalently bonded halogen atoms

Ab initio calculations have been performed on a series of complexes formed between halogen-containing molecules and ammonia to gain a deeper insight into the nature of halogen bonding. It appears that the dihalogen molecules form the strongest halogen-bonded complexes with ammonia, followed by HOX; the charge-transfer-type contribution has been demonstrated to dominate the halogen bonding in these complexes. For the complexes involving carbon-bound halogen molecules, our calculations clearly indicate that electrostatic interactions are mainly responsible for their binding energies. Whereas the halogen-bond strength is significantly enhanced by progressive fluorine substitution, the substitution of a hydrogen atom by a methyl group in the CH(3)X...NH(3) complex weakened the halogen bonding. Moreover, remote substituent effects have also been noted in the complexes of halobenzenes with different para substituents. The influence of the hybridization state of the carbon atom bonded to the halogen atom has also been examined and the results reveal that halogen-bond strengths decrease in the order HC triple bond CX > H(2)C=CHX approximately O=CHX approximately C(6)H(5)X > CH(3)X. In addition, several excellent linear correlations have been established between the interaction energies and both the amount of charge transfer and the electrostatic potentials corresponding to an electron density of 0.002 au along the R-X axis; these correlations provide good models with which to evaluate the electron-accepting abilities of the covalently bonded halogen atoms. Finally, some positively charged halogen-bonded systems have been investigated and the effect of the charge has been discussed.     

On the Controversial Fitting of Susceptibility Curves of Ferromagnetic CuII Cubanes: Insights from Theoretical Calculations

This paper reports a theoretical analysis of the electronic structure and magnetic properties of a tetranuclear Cu^II complex, [Cu₄(HL)₄], which has a 4+2 cubane‐like structure (H₃L=N,N′‐(2‐hydroxypropane‐1,3‐diyl)bis(acetylacetoneimine)). These theoretical calculations indicate a quintet (S=2) ground state; the energy‐level distribution of the magnetic states confirm Heisenberg behaviour and correspond to an S₄ spin–spin interaction model. The dominant interaction is the ferromagnetic coupling between the pseudo‐dimeric units (J₁=22.2 cm^?1), whilst a weak and ferromagnetic interaction is found within the pseudo‐dimeric units (J₂=1.4 cm^?1). The amplitude and sign of these interactions are consistent with the structure and arrangement of the magnetic Cu 3d orbitals; they accurately simulate the thermal dependence of magnetic susceptibility, but do not agree with the reported J values (J₁=38.4 cm^?1, J₂=?18.0 cm^?1) that result from the experimental fitting. This result is not an isolated case; many other polynuclear systems, in particular 4+2 Cu^II cubanes, have been reported in which the fitted magnetic terms are not consistent with the geometrical features of the system. In this context, theoretical evaluation can be considered as a valuable tool in the interpretation of the macroscopic behaviour, thus providing clues for a rational and directed design of new materials with specific properties.     

Pentacoordinate NiII Complexes: Preparation,Magnetic Measurements,and Ab Initio Calculations of the Magnetic Anisotropy Terms

Two novel mononuclear five‐coordinate nickel complexes with distorted square‐pyramidal geometries are presented. They result from association of a tridentate “half‐unit” ligand and 6,6′‐dimethyl‐2,2′‐bipyridine according to a stepwise process that highlights the advantage of coordination chemistry in isolating an unstable tridentate ligand by nickel chelation. Their zero‐field splittings (ZFS) were studied by means of magnetic data and state‐of‐the‐art ab initio calculations. Good agreement between the experimental and theoretical axial D parameters confirms that large single‐ion nickel anisotropies are accessible. The synthetic process can also yield dinuclear nickel complexes in which the nickel ions are hexacoordinate. This possibility is facilitated by the presence of phenoxo oxygen atoms in the tridentate ligand that can introduce a bridge between the two nickel ions. Two different double bridges are characterized, with the bridging oxygen atoms coming from each nickel ion or from the same nickel ion. This coordination change introduces a difference in the antiferromagnetic interaction parameter J. Although the magnetic data confirm the presence of single‐ion anisotropies in these complexes, these terms cannot be determined in a straightforward way from experiment due to the mismatch between the principal axes of the local anisotropies and the presence of intersite anisotropies.     

Local atomic and electronic structure in nanocrystalline Sn-doped anatase TiO2.

Tin-doped anatase TiO(2) nanopowders and nanoceramics with particle sizes between 12 and 30 nm are investigated by X-ray absorption fine-structure (EXAFS) and M?ssbauer spectroscopies. Furthermore, ab initio calculations based on the density functional theory are performed to analyze changes in the electronic structure due to Sn doping. The three approaches consistently show that Sn is dissolved on substitutional bulk sites with a slight increase of the bond lengths of the inner coordination shells. The Debye-Waller factors show that the nanocrystallites are highly ordered. There is no indication of defect states or bandgap changes with Sn doping.     

Ab initio and DFT study of the electronic structures and spectroscopic properties of pyrene ligands and their cyclometalated complexes

Two ligands 1‐diphenylphosphinopyrene (1‐PyP) ( L ₁ ), 1,6‐bis(diphenylphosphino)‐pyrene (1,6‐PyP) ( L ₂ ) and their cyclometalated complexes [Pt(dppm)(1‐PyP‐H)]⁺ ( 1 ), [Pt₂(dppm)₂(1,6‐PyP‐H₂)]²⁺ (dppm = bis(diphenylphosphino)methane ( 2 ), and [Pd(dppe)(1‐PyP‐H)⁺ (dppe = bis(diphenylphosphino)ethane) ( 3 ) are investigated theoretically to explore their electronic structures and spectroscopic properties. The ground‐ and excited‐state structures are optimized by the density functional theory (DFT) and single‐excitation configuration interaction method, respectively. At the time‐dependent DFT (TDDFT) and B3LYP level, the absorption and emission spectra in solution are obtained. As revealed from the calculations, the lowest‐energy absorptions of 1 and 3 are attributed to the mixing ligand‐to‐metal charge transfer (CT)/intraligand (IL)/ligand‐to‐ligand CT transitions, while that of 2 is attributed to the IL transition. The lowest‐energy phosphorescent emissions of the cyclometalated complexes are attributed to coming from the ³ILCT transitions. With the increase of the spin‐orbit coupling effect, the phosphorescence intensities and the emissions wavelength are correspondingly increased. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical studies of ground and excited electronic states in a series of heteroleptic iridium complexes using density functional theory

A series of new heteroleptic iridium(III) complexes [Ir(C^?N)₂(N^?N)]PF₆ ( 1 ‐ 6 ) (each with two cyclometalating C^?N ligands and one neutral N^?N ancillary ligand, where C^?N = 2‐phenylpyridine (ppy), 5‐methyl‐2‐(4‐fluoro)phenylpyridine (F‐mppy), and N^?N = 2,2′‐dipyridyl (bpy), 1,10‐phenanthroline (phen), 4,4′‐diphenyl‐2,2′‐dipyridy (dphphen) were found to have rich photophysical properties. Theoretical calculations are employed for studying the photophysical and electrochemical properties. All complexes are investigated using density functional theory. Excited singlet and triplet states are examined using time‐dependent density functional theory. The low‐lying excited‐state geometries are optimized at the ab initio configuration interaction singles level. Then, the excited‐state properties are investigated in detail, including absorption and emission properties, photoactivation processes. The excited state of complexes is complicated and contains triplet metal‐to‐ligand charge transfer, triplet ligand‐to‐ligand charge transfer simultaneously. Importantly, the absorption spectra and emission maxima can be tuned significantly by changing the N^?N ligands and C^?N ligands. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Ab initio and DFT studies on van der Waals trimers: the OCS.(CO2)2 complexes

Ab initio calculations [MP2, MP4SDTQ, and QCISD(T)] using different basis sets [6-31G(d,p), cc-pVXZ (X = D, T, Q), and aug-cc-pVDZ] and density functional theory [B3LYP/6-31G(d,p)] calculations were carried out to study the OCS.(CO2)2 van der Waals trimer. The DFT has proved inappropriate to the study of this type of systems where the dispersion forces are expected to play a relevant role. Three minima isomers (two noncyclic and one cyclic) were located and characterized. The most stable isomer exhibits a noncyclic barrel-like structure whose bond lengths, angles, rotational constants, and dipole moment agree quite well with the corresponding experimental values of the only structure observed in recent microwave spectroscopic studies. The energetic proximity of the three isomers, with stabilization energies of 1442, 1371, and 1307 cm-1, respectively, at the CBS-MP2/cc-pVXZ (X = D, T, Q) level, strongly suggests that the two unobserved structures should also be detected as in the case of the (CO2)3 trimer where both noncyclic and cyclic isomers have been reported to exist. The many-body symmetry-adapted perturbation theory is employed to analyze the nature of the interactions leading to the formation of the different structures. The three-body contributions are small and stabilizing for the two most stable structures and almost negligible for the cyclic isomer.     

On the DFT Ground State of Crystalline Bromine and Iodine

We report on an erroneous ground state within common density functional theory (DFT) methods for the solid elements bromine and iodine. Phonon computations at the GGA level for both molecular crystals yield imaginary vibrational modes, erroneously indicating dynamic instability—that fact alone could easily pass as a computational artefact, but these imaginary modes lead to energetically more favorable and dynamically stable structures, made up of infinite monoatomic chains. In contrast, meta‐GGA and hybrid functionals yield the correct energetic order for bromine, while for iodine, most global hybrids do not improve the GGA result significantly. The qualitatively correct answer, in both cases, is given by the long‐range corrected hybrid LC‐ωPBE, the Minnesota functionals M06L and M06, and by periodic Hartree–Fock and MP2 theory. This poor performance of economic DFT functionals should be kept in mind, for example, during global structure optimizations of systems with significant contributions from halogen bonds.     

Intermolecular Interactions in Weak Donor–Acceptor Complexes from Symmetry‐Adapted Perturbation and Coupled‐Cluster Theory: Tetracyanoethylene–Benzene and Tetracyanoethylene–p‐Xylene

The interactions in the complexes of tetracyanothylene (TCNE) with benzene and p‐xylene, often classified as weak electron donor–acceptor (EDA) complexes, are investigated by a range of quantum chemical methods including intermolecular perturbation theory at the DFT‐SAPT (symmetry‐adapted perturbation theory combined with density functional theory) level and explicitly correlated coupled‐cluster theory at the CCSD(T)‐F12 level. The DFT‐SAPT interaction energies for TCNE–benzene and TCNE–p‐xylene are estimated to be ?35.7 and ?44.9 kJ mol^?1, respectively, at the complete basis set limit. The best estimates for the CCSD(T) interaction energy are ?37.5 and ?46.0 kJ mol^?1, respectively. It is shown that the second‐order dispersion term provides the most important attractive contribution to the interaction energy, followed by the first‐order electrostatic term. The sum of second‐ and higher‐order induction and exchange–induction energies is found to provide nearly 40 % of the total interaction energy. After addition of vibrational, rigid‐rotor, and translational contributions, the computed internal energy changes on complex formation approach results from gas‐phase spectrophotometry at elevated temperatures within experimental uncertainties, while the corresponding entropy changes differ substantially.     

A theoretical study of the structure and bonding of UOX4 (X = F, Cl, Br, I) molecules: the importance of inverse trans influence.

During nitroxide-mediated polymerization (NMP) in the presence of a nitroxide R2(R1)NO*, the reversible formation of N-alkoxyamines [P-ON(R1)R2] reduces significantly the concentration of polymer radicals (P*) and their involvement in termination reactions. The control of the livingness and polydispersity of the resulting polymer depends strongly on the magnitude of the bond dissociation energy (BDE) of the C-ON(R1)R2 bond. In this study, theoretical BDEs of a large series of model N-alkoxyamines are calculated with the PM3 method. In order to provide a predictive tool, correlations between the calculated BDEs and the cleavage temperature (T(c)), and the dissociation rate constant (k(d)), of the N-alkoxyamines are established. The homolytic cleavage of the N-OC bond is also investigated at the B3P86/6-311++G(d,p)//B3LYP/6-31G(d), level. Furthermore, a natural bond orbital analysis is carried out for some N-alkoxyamines with a O-C-ON(R1)R2 fragment, and the strengthening of their C-ON(R1)R2 bond is interpreted in terms of stabilizing anomeric interactions.