Energy Isomers Research Articles

The Lowest-Energy Isomer of C2Si2H4 Is a Bridged Ring: Reinterpretation of the Spectroscopic Data Based on DFT and Coupled-Cluster Calculations

The lowest-energy isomer of C 2 Si 2 H 4 is determined by high-accuracy ab initio calculations to be the bridged four-membered ring 1,2-didehydro-1,3-disilabicyclo[1.1.0]butane (1), contrary to prior theoretical and experimental studies favoring the three-member ring silylsilacyclopropenylidene (2). These and eight other low-lying minima on the potential energy surface are characterized and ordered by energy using the CCSD(T) method with complete basis set extrapolation, and the resulting benchmark-quality set of relative isomer energies is used to evaluate the performance of several comparatively inexpensive approaches based on many-body perturbation theory and density functional theory (DFT). Double-hybrid DFT methods are found to provide an exceptional balance of accuracy and efficiency for energy-ordering isomers. Free energy profiles are developed to reason the relatively large abundance of isomer 2 observed in previous measurements. Infrared spectra and photolysis reaction mechanisms are modeled for isomers 1 and 2, providing additional insight about previously reported spectra and photoisomerization channels.

Towards a precise determination of the excitation energy of the Thorium nuclear isomer using a magnetic bottle spectrometer

229Th is the only known nucleus with an excited state that offers the possibility for a direct laser excitation using existing laser technology. Its excitation energy has been measured indirectly to be 7.8(5) eV (≈160 nm). The energy and lifetime of the isomeric state make it the presently only suitable candidate for a nuclear optical clock, the uncertainty of the excitation energy is, however, still too large to allow for a direct laser excitation in a Paul trap. Therefore, a major goal during the past years has been an improved energy determination. One possible approach is to measure the kinetic energy of electrons which are emitted in the internal conversion decay of the first isomeric state in 229Th. For this reason an electron spectrometer based on a magnetic bottle combined with electrical retarding fields has been built. Its design, as well as first test measurements are presented, which reveal a relative energy resolution of 3 % and thus enable to measure the electrons’ expected kinetic energy to better than 0.1 eV. This is sufficiently precise to specify a laser system able to drive the nuclear clock transition in 229Th.

Fingerprints of microscopic superfluidity in HHen+ clusters

ABSTRACTThe structures and the vibrational dynamics of the complexes HHe are investigated experimentally (via mass spectrometry (MS)) and at high levels of electronic-structure theory. The MS measurements reveal interesting trends about the stability of the starting members of the HHe family. The computations establish that the basically linear, strongly bound, symmetric triatomic molecular ion He(H)He, with an equilibrium H–He distance of 0.925 Å and about 2/3 but at least 1/2 of the positive charge on H, is the molecular core of all of the complexes. Definitive quantum-chemical results are obtained for HHe and HHe, including the proton affinity of He (computed to be cm−1 via the focal-point analysis (FPA) scheme), the FPA isomerisation energy between the two linear isomers of HHe ( cm−1), and the dissociation energy of the HHe HHe + He reaction, with an FPA estimate of cm−1. The structural isomers of the He-solvated complexes are discussed up to n=18. A useful notation, [k−l−m]-HHe, is introduced to characterise qualitatively the three possible belts around the He–H–He core in HHe (), where l denotes the number of He atoms in the central belt and denote the number of He atoms in the top and bottom belts. Capping He atoms attached to the belts can be indicated by sub- and superscripts. Several possible indicators of microscopic superfluidity are investigated: He evaporation energies, rotational constants, and vibrational fundamentals.

Exploration of large amplitude motions in the Ca+Ar2 complex

ABSTRACTWe present a theoretical study of the ground electronic state potential of the Ca+Ar2 complex and of its photoabsorption spectra, simulated at temperatures ranging between 20 and 220 K. These calculations exploit a Monte-Carlo (MC) method, based on a one-electron pseudo-potential approach. A pairwise additive potential fitted to coupled cluster ab initio points, is used to model the Ca+Ar2 complex. Our study shows that the most stable form of Ca+Ar2 is a bent C2v structure, whereas the linear isomer is located at around 90 ± 10 cm−1 above in energy. The analysis of the photoabsorption spectra establishes that a structural transition from bent Ca+Ar2 to linear ArCa+Ar occurs at T∼100 K. Trends in binding energies of both isomers, bond lengths and bond angles are also discussed. Molecular orbital overlaps provide an explanation for the order of stability between the bent and linear structures.

