First-Principle Characterization of Structural, Electronic, and Optical Properties of Tin-Halide Monomers.

The growing interest in tin-halide semiconductors for photovoltaic applications demands in-depth knowledge of the fundamental properties of their constituents, starting from the smallest monomers entering the initial stages of formation. In this first-principles work based on time-dependent density-functional theory, we investigate the structural, electronic, and optical properties of tin-halide molecules SnXn 2-n, with n = 1 , 2 , 3 , 4 ${n = 1,2,3,4}$ and X=Cl, Br, I, simulating these compounds in vacuo as well as in an implicit solvent. We find that structural properties are very sensitive to the halogen species while the charge distribution is also affected by stoichiometry. The ionicity of the Sn-X bond is confirmed by the Bader charge analysis albeit charge displacement plots point to more complex metal-halide coordination. Particular focus is posed on the neutral molecules SnX2, for which electronic and optical properties are discussed in detail. Band gaps and absorption onset decrease with increasing size of the halogen species, and despite general common features, each molecule displays peculiar optical signatures. Our results are elaborated in the context of experimental and theoretical literature, including the more widely studied lead-halide analogs, aiming to contribute with microscopic insight to a better understanding of tin-halide perovskites.

TAO-DFT with the Polarizable Continuum Model.

For the ground-state properties of gas-phase nanomolecules with multi-reference character, thermally assisted occupation (TAO) density functional theory (DFT) has recently been found to outperform the widely used Kohn-Sham DFT when traditional exchange-correlation energy functionals are employed. Aiming to explore solvation effects on the ground-state properties of nanomolecules with multi-reference character at a minimal computational cost, we combined TAO-DFT with the PCM (polarizable continuum model). In order to show its usefulness, TAO-DFT-based PCM (TAO-PCM) was used to predict the electronic properties of linear acenes in three different solvents (toluene, chlorobenzene, and water). According to TAO-PCM, in the presence of these solvents, the smaller acenes should have nonradical character, and the larger ones should have increasing polyradical character, revealing striking similarities to the past findings in the gas phase.

The second derivative of the electronic energy with respect to the compression scaling factor in the XP-PCM model: Theory and applications to compression response functions of atoms.

We present the analytical theory for the second derivative of the electronic energy with respect to the scaling factor of the compression cavity within the eXtreme pressure polarizable continuum model (XP-PCM) for the study of compressed atomic and molecular systems. The theory has been exploited to study compression response functions describing how the atomic/molecular properties are effected by an external pressure. The response functions considered include the atomic compressibility and the pressure coefficients of the ionization energy (IE) and electron affinity (EA). The theory has been validated by numerical application to compressed neon, argon, and krypton atoms.

Synthesis, spectroscopic, and computational studies on molecular charge-transfer complex of 2-((2-hydroxybenzylidene) amino)-2-(hydroxymethyl) propane-1, 3-diol with chloranilic acid: Potential antiviral activity simulation of CT-complex against SARS-CoV-2.

