A novel program for the search of global minimum structures of atomic clusters and molecules in the gas phase, AUTOMATON, is introduced in this work. This program involves the following: first, the generation of an initial population, using a simplified probabilistic cellular automaton method, which allows easy control of the adequate distribution of atoms in space; second, the fittest individuals are selected to evolve, through genetic operations (mating and mutations), until the best candidate for a global minimum surfaces.In addition, we propose a simple way to build the descendant structures by establishing a ranking of genes to be inherited. Thus, by means of a chemical formula checker procedure, genes are transferred to the offspring, ensuring that they always have the appropriate type and number of atoms. It is worth noting that a fraction of the fittest group is subject to mutation operations. This program also includes algorithms to identify duplicate structures: one based on geometric similarity and another on the similar distribution of atomic charges. The effectiveness of the program was evaluated in a group of 45 molecules, considering organic and organometallic compounds (benzene, cyclopentadienyl anion, and ferrocene), Zintl ion clusters [Sn 9−m−n Ge m Bi n ] (4−n)− (n = 1−4 and m = 0−(9−n)), star-shaped clusters (Li 7 E 5+ , E = BH, C, Si, Ge) and a variety of boron-based clusters. The global minimum and the lowest-energy isomers reported in the literature were found for all the cases considered in this article. These results successfully prove AUTOMATON's effectiveness on the identification of energetically preferred structures of a wide variety of chemical species.

A simple and chemically intuitive approach is used to design ptC-containing E–C clusters (E = Si–Pb).

The magnetic aromaticity of 6‐membered monoheterocycles containing Group 13 to 16 elements (C5H5X, where X = SiH, GeH, N, P, As, O+, S+, Se+) was assessed by using 2 magnetic descriptors: the π‐electron contribution to the out‐of‐plane component of the nucleus‐independent chemical shifts (NICSzz,π) and ring current strength. The results show that both descriptors lead to the same conclusion regarding magnetic aromaticity. However, they do not agree with the predictions obtained by isotropic NICS, which is a most commonly used method. Ring current strength and NICSπ predict that benzene is the most aromatic molecule of the series, with an only slightly less aromatic pyridine. Additionally, aromaticity decreases when going down in the same group of the periodic system. The only exception is the pyrylium cation, which is predicted as the least aromatic species of this series.

The present work is focused on the elucidation of the electronic structure, bonding nature and optical properties of a series of low symmetry (C2) coordination compounds of type [Ln(III)HAM](3+), where "Ln(III)" are the trivalent lanthanide ions: La(3+), Ce(3+), Eu(3+) and Lu(3+), while "HAM" is the neutral six-nitrogen donor macrocyclic ligand [C22N6H26]. This systematic study has been performed in the framework of the Relativistic Density Functional Theory (R-DFT) and also using a multi-reference approach via the Complete Active Space (CAS) wavefunction treatment with the aim of analyzing their ground state and excited state electronic structures as well as electronic correlation. Furthermore, the use of the energy decomposition scheme proposed by Morokuma-Ziegler and the electron localization function (ELF) allows us to characterize the bonding between the lanthanide ions and the macrocyclic ligand, obtaining as a result a dative-covalent interaction. Due to a great deal of lanthanide optical properties and their technological applications, the absorption spectra of this set of coordination compounds were calculated using the time-dependent density functional theory (TD-DFT), where the presence of the intense Ligand to Metal Charge Transfer (LMCT) bands in the ultraviolet and visible region and the inherent f-f electronic transitions in the Near-Infra Red (NIR) region for some lanthanide ions allow us to propose these systems as "single antenna molecules" with potential applications in NIR technologies.

We have theoretically evaluated the behavior of the Re3(μ-Cl)3Cl6 cluster under magnetic perturbation, and it clearly shows that the magnetic response within the Re3(μ-Cl)3 plane is highly diatropic in nature. An analysis of both the magnetically induced current density (MICD) and induced magnetic field (B(ind)) allows us to classify this cluster as doubly σ- and also π-aromatic on the magnetic criterion. These findings contradict the classical Re-Re double bond representation and favor a chemical bonding pattern that involves delocalized bonds.