Inorganic Molecular Electride Mg4 O3 : Structure, Bonding, and Nonlinear Optical Properties.

Growing demands of material science and, in particular, in the field of nonlinear optics (NLO) encourage us to look for stable highly polarizable molecules with excess diffuse electrons. An unusual class of compounds called electrides comply with these requirements. Many attempts have been made, yet only few electrides have been synthesized as solids and none of them as molecular species. In this paper, a new theoretically designed molecular species with electride characteristics is reported. The idea of this molecular electride comes from the formation of electride-like features in the MgO crystal with defect F-centers. The geometry of the investigated molecule can be described as a Mg4 O4 cube with one oxygen atom removed. In Mg4 O3 , two 3s electrons are pushed out from the inner area of the molecule forming a diffuse electride multicentered bond. Our calculations show that this electride-like cluster possesses a noticeably large first hyperpolarizability β=5733 au. At the same time, a complete cube Mg4 O4 and Mg4 O3 2+ without electride electron pair have much smaller β: 0 au and 741 au, respectively. This fact indicates the decisive role of the electride electron pair in NLO properties. Additionally, vertical detachment energies of isomers (VDE), excitation energies ΔE, polarizabilities α, and IR spectra were calculated. These properties, including β, are supposed to be observable experimentally and can serve as indirect evidence of the stable molecular electride formation.

The concept of laser-based conversion electron Mössbauer spectroscopy for a precise energy determination of 229mTh

$^{229}$Th is the only nucleus currently under investigation for the development of a nuclear optical clock (NOC) of ultra-high accuracy. The insufficient knowledge of the first nuclear excitation energy of $^{229}$Th has so far hindered direct nuclear laser spectroscopy of thorium ions and thus the development of a NOC. Here, a nuclear laser excitation scheme is detailed, which makes use of thorium atoms instead of ions. This concept, besides potentially leading to the first nuclear laser spectroscopy, would determine the isomeric energy to 40 $\mu$eV resolution, corresponding to 10 GHz, which is a $10^4$ times improvement compared to the current best energy constraint. This would determine the nuclear isomeric energy to a sufficient accuracy to allow for nuclear laser spectroscopy of individual thorium ions in a Paul trap and thus the development of a single-ion nuclear optical clock.

Thermal catalytic degradation of α-HBCD, β-HBCD and γ-HBCD over Fe3O4 micro/nanomaterial: Kinetic behavior, product analysis and mechanism hypothesis

The new persistent organic pollutant (POP), 1,2,5,6,9,10-hexabromocyclododecane (HBCD), has been widely detected in various environmental media and proved to be biotoxic. However, the research on catalytic degradation of HBCD is in its infancy. Herein, we examined the degradation of α-HBCD, β-HBCD and γ-HBCD, over Fe3O4 micro/nanomaterial at 200 °C. The pseudo-first-order kinetic rate constants were in the range of 0.04–0.15 min−1, with half-life values of 5–19 min. γ-HBCD is slightly less stable than β-HBCD, but both of them readily convert into α-HBCD, as consistent with the Gibbs free energies of isomers themselves. The four products containing pentabromocyclododecene, two isomers of tetrabromocyclododecene and 1,5,9-cyclododecatriene were detected by conventional GC–MS. Interestingly, a high-throughput non-target product detection were performed by ESI-FT-ICR-MS, where up to 59 types of intermediate products were determined. It is tentatively proposed that different types of bromine-removed products (C12H17Br5, C12H18Br4, C12H18, C12H19Br5, C12H24 and C12H19Br5O) and cyclododecane ring-opened products (C12H19Br7, C12H20Br6O and C12H20Br6) form via elimination reaction, nucleophilic substitution, hydrodebromination and addition reaction. Besides, most of the products that were detected contained oxygen. The average carbon oxidation state (OSc¯) of the products indicate that the oxidation reaction is the dominant reaction type. Deep oxidation products, such as small molecular organic acids (formic, acetic, propionic, and butyric acids) and gas-phase oxidation products (CO2 and CO) were further detected by ion chromatography and GC-FID, respectively. This study might provide an alternative technique for the low-cost treatment of HBCD waste.

Theoretical study of the complexes of dichlorobenzene isomers with argon. II. SAPT analysis of the intermolecular interaction

