Abinitiorelativistic effective potentials with spin–orbit operators. III. Rb through Xe

Structure of small yttrium monoxide clusters, chemical bonding, and photoionization: threshold photoionization and density functional theory investigations

The photoionization (PI) spectra of small gas-phase yttrium monoxide clusters, YnO (n = 1-8), are investigated, and the adiabatic ionization energies are determined. The stable structures are obtained from density functional theory (DFT) calculations. The ground state structures are further confirmed by the CCSD(T) method. The PI spectra are calculated for these stable structures and are compared with the experimental PI spectra. The ground-state structures of the neutral and cation clusters are experimentally assigned with confidence on the basis of a favourable agreement between the experimental and calculated PI spectra. New structures are proposed for Y2O, Y6O, and Y8O compared to the previous literature. Y2O is a linear molecule in the ground state that was previously proposed as a C2v bent molecule. The YnO clusters become 3-dimensional from n ≥ 3. The O atom stays outside, bridging a triangular face of yttrium clusters. Chemical bonding between the yttrium and oxygen atoms is mostly ionic. The excess charge on the oxygen atom is around 1.4e-, transferred from the yttrium atoms bonded with it. Yttrium atoms are mostly covalently bonded. However, for the bigger clusters, free charges of both polarities appear on yttrium atoms that are not bonded with oxygen, indicating ionic interactions. Frontier orbitals consist of mainly delocalized 4d electrons with some 5s contributions, forming Y-Y bonding interactions, but with little contribution and zero contribution from the oxygen orbitals, regardless of the cluster size. The lost electron of YnO+ mostly comes from the 5s orbitals of all Y atoms in the cluster up to size n = 4, and then from 4d-5s hybrid orbitals from n ≥ 5, with the d contribution increasing with size. This is contrary to the previous view in the literature that photoionization occurs from a localized 4d orbital.

A new generation of effective core potentials: Selected lanthanides and heavy elements

We construct correlation-consistent effective core potentials (ccECPs) for a selected set of heavy atoms and f elements that are currently of significant interest in materials and chemical applications, including Y, Zr, Nb, Rh, Ta, Re, Pt, Gd, and Tb. As is customary, ccECPs consist of spin-orbit (SO) averaged relativistic effective potential (AREP) and effective SO terms. For the AREP part, our constructions are carried out within a relativistic coupled-cluster framework while also taking into account objective function one-particle characteristics for improved convergence in optimizations. The transferability is adjusted using binding curves of hydride and oxide molecules. We address the difficulties encountered with f elements, such as the presence of large cores and multiple near-degeneracies of excited levels. For these elements, we construct ccECPs with core-valence partitioning that includes 4f subshell in the valence space. The developed ccECPs achieve an excellent balance between accuracy, size of the valence space, and transferability and are also suitable to be used in plane wave codes with reasonable energy cutoffs.

Understanding the surrounding effects on Raman optical activity signatures of a chiral cage system: Cryptophane-PP-111

Cryptophane molecules are cage-like structures consisting in two hemispheres, each made of three benzene rings. These hemispheres are bound together with three O(CH2)nOlinkers of various lengths giving rise to a plethora of cryptophane derivatives. Moreover, they are able to encapsulate neutral guests: CH2Cl2, CHCl3, …; and charged species: Cs+, Tl+, …. Finally, they exhibit chiroptical properties thanks to the anti arrangement of the linkers between the hemispheres. This work focuses on the Raman optical activity (ROA) signatures of Cryptophane-111 (n=1 for each linker). More specifically, we aim at simulating accurately its ROA spectra with and without a xenon atom inside its cavity. Experimental data (Buffeteau et al., 2017) have already demonstrated the effect of the encapsulation in the low-wavenumbers region. To generate the initial structures, we rely on the novel Conformer-Rotamer Ensemble Sampling Tool (CREST) program, developed by S. Grimme and co-workers. This is required due to the flexibility provided by the linkers. The CREST algorithm seems promising and has already been used to sample the potential energy surface (PES) of target systems before the simulation of their vibrational spectroscopies (Eikås et al., 2022). We observe large similarities between the two sets of conformers (one with and one without Xe encapsulated), demonstrating the robustness of the CREST algorithm. For corresponding structures, the presence of xenon pushed the two hemispheres slightly further apart. After optimization at the DFT level, only one unique conformer has a Boltzmann population ratio greater than 1%, pointing out the relative rigidity of the cage. Based on this unique conformer, our simulations are in good agreement with the experimental data. Regarding xenon encapsulation, the (experimental and theoretical) ROA signatures at low wavenumbers are impacted: slight shifts in wavenumbers are observed as well as a decrease in relative ROA intensity for bands around 150 cm-1. The wavenumber shifts were very well reproduced by our simulations, but the experimental decrease in the ROA intensity was unfortunately not reproduced.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory

