Solvent effects on nuclear shieldings and spin–spin couplings of hydrogen selenide

New pecJ-n (n = 1, 2) Basis Sets for Selenium Atom Purposed for the Calculations of NMR Spin-Spin Coupling Constants Involving Selenium

We present new compact pecJ-n (n = 1, 2) basis sets for the selenium atom developed for the quantum-chemical calculations of NMR spin-spin coupling constants (SSCCs) involving selenium nuclei. These basis sets were obtained at the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes (SOPPA(CCSD)) level with the property-energy consistent (PEC) method, which was introduced in our previous papers. The existing SSCC-oriented selenium basis sets are rather large in size, while the PEC method gives more compact basis sets that are capable of providing accuracy comparable to that reached using the property-oriented basis sets of larger sizes generated with a standard even-tempered technique. This is due to the fact that the PEC method is very different in its essence from the even-tempered approaches. It generates new exponents through the total optimization of angular spaces of trial basis sets with respect to the property under consideration and the total molecular energy. New basis sets were tested on the coupled cluster singles and doubles (CCSD) calculations of SSCCs involving selenium in the representative series of molecules, taking into account relativistic, solvent, and vibrational corrections. The comparison with the experiment showed that the accuracy of the results obtained with the pecJ-2 basis set is almost the same as that provided by a significantly larger basis set, aug-cc-pVTZ-J, while that achieved with a very compact pecJ-1 basis set is only slightly inferior to the accuracy provided by the former.

Redfield Propagation of Photoinduced Electron Transfer Reactions in Vacuum and Solution

Dynamical simulations of ultrafast electron transfer reactions are of utmost interest. To allow for energy dissipation directly into an external surrounding environment, a solvent coupling model has been deduced, implemented, and utilized to describe the photoinduced electron transfer dynamics within a model triad system herein. The model is based on Redfield theory, and the environment is represented by harmonic oscillators filled with bosonic quanta. To imitate real solvents, the oscillators have been equipped with frequencies and polarization lifetimes characteristic of the corresponding solvent. The population was found to transfer through the energetically lowest electron transfer route regardless of the medium. The condensed population transfer dynamics were observed to be highly dependent on the solvent parameters. In particular, an increase in the solvent coupling entailed a detainment in the population transfer from the initially prepared diabatic state and a promotion in the population transfer through the other electron transfer route. Two explanations based on the diagonal and off-diagonal matrix elements of the Kohn-Sham Fock matrix, respectively, have been provided. The lifetime of the populated partially charge-separated state was prolonged with increasing solvent polarity, and it was explained in terms of attractive interactions between the solvent's dipole moments and the fragments' charges. The high-frequency vibrational fine-structure in the correlation function was demonstrated to be important for the transfer dynamics, and the importance of dephasing effects in polar solvents was verified and precised to concern the optical polarization of the solvents.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

Quantum chemical calculations of 77 Se and 125 Te nuclear magnetic resonance spectral parameters and their structural applications

An accurate quantum chemical (QC) modeling of ⁷⁷ Se and ¹²⁵ Te nuclear magnetic resonance (NMR) spectra is deeply involved in the NMR structural assignment for selenium and tellurium compounds that are of utmost importance both in organic and inorganic chemistry nowadays due to their huge application potential in many fields, like biology, medicine, and metallurgy. The main interest of this review is focused on the progress in QC computations of ⁷⁷ Se and ¹²⁵ Te NMR chemical shifts and indirect spin-spin coupling constants involving these nuclei. Different computational methodologies that have been used to simulate the NMR spectra of selenium and tellurium compounds since the middle of the 1990s are discussed with a strong emphasis on their accuracy. A special accent is placed on the calculations resorting to the relativistic methodologies, because taking into account the relativistic effects appreciably influences the precision of NMR calculations of selenium and, especially, tellurium compounds. Stereochemical applications of quantum chemical calculations of ⁷⁷ Se and ¹²⁵ Te NMR parameters are discussed so as to exemplify the importance of integrated approach of experimental and computational NMR techniques.

Value of NMR Parameters and DFT Calculations for Quantum Information Processing Utilizing Phosphorus Heterocycles

