Semiempirical Approach to the Fukui Function Analysis of Uric Acid under Different pH Conditions.

Analytic Fukui functions calculated at a first-principles level are combined with experimental pKa values and the calculation of tautomerization energies to obtain the effective regioselectivity of uric acid toward electron-transfer reactions under different pH conditions. Second-order electron binding energies are also computed to determine which of the tautomers is more likely to participate in the electron transfer. A comparison of vertical and adiabatic proton detachment energies allows us to conclude that tautomerization is not mediating deprotonation and that two monoanionic species are of comparable relevance. The main difference between these monoanionic species is the ring that has been deprotonated. Both monoanionic species are produced from a single neutral tautomer and mainly produce a single dianionic tautomer. As a method for the analysis of systems affected by pH such as uric acid, we propose to plot condensed Fukui functions versus pH, allowing us to draw the effect of pH on the regioselectivity of electron transfer in a single image.

1H NMR global diatropicity in copper hydride complexes.

Understanding the magnetic response of electrons in nanoclusters is essential to interpret their NMR spectra thereby providing guidelines for their synthesis towards various target applications. Here, we consider two copper hydride clusters that have applications in hydrogen storage and release under standard temperature and pressure. Through Born-Oppenheimer molecular dynamics simulations, we study dynamics effects and their contributions to the NMR peaks. Finally, we examine the electrons' magnetic response to an applied external magnetic field using the gauge-including magnetically induced currents theory. Local diatropic currents are generated in both clusters but an interesting global diatropic current also appears. This diatropic current has contributions from three µ3-H hydrides and six Cu atoms that form a chain together with three S atoms from the closest ligands resulting in a higher shielding of these hydrides' 1H NMR response. This explains the unusual upfield chemical shift compared to the common downfield shift in similarly coordinated hydrides both observed in previous experimental reports.

Metallaborazines: To Be or Not To Be Delocalized.

In the current work, the analysis of the electronic delocalization of some metallacycles, based on borazine, was realized by employing magnetic criteria, such as the induced magnetic field and magnetically induced current densities, and electronic criteria, such as adaptative natural density partitioning and the analysis of molecular orbitals. The current metallaborazines were generated from isoelectronic substitutions. The main question is whether the electronic delocalization increases or decreases. The results showed that metal-N bonded borazines could be cataloged as delocalized compounds. On the other hand, the metal-B bonded borazines could be cataloged as nonaromatic (or weak aromatic) compounds based on the results of this analysis.

Cubic aromaticity in ligand-stabilized doped Au superatoms.

The magnetic response of valence electrons in doped gold-based M@Au8L8 q superatoms (M = Pd, Pt, Ag, Au, Cd, Hg, Ir, and Rh; L = PPh3; and q = 0, +1, +2) is studied by calculating the gauge including magnetically induced currents (GIMIC) in the framework of the auxiliary density functional theory. The studied systems include 24 different combinations of the dopant, total cluster charge, and cluster structure (cubic-like or oblate). The magnetically induced currents (both diatropic and paratropic) are shown to be sensitive to the atomic structure of clusters, the number of superatomic electrons, and the chemical nature of the dopant metal. Among the cubic-like structures, the strongest aromaticity is observed in Pd- and Pt-doped M@Au8L8 0 clusters. Interestingly, Pd- and Pt-doping increases the aromaticity as compared to a similar all-gold eight-electron system Au9L8 +1. With the recent implementation of the GIMIC in the deMon2k code, we investigated the aromaticity in the cubic and butterfly-like M@Au8 core structures, doped with a single M atom from periods 5 and 6 of groups IX-XII. Surprisingly, the doping with Pd and Pt in the cubic structure increases the aromaticity compared to the pure Au case not only near the central atom but encompassing the whole metallic core, following the aromatic trend Pd > Pt > Au. These doped (Pd, Pt)@Au8 nanoclusters show a closed shell 1S21P6 superatom electronic structure corresponding to the cubic aromaticity rule 6n + 2.

Magnetically induced currents and aromaticity in ligand-stabilized Au and AuPt superatoms.

Understanding magnetically induced currents (MICs) in aromatic or metallic nanostructures is crucial for interpreting local magnetic shielding and NMR data. Direct measurements of the induced currents have been successful only in a few planar molecules but their indirect effects are seen in NMR shifts of probe nuclei. Here, we have implemented a numerically efficient method to calculate gauge-including MICs in the formalism of auxiliary density functional theory. We analyze the currents in two experimentally synthesized gold-based, hydrogen-containing ligand-stabilized nanoclusters [HAu9(PPh3)8]2+ and [PtHAu8(PPh3)8]+. Both clusters have a similar octet configuration of Au(6s)-derived delocalized "superatomic" electrons. Surprisingly, Pt-doping in gold increases the diatropic response of the superatomic electrons to an external magnetic field and enhances the aromaticity of [PtHAu8(PPh3)8]+. This is manifested by a stronger shielding of the hydrogen proton in the metal core of the cluster as compared to [HAu9(PPh3)8]2+, causing a significant upfield shift in agreement with experimental proton NMR data measured for these two clusters. Our method allows the determination of local magnetic shielding properties for any component in large 3D nanostructures, opening the door for detailed interpretation of complex NMR spectra.

Analysis of the electronic delocalization in some isoelectronic analogues of B12 doped with beryllium and/or carbon.