The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol^-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.

State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives

Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au₂₅(SR)₁₈ ^- is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.

ωB97X-3c: A composite range-separated hybrid DFT method with a molecule-optimized polarized valence double-ζ basis set

A new composite density functional theory (DFT) method is presented. It is based on ωB97X-V as one of the best-performing density functionals for the GMTKN55 thermochemistry database and completes the family of "3c" methods toward range-separated hybrid DFT. This method is consistently available for all elements up to Rn (Z = 1-86). Its further key ingredients are a polarized valence double-ζ (vDZP) Gaussian basis set, which was fully optimized in molecular DFT calculations, in combination with large-core effective core potentials and a specially adapted D4 dispersion correction. Unlike most existing double-ζ atomic orbital sets, vDZP shows only small basis set superposition errors (BSSEs) and can compete with standard sets of triple-ζ quality. Small residual BSSE effects are efficiently absorbed by the D4 damping scheme, which overall eliminates the need for an explicit treatment or empirical corrections for BSSE. Thorough tests on a variety of thermochemistry benchmark sets show that the new composite method, dubbed ωB97X-3c, is on par with or even outperforms standard hybrid DFT methods in a quadruple-zeta basis set at a small fraction of the computational cost. Particular strengths of this method are the description of non-covalent interactions and barrier heights, for which it is among the best-performing density functionals overall.

A new generation of effective core potentials from correlated and spin-orbit calculations: selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference.

Tuning Photophysical Properties by p-Functional Groups in Zn(II) and Cd(II) Complexes with Piperonylic Acid

Aggregation between discrete molecules is an essential factor to prevent aggregation-caused quenching (ACQ). Indeed, functional groups capable of generating strong hydrogen bonds are likely to assemble and cause ACQ and photoinduced electron transfer processes. Thus, it is possible to compare absorption and emission properties by incorporating two ligands with a different bias toward intra- and intermolecular interactions that can induce a specific structural arrangement. In parallel, the π electron-donor or electron-withdrawing character of the functional groups could modify the Highest Ocuppied Molecular Orbital (HOMO)–Lowest Unocuppied Molecular Orbital (LUMO) energy gap. Reactions of M(OAc)₂·2H₂O (M = Zn(II) and Cd(II); OAc = acetate) with 1,3-benzodioxole-5-carboxylic acid (Piperonylic acid, HPip) and 4-acetylpyridine (4-Acpy) or isonicotinamide (Isn) resulted in the formation of four complexes. The elucidation of their crystal structure showed the formation of one paddle-wheel [Zn(μ-Pip)₂(4-Acpy)]₂ (1); a mixture of one dimer and two monomers [Zn(µ-Pip)(Pip)(Isn)₂]₂·2[Zn(Pip)₂(HPip)(Isn)]·2MeOH (2); and two dimers [Cd(μ-Pip)(Pip)(4-Acpy)₂]₂ (3) and [Cd(μ-Pip)(Pip)(Isn)₂]₂·MeOH (4). They exhibit bridged (1, µ₂-η¹:η¹), bridged, chelated and monodentated (2, µ₂-η¹:η¹, µ₁-η¹:η¹ and µ₁-η¹), or simultaneously bridged and chelated (3 and 4, µ₂-η²:η¹) coordination modes. Zn(II) centers accommodate coordination numbers 5 and 6, whereas Cd(II) presents coordination number 7. We have related their photophysical properties and fluorescence quantum yields with their geometric variations and interactions supported by TD-DFT calculations.

