Theoretical study of electronic structure of rhodium mononitride and interpretation of experimental spectra

Potential energy curves of 22 electronic states of RhN have been calculated by the complete active space second‐order perturbation theory method. The X¹Σ₀⁺ is assigned as the ground state, and the first excited state a³Π₀+ is 978 cm^?1 higher. The ¹Δ(I) and B¹Σ⁺ states are located at 9521 and 13,046 cm^?1 above the ground state, respectively. The B¹Σ⁺ state should be the excited state located 12,300 cm^?1 above the ground state in the experimental study. Moreover, two excited states, C¹Π and b³Σ⁺, are found 14,963 and 15,082 cm^?1 above the X¹Σ⁺ state, respectively. The transition C¹Π₁–X¹Σ₀⁺ may contribute to the experimentally observed bands headed at 15,071 cm^?1. There are two excited states, D¹Δ and E¹Σ⁺, situate at 20,715 and 23,145 cm^?1 above the X¹Σ⁺ state. The visible bands near 20,000 cm^?1 could be generated by the electronic transitions D¹Δ₂–a³Π₁ and E¹Σ⁺₀–X¹Σ⁺₀ because of the spin–orbit coupling effect. © 2013 Wiley Periodicals, Inc.     

129Xe NMR chemical shift in Xe@C60 calculated at experimental conditions: Essential role of the relativity,dynamics, and explicit solvent

The isotropic ¹²⁹Xe nuclear magnetic resonance (NMR) chemical shift (CS) in Xe@C₆₀ dissolved in liquid benzene was calculated by piecewise approximation to faithfully simulate the experimental conditions and to evaluate the role of different physical factors influencing the ¹²⁹Xe NMR CS. The ¹²⁹Xe shielding constant was obtained by averaging the ¹²⁹Xe nuclear magnetic shieldings calculated for snapshots obtained from the molecular dynamics trajectory of the Xe@C₆₀ system embedded in a periodic box of benzene molecules. Relativistic corrections were added at the Breit–Pauli perturbation theory (BPPT) level, included the solvent, and were dynamically averaged. It is demonstrated that the contribution of internal dynamics of the Xe@C₆₀ system represents about 8% of the total nonrelativistic NMR CS, whereas the effects of dynamical solvent add another 8%. The dynamically averaged relativistic effects contribute by 9% to the total calculated ¹²⁹Xe NMR CS. The final theoretical value of 172.7 ppm corresponds well to the experimental ¹²⁹Xe CS of 179.2 ppm and lies within the estimated errors of the model. The presented computational protocol serves as a prototype for calculations of ¹²⁹Xe NMR parameters in different Xe atom guest–host systems. © 2013 Wiley Periodicals, Inc.     

Comparative assessment of computational methods for the determination of solvation free energies in alcohol‐based molecules

The determination of differences in solvation free energies between related drug molecules remains an important challenge in computational drug optimization, when fast and accurate calculation of differences in binding free energy are required. In this study, we have evaluated the performance of five commonly used polarized continuum model (PCM) methodologies in the determination of solvation free energies for 53 typical alcohol and alkane small molecules. In addition, the performance of these PCM methods, of a thermodynamic integration (TI) protocol and of the Poisson–Boltzmann (PB) and generalized Born (GB) methods, were tested in the determination of solvation free energies changes for 28 common alkane‐alcohol transformations, by the substitution of an hydrogen atom for a hydroxyl substituent. The results show that the solvation model D (SMD) performs better among the PCM‐based approaches in estimating solvation free energies for alcohol molecules, and solvation free energy changes for alkane‐alcohol transformations, with an average error below 1 kcal/mol for both quantities. However, for the determination of solvation free energy changes on alkane‐alcohol transformation, PB and TI yielded better results. TI was particularly accurate in the treatment of hydroxyl groups additions to aromatic rings (0.53 kcal/mol), a common transformation when optimizing drug‐binding in computer‐aided drug design. © 2013 Wiley Periodicals, Inc.     

Synthesis,Structure, and Reactivity of Lewis Base Stabilized Plumbacyclopentadienylidenes

