Structures and Heats of Formation of Simple Alkaline Earth Metal Compounds: Fluorides, Chlorides, Oxides, and Hydroxides for Be, Mg, and Ca

Leaching behavior and kinetics of beryllium in beryllium-containing sludge (BCS)

Beryllium-containing sludge (BCS) is a byproduct of the physicochemical treatment of beryllium smelting wastewater. The pollutant element beryllium within BCS is highly unstable and extremely toxic, characterized by its small ionic radius and low charge density, resulting in a high risk of leaching and migration. This study is the first to investigate the leaching behavior, influencing mechanisms, and kinetic processes of beryllium in BCS under various environmental conditions. The results indicate that, under national standard conditions, beryllium exhibits a rapid leaching phase within the first 5 h, which then stabilizes after 10 h, with the total leached content significantly exceeding the leaching toxicity identification standards. Under mildly acidic (pH ≤ 5) or highly alkaline (pH = 14) conditions, beryllium demonstrates pronounced leaching and migration behaviors. Notably, in acidic conditions, the leaching rate exceeds 80% within 5 h. Combining the treatment process of beryllium-containing wastewater with analytical methods such as SEM, XPS, ToF-SIMS, and FTIR, it is revealed that due to the heterogeneous nature of BCS, the particle aggregates dissociate over time under acidic conditions. The particle surfaces become increasingly rough, leading to dissolution and the emergence of more reactive sites, resulting in a high proportion of beryllium leaching. Under these conditions, the gradual reaction of Be(OH)2 in BCS to form soluble Be2+ and its hydrolytic complexes is identified as the primary mechanism for extensive beryllium migration. The process encounters minimal diffusion resistance and is classified as reaction-controlled. In acidic conditions with pH = 4, the leaching rate of beryllium significantly increases with rising temperature. The leaching kinetics equation is [(1-x)-0.44]=e(18.26-53050RT)·t, with an apparent activation energy of 53.05 kJ mol-1.

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO₃, M(HCO₃)₂, and M(HCO₃)(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller-Peterson-Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO₂, and/or H₂O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M²⁺(aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange-correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. The best agreement with the FPD results is predicted for the τ-HCTH functional with an average unsigned error of 8.3 kcal/mol.

Electronically Excited Complex Formation in Magnesium Cluster-Halogen Atom Reactions

A near ultraviolet transition of Mg₂F has been observed in emission from the reaction between magnesium clusters, most likely Mg₃, and fluorine atoms. Because there is little evidence for upper-state internal excitation, the spectrum is assigned assuming that the upper state is quenched to its lowest vibrational levels. Two of possibly three ground-state vibrational frequencies, υ₁ = 516 ± 10 cm^-1 and υ₂ = 104 ± 10 cm^-1, have been established. Dispersed laser-induced fluorescence studies extrapolating on the observed chemiluminescence indicate an excited-state symmetric stretch frequency of order 370 ± 30 cm^-1. Electronic structure calculations at the CCSD(T)/CBS level predict that the ground state of Mg₂F has C_2v symmetry and can be described as an Mg₂⁺F^- ion pair with two Mg-F bonds. Like the MgF A-X transition that is largely a transition between Mg orbitals, the observed transition in Mg₂F is largely between orbitals on the magnesium dimer ion. The asymmetric C_∞v Mg₂⁺F^- complex is also a minimum and is predicted to be 6.7 kcal/mol higher in energy. Calculated structures for the Mg₂Cl isomers are also presented and used to further interpret the experimental results for the reaction of Mg clusters with Cl atoms. In contrast to Mg₂F, the ground state of Mg₂Cl is a linear C_∞v MgMgCl structure with the C_2v and D_∞h isomers of the MgClMg structure slightly higher in energy.

Binding and stability of MgO monomers on anatase TiO2(101)

In catalysis, MgO is often used to modify the acid-base properties of support oxides and to stabilize supported metal atoms and particles on oxides. In this study, we show how the sublimation of MgO powder can be used to deposit MgO monomers, hither on anatase TiO₂(101). A combination of x-ray electron spectroscopy, high-resolution scanning tunneling microscopy, and density functional theory is employed to gain insight into the MgO monomer binding, electronic and vibrational properties, and thermal stability. In the most stable configuration, the Mg and O of the MgO monomer bind to two surface oxygens and one undercoordinated surface titanium, respectively. The additional binding weakens the Mg-O monomer bond and makes Mg more ionic. The monomers are thermally stable up to 600 K, where the onset of diffusion into the TiO₂ bulk is observed. The monomeric MgO species on TiO₂(101) represent an ideal atomically precise system with modified acid-base properties and will be employed in our future catalytic studies.

Turning bulk materials into 0D, 1D and 2D metallic nanomaterials by selective aqueous corrosion

The realization of the facile and green synthesis of low-dimensional nanomaterials is critical not only for energy storage but also for catalysis. A selective aqueous corrosion strategy is presented here for obtaining low-dimensional metals, including nanoparticles, nanofibers and nanosheets, based on the dealloying of aqueous-favoring metal from its bulk alloy.

Pair natural orbital and canonical coupled cluster reaction enthalpies involving light to heavy alkali and alkaline earth metals: the importance of sub-valence correlation

In this work, we tested canonical and domain based pair natural orbital coupled cluster methods (CCSD(T) and DLPNO-CCSD(T), respectively) for a set of 32 ligand exchange and association/dissociation reaction enthalpies involving ionic complexes of Li, Be, Na, Mg, Ca, Sr, Ba and Pb(ii). Two strategies were investigated: in the former, only valence electrons were included in the correlation treatment, giving rise to the computationally very efficient FC (frozen core) approach; in the latter, all non-ECP electrons were included in the correlation treatment, giving rise to the AE (all electron) approach. Apart from reactions involving Li and Be, the FC approach resulted in non-homogeneous performance. The FC approach leads to very small errors (<2 kcal mol^-1) for some reactions of Na, Mg, Ca, Sr, Ba and Pb, while for a few reactions of Ca and Ba deviations up to 40 kcal mol^-1 have been obtained. Large errors are both due to artificial mixing of the core (sub-valence) orbitals of metals and the valence orbitals of oxygen and halogens in the molecular orbitals treated as core, and due to neglecting core-core and core-valence correlation effects. These large errors are reduced to a few kcal mol^-1 if the AE approach is used or the sub-valence orbitals of metals are included in the correlation treatment. On the technical side, the CCSD(T) and DLPNO-CCSD(T) results differ by a fraction of kcal mol^-1, indicating the latter method as the perfect choice when the CPU efficiency is essential. For completely black-box applications, as requested in catalysis or thermochemical calculations, we recommend the DLPNO-CCSD(T) method with all electrons that are not covered by effective core potentials included in the correlation treatment and correlation-consistent polarized core valence basis sets of cc-pwCVQZ(-PP) quality.

Ab initio study of the neutral and anionic alkali and alkaline earth hydroxides: Electronic structure and prospects for sympathetic cooling of OH. J Chem Phys 2017;146:194309. [PMID: 28527437 PMCID: PMC5438307 DOI: 10.1063/1.4983627]

We have performed a systematic ab initio study on alkali and alkaline earth hydroxide neutral (MOH) and anionic (MOH-) species where M = Li, Na, K, Rb, Cs or Be, Mg, Ca, Sr, Ba. The CCSD(T) method with extended basis sets and Dirac-Fock relativistic effective core potentials for the heavier atoms has been used to study their equilibrium geometries, interaction energies, electron affinities, electric dipole moment, and potential energy surfaces. All neutral and anionic species exhibit a linear shape with the exception of BeOH, BeOH-, and MgOH-, for which the equilibrium structure is found to be bent. Our analysis shows that the alkaline earth hydroxide anions are valence-bound whereas the alkali hydroxide anions are dipole bound. In the context of sympathetic cooling of OH- by collision with ultracold alkali and alkaline earth atoms, we investigate the 2D MOH- potential energy surfaces and the associative detachment reaction M + OH→- MOH + e-, which is the only energetically allowed reactive channel in the cold regime. We discuss the implication for the sympathetic cooling of OH- and conclude that Li and K are the best candidates for an ultracold buffer gas.

The nature of the chemical bond in BeO0,-, BeOBe+,0,-, and in their hydrogenated products HBeO0,-, BeOH, HBeOH, BeOBeH+,0,-, and HBeOBeH

The nature of the chemical bond in BeO^0,-, BeOBe^+,0,-, and in their hydrogenated products HBeO^0,-, BeOH, HBeOH, BeOBeH^+,0,-, and HBeOBeH has been studied through single and multi reference correlation methods. In all these species, excited and ionized atomic states participate in a resonant way making chemically possible molecules that have been termed hypervalent and explain also the "incompatible" geometrical structure of some species.

Laser cooling of BeCl and BeBr molecules in an ab initio method

In this study, the feasibility of laser-cooling of BeCl and BeBr molecules is studied using ab initio quantum chemistry. The potential energy curves for the X(2)Σ(+), A(2)Π, and 2(2)Π electronic states of BeCl and BeBr are plotted based on multi-reference configuration interaction plus Davidson corrections (MRCI + Q), and the spin-orbit coupling (SOC) effects are considered at the MRCI + Q level. The calculated spectroscopic parameters agree with the experimental data. Highly diagonally distributed Franck-Condon factors are determined for the A(2)Π(ν' = 0) ← X(2)Σ(+)(ν'' = 0) transition: f00(BeCl) = 0.947 and f00(BeBr) = 0.966. Moreover, the suitable radiative lifetimes τ of the A(2)Π(ν' = 0) state are determined for rapid laser cooling: τ(BeCl) = 18.38 ns and τ(BeBr) = 27.09 ns. The proposed cooling wavelengths of both BeCl and BeBr are within the ultraviolet region at λ00(BeCl) = 358.51 nm and λ00(BeBr) = 379.38 nm. Laser cooling schemes for BeCl and BeBr molecules are also developed in consideration of the SOC effects. These results indicate that the inclusion of SOC effects does not affect the judgment of the feasibility of laser cooling of BeCl and BeBr molecules, even for the given BeBr molecules in which the SOC effect is significant.

Multiconfiguration Pair-Density Functional Theory

We present a new theoretical framework, called Multiconfiguration Pair-Density Functional Theory (MC-PDFT), which combines multiconfigurational wave functions with a generalization of density functional theory (DFT). A multiconfigurational self-consistent-field (MCSCF) wave function with correct spin and space symmetry is used to compute the total electronic density, its gradient, the on-top pair density, and the kinetic and Coulomb contributions to the total electronic energy. We then use a functional of the total density, its gradient, and the on-top pair density to calculate the remaining part of the energy, which we call the on-top-density-functional energy in contrast to the exchange-correlation energy of Kohn-Sham DFT. Because the on-top pair density is an element of the two-particle density matrix, this goes beyond the Hohenberg-Kohn theorem that refers only to the one-particle density. To illustrate the theory, we obtain first approximations to the required new type of density functionals by translating conventional density functionals of the spin densities using a simple prescription, and we perform post-SCF density functional calculations using the total density, density gradient, and on-top pair density from the MCSCF calculations. Double counting of dynamic correlation or exchange does not occur because the MCSCF energy is not used. The theory is illustrated by applications to the bond energies and potential energy curves of H2, N2, F2, CaO, Cr2, and NiCl and the electronic excitation energies of Be, C, N, N(+), O, O(+), Sc(+), Mn, Co, Mo, Ru, N2, HCHO, C4H6, c-C5H6, and pyrazine. The method presented has a computational cost and scaling similar to MCSCF, but a quantitative accuracy, even with the present first approximations to the new types of density functionals, that is comparable to much more expensive multireference perturbation theory methods.

Experimental and theoretical characterization of the 2(2)A'-1(2)A' transition of BeOH/D

The hydroxides of Ca, Sr, and Ba are known to be linear molecules, while MgOH is quasilinear. High-level ab initio calculations for BeOH predict a bent equilibrium structure with a bond angle of 140.9°, indicating a significant contribution of covalency to the bonding. However, experimental confirmation of the bent structure is lacking. In the present study, we have used laser excitation techniques to observe the 2(2)A'-1(2)A' transition of BeOH/D in the energy range of 30300-32800 cm(-1). Rotationally resolved spectra were obtained, with sufficient resolution to reveal spin splittings for the electronically excited state. Two-color photoionization was used to determine an ionization energy of 66425(10) cm(-1). Ab initio calculations were used to guide the analysis of the spectroscopic data. Multireference configuration interaction calculations were used to construct potential energy surfaces for the 1(2)A', 2(2)A', and 1(2)A" states. The rovibronic eigenstates supported by these surfaces were determined using the Morse oscillator rigid bender internal dynamics Hamiltonian. The theoretical results were in sufficiently good agreement with the experimental data to permit unambiguous assignment. It was confirmed that the equilibrium geometry of the ground state is bent and that the barrier to linearity lies below the zero-point energies for both BeOH and BeOD.

New determination of the adiabatic ionization potential of the BaOH radical from laser photoionization-molecular beam experiments and ab initio calculations

The adiabatic ionization potential of the BaOH radical, as generated in a laser vaporization-supersonic expansion source has been determined by laser photoionization experiments to be (4.55 ± 0.03) eV. This value supports the three lowest out of seven previous experimental estimates, the former ranging from 4.35 to 4.62 eV. The present result is compared to ab initio calculations, as performed using both quantum chemistry at different levels of theory and density functional theory, and trying several effective core potentials and their accompanying basis sets for Ba. The most satisfactory agreement is obtained for either the adiabatic or vertical ionization potentials that derive from post-Hartree-Fock [MP2 and CCSD(T)] treatments of electron correlation, along with consideration of relativistic effects and extensive basis sets for Ba, in both BaOH and BaOH(+). Such conclusions extend to the results of related calculations on the Ba-OH dissociation energies of BaOH and BaOH(+), which were performed to help in calibrating the present computational study. Bonding in BaOH/BaOH(+), as well as possible sources of discrepancy with previous experimental determinations of the BaOH adiabatic ionization potential are discussed.