Perturbation Theory Treatment of Spin-Orbit Coupling, Part I: Double Perturbation Theory Based on a Single-Reference Initial Approximation

We develop a perturbation theory for solving the many-body Dirac equation within a given relativistic effective-core potential approximation. Starting from a scalar-relativistic unrestricted Hartree-Fock (SR UHF) solution, we carry out a double perturbation expansion in terms of spin-orbit coupling (SOC) and the electron fluctuation potential. Computationally convenient energy expressions are derived through fourth order in SOC, second order in the electron fluctuation potential, and a total of third order in the coupling between the two. Illustrative calculations on the halogen series of neutral and singly positive diatomic molecules show that the perturbation expansion is well-converged by taking into account only the leading (nonvanishing) term at each order of the electron fluctuation potential. Our perturbation theory approach provides a computationally attractive alternative to a two-component self-consistent field treatment of SOC. In addition, it includes coupling with the fluctuation potential through third order and can be extended (in principle) to multireference calculations, when necessary for both closed- and open-shell cases, using quasi-degenerate perturbation theory.

Steric and Electronic Effects on the Structure and Photophysical Properties of Hg(II) Complexes

Since many factors influence the coordination around a metal center, steric and electronic effects of the ligands mainly determine the connectivity and, thus, the final arrangement. This is emphasized on Hg(II) centers, which have a zero point stabilization energy and, thus, a flexible coordination environment. Therefore, the unrestricted Hg(II) geometry facilitates the predominance of the ligands during the structural inception. Herein, we synthesized and characterized a series of six Hg(II) complexes with general formula (Hg(Pip)₂(dPy)) (Pip = piperonylate, dPy = 3-phenylpyridine (3-phpy) (1), 4-phenylpyridine (4-phpy) (2), 2,2′-bipyridine (2,2′-bipy) (3), 1,10-phenanthroline (1,10-phen) (4), 2,2′:6′,2′-terpyridine (terpy) (5), or di(2-picolyl)amine (dpa) (6)). The elucidation of their crystal structures revealed the arrangement of three monomers (3, 5, and 6), one dimer (4), and two coordination polymers (1 and 2) depending on the steric requirements of the dPy and predominance of the ligands. Besides, the study of their photophysical properties in solution supported by TD-DFT calculations enabled us to understand their electronic effects and the influence of the structural arrangement on them.

Six Hg(II) complexes containing 1,3-benzodioxole-5-carboxylic acid and pyridine derivatives bearing N-, N^N′-, and N^N^N-donor sites have been characterized and their crystal structure elucidated. Steric crowding defined the predominance of the linkers in the structural arrangement, varying from monomers to coordination polymers. Their photophysical properties were analyzed through experimental and TD-DFT calculations, identifying the character and regions involved on each electronic transition. In addition, seclusion of the HOMO from the LUMO defined them as appropriate to selective band gap tuning.

A Heteropolynuclear Pt-Ag System Having Cycloplatinated Rollover Bipyridyl Units

