Contracted Gaussian-type basis functions revisited II. Atoms Na through Ar

Mass-selected ion-molecule cluster beam apparatus for ultrafast photofragmentation studies

We describe an apparatus for investigating the excited-state dissociation dynamics of mass-selected ion-molecule clusters by mass-resolving and detecting photofragment-ions and neutrals, in coincidence, using an ultrafast laser operating at high repetition rates. The apparatus comprises a source that generates ion-molecule clusters, a time-of-flight spectrometer, and a mass filter that selects the desired anions, and a linear-plus-quadratic reflectron mass spectrometer that discriminates the fragment anions after the femtosecond laser excites the clusters. The fragment neutrals and anions are then captured by two channeltron detectors. The apparatus performance is tested by measuring the photofragments: I-, CF3I-, and neutrals from photoexcitation of the ion-molecule cluster CF3I·I- using femtosecond UV laser pulses with a wavelength of 266 nm. The experimental results are compared with our ground state and excited state electronic structure calculations as well as the existing results and calculations, with particular attention to the generation mechanism of the anion fragments and dissociation channels of the ion-molecule cluster CF3I·I- in the charge-transfer excited state.

Relativistic nonorthogonal configuration interaction: application to L2,3-edge X-ray spectroscopy

In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF₆, SiCl₄ and [FeCl₆]^3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L_2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.

Generation, contraction, and polarisation of Gaussian basis sets for atomic and molecular calculations using the generator coordinate method with polynomial discretisation: atoms from Na through Cl

The polynomial Generator Coordinate Hartree-Fock Gaussian basis sets, pGCHF, for the atoms Na, Mg, Al, Si, P, S, and Cl were generated using the generator coordinate method based on polynomial integral expansion to discretise the Griffin-Wheeler-Hartree-Fock equations. The pGCHF basis sets were contracted with the CONTRACT program based on the Davidson contraction model through which a set of 9s8p functions for the atoms Na through Cl were obtained. Polarisation exponents generated using the POLARIZATION program were added to the contracted pGCHF Gaussian basis sets. Molecular calculations at the DFT level of theory showed that the pGCHF basis sets can be used to calculate the atomisation energy with accuracy comparable to the well-established pcseg-3, def2-QZVP, and Sapporo-QZP basis sets; also, the complete basis set (CBS) limit estimate was obtained with the pcseg-3/pcseg-4 basis sets.

Analysis of the relative stability of lithium halide crystal structures: Density functional theory and classical models

All lithium halides exist in the rock salt crystal structure under ambient conditions. In contrast, common lithium halide classical force fields more often predict wurtzite as the stable structure. This failure of classical models severely limits their range of application in molecular simulations of crystal nucleation and growth. Employing high accuracy density functional theory (DFT) together with classical models, we examine the relative stability of seven candidate crystal structures for lithium halides. We give a detailed examination of the influence of DFT inputs, including the exchange-correlation functional, basis set, and dispersion correction. We show that a high-accuracy basis set, along with an accurate description of dispersion, is necessary to ensure prediction of the correct rock salt structure, with lattice energies in good agreement with the experiment. We also find excellent agreement between the DFT-calculated rock salt lattice parameters and experiment when using the TMTPSS-rVV10 exchange-correlation functional and a large basis set. Detailed analysis shows that dispersion interactions play a key role in the stability of rock salt over closely competing structures. Hartree-Fock calculations, where dispersion interactions are absent, predict the rock salt structure only for LiF, while LiCl, LiBr, and LiI are more stable as wurtzite crystals, consistent with radius ratio rules. Anion-anion second shell dispersion interactions overcome the radius ratio rules to tip the structural balance to rock salt. We show that classical models can be made qualitatively correct in their structural predictions by simply scaling up the pairwise additive dispersion terms, indicating a pathway toward better lithium halide force fields.

Unifying General and Segmented Contracted Basis Sets. Segmented Polarization Consistent Basis Sets

We propose a method, denoted P-orthogonalization, for converting a general contracted basis set to a computationally more efficient segmented contracted basis set, while inheriting the full accuracy of the general contracted basis set. The procedure can be used for any general contracted basis set to remove the redundancies between general contracted functions in terms of primitive functions. The P-orthogonalization procedure is used to construct a segmented contracted version of the polarization consistent basis sets, which are optimized for density functional theory calculations. Benchmark calculations show that the new pcs-n basis sets provide uniform error control of the basis set incompleteness for molecular systems composed of atoms from the first three rows in the periodic table (H-Kr) and for different exchange-correlation functionals. The basis set errors at a given zeta quality level are lower than other existing basis sets, and the pcs-n basis sets are furthermore shown to be among the computationally most efficient. The pcs-n basis sets are available in qualities ranging from (unpolarized) double-zeta to pentuple zeta quality and should therefore be well suited for both routine and benchmark calculations using density functional theory methods in general.

The nature of the interaction of dimethylselenide with IIIA group element compounds

