A theoretical study on the stability difference of the borane and carborane C2Bn−2Hn (5 ≤ n ≤ 7) clusters

In order to study the electronic structure and structural stability of borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters, especially the stability difference between the borane and carborane C₂B₃H₅. The frontier orbital energy levels of the borane and carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters are calculated at CCSD(T)/aug‐cc‐pVXZ//B3LYP/def2‐TZVPP level. The results are further analyzed by qualitative frontier orbital method based on the cap–ring interaction. The results reveal that: (1) the larger E_gap(HOMO‐LUMO energy gap) of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than borane (5 ≤ n ≤ 7) clusters originates from the more effective cap–ring orbital overlap of carborane C₂B_n?2H_n (5 ≤ n ≤ 7) clusters than that of borane (5 ≤ n ≤ 7) clusters; (2) the smallest E_gap of the borane results from the highest energy level of the ring symmetry‐adapted linear combination orbital of cluster; and (3) the largest E_gap of the carborane C₂B₃H₅ is induced by the most effective cap–ring orbital interaction of C₂B₃H₅ cluster. © 2014 Wiley Periodicals, Inc.     

Theoretical prediction of the optical rotation of chiral molecules in ordered media: A computational study of (Ra)‐1,3‐dimethylallene, (2R)‐2‐methyloxirane,and (2R)‐N‐methyloxaziridine

A theoretical procedure has been developed and implemented to calculate the optical rotation of chiral molecules in ordered phase via origin‐independent diagonal components , of the optical activity tensor and origin‐independent components , for , of the mixed electric dipole‐electric quadrupole polarizability. Origin independence was achieved by referring these tensors to the principal axis system of the electric dipole dynamic polarizability at the same laser frequency ω. The approach has been applied, allowing for alternative quantum mechanical methods based on different gauges, to estimate near Hartree–Fock values for three chiral molecules, (2R)‐N‐methyloxaziridine C₂NOH₅, (2R)‐2‐methyloxirane (also referred to as propylene oxide) C₃OH₆, and ( )‐1,3‐dimethylallene C₅H₈, at two frequencies. The theoretical predictions can be useful for an attempt at measuring correspondent experimental values in crystal phase. © 2015 Wiley Periodicals, Inc.     

General build up of basis and matrix in the diagonalization approach. Determination of Kramers configuration state functions

Algorithms to build the basis and matrix representation to obtain the Kramers configuration space functions (KCSFs) via diagonalization will be formally generalized to an arbitrary number of unpaired (open shell) fermions. Effective build up of the matrix representation will be outlined (including threading and graphical processing unit parallelism) to subsequently obtain the KCSFs via calling external/numerical library routines for diagonalization. The effective build up of the matrix representation relays on a binary tree search algorithm to allow evaluation the action on a given basis vector. The binary tree search avoids the treatment of zero matrix elements which leads to an exponential acceleration. The implementation ( basis creation, matrix representation, and matrix diagonalization) will be done in an all in core and all at once manner, hence the available core memory sets the physical limits in practical applications. Memory limitations, sparsity of the matrix, general case of n fermions in m spinors, and the application of KCSFs will be put into further perspective.     

(Rg = He ∼ Rn,n = 1–4): In quest of the potential trapping ability of the aromatic ring

A new series of divalent boron‐rare gas cations (Rg = He ∼ Rn, n = 1–4) have been predicted theoretically at the B3LYP, MP2, and CCSD(T) levels to present the structures, stability, charge distributions, bond natures, and aromaticity. The Rg B bond energies are quite large for heavy rare gases and increase with the size of the Rg atom. Because of steric hindrance new Rg atoms introduced to the B₄ ring will weaken the Rg B bond. Thus in the Rg B bond has the largest binding energy 90–100 kcal/mol. p‐ has a slightly shorter Rg B bond length and a larger bond energy than o‐ . NBO and AIM analyses indicate that for the heavy Rg atoms Ar ∼ Rn the B Rg bonds have character of typical covalent bonds. The energy decomposition analysis shows that the σ‐donation from rare gases to the boron ring is the major contribution to the Rg B bonding. Adaptive natural density partitioning and nuclear‐independent chemical shift analyses suggest that both and have obvious aromaticity.     