The synthesis and photophysical properties of the heteropolynuclear Pt-Ag complex having cyclometalated rollover bipyridine ligands (bpy*) and bridging pyrazolato ligands are reported. The Pt₂Ag₂ complex was synthesized by two step reactions from a Pt(II) complex precursor having the rollover bpy* ligand, [Pt(bpy*)(dmso)Cl], with 3,5-dimethylpyrazole (Me₂pzH) and a subsequent replacement of NH protons in the Me₂pzH moieties with the Ag(I) ion. The Pt₂Ag₂ complex exists as a mixture of U- and Z-shaped isomers in solution, whose structures were clearly determined by single-crystal X-ray structural analyses. NMR studies using the single crystals revealed rapid isomerization of the Pt₂Ag₂ complexes in solution, although the Pt₂Ag₂ structures were supported effectively by the multiple metal-metal interactions. Furthermore, the Pt₂Ag₂ framework can capture a Ag(I) ion during the U-Z isomerization to afford a Pt₂Ag₃ core with the formation of Pt → Ag dative bonds. The Pt₂Ag₃ complex showed further aggregation to form a dimer structure in the presence of coordinating solvent via the crystallization process. The formation of Pt → Ag dative bonds significantly changes the emission energy from the Pt₂Ag₂ complex, while the emission spectra of U- and Z-isomers of Pt₂Ag₂ complex almost coincide with each other and their emissive properties are very similar. The density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations revealed the effects of additional Ag(I) ion on the photophysical properties of the heteropolynuclear Pt-Ag complexes bearing the rollover bpy* ligands.

Impact of the Carbene Derivative Charge on the Decomposition Rates of Hoveyda-Grubbs-like Metathesis Catalysts

Hoveyda–Grubbs metathesis catalysts undergo a relatively fast decomposition in the presence of olefins. Using a computational density functional theory approach, we show that positively charged derivatives of N-heterocyclic carbenes have little impact on the degradation/deactivation rates of such catalysts with respect to neutral carbenes. On the other hand, the hypothetical anionic Hoveyda–Grubbs-like catalysts are predicted to less likely undergo degradation in the presence of the olefin, while being as active as standard, neutral Hoveyda–Grubbs catalysts.

Indium(iii) {2}-metallacryptates assembled from 2,6-dipicolinoyl-bis(N,N-diethylthiourea)

{2}-Metallacryptates are constructed from the self-assembly of 2,6-dipicolinoylbis(N,N-diethylthiourea) ligands and mixtures of indium(iii) chloride and chlorides of monovalent cations.

The nature of the T=T double bond (T = B, Al, Ga, In) in dialumene and its derivatives: topological study of the electron localization function (ELF)

The local electronic structure of the Al=Al bond was studied in dialumene and derivatives of dialumene in which the Al atoms were substituted by B, Ga, or In atoms. DFT calculations were performed using the B3LYP, B3PW91, PBE0, M06-L, and M06-2X functionals. Topological analysis of the electron localization function described the covalent bonds mentioned above using the disynaptic basins V_i=1,2(B,B), V_i=1,2(Al,Al), V(Ga,Ga), and V_i=1,2(In,In). The basin populations were smaller than 4 e, as expected for a double bond: B=B 2.97 e, Al=Al 3.44–3.5 e, Ga=Ga 3.58 e, and In=In 3.86 e. The Al=Al, Ga=Ga, and In=In bonds were found to be intermediate in character between single and double bonds. Topological analysis of the ρ(r) field for dialumene showed a non-nuclear attractor along the Al=Al bond, with a pseudoatom basin population of 0.937 e. NBO analysis suggested that a double bond occurred only in the molecules containing Al, Ga, or In atoms. The character of the Ga=Ga bond was observed to be strongly dependent on the effective core potential used in the calculations.

Potential Energy Surface of Group 11 Trimers (Cu, Ag, Au): Bond Angle Isomerism in Au3

Potential energy surfaces for the group 11 trimers were generated at various levels of coupled-cluster theory to examine the effects of Jahn-Teller distortions. Our calculations show that the lowest-energy conformer for Cu₃, Ag₃, and Au₃ is the ²B₂ (∼65° isomer) without spin-orbit corrections. Spin-orbit corrections have negligible contributions to the relative energies for the angle dependence of the potential energy surfaces for Cu₃ and Ag₃. The inclusion of spin-orbit corrections for Au₃ makes the ²B₂ (∼65°) and ²A₁ (∼55°) states approximately degenerate. A novel ²B₂ isomer of Au₃ at an obtuse angle of ∼125° was also characterized, providing evidence for bond angle isomerism on the same ²B₂ potential energy surface. Spin-orbit corrections increase the barrier height between the ²B₂ (65°) and ²B₂ (125°) bond angle isomers of Au₃. The calculated symmetric stretch vibrational frequencies are in good agreement with the available experimental values. All frequencies calculated for the Au₃ ²B₂ (∼125°) state are real, and there is at least one bound bending vibration for this state. Jahn-Teller parameters are also derived for each trimer.