Quantum computing is the field of science that uses quantum-mechanical phenomena, such as superposition and entanglement, to perform operations on data. The fundamental information unit used in quantum computing is the quantum bit or qubit. It is well-known that quantum computers could theoretically be able to solve problems much more quickly than any classical computers. Currently, the first and still the most successful implementations of quantum information processing (QIP) have been based on nuclear spins in liquids. However, molecules that enable many qubits NMR QIP implementations should meet some conditions: have large chemical shifts and be appropriately dispersed for qubit addressability, appreciable spin-spin coupling between any pair of spins, and a long relaxation time. In this line, benzyldene-2,3-dihydro-1H-[1,3]diphosphole (BDF) derivatives have been theoretically tested for maximizing large chemical shifts, spin-spin coupling, and minimizing the hyperfine coupling constant. Thus, the structures were optimized at the B3LYP/6-311G(d,p) level and showed a significant similarity with the experimental geometrical parameters. The NMR spectroscopic parameters (δ and J) were calculated with six different DFT functionals. The τ-HCTH/6-31G(2d) level is in better agreement with the experimental data of ³¹P and ¹³C chemical shifts, while PCM-B3LYP/cc-pVDZ level shows a decrease on deviation between calculated and experimental values for P-P and P-C SSCC. The surface response technique was employed to rationalize how the hyperfine constant varies with the chemical shifts and coupling constants values. From our findings, BDF-NO₂ was the best candidate for NMR quantum computations (NMR-QC) among the studied series.

Solvent Dependence of (14)N Nuclear Magnetic Resonance Chemical Shielding Constants as a Test of the Accuracy of the Computed Polarization of Solute Electron Densities by the Solvent

Although continuum solvation models have now been shown to provide good quantitative accuracy for calculating free energies of solvation, questions remain about the accuracy of the perturbed solute electron densities and properties computed from them. Here we examine those questions by applying the SM8, SM8AD, SMD, and IEF-PCM continuum solvation models in combination with the M06-L density functional to compute the (14)N magnetic resonance nuclear shieldings of CH3CN, CH3NO2, CH3NCS, and CH3ONO2 in multiple solvents, and we analyze the dependence of the chemical shifts on solvent dielectric constant. We examine the dependence of the computed chemical shifts on the definition of the molecular cavity (both united-atom models and models based on superposed individual atomic spheres) and three kinds of treatments of the electrostatics, namely the generalized Born approximation with the Coulomb field approximation, the generalized Born model with asymmetric descreening, and models based on approximate numerical solution schemes for the nonhomogeneous Poisson equation. Our most systematic analyses are based on the computation of relative (14)N chemical shifts in a series of solvents, and we compare calculated shielding constants relative to those in CCl4 for various solvation models and density functionals. While differences in the overall results are found to be reasonably small for different solvation models and functionals, the SMx models SM8, and SM8AD, using the same cavity definitions (which for these models means the same atomic radii) as those employed for the calculation of free energies of solvation, exhibit the best agreement with experiment for every functional tested. This suggests that in addition to predicting accurate free energies of solvation, the SM8 and SM8AD generalized Born models also describe the solute polarization in a manner reasonably consistent with experimental (14)N nuclear magnetic resonance spectroscopy. Models based on the nonhomogeneous Poisson equation show slightly reduced accuracy. Scaling the intrinsic Coulomb radii to larger values (as has sometimes been suggested in the past) does not uniformly improve the results for any kind of solvent model; furthermore it uniformly degrades the results for generalized Born models. Use of a basis set that increases the outlying charge diminishes the accuracy of continuum models that solve the nonhomogeneous Poisson equation, which we ascribe to the inability of the numerical schemes for approximately solving the nonhomogeneous Poisson equation to fully account for the effects of electronic charge outside the solute cavity.

Divinyl selenide: conformational study and stereochemical behavior of its 77Se-1H spin-spin coupling constants

Theoretical energy-based conformational analysis of divinyl selenide performed at the MP2/6-311G** level is substantiated by the second-order polarization propagator approach (SOPPA) calculations and experimental measurements of its (77)Se-(1)H spin-spin coupling constants, demonstrating marked stereochemical behavior in respect of the internal rotation of both vinyl groups around the Se-C bonds. Based on these data, divinyl selenide is shown to exist in an equilibrium mixture of three nonplanar conformers: one the preferred syn-s-cis-s-trans and two minor anti-s-trans-s-trans and syn-s-trans-s-trans forms.

Dielectric permittivity and temperature effects on spin-spin couplings studied on acetonitrile

