Solvent effects on acid-base complexes. What is more important: A macroscopic reaction field or solute-solvent interactions?

Reference interaction site model self-consistent field with constrained spatial electron density approach for nuclear magnetic shielding in solution

We propose a new hybrid approach combining quantum chemistry and statistical mechanics of liquids for calculating the nuclear magnetic resonance (NMR) chemical shifts of solvated molecules. Based on the reference interaction site model self-consistent field with constrained spatial electron density distribution (RISM–SCF–cSED) method, the electronic structure of molecules in solution is obtained, and the expression for the nuclear magnetic shielding tensor is derived as the second-order derivative of the Helmholtz energy of the solution system. We implemented a method for calculating chemical shifts and applied it to an adenine molecule in water, where hydrogen bonding plays a crucial role in electronic and solvation structures. We also performed the calculations of 17O chemical shifts, which showed remarkable solvent dependence. While converged results could not be sometimes obtained using the conventional method, in the present framework with RISM–SCF–cSED, an adequate representation of electron density is guaranteed, making it possible to obtain an NMR shielding constant stably. This introduction of cSED is key to extending the method’s applicability to obtain the chemical shift of various chemical species. The present demonstration illustrates our approach’s superiority in terms of numerical robustness and accuracy.

Application of TD-DFT Theory to Studying Porphyrinoid-Based Photosensitizers for Photodynamic Therapy: A Review

An important focus for innovation in photodynamic therapy (PDT) is theoretical investigations. They employ mostly methods based on Time-Dependent Density Functional Theory (TD-DFT) to study the photochemical properties of photosensitizers. In the current article we review the existing state-of-the-art TD-DFT methods (and beyond) which are employed to study the properties of porphyrinoid-based systems. The review is organized in such a way that each paragraph is devoted to a separate aspect of the PDT mechanism, e.g., correct prediction of the absorption spectra, determination of the singlet-triplet intersystem crossing, and interaction with molecular oxygen. Aspects of the calculation schemes are discussed, such as the choice of the most suitable functional and inclusion of a solvent. Finally, quantitative structure-activity relationship (QSAR) methods used to explore the photochemistry of porphyrinoid-based systems are discussed.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.

NMR Properties of the Cyanide Anion, a Quasisymmetric Two-Faced Hydrogen Bonding Acceptor

The isotopically enriched cyanide anion, (13C≡15N)−, has a great potential as the NMR probe of non-covalent interactions. However, hydrogen cyanide is highly toxic and can decompose explosively. It is therefore desirable to be able to theoretically estimate any valuable results of certain experiments in advance in order to carry out experimental studies only for the most suitable molecular systems. We report the effect of hydrogen bonding on NMR properties of 15N≡13CH···X and 13C≡15NH···X hydrogen bonding complexes in solution, where X = 19F, 15N, and O=31P, calculated at the ωB97XD/def2tzvp and the polarizable continuum model (PCM) approximations. In many cases, the isotropic 13C and 15N chemical shieldings of the cyanide anion are not the most informative NMR properties of such complexes. Instead, the anisotropy of these chemical shieldings and the values of scalar coupling constants, including those across hydrogen bonds, can be used to characterize the geometry of such complexes in solids and solutions. 1J(15N13C) strongly correlates with the length of the N≡C bond.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

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.

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.

Editorial to the Special Issue "Gulliver in the Country of Lilliput: An Interplay of Noncovalent Interactions"

Noncovalent interactions allow our world to exist [...].

Sterically Facilitated Intramolecular Nucleophilic NMe2 Group Substitution in the Synthesis of Fused Isoxazoles: Theoretical Study

The influence of steric repulsion between the NMe₂ group and a second ortho-(peri-)substituent in the series of 1-dimethylaminonaphthalene and N,N-dimethylanilene ortho-oximes on the ease of the NMe₂ group’s intramolecular nucleophilic substitution is studied. Possible reaction intermediates for three mechanisms are calculated (ωB97xd/def-2-TZVP), and their free Gibbs energies are compared to model reaction profiles. Supporting experiments have proved the absence of studied reactivity in the case of simple 2-dimethylaminobenzaldoxime, which allowed us to establish reactivity limits. The significant facilitation of NMe₂ group displacement in the presence of bulky substituents is demonstrated. The possibility of fused isoxazoles synthesis via the intramolecular nucleophilic substitution of a protonated NMe₂ group in the aniline and naphthalene series is predicted.