An implementation of spin-orbit coupling for band structure calculations with Gaussian basis sets: Two-dimensional topological crystals of Sb and Bi

We present an implementation of spin-orbit coupling (SOC) for density functional theory band structure calculations that makes use of Gaussian basis sets. It is based on the explicit evaluation of SOC matrix elements, both the radial and angular parts. For all-electron basis sets, where the full nodal structure is present in the basis elements, the results are in good agreement with well-established implementations such as VASP. For more practical pseudopotential basis sets, which lack nodal structure, an ad-hoc increase of the effective nuclear potential helps to capture all relevant band structure variations induced by SOC. In this work, the non-relativistic or scalar-relativistic Kohn-Sham Hamiltonian is obtained from the CRYSTAL code and the SOC term is added a posteriori. As an example, we apply this method to the Bi(111) monolayer, a paradigmatic 2D topological insulator, and to mono- and multilayer Sb(111) (also known as antimonene), the former being a trivial semiconductor and the latter a topological semimetal featuring topologically protected surface states.

Efficient Method for Calculating Effective Core Potential Integrals

Effective core potential (ECP) integrals are among the most difficult one-electron integrals to calculate due to the projection operators. The radial part of these operators may include r⁰, r^-1, and r^-2 terms. For the r⁰ terms, we exploit a simple analytic expression for the fundamental projected integral to derive new recurrence relations and upper bounds for ECP integrals. For the r^-1 and r^-2 terms, we present a reconstruction method that replaces these terms by a sum of r⁰ terms and show that the resulting errors are chemically insignificant for a range of molecular properties. The new algorithm is available in Q-Chem 5.0 and is significantly faster than the ECP implementations in Q-Chem 4.4, GAMESS (US) and Dalton 2016.

First-principles modeling of hafnia-based nanotubes

Hybrid density functional theory calculations were performed for the first time on structure, stability, phonon frequencies, and thermodynamic functions of hafnia-based single-wall nanotubes. The nanotubes were rolled up from the thin free layers of cubic and tetragonal phases of HfO₂ . It was shown that the most stable HfO₂ single-wall nanotubes can be obtained from hexagonal (111) layer of the cubic phase. Phonon frequencies have been calculated for different HfO₂ nanolayers and nanotubes to prove the local stability and to find the thermal contributions to their thermodynamic functions. The role of phonons in stability of nanotubes seems to be negligible for the internal energy and noticeable for the Helmholtz free energy. Zone folding approach has been applied to estimate the connection between phonon modes of the layer and nanotubes and to approximate the nanotube thermodynamic properties. It is found that the zone-folding approximation is sufficiently accurate for heat capacity, but less accurate for entropy. The comparison has been done between the properties of TiO₂ , ZrO₂ , and HfO₂ . © 2017 Wiley Periodicals, Inc.

New insight on cubic–tetragonal–monoclinic phase transitions in ZrO2: ab initio study and symmetry analysis

A group-theory analysis of temperature-induced phase transitions in ZrO2 has been performed in the framework of the group–subgroup relationship tree (Bärnighausen tree) with the computer tools of the Bilbao Crystallographic Server. The transition paths including symmetry-allowed intermediate phases have been established. The active irreducible representations corresponding to soft phonon modes and spontaneous deformation strains responsible for the phase transitions have been determined. The phonon mode frequencies at the symmetry points of the Brillouin zones of cubic, tetragonal and monoclinic phases have been calculated using the ab initio density functional theory method. As a result, the soft modes and their symmetries have been revealed, which are in a complete agreement with the group-theoretical predictions.

