Spin–orbit effects on magnetically induced current densities in the () clusters

This study reports the spin–orbit effects on the aromaticity of the , , , , , and anionic clusters via the magnetically induced current‐density method. All‐electron density functional theory (DFT) calculations were carried out using the four‐component Dirac‐Coulomb (DC) hamiltonian, including scalar and spin–orbit relativistic effects. The magnetic index of aromaticity was calculated by numerical integration over the current flow between two atoms in the pentagonal ring. These values were compared to the spin‐free values (spin–orbit coupling switched off), in order to assess the spin–orbit effect on aromaticity. It was found that in the heavy anions, and , there is a significant influence of the spin–orbit coupling. © 2018 Wiley Periodicals, Inc.     

Accurate,robust, and reliable calculations of Poisson–Boltzmann binding energies

Poisson–Boltzmann (PB) model is one of the most popular implicit solvent models in biophysical modeling and computation. The ability of providing accurate and reliable PB estimation of electrostatic solvation free energy, , and binding free energy, , is important to computational biophysics and biochemistry. In this work, we investigate the grid dependence of our PB solver (MIBPB) with solvent excluded surfaces for estimating both electrostatic solvation free energies and electrostatic binding free energies. It is found that the relative absolute error of obtained at the grid spacing of 1.0 Å compared to at 0.2 Å averaged over 153 molecules is less than 0.2%. Our results indicate that the use of grid spacing 0.6 Å ensures accuracy and reliability in calculation. In fact, the grid spacing of 1.1 Å appears to deliver adequate accuracy for high throughput screening. © 2017 Wiley Periodicals, Inc.     

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.     

Nuclear quantum effects induce metallization of dense solid molecular hydrogen

We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band‐gap energies of the , Pc, and structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab initio path‐integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero‐point energy and temperature‐induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band‐gaps substantially, leading to metallization of the and Pc phases via band overlap; the effect on the band‐gap of the structure is less pronounced. Our combined DMC‐PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the and Pc phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band‐gap energy and vibron modes. This strong coupling induces a band‐gap reduction of more than 2.46 eV in high‐pressure solid molecular hydrogen. Comparing our DMC‐PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen. © 2017 Wiley Periodicals, Inc.     

Ab initio potential energy surface and vibration–rotation energy levels of germanium dicarbide,GeC2

The accurate ground‐state potential energy surface of germanium dicarbide, GeC₂, has been determined from ab initio calculations using the coupled‐cluster approach. The core–electron correlation, higher‐order valence‐electron correlation, and scalar relativistic effects were taken into account. The potential energy surface of GeC₂ was shown to be extraordinarily flat near the T‐shaped equilibrium configuration. The potential energy barrier to the linear CCGe configuration was predicted to be 1218 cm⁻¹. The vibration–rotation energy levels of some GeC₂ isotopologues were calculated using a variational method. The vibrational bending mode ν₃ was found to be highly anharmonic, with the fundamental wavenumber being only 58 cm⁻¹. Vibrational progressions due to this mode were predicted for the , , and states of GeC₂. © 2018 Wiley Periodicals, Inc.     

Structure and electrochemical properties for complexes of nitrocompounds with inorganic ions: A theoretical approach

Reduction and oxidation (redox) reactions are widely used for removal of nitrocompounds from contaminated soil and water. Structures and redox properties for complexes of nitrocompounds, such as 2,4,6‐trinitrotoluene (TNT), 2,4‐dinitrotoluene (DNT), 2,4‐dinitroanisole (DNAN), and 5‐nitro‐2,4‐dihydro‐3H?1,2,4‐triazol‐3‐one (NTO), with common inorganic ions (Na⁺, Cl^?, ) were investigated at the SMD(Pauling)/PCM(Pauling)/MPWB1K/TZVP level of theory. Atoms in molecules (AIM) theory was applied to analyze the topological properties of the bond critical points involved in the interactions between the nitrocompounds and the ions. Topological analyses show that intermolecular interactions of the types O(N)…Na⁺, C‐H…Cl^?( ), and C…Cl^?( ) may be discussed as noncovalent closed‐shell interactions, while N‐H···Cl^?( ) hydrogen bonds are partially covalent in nature. Complexation causes significant decrease of redox activity of the nitrocompounds. Analysis of the reduction potentials of the complexes obtained through application of the Pourbaix diagram of an iron/water system revealed that sodium complexes of NTO might be reduced by metallic iron. © 2016 Wiley Periodicals, Inc.     