Dielectric permittivity (epsilon) and temperature effects on indirect spin-spin coupling constants were studied using acetonitrile as a probe molecule. Experiments were accompanied by hybrid DFT (density functional theory) studies, where the solvent was modeled using the polarization continuum model. Owing to its numerous types of J-couplings, acetonitrile is a very convenient molecule against which various basis sets can be tested or the best basis set can be selected for a given study. The results show reasonable agreement between calculated and experimental values. According to our data, scalar spin-spin coupling constants undergo substantial shifts at lower values of the dielectric constant. Thus J-coupling values are not transferable between measurements made at differing epsilon-conditions, and the assumption of the epsilon-independence of the J-coupling can lead to crucial mistakes in experiments using low-epsilon media. Dielectric permittivity also causes small geometric fluctuations within the molecule, which themselves can affect J-coupling values. Examinations of the results computed with frozen and relaxed geometries show that geometry mediation mostly affects the spin-dipole term of the J-coupling; hence, for accurate evaluation of the latter, frozen geometries are not acceptable. Another interesting fact revealed is the connection between the solvent dielectric properties and the temperature-dependence slopes of J-couplings in corresponding media.

Theory and computation of nuclear magnetic resonance parameters

The art of quantum chemical electronic structure calculation has over the last 15 years reached a point where systematic computational studies of magnetic response properties have become a routine procedure for molecular systems. One of their most prominent areas of application are the spectral parameters of nuclear magnetic resonance (NMR) spectroscopy, due to the immense importance of this experimental method in many scientific disciplines. This article attempts to give an overview on the theory and state-of-the-art of the practical computations in the field, in terms of the size of systems that can be treated, the accuracy that can be expected, and the various factors that would influence the agreement of even the most accurate imaginable electronic structure calculation with experiment. These factors include relativistic effects, thermal effects, as well as solvation/environmental influences, where my group has been active. The dependence of the NMR spectra on external magnetic and optical fields is also briefly touched on.

Calculation of nuclear magnetic resonance shieldings using frozen-density embedding

We have extended the frozen-density embedding (FDE) scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] to include external magnetic fields and applied this extension to the nonrelativistic calculation of nuclear magnetic resonance (NMR) shieldings. This leads to a formulation in which the electron density and the induced current are calculated separately for the individual subsystems. If the current dependence of the exchange-correlation functional and of the nonadditive kinetic-energy functional are neglected, the induced currents in the subsystems are not coupled and each of them can be determined without knowledge of the induced current in the other subsystem. This allows the calculation of the NMR shielding as a sum of contributions of the individual subsystems. As a test application, we have calculated the solvent shifts of the nitrogen shielding of acetonitrile for different solvents using small geometry-optimized clusters consisting of acetonitrile and one solvent molecule. By comparing to the solvent shifts obtained from supermolecular calculations we assess the accuracy of the solvent shifts obtained from FDE calculations. We find a good agreement between supermolecular and FDE calculations for different solvents. In most cases it is possible to neglect the contribution of the induced current in the solvent subsystem to the NMR shielding, but it has to be considered for aromatic solvents. We demonstrate that FDE can describe the effect of induced currents in the environment accurately.

Comparison of frozen-density embedding and discrete reaction field solvent models for molecular properties

We investigate the performance of two discrete solvent models in connection with density functional theory (DFT) for the calculation of molecular properties. In our comparison we include the discrete reaction field (DRF) model, a combined quantum mechanics and molecular mechanics (QM/MM) model using a polarizable force field, and the frozen-density embedding (FDE) scheme. We employ these solvent models for ground state properties (dipole and quadrupole moments) and response properties (electronic excitation energies and frequency-dependent polarizabilities) of a water molecule in the liquid phase. It is found that both solvent models agree for ground state properties, while there are significant differences in the description of response properties. The origin of these differences is analyzed in detail and it is found that they are mainly caused by a different description of the ground state molecular orbitals of the solute. In addition, for the calculation of the polarizabilities, the inclusion of the response of the solvent to the polarization of the solute becomes important. This effect is included in the DRF model, but is missing in the FDE scheme. A way of including it in FDE calculations of the polarizabilities using finite field calculations is demonstrated.

Calculation of binary magnetic properties and potential energy curve in xenon dimer: Second virial coefficient of 129Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order Møller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.

Influence of the correlation, aggregation, and solvation on ab initio computed Li-C, Li-N, and Li-Li NMR coupling constants

The (1)J and (3)J(C-Li), (1)J(N-Li), and (2)J(Li-Li) NMR coupling constants have been calculated for various homogeneous and heterogeneous aggregates of methyllithium and lithium dimethylamide at the HF and MP2 levels of calculation. Ethereal solvation has also been taken into account either through a continuum model or through the explicit introduction of Me(2)O molecules. The results obtained are in good general agreement with the experimental data available for methyllithium itself or model alkyllithiums and supports the empirical rule proposed by Bauer, Winchester, and Schleyer to evaluate (1)J(C-Li) provided that calculations include solvent and/or aggregation effects.