The interaction of argon with dichlorobenzene isomers (DCB-Ar) has been analyzed with the help of the symmetry-adapted perturbation theory based on the density functional description of monomer properties (DFT-SAPT). The global potential energy surface (PES) of these complexes determined from the DFT-SAPT interaction energy (Eint) values has been compared to the CCSD(T) (coupled cluster method including single and double excitations with perturbative triple excitations) PES reported in the companion Paper I [J. Makarewicz and L. Shirkov, J. Chem. Phys. 150, 074301 (2019)]. The equilibrium structures and the binding energies found using DFT-SAPT and CCSD(T) methods combined with adequate basis sets are in good agreement. Besides DCB-Ar, we confirmed that DFT-SAPT gives accurate values of these quantities for other complexes containing an aromatic molecule and Ar. However, DFT-SAPT PES of DCB-Ar is flatter than the corresponding CCSD(T) one. As a result, the intermolecular vibrational energies are systematically underestimated. The analytical form of the important interrelations between SAPT components of Eint, established previously by us [J. Makarewicz and L. Shirkov, J. Chem. Phys. 144, 204115 (2016)], has been approved for the DCB-Ar complexes. Simplified SAPT models based on these relations have been employed to explain physical reasons for differences in the structures and the binding energies of DCB-Ar isomers. It is shown that the equilibrium distance of Ar to DCB plane and the binding energy are determined mainly by dispersion energy. The shift of Ar toward Cl is caused by both exchange and dispersion terms.

Relative stability of the cis and trans isomers of octahedral cobalt(iii) complexes of the type [Co(N2O2)X2] with tetradentate salen type Schiff bases: a combined theoretical and experimental study

Energies of cis and trans isomers of a series of octahedral cobalt(iii) complexes with salen type Schiff bases have been estimated by DFT calculations.

Masses and Beta-decay Studies of Neutron-rich Nuclei using the X-array and Gammasphere

Properties of neutron-rich nuclei in the A˜160 region are important for achieving a better understanding of the nuclear structure in this region where little is known owing to diffculties in the production of these nuclei at the present nuclear physics facilities. These properties are essential ingredients in the interpretation of the rareearth peak at A˜160 in the r process abundance distribution, since theoretical models are sensitive to nuclear structure input. Predicated on these ideas, we have initiated a new experimental program at Argonne National Laboratory. During the first experiment, beams from the Californium Rare Isotope Breeder Upgrade radioactive beam facility were used in conjunction with the SATURN decay station and the X-array. We focused initially on several odd-odd nuclei, where β decays of both the ground state and an excited isomer were investigated. Because of the spin difference, a variety of structures in the daughter nuclei were selectively populated and characterized based on their decay properties. Mass measurements using the Canadian Penning Trap aimed at establishing the excitation energy of the β-decaying isomers were also carried out. Evidence was found for a change in the single-particle structure, which in turn results in the formation of a sizable N=98 sub-shell gap at large deformation. Results from the first experimental campaign using the newly-commissioned β-decay station at Gammasphere are also presented.

Models of the iron-only hydrogenase enzyme: structure, electrochemistry and catalytic activity of Fe2(CO)3(μ-dithiolate)(μ,κ1,κ2-triphos).