An innovative charge-transfer complex between the Schiff base 2-((2-hydroxybenzylidene) amino)-2-(hydroxymethyl) propane-1,3-diol [SAL-THAM] and the π-acceptor, chloranilic acid (CLA) within the mole ratio (1:1) was synthesized and characterized aiming to investigate its electronic transition spectra in acetonitrile (ACN), methanol (MeOH) and ethanol (EtOH) solutions. Applying Job`s method in the three solvents supported the 1:1 (CLA: SAL-THAM) mole ratio complex formation. The formation of stable CT- complex was shown by the highest values of charge-transfer complex formation constants, KCT, calculated using minimum-maximum absorbance method, with the sequence, acetonitrile > ethanol > methanol DFT study on the synthesized CT complex was applied based on the B3LYP method to evaluate the optimized structure and extract geometrical and reactivity parameters. Based on TD-DFT theory, the electronic properties, 1H and 13C NMR, IR, and UV-Vis spectra of the studied system in different solvents showing good agreement with the experimental studies. MEP map described the possibility of hydrogen bonding and charge transfer in the studied system. Finally, a computational approach for screening the antiviral activity of CT - complex towards SARS-CoV-2 coronavirus protease via molecular docking simulation was conducted and confirmed with molecular dynamic (MD) simulation.

Density functional theory studies of the antioxidants-a review.

The following review article attempts to compare the antioxidant activity of the compounds. For this purpose, density functional theory/Becke three-parameter Lee-Yang-Parr (DFT/B3LYP) methodology was carried out instead of using pharmacological methodologies because of economic benefits and high accuracy. This methodology filtrates the compounds with the lowest antioxidant activity. At first, the Koopmans' theorem was carried out to calculate some descriptors to compare antioxidants. The energy of the highest occupied molecular orbitals (HOMO) was accepted as the best indicator, and then some studies confirmed that the highest occupied molecular orbital/lowest unoccupied molecular orbital (HOMO-LUMO) energy gap is the more precise descriptor. Although it would be better to compare spin density distribution (SDD) on the oxygen of the corresponding radical in the polarizable continuum model (PCM) to evaluate their capability to chain reaction inhibition. Next, it was mentioned that in the multi-target directed ligands (MTDLs), the antioxidant is connected to other moieties in para positions to create better antioxidants or novel hybrid compounds. Indeed, SDD was introduced as a descriptor for MTDL antioxidant effectiveness. Then, the relation between antioxidants and aromaticity was investigated. The more the aromaticity of an antioxidant, the more stable the corresponding radical is. Subsequently, in preferred antioxidant activity, it was defined that the hydrogen atom transfer (HAT) mechanism is more favored in metabolism phase I. It has been seen that the solvent model can change the antioxidant mechanism. Therefore, the solvent model is more important than the chemical structure of antioxidants, and an ideal antioxidant should be evaluated in PCM for pharmacological evaluations.

Ab initiomodelling of interfacial electrochemical properties: beyond implicit solvation limitations.

First-principles calculations are an important tool to investigate the complex processes occurring at solid/liquid interfaces which are at the heart of modern technologies. Currently, capturing the whole electrochemical environment at an interface, including the applied potential and solvation, still remains challenging as it necessitates to couple different approaches whose interactions are not fully understood. In this work, a grand canonical density functional theory approach is coupled with solvation models to investigate the electrochemical interfaces under applied potential. We show that a parametrized polarizable continuum model (PCM) which represent solvation in a mean field approach by a continuous polarizable media, possesses catastrophic limitations for the modelling of ionic and charged interfaces. We reveal the origin of PCM instabilities under chemical or electrochemical strong oxidation to be the consequence of a phase transition in the surface Li electronic structure. Thus, PCM undergoes an unphysical response to this phase transition by penetrating within the atomic radius of surface Li atoms. To recover a physical response, an explicit first solvation shell has to be included in addition to the PCM in order to properly describe the electrochemistry of the interface. The Fukui functions show that the first solvation shell becomes involved in the redox process as solvent electron doublet is transferred to the acidic Li+. If another explicit solvent layer is added, the interface electrochemical properties become independent of the PCM parameters: in particular capacitance can then be computed from a parameter-free electrochemical approach. This is an important conclusion as the experimental electrochemical capacitance are not easily found and thus the parametrization of the PCM for electrochemical interface can be difficult. This approach can easily be applied to investigate electrochemical properties at the atomic scale and generalized to any electrochemical device for which interfaces play a crucial role.

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 31P 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 31P 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.

Enhanced Biodegradation/Photodegradation of Organophosphorus Fire Retardant Using an Integrated Method of Modified Pharmacophore Model with Molecular Dynamics and Polarizable Continuum Model.

A comprehensive 3D-quantitative structure-activity relationship (QSAR) pharmacophore model was constructed using the values of comprehensive biodegradation/photodegradation effects of 17 organophosphorus flame retardants (OPFRs) evaluated by a normalization method to modify OPFRs with high biodegradation/photodegradation, taking tris(chloro-isopropyl) phosphate (TCPP), tris(2-chloroethyl) phosphate (TCEP) and tris(1-chloro-2-propyl) phosphate (TCIPP)-which occur frequently in the environment, and are the most difficult to degrade as target molecules. OPFR-derivative molecules TCPP-OH shows the highest improvement in biodegradation and photodegradation (55.48% and 46.37%, respectively). On simulating the biodegradation path and photodegradation path, it is found that the energy barrier of TCPP-OH for phosphate bond cleavage is reduced by 15.73% and 52.52% compared to TCPP after modification, respectively. Finally, in order to further significantly improve its biodegradability and photodegradation, the efficiency enhancement in the biodegradation and photodegradation of TCPP-OH are analyzed under the simulated environment by molecular dynamics and polarizable continuum model, respectively. The results of molecular dynamics show that the biodegradation efficiency of the TCPP-OH increased by 75.52% compared to TCPP. The UV spectral transition energy (4.07 eV) of TCPP-OH under the influence of hydrogen peroxide solvation effect is 44.23% lower than the actual transition energy (7.29 eV) of TCPP.

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 CDF3/CDClF2. With some limitations, they can be applied to the chemically similar CHCl3 and CH2Cl2. 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.

The Role of Computational Chemistry in the Experimental Determination of the Dipole Moment of Molecules in Solution.

The methods for the experimental determination of electric dipole moment of molecules in solution from measurements of dielectric permittivity and refractive index are traditionally based on the classical Onsager model. In this model the molecular solute is approximated as a simple polarizable point dipole inside a spherical or ellipsoidal cavity of a dielectric medium representing the solvent. However, the inadequacies of the model resulting from the assumption of a simple shape of the cavity, for the evaluation of the cavity field effect, and from the uncertainty of the polarizability of the molecular solute influences the results and hampers the comparison with the electric dipole moments computed from quantum chemical solvation models. In this article we propose a new method for the experimental determination of the electric dipole moment in solution in which information from the Polarizable Continuum Model calculations are used in place of the Onsager model. The new method overcomes the limitations of this latter model regarding both the cavity field effect and the polarizability of the molecular solutes, and thus allows a coherent comparison between experimental and computed dipole moments of solvated molecules. © 2019 Wiley Periodicals, Inc.

Helical Multi-walled Carbon Nanotubes as an Efficient Material for the Dispersive Solid-Phase Extraction of Low and High Molecular Weight Polycyclic Aromatic Hydrocarbons from Water Samples: Theoretical Study.

The differences in effectiveness of multi-walled carbon nanotubes (MWCNTs) as the dispersive solid-phase extraction (dSPE) sorbent for the selective extraction of polycyclic aromatic hydrocarbons (PAHs) were explained on the basis of theoretical study. It was observed that for low molecular weight PAHs, the recoveries using non-helical and helical MWCNTs were similar. In contrary, for PAHs containing five or more aromatic rings, the extraction efficiency was higher using HMWCNTs than for non-helical ones. Principle component analysis (PCA) as well as providing structural parameters and interaction energies for adsorption processes (PAH + CNT → PAH-CNT) have been used for this purpose. All the PAH + CNT → PAH-CNT adsorption processes considered were found to be thermodynamically favorable. However, the adsorption energies (Eads) for PAHs and the helical carbon nanotube surface estimated for the B(a)P-HCNT and I(1,2,3-cd)P-HCNT are substantially less negative than those observed for PAH molecules interacting with the non-helical CNT. Namely, the Eads calculated in simulated aqueous environment for the B(a)P-MWCNT(6,2) and I(1,2,3-cd)P-MWCNT(6,2) were respectively - 43.32 and - 59.98 kcal/mol, while values of only - 7.75 kcal/mol (B(a)P-HCNT) and - 9.13 kcal/mol (I(1,2,3-cd)P-HCNT) were found for the corresponding PAH-HCNT systems. Therefore, we conclude that the replacement of MWCNTs with HCNTs leads to PAH-HCNT systems in which the interaction energies are much smaller than those estimated for the corresponding PAH-MWCNT systems. HMWCNTs are therefore recommended as the dSPE sorbent phase for the extraction of both low and high molecular weight PAHs from water samples.

Dielectric Screening Meets Optimally Tuned Density Functionals.

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Computational study of Simultaneous synthesis of optically active (RS)-1,2,4-butanetriol trinitrate (BTTN).

In respective water or ethanol polarizable continuum cavity environments, simultaneous aldol condensation was performed using density functional theory (DFT) computational method to model the synthesis of optically active (RS)-1,2,4-butanetriol trinitrate (BTTN). The results of reaction energy barrier analysis suggested feasible routes with lower activation energies to obtain either the (R)- or (S)-configuration product in ethanolic solution. In addition, local analysis of average inter-particulate distances of reaction species revealed that a stronger inter-particulate interaction accompanied a shorter average distance in the ethanol system. The stabilization effect also indicated that related syntheses would be able to proceed in ethanol. Furthermore, relative to the production of (R)-BTTN, a lower overall energy of 425.3 kJ/mol was required for the synthesis of (S)-BTTN. Through analysis of the effects of temperature on the reaction rates of individual parallel stages of (R)- and (S)-species synthesis, it was simple to adjust the reaction temperature accordingly to differentiate between relative rates in order to obtain a product of a specific configuration. Graphical abstract ᅟ.

Three-body expansion of the fragment molecular orbital method combined with density-functional tight-binding.

The three-body fragment molecular orbital (FMO3) method is formulated for density-functional tight-binding (DFTB). The energy, analytic gradient, and Hessian are derived in the gas phase, and the energy and analytic gradient are also derived for polarizable continuum model. The accuracy of FMO3-DFTB is evaluated for five proteins, sodium cation in explicit solvent, and three isomers of polyalanine. It is shown that FMO3-DFTB is considerably more accurate than FMO2-DFTB. Molecular dynamics simulations for sodium cation in water are performed for 100 ps, yielding radial distribution functions and coordination numbers. © 2017 Wiley Periodicals, Inc.

Comparative simulation study of chemical synthesis of functional DADNE material.

Amorphous molecular simulation to model the reaction species in the synthesis of chemically inert and energetic 1,1-diamino-2,2-dinitroethene (DADNE) explosive material was performed in this work. Nitromethane was selected as the starting reactant to undergo halogenation, nitration, deprotonation, intermolecular condensation, and dehydration to produce the target DADNE product. The Materials Studio (MS) forcite program allowed fast energy calculations and reliable geometric optimization of all aqueous molecular reaction systems (0.1-0.5 M) at 283 K and 298 K. The MS forcite-computed and Gaussian polarizable continuum model (PCM)-computed results were analyzed and compared in order to explore feasible reaction pathways under suitable conditions for the synthesis of DADNE. Through theoretical simulation, the findings revealed that synthesis was possible, and a total energy barrier of 449.6 kJ mol-1 needed to be overcome in order to carry out the reaction according to MS calculation of the energy barriers at each stage at 283 K, as shown by the reaction profiles. Local analysis of intermolecular interaction, together with calculation of the stabilization energy of each reaction system, provided information that can be used as a reference regarding molecular integrated stability. Graphical Abstract Materials Studio software has been suggested for the computation and simulation of DADNE synthesis.

A new extension of the polarizable continuum model: Toward a quantum chemical description of chemical reactions at extreme high pressure.

A quantum chemical method for studying potential energy surfaces of reactive molecular systems at extreme high pressures is presented. The method is an extension of the standard Polarizable Continuum Model that is usually used for Quantum Chemical study of chemical reactions at a standard condition of pressure. The physical basis of the method and the corresponding computational protocol are described in necessary detail, and an application of the method to the dimerization of cyclopentadiene (up to 20 GPa) is reported.

Calculation of Electrochemical Reorganization Energies for Redox Molecules at Self-Assembled Monolayer Modified Electrodes.

Electrochemical electron transfer reactions play an important role in energy conversion processes with many technological applications. Electrodes modified by self-assembled monolayers (SAMs) exhibit reduced double layer effects and are used in molecular electronics. An important quantity for calculating the electron transfer rate constant is the reorganization energy, which is associated with changes in the solute geometry and the environment. In this Letter, an approach for calculating the electrochemical reorganization energy for a redox molecule attached to or near a SAM modified electrode is presented. This integral equations formalism polarizable continuum model (IEF-PCM) approach accounts for the detailed electronic structure of the molecule, as well as the contributions from the electrode, SAM, and electronic and inertial solvent responses. The calculated total reorganization energies are in good agreement with experimental data for a series of metal complexes in aqueous solution. This approach will be useful for calculating electron transfer rate constants for molecular electrocatalysts.

Experimental and theoretical study of the cation binding properties of macrocyclic dehydrodibenzopyrido[15]annulenes.

The UV/Vis titration measurements, vibrational and NMR spectroscopy of isomeric dehydrodibenzopyrido[15]annulenes (DBPA) 1 and 2 clearly show that under proper conditions these macrocycles can achieve fast, quantitative and unselective binding of metal ions. The macrocycle 1 is an example of a hindered amine 2,6-bis(R)pyridine and its isomer 2 of a non-hindered amine 3,5-bis(R)pyridine. The protonation stoichiometry for both 1 and 2 was assumed to be DBPA:H(+)=1:1 and the formation constants logK=4.77±0.02 for 1, and logK=6.78±0.08 for 2 were obtained that well agree with those obtained under similar conditions for a macrocycle containing bipyridine units. The protonation of 2 gave the estimated stoichiometry of 2:H(+)=1:1 while the stoichiometric protonation of macrocycle 1 could not be achieved and the lower stability of the ion pair containing 1H(+) is most likely due to the inaccessibility of the nitrogen atom of 1 to the counterions and solvent molecules. The structures and electronic absorption spectra of 1 and 2, as well as the structures and spectra of 1H(+) and 2H(+), i.e. the species formed by protonation of the pyridine nitrogen, were calculated with the time-dependent DFT method with a B3LYP functional and a 6-31+G(d) basis set. The solvent effects were incorporated by means of the polarizable continuum model (PCM). The agreement of the calculated absorption data for the parent and protonated species with the observed spectra is rather satisfactory. Vibrational IR and Raman spectra of 1, 2, 1H(+) and 2H(+) in vacuo were calculated at the B3LYP/cc-pVTZ level of theory. Macrocycles 1 and 2, and their products protonated by trifluoromethanesulfonic acid (1HOTf and 2HOTf) were also characterized by temperature-dependent FTIR technique known as two dimensional IR correlation analysis. Quite large difference in degradation temperature between macrocycle 1 and 2 and their protonated complexes was measured, indicating that inclusion of proton leads to significant thermal stabilization of dehydroannulene ring.

Solvent effect on the vibrational spectrum of Michler's ketone. Experimental and theoretical investigations.

We examined solvent effect on the IR and Raman spectra of MK in several solvents of different polarity and proticity, for understanding of intermolecular interactions, focusing on solvent effect in detail. It has been found that change of solvent polarity has an ambiguous influence on solvatochromism of MK. We have observed that not only vibrations of carbonyl group are affected by the solvent polarity, but also mode ν(CN) and ν(CC) in IR and Raman spectra of MK. Experimental investigations have been supported by the quantum-mechanical computations to gain more insight into the solvatochromic behavior of Michler's ketone. Calculations have been carried using Kohn-Sham formulation of Density Functional Theory (DFT) and the Polarizable Continuum Model (PCM) was employed to account for solute solvent interactions.