Carbon–sulfur chains: A high‐resolution infrared and quantum‐chemical study of C3S and SC7S

In the course of a 5 μm high‐resolution infrared study of laser ablation products from carbon–sulfur targets, the ν₁ vibrational mode region of linear C₃S has been studied continuously from 2046 to 2065 cm^?1. Besides the prominent vibrational fundamental, the region was found to feature the , and even hot bands, the latter two of which were observed for the first time. Owing to the high signal‐to‐noise ratio obtained, the ν₁ mode of S could also be observed in natural abundance for the first time at high spectral resolution in the infrared. At 2061 cm^?1, hidden inside the branch of the C₃S ν₁ fundamental mode, a weak new band is observed which exhibits very tight line spacing and stems from a heavy both carbon and sulfur containing carrier. On the basis of high‐level quantum‐chemical calculations of selected carbon–sulfur chains and other carbon‐rich cumulenes, this feature is attributed to the ν₅ vibrational fundamental of linear SC₇S, which stands for the first gas‐phase spectroscopic detection of this long cumulenic chain.     

Quantum control of electron‐proton symmetry breaking in dissociative ionization of H2 by intense laser pulses

Forward and backward electron/proton ionization/dissociation spectra from one‐dimensional non‐Born‐Oppenheimer H₂ molecule exposed to ultrashort intense laser pulses ( W/cm², λ = 800 nm) have been computed by numerically solving the time‐dependent Schrödinger equation. The resulting above‐threshold ionization and above‐threshold dissociation spectra exhibit the characteristic forward‐backward asymmetry and sensitivity to the carrier‐envelope phase (CEP), particularly for high energies. A general framework for understanding CEP effects in the asymmetry of dissociative ionization of H₂ has been established. It is found that the symmetry breaking of electron‐proton distribution with π periodic modulation occurs for all CEPs except for ( integer) and the largest asymmetry coming from the CEP of . At least one of the electron and proton distributions is asymmetric when measured simultaneously. Inspection of the nuclear and electron wave packet dynamics provides further information about the relative contribution of the gerade and ungerade states of to the dissociation channel and the time delay of electrons in asymmetric ionization. © 2014 Wiley Periodicals, Inc.     

A theoretical study on cyclacenes: Analytical tight‐binding approach

We present a theoretical study of cyclacene molecules performed at tight‐binding level. The orbital energies and eigenvectors have been analytically computed, and exact expressions for the axial component of the total position spread and polarizability tensors have been obtained. In absence of dimerization, the system has a D_nh symmetry, where n is the number of hexagonal units. The energy bands present no gap at the Fermi level, and to this fact it corresponds a diverging (per‐electron) polarizability for in the direction of the system symmetry axis. The two (degenerate) components of the polarizability on the σ_h symmetry plane, conversely, remain finite for . The total position spread tensor presents a qualitatively different behavior, since all the three components of the position spread per electron remain finite for . The results are analyzed and discussed for both axial and planar components separately as these are affected differently with respect to the increasing system size. Both dipole polarizability and total position spread have been computed using an ab initio approach for the smallest systems, to compare the analytical tight‐binding expressions with a higher‐level theory.     

Position and momentum information‐theoretic measures of the pseudoharmonic potential

In this study, the information‐theoretic measures in both the position and momentum spaces for the pseudoharmonic potential using Fisher information, Shannon entropy, Renyi entropy, Tsallis entropy, and Onicescu information energy are investigated analytically and numerically. The results obtained are applied to some diatomic molecules. The Renyi and Tsallis entropies are analytically obtained in position space using Srivastava–Niukkanen linearization formula in terms of the Lauricella hypergeometric function. Also, they are obtained in the momentum space in terms of the multivariate Bell polynomials of Combinatorics. We observed that the Fisher information increases with n in both the position and momentum spaces, but decreases with for all the diatomic molecules considered. The Shannon entropy also increases with increasing n in the position space and decreases with increasing . The variations of the Renyi and Tsallis entropies with are also discussed. The exact and numerical values of the Onicescu information energy are also obtained, after which the ratio of information‐theoretic impetuses to lengths for Fisher, Shannon, and Renyi are obtained. © 2015 Wiley Periodicals, Inc.     