Computational study of basis set and electron correlation effects on anapole magnetizabilities of chiral molecules

In the presence of a static, nonhomogeneous magnetic field, represented by the axial vector at the origin of the coordinate system and by the polar vector , assumed to be spatially uniform, the chiral molecules investigated in this paper carry an orbital electronic anapole, described by the polar vector . The electronic interaction energy of these molecules in nonordered media is a cross term, coupling and via , one third of the trace of the anapole magnetizability a_αβ tensor, that is, . Both and W^BC have opposite sign in the two enantiomeric forms, a fact quite remarkable from the conceptual point of view. The magnitude of predicted in the present computational investigation for five chiral molecules is very small and significantly biased by electron correlation contributions, estimated at the density functional level via three different functionals. © 2016 Wiley Periodicals, Inc.     

Exploring possible reaction pathways for the o‐atom transfer reactions to unsaturated substrates catalyzed by a [Ni‐NO2] ↔ [Ni‐NO] redox couple using DFT methods

The (nitro)(N‐methyldithiocarbamato)(trimethylphospane)nickel(II), [Ni(NO₂)(S₂CNHMe)(PMe₃)] complex catalyses efficiently the O‐atom transfer reactions to CO and acetylene. Energetically feasible sequence of elementary steps involved in the catalytic cycle of the air oxidation of CO and acetylene are proposed promoted by the Ni(NO₂)(S₂CNHMe)(PMe₃)] ↔ Ni(NO₂)(S₂CNHMe)(PMe₃) redox couple using DFT methods both in vacuum and dichloromethane solutions. The catalytic air oxidation of HC≡CH involves formation of a five‐member metallacycle intermediate, via a [3 + 2] cyclo‐addition reaction of HC≡CH to the Ni‐N = O moiety of the Ni(NO₂)(S₂CNHMe)(PMe₃)] complex, followed by a β H‐atom migration toward the C_α carbon atom of the coordinated acetylene and release of the oxidation product (ketene). The geometric and energetic reaction profile for the reversible [Ni( ‐NO₂)(S₂CNHMe)(PMe₃)] [Ni( ‐ONO)(S₂CNHMe)(PMe₃)] linkage isomerization has also been modeled by DFT calculations. © 2017 Wiley Periodicals, Inc.     

Development of spin‐orbit coupling for stochastic configuration interaction techniques

To perform spin‐orbit coupling calculations on atoms and molecules, good zeroth‐order wavefunctions are necessary. Here, we present the software development of the Monte Carlo Configuration Interaction (MCCI) method, to enable calculation of such properties, where MCCI iteratively constructs a multireference wavefunction using a stochastic procedure. In this initial work, we aim to establish the efficacy of this technique in predicting the splitting of otherwise degenerate energy levels on a range of atoms and small diatomic molecules. It is hoped that this work will subsequently act as a gateway toward using this method to investigate singlet‐triplet interactions in larger multireference molecules. We show that MCCI can generate very good results using highly compact wavefunctions compared to other techniques, with no prior knowledge of important orbitals. Higher‐order relativistic effects are neglected and spin‐orbit coupling effects are incorporated using first‐order degenerate perturbation theory with the Breit‐Pauli Hamiltonian and effective nuclear charges in the one‐electron operator. Results are obtained and presented for B, C, O, F, Si, S, and Cl atoms and OH, CN, NO, and C₂ diatomic radicals including spin‐orbit coupling constants and the relative splitting of the lowest energy degenerate state for each species. Convergence of MCCI to the full configuration interaction result is demonstrated on the multireference problem of stretched OH. We also present results from the singlet‐triplet interaction between the and both the and states of the O₂ molecule. © 2017 Wiley Periodicals, Inc.     

Calculation of photoionization differential cross sections using complex Gauss‐type orbitals

Accurate theoretical calculation of photoelectron angular distributions for general molecules is becoming an important tool to image various chemical reactions in real time. We show in this article that not only photoionization total cross sections but also photoelectron angular distributions can be accurately calculated using complex Gauss‐type orbital (cGTO) basis functions. Our method can be easily combined with existing quantum chemistry techniques including electron correlation effects, and applied to various molecules. The so‐called two‐potential formula is applied to represent the transition dipole moment from an initial bound state to a final continuum state in the molecular coordinate frame. The two required continuum functions, the zeroth‐order final continuum state and the first‐order wave function induced by the photon field, have been variationally obtained using the complex basis function method with a mixture of appropriate cGTOs and conventional real Gauss‐type orbitals (GTOs) to represent the continuum orbitals as well as the remaining bound orbitals. The complex orbital exponents of the cGTOs are optimized by fitting to the outgoing Coulomb functions. The efficiency of the current method is demonstrated through the calculations of the asymmetry parameters and molecular‐frame photoelectron angular distributions of and . In the calculations of , the static exchange and random phase approximations are employed, and the dependence of the results on the basis functions is discussed. © 2017 Wiley Periodicals, Inc.     