Plumbacyclopentadienylidenes, in which the lead atoms have divalent states and are coordinated by THF, pyridine and N‐heterocyclic carbene, were synthesized and characterized. The THF‐ and pyridine‐stabilized compounds can be regarded as rare examples of hypervalent 10‐X‐4 species. The equilibrium between the THF adduct and the free plumbacyclopentadienylidene was evidenced by spectroscopic analysis and theoretical calculations. The THF adduct in benzene converted into a plumbylene dimer, where one of the lead centers is coordinated by THF and the other lead atom is coordinated by a divalent lead atom, the dimer gradually decomposing into spiroplumbole. The THF adduct unexpectedly reacted with trifluoroborane and trichlorogallane to afford fluoroborole and chlorogallole, which are the first examples of non‐annulated fluoroborole and gallole, respectively.     

Interactions of delta shock waves for the relativistic Chaplygin Euler equations with split delta functions

In this article, we are concerned with the interactions of delta shock waves with contact discontinuities for the relativistic Euler equations for Chaplygin gas by using split delta functions method. The solutions are obtained constructively and globally when the initial data consists of three piecewise constant states. The global structure and large time‐asymptotic behaviors of the solutions are analyzed case by case. During the process of the interaction, the strengths of delta shock waves are computed completely. Moreover, it can be found that the Riemann solutions are stable for such small perturbations with special initial data by letting perturbed parameter ε tends to zero. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling of Solute-Solvent Interactions Using an External Electric Field—From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute–solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.     

1,3,5-Triaza-7-Phosphaadamantane (PTA) as a 31P NMR Probe for Organometallic Transition Metal Complexes in Solution

Due to the rigid structure of 1,3,5-triaza-7-phosphaadamantane (PTA), its ³¹P chemical shift solely depends on non-covalent interactions in which the molecule is involved. The maximum range of change caused by the most common of these, hydrogen bonding, is only 6 ppm, because the active site is one of the PTA nitrogen atoms. In contrast, when the PTA phosphorus atom is coordinated to a metal, the range of change exceeds 100 ppm. This feature can be used to support or reject specific structural models of organometallic transition metal complexes in solution by comparing the experimental and Density Functional Theory (DFT) calculated values of this ³¹P chemical shift. This approach has been tested on a variety of the metals of groups 8–12 and molecular structures. General recommendations for appropriate basis sets are reported.     

Experimental and four‐component relativistic DFT studies of tungsten carbonyl complexes

We present a theoretical and experimental study of the structure and nuclear magnetic resonance (NMR) parameters of the pentacarbonyltungsten complexes of η¹‐2‐(trimethylstannyl)‐4,5‐dimethylphosphinine, η²‐norbornene, and imidazolidine‐2‐thione. The three complexes have a pseudo‐octahedral molecular structure with the six ligands bonded to the tungsten atom. The η¹‐2‐(trimethylstannyl)‐4,5‐dimethylphosphinine‐pentacarbonyl tungsten complex was synthesized for the first time. For all compounds, we present four‐component relativistic calculations of the NMR parameters at the Dirac–Kohn–Sham density functional level of theory using hybrid functionals. These large‐scale relativistic calculations of NMR chemical shifts and spin–spin coupling constants were compared with available experimental data, either taken from the literature or measured in this work. The inclusion of solvent effects modeled using a conductor‐like screening model was found to improve agreement between the calculated and experimental NMR parameters, and our best estimates for the NMR parameters are generally in good agreement with available experimental results. The present work demonstrates that four‐component relativistic theory has reached a level of maturity that makes it a convenient and accurate tool for modeling and understanding chemical shifts and indirect spin–spin coupling constants of organometallic compounds containing heavy elements, for which conventional non‐relativistic theory breaks down. Copyright © 2015 John Wiley & Sons, Ltd.     

A generalised model of hydrogen diffusion in metals with multiple trap types

Continuum modelling of hydrogen diffusion in metals, which accounts for both trapping and an imposed force field, is revisited. A generalised model of hydrogen diffusion and trapping is developed as a continuous interpretation of the discrete random-walk theory. A system of nonlinear equations describing the phenomenon of diffusion with multiple types of traps is derived without the assumption of a local equilibrium among hydrogen populations in dissimilar positions. Lattice-trap interchange kinetics can degenerate into local equilibrium as a limit case. Moreover, certain terms in general equations may be negligible in specific situations. By removing these terms, known particularised models of hydrogen diffusion and trapping are recovered. Determining the terms, which are disregarded in reduced models, enables a straightforward assessment of the applicability of these models. The advantages and limitations of particularised models applied to hydrogen embrittlement analyses are discussed.