Full–dimensional quantum molecular dynamics calculations of the rovibrationally mediated X → 2 transition of nitrous oxide

A full dimensional time‐dependent quantum wavepacket approach is used to study the photodissociation dynamics of nitrous oxide for the X → 2 bound–bound transition based on new highly accurate potential energy and transition dipole moment surfaces. The computed 2 absorption spectra at room temperature are characterized by sharp vibrational structures that contribute slightly to the diffuse vibrational structures around the maximum peak at 180 nm of the first ultraviolet absorption band (from the contribution of 2 , 1 , and 2 states) of N₂O. Transitions from different initial rovibrational states reveal that the sharp structures arise mainly from N₂? O bending vibrations, whereas, at higher temperatures, the N₂? O and N? NO stretching vibrations are responsible for enhancing the intensity of the structures. At absorption wavelengths 166 nm and 179 nm, vibrational quantum state distributions of N₂ product fragments decrease monotonically with increasing vibrational quantum number v = 0, 1, 2. At 166 nm, rotational quantum state distributions of N₂ at fixed v = 0 and v = 1 display multimodal profiles with maximum peaks at j = 77 and j = 75, respectively, whereas, the distributions at the 179 nm absorption wavelength display bimodal profiles with maximum peaks at j = 73 and j = 71, respectively. Accordingly, the presence of rotationally hot N₂ from previous experimental and theoretical works in the first band strongly implies a significant influence of the 2 state in determining the final dissociation pathway of N₂ + O. © 2016 Wiley Periodicals, Inc.     

Einstein A Coefficients,Oscillator Strengths,and Absolute Band Strengths for the First Negative System,and Franck-Condon Factors for the Second Negative System of

The Franck-Condon factors and r centroids for the first negative and the second negative band systems of the molecule, based on the Rydberg-Klein-Rees potential energy curves, have been computed. The variation of the electronic transition moment with the internuclear separation has been studied for the first negative bands and the Einstein A coefficients, the oscillator strengths, and the absolute band strengths for this system have been calculated by adopting the recent experimental data on the lifetimes of the various levels of the excited B state. On a calculé pour la molécule les facteurs de Franck-Condon et les r centroïds pour les deux premiers systèmes de bandes négatives avec des courbes de potential de Rydberg-Klein-Rees. On a étudié la variation du moment de transition éléctronique avec la séparation internucléaire pour les premières bandes négatives. Les constantes A d'Einstein, les forces d'oscillateur et les intensités absolues des bandes ont été calculé des données expérimentales récentes sur les durées de vie des niveaux différents de l'état excité B . Die Franck-Condon-Faktoren und die r-Zentroide der zwei ersten negativen Banden-systeme des -Moleküls, wurden mit den Rydberg-Klein-kees Potentialkurven berechnet. Die Variation des elektronischen Übergangsmoment mit dem Kernabstand würde für die ersten negativen Bänder studiert. Die Einsteinschen A-Koeffizienten, die Oszillatorstärken and die absoluten Bandenstärken wurden mit die neuen experimentellen Tatsachen über die Lebensdauern der verschiedenen Niveaus des angeregten B Zustands berechnet.     

Energies for cyclic and acyclic aggregations of adamantane and diamantane units sharing vertices,edges, or six‐membered rings