A series of diiron bis(2-diphenylphosphinoethyl)phenylphosphine (triphos) complexes Fe2(CO)3(μ-dithiolate)(μ,κ1,κ2-triphos) (1-4) [dithiolate = 1 pdt; 2 edt; 3 adt (R = Bz), 4 (SMe)2] have been prepared and investigated as biomimics of the diiron site of [FeFe]-hydrogenases. The triphos ligand bridges the diiron vector whilst also chelating to one iron and 1-3 exist as a mixture of basal-basal-apical (bba) and basal-basal-basal (bbb) isomers which differ in the mode of chelation. In solution the bba and bbb forms do not interconvert on the NMR time scale, but the bba isomers are fluxional, and at low temperature four forms of 1bba are seen as the conformations for the pdt ring and triphos methylene groups are frozen. Crystallographic studies have established bba (pdt) and bbb (adt) ground state conformations and in both there is a significant deviation away from the expected eclipsed conformation (Lap-Fe-Fe-Lap torsion angle 0°) by 49.4 and 24.9° respectively, suggesting that introduction of triphos leads to significant strain and DFT calculations have been used to understand the relative energies of isomers. The electron rich nature of the diiron centre in 1-4 would suggest rapid protonation, but while bridging hydride complexes such as [Fe2(CO)3(μ-pdt)(μ,κ1,κ2-triphos)(μ-H)][BF4] (1H+) can be formed the process is slow. This behavior is likely a result of the high energy barrier in forming the initial (not observed) terminal hydride which requires a significant conformational change in triphos coordination. CV studies show that all starting compounds oxidize at low potentials and the addition of [Cp2Fe][PF6] to 1 affords [Fe2(CO)3(μ-pdt)(μ,κ1,κ2-triphos)][PF6] (1+) which has been characterised by IR spectroscopy. DFT studies suggest a ground state for 1+ with a partially rotated Fe(CO)2P moiety that yields a weak semi-bridging carbonyl with the adjacent Fe(CO)P2 group. No reduction peaks are seen for 1-4 within the solvent window but 1H+ undergoes reduction at -1.7 V. All complexes act as proton-reduction catalysts in the presence of HBF4·Et2O. For 1, three separate processes are observed and their dependence on acid concentration has been probed, and a mechanistic scheme is proposed based on formation via a CECE process of 1(μ-H)H which can either slowly release H2 or undergo further reduction. Relative contributions of the three processes to the total current were found to be highly dependent upon the background electrolyte, being attributed to their relative abilities to facilitate proton transfer processes. While 2 and 4 show similar proton reduction behaviour, the adt complex 3 is quite different being attributed to facile protonation of nitrogen which is followed by addition of a second proton at the diiron centre.

The effect of two types of dibenzoannulation of pentalene on molecular energies and magnetically induced currents.

The effect of two types of dibenzo-fusion of pentalene in the singlet and triplet states on its molecular energies and magnetically induced ring currents was examined via density functional calculations. The isomerization energy decomposition analysis (IEDA) together with the calculated aromaticity indices (NICS(1)zz, HOMA and FLUπ), estimation of resonance energies (RE) and extra cyclic resonance energies (ECRE) via the NBO method, and the NICS-XY-scans revealed that the π-electronic system is the most important factor controlling the molecular energies. The [a,f] topology features greater delocalization, which results in two opposing effects: larger ECRE, but weaker π-bonding. The latter is mainly responsible for the higher energy of [a,f]-dibenzopentalene (DBP) (ΔEiso = 21.7 kcal mol-1), with the other effects being σ-orbital and electrostatic interactions. The reversal of energetic stability in the triplet states (ΔEiso = -10.8 kcal mol-1) mainly comes from the reduced Pauli repulsion in [a,f]-DBP, which stabilizes the unpaired spin density over the central trimethylenemethane subunit vs. the central pentalene subunit in [a,e]-DBP. Although the [a,e] topology only reduces the diatropic and paratropic currents of the elementary subunits, benzene and pentalene, the [a,f] topology also creates strong global paratropicity involving the benzene rings. Both DBP isomers are characterized by global and smaller semi-global and local diatropic currents in the triplet state.

Performance of Density-Functional Tight-Binding in Comparison to Ab Initio and First-Principles Methods for Isomer Geometries and Energies of Glucose Epimers in Vacuo and Solution.

Density functional theory (DFT) is a widely used methodology for the computation of molecular and electronic structure, and we confirm that B3LYP and the high-level ab initio G3B3 method are in excellent agreement for the lowest-energy isomers of the 16 glucose epimers. Density-functional tight-binding (DFTB) is an approximate version of DFT with typically comparable accuracy that is 2 to 3 orders of magnitude faster, therefore generally very suitable for processing large numbers of complex structures. Conformational isomerism in sugars is well known to give rise to a large number of isomer structures. On the basis of a comprehensive study of glucose epimers in vacuo and aqueous solution, we found that the performance of DFTB is on par to B3LYP in terms of geometrical parameters excluding hydrogen bonds and isomer energies. However, DFTB underestimates both hydrogen bonding interactions as well as torsional barriers associated with rotations of the hydroxy groups, resulting in a counterintuitive overemphasis of hydrogen bonding in both gas phase as well as in water. Although the associated root mean squared deviation from B3LYP within epimer isomer groups is only on the order of 1 kcal/mol, this deviation affects the correct assignment of major isomer ordering, which span less than 10 kcal/mol. Both second- as well as third-order DFTB methods are exhibiting similar deviations from B3LYP. Even after the inclusion of empirical dispersion corrections in vacuum, these deviations remain for a large majority of isomer energies and geometries when compared to dispersion-corrected B3LYP.

