On the JAP Method for the Indirect Determination of Delocalized Currents from Experimental Chemical Shifts

The JAP model (after Jirásek, Anderson, and Peeks) to retrieve global current strengths from experimental ¹ H chemical shifts has been tested with DFT computations. Both global and local tropicities are correctly predicted in most cases and the quantitative agreement is overall fair. An extension of the model is found to give improvement in an exemplary critical case, where the global delocalized current is negligible and the current density map is dominated by local currents.

Electronic structure, stability, and aromaticity of M2B6 (M = Mg, Ca, Sr, and Ba): an interplay between spin pairing and electron delocalization

It has been shown in previous studies that the Be₂B₆ complex exhibits a triplet ground state with double aromaticity. In this work, the stability, electronic structure, and aromaticity of the homologous series M₂B₆ (M = Mg, Ca, Sr and Ba) were examined and compared to those of Be₂B₆. At the CCSD(T)/def2-TZVP//B3LYP/def2-TZVP level of theory, the target molecules were found to be more stable in the singlet than in the triplet spin state. Magnetically induced current densities and multicentre delocalization index (MCI) were employed to assess the aromatic character of the studied complexes. Both employed methods agree that M₂B₆ (M = Mg, Ca, Sr and Ba) are π aromatic and σ nonaromatic in the singlet ground state, and double aromatic in the triplet state. It was demonstrated that the electron counting rules of aromaticity cannot be used to correctly predict the aromaticity and relative stability of the examined molecules in different spin states.

Spatial and electronic structures of BeB8 and MgB8. How far the analogy goes?

Doping of boron clusters with Be and its heavier alkaline-earth congener, Mg usually leads to complexes of different geometry and electronic structure. In this work we showed that both neutral BeB 8 and MgB 8 exhibit a singlet ground state umbrella-like form. In addition, the stability, electronic structure, and aromaticity of the target molecules were compared. The magnetically induced current densities showed that BeB 8 and MgB 8 are double aromatic systems: π and σ electrons induce strong diatropic currents. The current densities induced in the studied complexes are of very similar intensity, but with a different spatial distribution. The found differences between the current density patterns in BeB 8 and MgB 8 arise from the very nature of the bonding interactions between the M atom and B 8 fragment, as demonstrated through the energy decomposition analysis (EDA).

Prediction of beryllium clusters (Ben; n = 3-25) from first principles

Evolutionary searches using the USPEX method (Universal Structure Predictor: Evolutionary Xtallography) combined with density functional theory (DFT) calculations were performed to obtain the global minimum structures of beryllium (Be_n, n = 3-25) clusters. The thermodynamic stability, optoelectronic and photocatalytic properties as well as the nature of bonding are considered for the most stable clusters. It is found that the cluster with n = 15 is the transition point at which the configurations change from 3D hollow cages to filled cage structures (with an interior atom appearing in the structure). All the ground state structures are energetically favorable with negative binding energies, suggesting good synthetic feasibility for these structures. The calculated relative stabilities and electronic structure show that the Be₄, Be₁₀ and, Be₁₇ clusters are the most stable structures and can be considered as superatoms. The electron configurations of Be₄, Be₁₀ and Be₁₇ clusters with 8, 20 and 34 electrons are identified as 1S² 1P⁶, 1S² 1P⁶ 1D¹⁰ 2S², 1S² 1P⁶ 1D¹⁰ 2S² 1F¹⁴, respectively. Theoretical simulations determined that all the ground state structures exhibit excellent thermal stability, where the upper-limit temperature that the structures can tolerate is 900 K. During AIMD simulation of O₂ adsorption onto the Be₁₇ cluster an interesting phenomenon was happening in which the pristine Be₁₇ cluster becomes a new stable Be₁₇O₁₆ cluster. Based on ELF (electron localization function) analysis, it can be concluded that the Be-Be bonds in the small clusters are primarily of van der Waals type, while for the larger clusters, the bonds are of metallic nature. The Be_n clusters show very strong absorption in the UV and visible regions with absorption coefficients larger than 10⁵ cm^-1, which suggests a wide range of potential advanced optoelectronics applications. The Be₁₇ cluster has a suitable band alignment in the visible-light excitation region which will produce enhanced photocatalytic activities (making it a promising material for water splitting).

Relating nucleus independent chemical shifts with integrated current density strengths

Indices based on the nucleus independent chemical shift (NICS) are the most frequently used in analysis of magnetic aromaticity. The magnetically induced current density, on the other hand, is a key concept in defining magnetic aromaticity. The integrated current strength (current strength susceptibility) was found to be a very useful tool in aromaticity studies. There is widely accepted notion that the properly chosen NICS-based index can provide information on the current density strength and direction in a molecule of interest. In this work, a detailed numerical testing of the relationship between the integrated bond current strength and the most employed NICS indices was performed for a set of 43 monocyclic aromatic molecules. Based on the statistical data analysis, the relationship between the bond current strength and its π and σ electron components, on one side, and the isotropic NICS (NICSiso and NICSπ,iso) and zz-component of the NICS tensor (NICSzz and NICSπ,zz), on the other side, was examined. It was found that between the NICSπ,zz(1) and π-electron bond current strenghts there is very good linear correlation. Quite surprisingly, it was revealed that the NICSiso(1) and NICSzz(1) are not correlated with the π electron bond current strengths. On the other hand, a reasonably good linear correlation was found between the NICSzz(1) and total bond current strengths.

Program Package for the Calculation of Origin-Independent Electron Current Density and Derived Magnetic Properties in Molecular Systems

We present SYSMOIC, a program package for the calculation of the origin-independent current density induced at first order by an external magnetic field in planar and nonplanar molecular systems. Origin independence is obtained adopting the continuous transformation of the origin of the current density method, implemented at both density functional theory (DFT) and Hartree-Fock (HF) levels. Expansion coefficients for perturbed and unperturbed molecular orbitals, over basis sets containing up to m-type Gaussian functions, can be calculated by the package itself or obtained from a Gaussian calculation. A number of different functionalities presented so far in the literature that are connected to the induced current, such as current density maps for any orientation of the inducing magnetic field, net bond current strengths, stagnation graphs, magnetic shielding densities, vorticities, and anisotropies, are now made available all together in a single multiplatform package installation.

Singlet and triplet states of the sandwich-type Be2B6 and Be2B7+ clusters. A test for the electron counting rules of aromaticity

The studied complexes exhibit double aromaticity in their triplet states in line with the predictions of Hückel and Baird's rules.

Be2B6 and Be2B7+: two double aromatic inverse sandwich complexes with spin-triplet ground state

Be₂B₆ and Be₂B₇⁺ clusters adopt interesting inverse sandwich structures with double σ/π aromaticity, and the former possesses the smallest monocyclic boron ring motif.