Diamondoids are hydrocarbons having a carbon scaffold comprised from polymer‐like composites of adamantane cages. This article describes computed total energies and “SWB‐tension” energies (often referred to as “strain” energies) for species having n adamantane or diamantane units sharing pairwise: one carbon atom (spiro‐[n]adamantane or spiro‐[n]diamantane); one C? C bond (one‐bond‐sharing‐[n]adamantane or one‐bond‐sharing‐[n]diamantane); or one chair‐shaped hexagon of carbon atoms (1234‐helical‐cata‐[n]diamantanes). Each of the five investigated polymer‐like types is considered either as an acyclic or a cyclic chain of adamantane‐ or diamantane‐unit cages. With increasing n values, SWB‐tension energies for acyclic aggregates are found to increase linearly, while the net SWB‐tension energies of cyclic aggregates often go thru a minimum at a suitable value of . In all five cases, a limiting common energy per unit ( ) is found to be approached by both cyclic and acyclic chains as , as revealed from plots of versus 1/n for acyclic chains and of versus 1/n² for cyclic chains. © 2015 Wiley Periodicals, Inc.     

Structural,magnetic, and vibrational properties of stoichiometric clusters of CrN

We performed density‐functional‐theoretic calculations to investigate the structural, magnetic and vibrational properties of the stoichiometric clusters (CrN)_n ( ). We show that the building block of the ground‐state structures of these clusters is a square CrNCrN unit; the only exception with n > 2 occurs for (CrN)₃, but this cluster has an isomer not very far in energy from the ground state consisting of a pair of CrNCrN squares sharing a CrN bond. In the smaller CrN, (CrN)₂, and (CrN)₃ clusters the magnetic moments of the N atoms are non‐negligible and antiparallel to those of the Cr atoms, but for the larger species (CrN)₄, (CrN)₅, (CrN)₆, and (CrN)₉ the cluster magnetic moments are almost entirely due to the Cr atoms. Lack of imaginary vibrational frequencies in the predicted ground‐state structures of (CrN)_n ( ) confirms that they are mechanically stable equilibrium states. © 2015 Wiley Periodicals, Inc.     

Lanthanide metals in the boron cages: Computational prediction of M@Bn (M = Eu,Gd; n = 38, 40)

Endohedral metalloborofullerenes (EMBFs) are novel boron analogues of the famous endohedral metallofullerenes (EMFs). Many EMBFs have been proposed by theoretical calculations thus far. However, in sharp contrast to EMFs, which trap most of the lanthanides with f electrons inside the cages, the corresponding lanthanide‐based EMBFs have never been reported. In this work, the encapsulation of Eu and Gd in the B₃₈ and B₄₀ fullerenes was studied by means of density functional theory calculations. Our results revealed that Gd@B₃₈(⁹A), Eu@B₄₀(⁸B₂), and Gd@B₄₀(⁷A″) all favor the endohedral configuration, and the electronic structures can be described as Gd³⁺@ , Eu²⁺@ , and Gd³⁺@ with jailed f electron spins. The large binding energies and sizable HOMO–LUMO gaps suggest that they may be achieved experimentally. They feature σ and π double aromaticity, and their excellent stabilities were confirmed by the Born–Oppenheimer molecular dynamics simulations. Finally, the infrared and UV/vis spectra were simulated to assist experimental characterization.     

The nature of the chemical bond in borazine (B3N3H6), boroxine (B3O3H3), carborazine (B2N2C2H6), and related species

The bonding problem in borazine (B₃N₃H₆), boroxine (B₃O₃H₃), and carborazine (B₂N₂C₂H₆) is successfully addressed through the consideration of the excited states of the constituent fragments, namely BH( ), NH( ), and CH( ). We propose the participation of resonant structures for all three species that help to explain the experimental findings. A discussion on the chemical pattern of the parental molecule benzene (C₆H₆) helps to make coherent the whole bonding analysis on the titled species.     

Design of neutral organic superacids using fulvene derivatives with di‐enol substituent

A new class of superacids was designed using enolic derivatives of fulvene. After deprotonation, bond rearrangement leads in a stable conjugate base with an aromatic cyclopentadienyl ring and a carboxyl group. The gas phase enthalpies ( ) of the deprotonation, as an index of acidity, were calculated employing the B3LYP method and 6‐311++G(d,p) and aug‐cc‐PVDZ basis sets. The acidity of these compounds without any electron withdrawing groups was more than H₂SO₄ in gas phase. The acidity increased by substituting electron withdrawing groups (? F, ? CN, and ?O) into the molecules so that we could achieve a cyano derivative of fulvene with = 250 kcal/mol.     