BNg3F3: the first three noble gas atoms inserted into mono-centric neutral compounds - a theoretical study

Following the study of HXeOXeH and HXeCCXeH, in which two Xe atoms were inserted into H2O and C2H2 theoretically and experimentally, the structures and stability of BNg3F3 (Ng = Ar, Kr and Xe), in which three Ng atoms are inserted into BF3, have been explored theoretically using DFT and ab initio calculations. It is shown that BNg3F3 (Ng = Ar, Kr and Xe) with D3h symmetry are local minima with short B-Ng bond lengths of 1.966, 2.027 and 2.214 Å at the CCSD(T)/aug-cc-pVTZ/LJ18 level, which are close to their covalent limits. Note that although BNg3F3 (Ng = Kr and Xe) are energetically higher than the dissociation products 3Ng + BF3, they are still kinetically stable as metastable species with protecting barriers of 13.38 and 17.99 kcal mol(-1) for BKr3F3 and BXe3F3. Moreover, BKr3F3, the tri-Kr-inserted compound, even has comparable kinetic stability to HXeOXeH and HXeOXeF. In addition, upon the formation of BNg3F3, there is a large amount of charge transferred from B to Ng of at least 0.619 e. The calculated Wiberg Bond Indices (WBI) suggest that B-Ng bonds are naturally singly bonded; the large vibrational frequencies of B-Ng and Ng-F stretching modes and the negative Laplacian electron density of B-Ng bonds confirm further that BNg3F3 are stiff molecules with covalent B-Ng bonds. It should be noted that three Ng atoms inserted into mono-centric neutral molecules have not been reported previously. We hope that the present theoretical study may provide important evidence for the experimental synthesis of BNg3F3.

A Noble-Gas Hydride in a Nitrogen Medium: Structure, Spectroscopy, and Intermolecular Vibrations of HXeBr@(N2)22

Noble-gas hydrides have been extensively studied in noble gas matrices. However, little is known on their stability and properties in molecular hosts. Here, HXeBr in the N2 environment is modeled at the B3LYP-D level of theory in a complete single shell of 22 N2 molecules. The system is compared to similar models of HXeBr in CO2 and Xe clusters. The optimized structure of (HXeBr)@(N2)22 is of low symmetry and is highly anisotropic. None of the N2 molecules are freely rotating, and the host molecules are not symmetrically positioned with respect to the HXeBr axis. The axes of the N2 molecules are nonuniformly distributed. The computed anharmonic H-Xe stretching frequency of HXeBr in the N2 cluster is in good accord with the experimental value. The soft-mode frequencies of the cluster including both intermolecular vibrations and librations, have a broad distribution that ranges from 8.7 to 107 cm(-1). It is expected that these findings and specifically, the single-shell model, may shed light also on the local structure and vibrations of other impurities in a molecular media.

A theoretical study on the reaction of ozone with aqueous iodide

Atmospheric iodine chemistry plays a key role in tropospheric ozone catalytic destruction, new particle formation, and as one of the possible sinks of gaseous polar elemental mercury. Moreover, it has been recently proposed that reaction of ozone with iodide on the sea surface could be the major contributor to the chemical loss of atmospheric ozone. However, the mechanism of the reaction between aqueous iodide and ozone is not well known. The aim of this paper is to improve the understanding of such a mechanism. In this paper, an ab initio study of the reaction of aqueous iodide and ozone is presented, evaluating thermodynamic data of the different reactions proposed in previous experimental studies. In addition, the structures, energetics and possible evolution of the key IOOO(-) intermediate are discussed for the first time.

Size-dependent properties of transition metal clusters: from molecules to crystals and surfaces--computational studies with the program ParaGauss