Electric field effect on 31P NMR magnetic shielding

Magnetic shielding depends on molecular structure and noncovalent interactions. This study shows that it is also measurably dependent on the electric field generated by surrounding molecules. This effect has been observed explicitly for ³¹P nucleus using the adduct under field approach. The results obtained indicate that the field strength experienced by molecules in crystals consisting of molecules with large dipole moments is similar to that in polar solvents. Therefore, magnetic shielding should explicitly depend on solvent polarity. It is important to note that this effect cannot be reproduced correctly within the polarizable continuum model approach.

Inter- vs. Intramolecular Hydrogen Bond Patterns and Proton Dynamics in Nitrophthalic Acid Associates

Noncovalent interactions are among the main tools of molecular engineering. Rational molecular design requires knowledge about a result of interplay between given structural moieties within a given phase state. We herein report a study of intra- and intermolecular interactions of 3-nitrophthalic and 4-nitrophthalic acids in the gas, liquid, and solid phases. A combination of the Infrared, Raman, Nuclear Magnetic Resonance, and Incoherent Inelastic Neutron Scattering spectroscopies and the Car-Parrinello Molecular Dynamics and Density Functional Theory calculations was used. This integrated approach made it possible to assess the balance of repulsive and attractive intramolecular interactions between adjacent carboxyl groups as well as to study the dependence of this balance on steric confinement and the effect of this balance on intermolecular interactions of the carboxyl groups.

For Whom a Puddle Is the Sea? Adsorption of Organic Guests on Hydrated MCM-41 Silica

Thermal and hydration effects on the mobility of compact and branched organic molecules and a bulky pharmaceutical substance loaded in submonolayer amounts onto mesoporous silica have been elucidated using ¹H and ³¹P solid-state NMR. In all cases, the ambient hydration has a stronger effect than an increase in temperature to 370 K for water-free silica. The effect of hydration depends on the guest and ranges from complete solvation to a silica-water-guest sandwich structure to a silica-guest/silica-water pattern. The mobility of the guests under different conditions has been described. The specific structure of the MCM-41 surface allows one to study very slow surface diffusion, a diffusivity of about 10^-15-10^-16 m²/s. The data reported are relevant to any nonfunctionalized silica, while the method used is applicable to any phosphor-containing guest on any host.

Gas Phase Computational Study of Diclofenac Adsorption on Chitosan Materials

Environmental pollution with non-steroidal anti-inflammatory drugs and their metabolites exposes living organisms on their long-lasting, damaging influence. Hence, the ways of non-steroidal anti-inflammatory drugs (NSAIDs) removal from soils and wastewater is sought for. Among the potential adsorbents, biopolymers are employed for their good availability, biodegradability and low costs. The first available theoretical modeling study of the interactions of diclofenac with models of pristine chitosan and its modified chains is presented here. Supermolecular interaction energy in chitosan:drug complexes is compared with the the mutual attraction of the chitosan dimers. Supermolecular interaction energy for the chitosan-diclofenac complexes is significantly lower than the mutual interaction between two chitosan chains, suggesting that the diclofenac molecule will encounter problems when penetrating into the chitosan material. However, its surface adsorption is feasible due to a large number of hydrogen bond donors and acceptors both in biopolymer and in diclofenac. Modification of chitosan material introducing long-distanced amino groups significantly influences the intramolecular interactions within a single polymer chain, thus blocking the access of diclofenac to the biopolymer backbone. The strongest attraction between two chitosan chains with two long-distanced amino groups can exceed 120 kcal/mol, while the modified chitosan:diclofenac interaction remains of the order of 20 to 40 kcal/mol.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.