Additivity in kramers pairs symmetry: Multiplets with up to four unpaired electrons

Many fermions Kramers pairs formalism is considered from the prospective of the sum of individual single fermion time‐reversal operators. The obtained many fermions “pseudo Kramers pairs operator” ( ), as well as its square ( ), have formally the same structure as the many fermion spin operators and . Nevertheless, the shape of eigenfunctions with respect to and is different. Herein all Kramers adapted eigenfunctions of for cases of up to four unpaired fermions are compiled, and their properties with respect to further advocated. It will be shown that degeneracy of the multiplets recovers the proper behavior with respect to Pascal's triangle. A projection operator for obtaining the “high spin” Kramers adapted eigenfunctions is suggested. Noncommutation of with spin and angular momentum operators and degeneracy is discussed at last. © 2016 Wiley Periodicals, Inc.     

Explicitly correlated variational estimates of the energy levels of negative hydrogen ion under spatial confinement

Comprehensive investigations on the structural modifications of negative hydrogen ion within an impenetrable spherical domain has been performed in the framework of Ritz variational method. Electron correlation plays a major role in the formation of H^– ion. The Hylleraas‐type basis set expansion of wave function considered here incorporates the effect of electron correlation in an explicit manner. Energy values of and 1sn states of H^– ion within confined domain have been calculated. Although the singly excited states do not exist for a “free” H^– ion, well converged energy values of such states have been found within a wide range of confinement radius. The thermodynamic pressure felt by the ion inside the sphere is also estimated. The general trend shows successive destabilization of the excited energy levels with increase of pressure. The contribution of angular correlation in the energy values have been estimated. Evolution of and energy levels of H^– ion as quasi‐bound states are being reported.     

Rovibrational spectra of bounded diatomic molecules

Spectra of a bounded diatomic molecule is studied numerically. Shifted Deng–Fan oscillator potential has been used to model the molecule. The accurate five‐point finite difference method has been used to solve the Schrödinger equation for rovibrational motion of the molecule. The energies of the bound states as well as free states of the molecule have been calculated. In addition, radial matrix elements like , n = 1, 2, and 3 have been calculated. These have been used to calculate the ‐pole static polarizabilities. The variation of bound state energies, matrix elements and ‐pole static polarizabilities with the boundary radius has also been studied. The Stark effect in case of this bounded system has also been investigated.     

Quantum chemical prediction for complex organic molecules

Numerous types of quantum chemical calculations and protocols have been successfully applied to computing of small, uncomplicated organic molecules. Here, we argue for the need to shift attention to more challenging molecules that are marked by an interplay of complicating factors such as conformational, tautomeric, steric, and other effects. The challenge is not in choosing the right quantum chemical method and solvation model but in combining the existing methods to simultaneously and accurately describe the breadth of chemical and physical phenomena that give rise to the experimentally observed . The complexity of the phenomena that must be considered begs for the need for a greater automation of prediction workflows. We review our experience with these challenges and outline paths for future progress in the direction of tackling prediction of complex organic molecules.     

Geometric and electronic structures of silicon fluorides (N = 4 – 6) and potential energy surfaces for dissociation reactions → SiF4 + F– and → + F–

The geometric and electronic structures of a series of silicon fluorides (n = 4 ? 6) were computationally studied with the aid of density functional theory (DFT) method with B3LYP and M06‐2X functionals and coupled cluster (CCSD and CCSD(T)) methods with 6‐311++G(d,p) basis set. The nature of the Si‐F bonds in these compounds was analyzed in the framework of the natural bond orbital theory and natural resonance theory. Energy characteristics (heats of reactions and energy barriers) of the dissociation reactions → SiF₄ + F^– and → + F^– were calculated using the DFT and CCSD methods. The potential energy surface of elimination of a fluoride anion from has a specific topology with valley‐ridge inflection points corresponding to bifurcations of the minimal energy reaction path. © 2016 Wiley Periodicals, Inc.