In the so-called scalable regime the size-dependent behavior of the physical and chemical properties of transition metal clusters is described by scaling relationships. For most quantities this scalable regime is reached for cluster sizes between a few tens and a few hundreds of atoms, hence for systems for which an accurate treatment by density functional theory is still feasible. Thus, by invoking scaling relations one is able to obtain properties of very large nanoparticles and even the bulk limit from the results of a series of smaller cluster models. In this invited review we illustrate this strategy by exploiting results from computational studies that mostly were carried out with the density functional theory software ParaGauss. We address mainly the size-dependent behavior of the properties of transition metal clusters. To this end, we first present benchmark studies probing various approximations that are used in such density functional calculations. Subsequently we show how physical insight may be gained by exploring less understood types of systems. These applications range from bare clusters to nanoislands and nanoalloys to adsorption complexes.

Interpretation of the vacuum ultraviolet photoabsorption spectrum of iodobenzene by ab initio computations

Identification of many Rydberg states in iodobenzene, especially from the first and fourth ionization energies (IE1 and IE4, X(2)B1 and C(2)B1), has become possible using a new ultraviolet (UV) and vacuum-ultraviolet (VUV) absorption spectrum, in the region 29 000-87 000 cm(-1) (3.60-10.79 eV), measured at room temperature with synchrotron radiation. A few Rydberg states based on IE2 (A(2)A2) were found, but those based on IE3 (B(2)B2) are undetectable. The almost complete absence of observable Rydberg states relating to IE2 and IE3 (A(2)A2 and B(2)B2, respectively) is attributed to them being coupled to the near-continuum, high-energy region of Rydberg series converging on IE1. Theoretical studies of the UV and VUV spectra used both time-dependent density functional (TDDFT) and multi-reference multi-root doubles and singles-configuration interaction methods. The theoretical adiabatic excitation energies, and their corresponding vibrational profiles, gave a satisfactory interpretation of the experimental results. The calculations indicate that the UV onset contains both 1(1)B1 and 1(1)B2 states with very low oscillator strength, while the 2(1)B1 state was found to lie under the lowest ππ(∗) 1(1)A1 state. All three of these (1)B1 and (1)B2 states are excitations into low-lying σ(∗) orbitals. The strongest VUV band near 7 eV contains two very strong ππ(∗) valence states, together with other weak contributors. The lowest Rydberg 4b16s state (3(1)B1) is very evident as a sharp multiplet near 6 eV; its position and vibrational structure are well reproduced by the TDDFT results.

Energetic stability and photocatalytic activity of SrTiO3 nanowires: ab initio simulations

Ab initio simulations have been performed to describe, for the first time, energetic stability and photocatalytic activity of SrTiO₃ nanowires.

Ab initio study of the RbSr electronic structure: potential energy curves, transition dipole moments, and permanent electric dipole moments

Excited states and the ground state of the diatomic molecule RbSr were calculated by post Hartree-Fock molecular orbital theory up to 22 000 cm(-1). We applied a multireference configuration interaction calculation based on multiconfigurational self-consistent field wave functions. Both methods made use of effective core potentials and core polarization potentials. Potential energy curves, transition dipole moments, and permanent electric dipole moments were determined for RbSr and could be compared with other recent calculations. We found a good agreement with experimental spectra, which have been obtained recently by helium nanodroplet isolation spectroscopy. For the lowest two asymptotes (Rb (5s (2)S) + Sr (5s4d (3)P°) and Rb (5p (2)P°) + Sr (5s(2) (1)S)), which exhibit a significant spin-orbit coupling, we included relativistic effects by two approaches, one applying the Breit-Pauli Hamiltonian to the multireference configuration interaction wave functions, the other combining a spin-orbit Hamiltonian and multireference configuration interaction potential energy curves. Using the results for the relativistic potential energy curves that correspond to the Rb (5s (2)S) + Sr (5s4d (3)P°) asymptote, we have simulated dispersed fluorescence spectra as they were recently measured in our lab. The comparison with experimental data allows to benchmark both methods and demonstrate that spin-orbit coupling has to be included for the lowest states of RbSr.