(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.     

Low‐energy structures of benzene clusters with a novel accurate potential surface

The benzene‐benzene (Bz‐Bz) interaction is present in several chemical systems and it is known to be crucial in understanding the specificity of important biological phenomena. In this work, we propose a novel Bz‐Bz analytical potential energy surface which is fine‐tuned on accurate ab initio calculations in order to improve its reliability. Once the Bz‐Bz interaction is modeled, an analytical function for the energy of the clusters may be obtained by summing up over all pair potentials. We apply an evolutionary algorithm (EA) to discover the lowest‐energy structures of clusters (for ), and the results are compared with previous global optimization studies where different potential functions were employed. Besides the global minimum, the EA also gives the structures of other low‐lying isomers ranked by the corresponding energy. Additional ab initio calculations are carried out for the low‐lying isomers of and clusters, and the global minimum is confirmed as the most stable structure for both sizes. Finally, a detailed analysis of the low‐energy isomers of the n = 13 and 19 magic‐number clusters is performed. The two lowest‐energy isomers show S₆ and C₃ symmetry, respectively, which is compatible with the experimental results available in the literature. The structures reported here are all non‐symmetric, showing two central Bz molecules surrounded by 12 nearest‐neighbor monomers in the case of the five lowest‐energy structures. © 2015 Wiley Periodicals, Inc.     

Polygonal current models for polycyclic aromatic hydrocarbons and graphene sheets of various shapes

Assuming that graphene is an “infinite alternant” polycyclic aromatic hydrocarbon resulting from tessellation of a surface by only six‐membered carbon rings, planar fragments of various size and shape (hexagon, triangle, rectangle, and rhombus) have been considered to investigate their response to a magnetic field applied perpendicularly. Allowing for simple polygonal current models, the diatropicity of a series of polycyclic textures has been reliably determined by comparing quantitative indicators, the π‐electron contribution to I^B, the magnetic field‐induced current susceptibility of the peripheral circuit, to and to , respectively the out‐of‐plane components of the magnetizability tensor and of the magnetic shielding tensor at the center of mass. Extended numerical tests and the analysis based on the polygonal model demonstrate that (i) and yield inadequate and sometimes erroneous measures of diatropicity, as they are heavily flawed by spurious geometrical factors, (ii) I^B values computed by simple polygonal models are valid quantitative indicators of aromaticity on the magnetic criterion, preferable to others presently available, whenever current susceptibility cannot be calculated ab initio as a flux integral, (iii) the hexagonal shape is the most effective to maximize the strength of π‐electron currents over the molecular perimeter, (iv) the edge current strength of triangular and rhombic graphene fragments is usually much smaller than that of hexagonal ones, (v) doping by boron and nitrogen nuclei can regulate and even inhibit peripheral ring currents, (vi) only for very large rectangular fragments can substantial current strengths be expected. © 2017 Wiley Periodicals, Inc.     

The solvent effect on two competing reaction mechanisms involving hypervalent iodine reagents (λ3‐iodanes): Facing the limit of the stationary quantum chemical approach

Trifluoromethylation of acetonitrile with 3,3‐dimethyl‐1‐(trifluoromethyl)?1λ³,2‐ benziodoxol is assumed to occur via reductive elimination (RE) of the electrophilic CF₃‐ligand and MeCN bound to the hypervalent iodine. Computations in gas phase showed that the reaction might also occur via an S_N2 mechanism. There is a substantial solvent effect present for both reaction mechanisms, and their energies of activation are very sensitive toward the solvent model used (implicit, microsolvation, and cluster‐continuum). With polarizable continuum model‐based methods, the S_N2 mechanism becomes less favorable. Applying the cluster‐continuum model, using a shell of solvent molecules derived from ab initio molecular dynamics (AIMD) simulations, the gap between the two activation barriers ( ) is lowered to a few kcal mol^?1 and also shows that the activation entropies ( ) and volumes ( ) for the two mechanisms differ substantially. A quantitative assessment of will therefore only be possible using AIMD. A natural bond orbital‐analysis gives further insight into the activation of the CF₃‐reagent by protonation. © 2014 Wiley Periodicals, Inc.     

Force constant decomposition for penta‐coordinated XH3Cl2‐ (X = C,Si, Ge) structures

Based on the energy decomposition analysis of an interacting system, we propose a method for force constant decomposition analysis with respect to the specific normal coordinate. Using the presented method, we examined the penta‐coordinated system (X = C, Si, Ge), which possesses a three‐center four‐electron bond. The origin of the difference in the stability of the penta‐coordinated D_3h structures was clearly shown to be the effect of electron delocalization–polarization term. © 2018 Wiley Periodicals, Inc.     

Molecules under external electric field: On the changes in the electronic structure and validity limits of the theoretical predictions

A detailed analysis of the electronic structure of the ground and first excited spin state of three diatomic molecules ( and ) under static applied electric field is performed at CCSD(T), DFT, MRCI and MRCI(Q) levels of theory. Our findings have revealed that by boosting the applied field one induces changes in the occupation numbers of molecular orbitals, giving rise to changes in the equilibrium geometry and in the HOMO–LUMO energy gap. Specifically, singlet to triplet spin transition can be induced by increasing the applied electric field beyond a critical value. Accordingly, affecting the accuracy of the widely used expression of energy expanded in Taylor series with respect to the applied electric field. © 2018 Wiley Periodicals, Inc.     

Pentaatomic planar tetracoordinate silicon with 14 valence electrons: A large‐scale global search of (n + m = 4; q = 0, ±1, −2; X,Y = main group elements from H to Br)

Designing and characterizing the compounds with exotic structures and bonding that seemingly contrast the traditional chemical rules are a never‐ending goal. Although the silicon chemistry is dominated by the tetrahedral picture, many examples with the planar tetracoordinate‐Si skeletons have been discovered, among which simple species usually contain the 17/18 valence electrons. In this work, we report hitherto the most extensive structural search for the pentaatomic ptSi with 14 valence electrons, that is, (n + m = 4; q = 0, ±1, ?2; X, Y = main group elements from H to Br). For 129 studied systems, 50 systems have the ptSi structure as the local minimum. Promisingly, nine systems, that is, , HSiY₃ (Y = Al/Ga), Ca₃SiAl^?, Mg₄Si^2?, C₂LiSi, Si₃Y₂ (Y = Li/Na/K), each have the global minimum ptSi. The former six systems represent the first prediction. Interestingly, in HSiY₃ (Y = Al/Ga), the H‐atom is only bonded to the ptSi‐center via a localized 2c–2e σ bond. This sharply contradicts the known pentaatomic planar‐centered systems, in which the ligands are actively involved in the ligand–ligand bonding besides being bonded to the planar center. Therefore, we proposed here that to generalize the 14e‐ptSi, two strategies can be applied as (1) introducing the alkaline/alkaline‐earth elements and (2) breaking the peripheral bonding. In light of the very limited global ptSi examples, the presently designed six systems with 14e are expected to enrich the exotic ptSi chemistry and welcome future laboratory confirmation. © 2014 Wiley Periodicals, Inc.     

Effects of zero point vibration on the reaction dynamics of water dimer cations following ionization

Reactions of water dimer cation following ionization have been investigated by means of a direct ab initio molecular dynamics method. In particular, the effects of zero point vibration and zero point energy (ZPE) on the reaction mechanism were considered in this work. Trajectories were run on two electronic potential energy surfaces (PESs) of : ground state (²A″‐like state) and the first excited state (²A′ ‐ like state). All trajectories on the ground‐state PES lead to the proton‐transferred product: H₂O⁺(Wd)‐H₂O(Wa) → OH(Wd)‐H₃O⁺(Wa), where Wd and Wa refer to the proton donor and acceptor water molecules, respectively. Time of proton transfer (PT) varied widely from 15 to 40 fs (average time of PT = 30.9 fs). The trajectories on the excited‐state PES gave two products: an intermediate complex with a face‐to‐face structure (H₂O‐OH₂)⁺ and a PT product. However, the proton was transferred to the opposite direction, and the reverse PT was found on the excited‐state PES: H₂O(Wd)‐H₂O⁺ (Wa) → H₃O⁺(Wd)‐OH(Wa). This difference occurred because the ionizing water molecule in the dimer switched between the ground and excited states. The reaction mechanism of and the effects of ZPE are discussed on the basis of the results. © 2017 Wiley Periodicals, Inc.