The first systematic theoretical study of the nature of intermolecular bonding of dimethylselenide as donor and IIIA group element halides as acceptors was made with the help of the approach of Quantum Theory of Atoms in Molecules. Density Functional Theory with "old" Sapporo triple-ζ basis sets was used to calculate geometry, thermodynamics, and wave function of Me2Se···AX3 complexes. The analysis of the electron density distribution and the Laplacian of the electron density allowed us to reveal and explain the tendencies in the influence of the central atom (A = B, Al, Ga, In) and halogen (X = F, Cl, Br, I) on the nature of Se···A bonding. Significant changes in properties of the selenium lone pair upon complexation were described by means of the analysis of the Laplacian of the charge density. Charge transfer characteristics and the contributions to it from electron localization and delocalization were analyzed in terms of localization and delocalization indexes. Common features of the complexation and differences in the nature of bonding were revealed. Performed analysis evidenced that gallium and indium halide complexes can be attributed to charge transfer-driven complexes; aluminum halides complexes seem to be mainly of an electrostatic nature. The nature of bonding in different boron halides essentially varies; these complexes are stabilized mainly by covalent Se···B interaction. In all the complexes under study covalence of the Se···A interaction is rather high.

Co(III) complexes with N2(SO)2-type equatorial planar ligands similar to the active center of nitrile hydratase: role of the sulfenate group in the enzymatic reaction

In order to gain an understanding of the role of the sulfenyl group of nitrile hydratase (NHase), a new Co(III) complex with a sulfenyl-type ligand (LC=O:N2(SO)2), Na[CoIII(LC=O:N2(SO)2)(tBuNC)2] (2), was synthesized. The compound includes two amide groups, two sulfenate sulfurs in the equatorial plane, and two tBuNC molecules in the axial positions. Characterization of the compound was performed by UV-vis spectroscopic, IR spectral, thermogravimetric (TG), and X-ray structure analytical methods. The results are discussed in the context of Co(III) complexes containing the corresponding sulfur-type (LC=O:N2S2) (1) and sulfinyl-type ligands (LC=O:N2(SO2)2) (3). Complex 2 crystallized with the formula Na[CoIII(LC=O:N2(SO)2)(tBuNC)2].urea.2H2O.0.5EtOH. The X-ray structure revealed that the Co(III) complex has an octahedral geometry with Co-S=av. 2.221 A, Co-N=av. 1.998 A, and Co-C=av. 1.87 A. The sulfenyl oxygen and amidate carbonyl oxygen are linked to urea, water, EtOH, and Na+ and participate in a hydrogen-bond and an electrostatic interaction. IR and TG measurements demonstrated that the coordination strength of tBuNC to the Co atom increases as follows: 1<2<3. Complex 2 has almost the same stability as 3 in all solutions tested, although 1 exhibits a release of axial ligands in nonaqueous solutions. DFT calculations for 1, 2, and 3 demonstrated that Milliken atomic charges of the Co(III) centers are +1.466, +1.536, and +1.542, respectively, indicating that the extent of oxidation of the sulfur atoms increases the Lewis acidity of the Co(III) centers. Interestingly, the solution-state IR spectrum of 2 exhibits a solvent-dependent S-O stretching frequency. The frequency decreases with an increase in the electrophilicity (acceptor number) of the solvent. This solvent dependence was not observed for 3, which has a sulfinate (SO2) group, suggesting that the sulfenyl oxygen atom has nucleophilic character and promotes strong binding of the tBuNC molecule to lower the reaction barrier. These findings may suggest that the sulfenate oxygen in native NHase acts as a base (proton acceptor) and contributes to the activation of a water molecule and/or nitrile molecule.

Quality of contracted Gaussian-type function basis sets

The valence quality of contracted (C) Gaussian-type function (GTF) basis sets in molecular calculations is discussed for the first- through fourth-row atoms. The split-valence basis sets derived from minimal-type CGTF sets are compared with those derived from primitive (P) GTF sets. Using F, Cl, Br, and I atoms and their homonuclear diatomics as test species, we find that the split-valence CGTF sets have almost the same quality as PGTF sets with larger s and p expansion terms: for example, the (53/5), (533/53), (5333/533/5), and (53 333/5333/53) CGTF sets correspond approximately to the [9/5], [15/9], [19/15/5], and [22/18/7] PGTF sets for the first- to fourth-row atoms, respectively, where the slash separates the s, p, and d symmetries. For the main group atoms of the four rows, we recommend using the above-mentioned CGTFs or larger.

Compact and efficient basis sets of s- and p-block elements for model core potential method

We propose compact and efficient valence-function sets for s- and p-block elements from Li to Rn to appropriately describe valence correlation in model core potential (MCP) calculations. The basis sets are generated by a combination of split MCP valence orbitals and correlating contracted Gaussian-type functions in a segmented form. We provide three types of basis sets. They are referred to as MCP-dzp, MCP-tzp, and MCP-qzp, since they have the quality comparable with all-electron correlation consistent basis sets, cc-pVDZ, cc-pVTZ, and cc-pVQZ, respectively, for lighter atoms. MCP calculations with the present basis sets give atomic correlation energies in good agreement with all-electron calculations. The present MCP basis sets systematically improve physical properties in atomic and molecular systems in a series of MCP-dzp, MCP-tzp, and MCP-qzp. Ionization potentials and electron affinities of halogen atoms as well as molecular spectroscopic constants calculated by the best MCP set are in good agreement with experimental values.