In the current work, a new family of isoelectronic analogues to B12 is reported. The construction of this family was performed through the isoelectronic substitution principle to generate species such as B11C+, B11Be-, B10BeC, B10C22+, B10Be22- B9Be2C-, and B9BeC2+. The search for the global minimum was realized by utilizing genetic algorithms, while the induced magnetic field, electronic localization function, magnetic current densities, and multicenter aromaticity criteria were calculated to understand their electronic delocalization. Our results show that, in general, C atoms avoid hypercoordination, whereas we have found species with Be atoms located in hypercoordinated positions that are relatively stable. Our analysis of aromaticity indicates that B12 has double σ and π disk aromaticity. Mono, double or triple substitution of B by C+ or Be- reduces somewhat the aromaticity of the clusters, but less in the case of Be- substitution.

Calculation of the EPR g-tensor from auxiliary density functional theory.

The working equations for the calculation of the electron paramagnetic resonance (EPR) g-tensor within the framework of the auxiliary density functional theory (ADFT) are presented. The scheme known as gauge including atomic orbitals (GIAOs) is employed to treat the gauge origin problem. This ADFT-GIAO formulation possesses an inherent high computational performance, allowing for the calculation of the EPR g-tensor of molecules containing some hundreds of atoms in reasonable computational time employing moderate computational resources. The effect of the use of a gauge independent auxiliary density on the quality of the g-tensor calculation for the evaluation of the exchange-correlation contribution is analyzed in this work. The best agreement with the experiment is obtained with the BLYP functional (Becke 1988 exchange and Lee-Yang-Parr correlation) in combination with a double-ζ basis set, in particular aug-cc-pVDZ. Furthermore, models of endohedral fullerenes N@Cn, with n = {60, 70, 100, 180, 240}, were used for benchmarking its computational performance.

Efficient calculation of nuclear spin-rotation constants from auxiliary density functional theory.

The computation of the spin-rotation tensor within the framework of auxiliary density functional theory (ADFT) in combination with the gauge including atomic orbital (GIAO) scheme, to treat the gauge origin problem, is presented. For the spin-rotation tensor, the calculation of the magnetic shielding tensor represents the most demanding computational task. Employing the ADFT-GIAO methodology, the central processing unit time for the magnetic shielding tensor calculation can be dramatically reduced. In this work, the quality of spin-rotation constants obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good and very similar to the ones obtained by the coupled-cluster single-doubles-perturbative triples-GIAO methodology. With the improved computational performance achieved, the computation of the spin-rotation tensors of large systems or along Born-Oppenheimer molecular dynamics trajectories becomes feasible in reasonable times. Three models of carbon fullerenes containing hundreds of atoms and thousands of basis functions are used for benchmarking the performance. Furthermore, a theoretical study of temperature effects on the structure and spin-rotation tensor of the H(12)C-(12)CH-DF complex is presented. Here, the temperature dependency of the spin-rotation tensor of the fluorine nucleus can be used to identify experimentally the so far unknown bent isomer of this complex. To the best of our knowledge this is the first time that temperature effects on the spin-rotation tensor are investigated.

Efficient calculation of the rotational g tensor from auxiliary density functional theory.

The computation of the rotational g tensor with the recently developed auxiliary density functional theory (ADFT) gauge including atomic orbital (GIAO) methodology is presented. For the rotational g tensor, the calculation of the magnetizability tensor represents the most demanding computational task. With the ADFT-GIAO methodology, the CPU time for the magnetizability tensor calculation can be dramatically reduced. Therefore, it seems most desirable to employ the ADFT-GIAO methodology also for the computation of the rotational g tensor. In this work, the quality of rotational g tensors obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good. Furthermore, we also show that the ADFT-GIAO g tensor calculation is applicable to large systems like carbon nanotube models containing hundreds of atom and thousands of basis functions.

Magnetizability tensors from auxiliary density functional theory.

The working equations for the calculation of the magnetizability tensor in the framework of auxiliary density functional theory with gauge including atomic orbitals (ADFT-GIAO) are derived. Unlike in the corresponding conventional density functional theory implementations the numerical integration of the GIAOs is avoided in ADFT-GIAO. Our validation shows that this simplification has no effect on the accuracy of the methodology. As a result, a reliable and efficient implementation for the calculation of magnetizabilities of systems with more than 1000 atoms and 14,000 basis functions is presented.

NMR shielding tensors from auxiliary density functional theory.

The working equations for the calculation of NMR shielding tensors in the framework of auxiliary density functional theory are derived. It is shown that in this approach the numerical integration over gauge-including atomic orbitals can be avoided without the loss of accuracy. New integral recurrence relations for the required analytic electric-field-type integrals are derived. The computational performance of the resulting formalism permits shielding tensor calculations of systems with more than 1000 atoms and 15,000 basis functions.

A MinMax self-consistent-field approach for auxiliary density functional theory.

A MinMax self-consistent-field (SCF) approach is derived in the framework of auxiliary density functional theory. It is shown that the SCF convergence can be guided by the fitting coefficients that arise from the variational fitting of the Coulomb potential. An in-core direct inversion of the iterative subspace (DIIS) algorithm is presented. Due to its reduced memory demand this new in-core DIIS method can be applied without overhead to very large systems with tens of thousands of basis and auxiliary functions. Due to the new DIIS error definition systems with fractional occupation numbers can be treated, too.

Density functional theory optimized basis sets for gradient corrected functionals: 3d transition metal systems.

Density functional theory optimized basis sets for gradient corrected functionals for 3d transition metal atoms are presented. Double zeta valence polarization and triple zeta valence polarization basis sets are optimized with the PW86 functional. The performance of the newly optimized basis sets is tested in atomic and molecular calculations. Excitation energies of 3d transition metal atoms, as well as electronic configurations, structural parameters, dissociation energies, and harmonic vibrational frequencies of a large number of molecules containing 3d transition metal elements, are presented. The obtained results are compared with available experimental data as well as with other theoretical data from the literature.