Performance of DFT for C60 Isomerization Energies: A Noticeable Exception to Jacob's Ladder.

The ability to accurately calculate relative energies of fullerenes is important in many areas of computational nanotechnology. Because of the large size of fullerenes, their relative energies cannot normally be calculated by means of high-level ab initio procedures, and therefore, density functional theory (DFT) represents a cost-effective alternative. In an extensive benchmark study, we calculate the electronic energies of eight C60 isomers by means of the high-level G4(MP2) composite procedure. G4(MP2) isomerization energies span a wide range between 307.5 and 1074.0 kJ mol-1. We use these benchmark data to assess the performance of DFT, double-hybrid DFT (DHDFT), and MP2-based ab initio methods. Surprisingly, functionals from the second and third rungs of Jacob's Ladder (i.e., GGA and meta-GGA functionals) significantly and systematically outperform hybrid and hybrid-meta-GGA functionals, which occupy higher rungs of Jacob's Ladder. In addition, DHDFT functionals do not offer a substantial improvement over meta-GGA functionals, with respect to isomerization energies. Overall, the best performing functionals with mean absolute deviations (MADs) below 15.0 kJ mol-1 are (MADs given in parentheses) the GGA N12 (14.7); meta-GGAs M06-L (10.6), M11-L (10.8), MN15-L (11.9), and TPSS-D3BJ (12.8); and the DHDFT functionals B2T-PLYP (9.3), mPW2-PLYP (9.8), B2K-PLYP (12.1), and B2GP-PLYP (12.3 kJ mol-1). In light of these results, we recommend the use of meta-GGA functionals for the calculation of fullerene isomerization energies. Finally, we show that inclusion of very small percentages of exact Hartree-Fock exchange (3-5%) slightly improves the performance of the GGA and meta-GGA functionals. However, their performance rapidly deteriorates with the inclusion of larger percentages of exact Hartree-Fock exchange.

Theoretical Study of Isomers of Doped Clusters M2Al42 with Light Element Dopants M inside and on the Surface of the Aluminum Cage

Structures, energies, and spectroscopic characteristics of the isomers of the family of doubly doped М2Al42 clusters with dopants M from the first two periods and H, Cu, and Zn atoms located in different positions at the surface and in the inner cavity of the aluminum cage have been calculated by the density functional theory method. The effect of the dopant nature on the relative energies of isomers and on the energies of their dissociation along the channels М2Al42 → 2М + Al42 and М2Al42 + 2Al → 2М+ Al44. The results are compared with the results of previous DFT calculations of endohedral (M@Al12) and exohedral (Al@)MAl11) isomers of the simpler doped clusters MAl12 with the same dopants M. The influence of the aluminum cage size on the relative energy stability of the surface and and internal sopant positions is considered.

Acylhydrazone Switches: E/Z Stability Reversed by Introduction of Hydrogen Bonds

We report a series of acylhydrazone‐based switches, of which double bond rotation can be stimulated by a pH response. The switches prefer Z rather than E configuration in the ground state because of the intramolecular hydrogen bonding, which was analyzed both thermodynamically and kinetically. Solvation studies show a linear relationship between the isomerization energy in Z/E states (ΔGE/Zo) of the switches and the hydrogen bond basicity scales (β and DN) of the solvent. This solvent effect was also found to be relative to the skeleton of compound, the less shield of hydrogen bond domain the more sensitivity to the solvent environment. By introducing multiple intramolecular hydrogen bonds, the Z isomer completely predominates for compound 4, offering an advantage over a full E/Z isomerization switching.

Analysis of Stability and (Anti)aromaticity of BN‐Dibenzo[a,e]pentalenes

Relatively scarce literature data on BN/CC isosterism in 4nπ‐electronic systems have prompted us to investigate theoretically the influence of BN pair position within the central butadiene fragment of dibenzo[a,e]pentalene on two fundamental molecular properties: stability and (anti)aromaticity. It was found that stability and aromaticity follow the same trend only for BN‐orientational isomers. The source of different stability of other isomers was examined first by an analysis of bond types and their dissociation energies and then by isomerization energy decomposition analysis and was explained in terms of classical electrostatic interactions, quantum‐mechanical orbital interactions, structural changes and electronic changes (transition from charge‐separated π‐system to the neutral one). (Anti)aromaticity was investigated by using three kinds of indices, HOMA, FLUπ, and NICS(1)zz, which indicated that delocalization at the central pentalene motif is almost unaffected by various BN arrangements, molecular perimeter is slightly affected, while the most affected subunits are five‐membered rings and benzene‐fused five‐membered rings containing only one heteroatom.