Since the synthesis of the first molecular cluster [Re6(μ3-Q8)X6](4-), the substitutional lability of the terminal ligands prompted new developments in their chemistry, making these molecular clusters a reasonable point of departure for building new materials. The development of novel inorganic materials of technological interest certainly requires an understanding of the electronic structure, bonding, spectroscopy, photophysical and structural properties of these clusters. Taking into account the potential applications in material sciences and the lack of systematization in the study of these kinds of clusters, the proposal of the present work is to perform a detailed theoretical study of the [Re6(μ3-Q8)X6](4-) (Q = S(2-), Se(2-), Te(2-); X = F(-), Cl(-), Br(-), I(-), CN(-), NC(-), SCN(-), NCS(-), OCN(-), NCO(-)) clusters based on the detailed description of the electronic structure of these complexes and the bonding nature between the [Re6(μ3-Q8)](2+) core and several donor-acceptor peripheral ligands. All this work was developed on the framework of the relativistic density functional theory, in which relativistic effects were incorporated by means of a two-component Hamiltonian with the zeroth-order regular approximation. To describe the relative stability of these complexes, we employed the global descriptors of chemical hardness and softness introduced by Pearson. Moreover, an analysis of bonding energetics was performed by combining a fragment approach to the molecular structure with the decomposition of the total bonding energy according to the Morokuma-Ziegler energy partitioning scheme. After an analysis of these results, we found in all cases an extensive ionic character in the bonding between the core and each peripheral ligand. The interaction between the halide ligand and the core gives about 75% ionic character, whereas the other ligands show a more covalent interaction due to effective synergic mechanisms. We conclude that the most stable clusters are those that present the stronger σ-donor terminal ligands, whereas the cluster stability starts to decrease when the π-acceptor effect will be stronger; this fact is directly related to the terminal ligand lability and the strong electrophilic character of the [Re6(μ3-Q8)](2+) core.

The high interest in lanthanide chemistry, and particularly in their luminescence, has been encouraged by the need of understanding the lanthanide chemical coordination and how the design of new luminescent materials can be affected by this. This work is focused on the understanding of the electronic structure, bonding nature, and optical properties of a set of lanthanide hexaaza macrocyclic complexes, which can lead to potential optical applications. Here we found that the DFT ground state of the open-shell complexes are mainly characterized by the manifold of low lying f states, having small HOMO-LUMO energy gaps. The results obtained from the wave function theory calculations (SO-RASSI) put on evidence the multiconfigurational character of their ground state and it is observed that the large spin-orbit coupling and the weak crystal field produce a strong mix of the ground and the excited states. The electron localization function (ELF) and the energy decomposition analysis (EDA) support the idea of a dative interaction between the macrocyclic ligand and the lanthanide center for all the studied systems; noting that, this interaction has a covalent character, where the d-orbital participation is evidenced from NBO analysis, leaving the f shell completely noninteracting in the chemical bonding. From the optical part we observed in all cases the characteristic intraligand (IL) (π-π*) and ligand to metal charge-transfer (LMCT) bands that are present in the ultraviolet and visible regions, and for the open-shell complexes we found the inherent f-f electronic transitions on the visible and near-infrared region.

In this study, we report about the relativistic effects on the aromaticity of borazine, B3N3H6, and their halogenated derivatives (B3N3F6, B3N3Cl6, B3N3Br6, B3N3I6, and B3N3At6), via the magnetically‐induced current density method. All‐electron density functional theory calculations were carried out using the four‐component Dirac‐Coulomb hamiltonian, including scalar and spin‐orbit relativistic effects. Ring current strengths were obtained by numerical integration over the current flow. These values were compared to the spin‐free (scalar relativistic) and nonrelativistic values, to assess the corresponding contributions to aromaticity. It was found that in B3N3I6 and B3N3At6 there exists a significant spin‐orbit influence, in line with the expected relativistic effects associated to the heavy elements, iodine, and astatine. Borazine, B3N3H6, is known to be slightly aromatic, but much less aromatic than benzene. The application of an external magnetic field induces a current density, that accounts for the delocalization and mobility of electrons in a molecule. Using this theoretical method, the aromaticity of the derivatives B3N3X6(X = H, F, Cl, Br, I, At) was investigated. The inclusion of heavy elements requires a relativistic treatment where the spin‐orbit coupling must be included. The figure shows the three‐dimensional pathways of the current density flow in B3N3At6.