Assessing the performances of different continuum solvation models for the calculation of hydration energies of molecules, polymers and surfaces: a comparison between the SMD, VASPsol and FDPB models

TDDFT study on a fluorescent probe for distinguishing analogous thiols based on smiles rearrangement

The complete excited-state sensing mechanism of a fluorescent probe capable of distinguishing cysteine/homocysteine and glutathione from analogous biological thiols has been investigated. Using a TDDFT method, the nature of the fluorescence differences in the detection of thiols by the probe has been explained at the molecular level. Calculation results imply that the probe undergoes photoinduced electron transfer (PET) from the fluorophore to the nitrobenzooxadiazole (NBD)-based acceptor in the excited state. In the presence of a thiol, the NBD moiety is cleaved and the red fluorescence emission of the fluorophore is enhanced through inhibition of the PET process. The sulfur-substituted NBD-thiol product is predisposed to undergo excited-state torsion, leading to fluorescence quenching. However, for cysteine and homocysteine, their appropriate distances lead to Smiles rearrangements with relatively low activation energies (26.60 kJ/mol and 42.94 kJ/mol, respectively) and the emission of a distinct green fluorescence at ambient temperature. It has been theoretically confirmed that the distance between two reactive sites, such as sulfhydryl and amino moieties, can be used to distinguish different thiols, thus providing rational support for the control of fluorescence activity and the design of probe molecules.

Towards a transferable nonelectrostatic model for continuum solvation: The electrostatic and nonelectrostatic energy correction model

In this work, we introduce an electrostatic and non-electrostatic (ENE) correction to the solvation energy based on the Solvent-Accessible Surface Area (SASA) of the solute and the solvent static dielectric constant. The proposed correction was developed for neutral solutes in non-aqueous solvents, considering three different implicit solvation models based on a Self-Consistent Reaction Field treatment of solute-solvent mutual polarization using an Apparent Surface Charge formalism, namely the Integral Equation Formalism of the Polarizable Continuum Model using a continuous surface charge scheme (PCM), the Solvation Model based on solute electron density (SMD), and the generalized Finite-Difference Poisson-Boltzmann (FDPB) model. The proposed correction was parametrized on a diverse training set of 4980 solvation data from the Solv@tum database of experimental solvation energies, and validated on the non-aqueous subset of the MNSOL database comprising 2140 solvation energies. The performances of the proposed ENE models with minimal and extended parameters formulations have been analyzed and the latter variant has been further compared to the widely used Cavity, Dispersion, and Solvent structural effects (CDS) non-electrostatic model originally developed for the SMx family of implicit solvation models. Overall, a very good agreement between the computed solvation energies with the ENE correction and the reference experimental data has been found on both the training and test sets for all continuum solvation models considered. Furthermore, results for the ENE correction are on par with the reference CDS non-electrostatic model for both SMD and FDPB electrostatics, but with the advantage of using a lower number of parameters and thus an improved transferability between different electrostatics treatments.

Theoretical Studies of a Silica Functionalized Acrylamide for Calcium Scale Inhibition

The calcium carbonate (CaCO₃) scale is one of the most common oilfield scales and oil and gas production bane. CaCO₃ scale can lead to a sudden halt in production or, worst-case scenario, accidents; therefore, CaCO₃ scale formation prevention is essential for the oil and gas industry. Scale inhibitors are chemicals that can mitigate this problem. We used two popular theoretical techniques in this study: Density Functional Theory (DFT) and Ab Initio Molecular Dynamics (AIMD). The objective was to investigate the inhibitory abilities of mixed oligomers, specifically acrylamide functionalized silica (AM-Silica). DFT studies indicate that Ca²⁺ does not bind readily to acryl acid and acrylamide; however, it has a good binding affinity with PAM and Silica functionalized PAM. The highest binding affinity occurs in the silica region and not the -CONH functional groups. AIMD calculations corroborate the DFT studies, as observed from the MD trajectory that Ca²⁺ binds to PAM-Silica by forming bonds with silicon; however, Ca²⁺ initially forms a bond with silicon in the presence of water molecules. This bonding does not last long, and it subsequently bonds with the oxygen atoms present in the water molecule. PAM-Silica is a suitable calcium scale inhibitor because of its high binding affinity with Ca²⁺. Theoretical studies (DFT and AIMD) have provided atomic insights on how AM-Silica could be used as an efficient scale inhibitor.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Generalizing Continuum Solvation in Crystal to Nonaqueous Solvents: Implementation, Parametrization, and Application to Molecules and Surfaces

We present an extension of a generalized finite-difference Poisson-Boltzmann (FDPB) continuum solvation model based on a self-consistent reaction field treatment to nonaqueous solvents. Implementation and reparametrization of the cavitation, dispersion, and structural (CDS) effects nonelectrostatic model are presented in CRYSTAL, with applications to both finite and infinite periodic systems. For neutral finite systems, computed errors with respect to available experimental data on free energies of solvation of 2523 solutes in 91 solvents, as well as 144 transfer energies from water to 14 organic solvents are on par with the reference SM12 solvation model for which the CDS parameters have been developed. Calculations performed on a TiO₂ anatase surface and compared to VASPsol data revealed an overall very good agreement of computed solvation energies, surface energies, as well as band structure changes upon solvation in three different solvents, validating the general applicability of the reparametrized FDPB approach to neutral nonperiodic and periodic solutes in aqueous and nonaqueous solvents. For ionic species, while the reparametrized CDS model led to large errors on free energies of solvation of anions, addition of a corrective term based on Abraham's acidity of the solvent significantly improved the accuracy of the proposed continuum solvation model, leading to errors on aqueous pK_a of a test set of 83 solutes divided by a factor of 4 compared to the reference solvation model based on density (SMD). Overall, therefore, these encouraging results demonstrate that the generalized FDPB continuum solvation model can be applied to a broad range of solutes in various solvents, ranging from finite neutral or charged solutes to extended periodic surfaces.

Inspecting the structural characteristics of chiral drug penicillamine under different pH conditions using Raman optical activity spectroscopy and DFT calculations

The investigation of the structural characteristics of chiral drugs in physiological environments is a challenging research topic, which may lead to a better understanding of how the drugs work. Raman optical activity (ROA) spectroscopy in combination with density functional theory (DFT) calculations was exploited to inspect the structural changes in penicillamine under different acid-base states in aqueous solutions. The B3LYP/aug-cc-PVDZ method was employed and the implicit solvation model density (SMD) was considered for describing the solvation effect in H₂O. The conformations of penicillamine varied with pH, but penicillamine was liable to stabilize in the form of the P_C conformation (the sulfur atom is in a trans orientation with respect to carboxylate) in most cases for both D- and L-isomers. The relationship between the conformations of penicillamine and the ROA peaks, as well as peak assignments, were comprehensively studied and elucidated. In the fingerprint region, two ROA couplets and one ROA triplet with different patterns were recognized. The intensity, sign and frequency of the corresponding peaks also changed with varying pH. Deuteration was carried out to identify the vibrational modes, and the ROA peaks of the deuterated amino group in particular are sensitive to change in the ambient environment. The results are expected not only to serve as a reference for the interpretation of the ROA spectra of penicillamine and other chiral drugs with analogous structures but also to evaluate the structural changes of chiral molecules in physiological environments, which will form the basis of further exploration of the effects of structural characteristics on the pharmacological and toxicological properties of chiral drugs.