Computational characterization of the herbicide metolachlor and its mono-hydroxylated photodegradation products

Electronic structure calculations have been performed for the characterization of the different conformers of the herbicide metolachlor (MC), (2-chloro-N-(2-methyl-6-ethylphenyl)-N-(2-methoxy-1-methylethyl) acetamide), the radical structures obtained through the photo chemical breaking of the C–Cl bond and the most important, experimentally identified mono-hydroxylated photo-fragmentation products. Only the S-metolachlor enantiomer that shows herbicide activity was studied. The goal of the study has been twofold. The first target has been the conformational analysis of the neutral molecule and the corresponding radical. Ten different conformers were identified for the S-metolachlor enantiomer, all of them exhibiting a near perpendicular disposition of the amidyl and the phenyl planes. They were classified into two families of minima depending on the orientation of the carbonyl group (either cis or trans with respect to the phenyl ring). All the minima in the S-cis family were found to lie higher than those in the S-trans family. The free energy of isomerization between the lowest minima of each family was found to be $$ \Delta G_{298}^{0} $$ = 6.2 kcal/mol, while the isomerization barrier around the amide C–N bond causing atropisomerism was calculated to be 34.9 kcal/mol (in good agreement with the experimental energy of activation, 36.9 ± 3 kcal/mol) [45]. The second part of the work is devoted to the computational characterization of the most important isomeric forms of the metolachlor radical and the major photodegradation mono-hydroxylated metabolites formed through the coupling of the metolachlor radical with OH. Reaction enthalpies and free energies for the various hydroxylation pathways have been calculated. Based on their order, the relative significance of the various hydroxylation channels is justified and the preference in the phenyl ring hydroxylation is established. The lowest exothermicity value determined, − 5.8 kcal/mol at the M06/6-31 + G(d,p) level, corresponds to phenyl group hydroxylation at the p-position with regard to C–N bond, in excellent agreement with the experimental findings.

Revised M06 density functional for main-group and transition-metal chemistry

We present a hybrid metageneralized-gradient-approximation functional, revM06, which is based on adding Hartree-Fock exchange to the revM06-L functional form. Compared with the original M06 suite of density functionals, the resulting revM06 functional has significantly improved across-the-board accuracy for both main-group and transition-metal chemistry. The revM06 functional improves on the M06-2X functional for main-group and transition-metal bond energies, atomic excitation energies, isomerization energies of large molecules, molecular structures, and both weakly and strongly correlated atomic and molecular data, and it shows a clear improvement over M06 and M06-2X for noncovalent interactions, including smoother potential curves for rare-gas dimers. The revM06 functional also predicts more accurate results than M06 and M06-2X for most of the outside-the-training-set test sets examined in this study. Therefore, the revM06 functional is well-suited for a broad range of chemical applications for both main-group and transition-metal elements.

Solvation of barium atoms and singly charged cations in acetonitrile clusters.

The size distributions of neutral and cationic Ba x (CH3CN) n (x = 0, +1; n ≤ 7) clusters, as produced by a standard laser vaporization-supersonic expansion pick-up source, were determined from molecular beam experiments. The size distribution for cations is in the range of n = 1-7, whereas only the n = 1 complex is observed for neutral clusters, and these two features are unaffected by the variables controlling the performance of the cluster source. The distinct behavior is compatible with the expected charge-dipole interactions in the ionic species, which are stronger than the dipole induced-dipole interactions at play in neutral clusters, and it is corroborated by the relative magnitude of the theoretical successive binding energies (SBEs) for the lowest-lying isomers of cationic and neutral clusters with n = 1-5, as computed at the density functional theory level. The theoretical results also allow for the rationalization of the bimodal Ba+(CH3CN)1-7 size distribution, featuring an apparent minimum at n = 3, in terms of chiefly 6s-5d σ hybridization of the Ba+ ions, which ultimately leads to a relatively small third SBE for the Ba+(CH3CN)3 complex, as compared to those for n = 1, 2, and 4. Additional Born-Oppenheimer molecular dynamics simulations on the Ba+(CH3CN)2-4 clusters suggest that all of the ligands are coordinated to the Ba+ ion and prevent considering completion of the first solvent shell as responsible for the bimodal size distribution.