Electronic polarization in liquid acetonitrile: A sequential Monte Carlo/quantum mechanics investigation

Exploring the molecular solvatochromism, stability, reactivity, and non-linear optical response of resveratrol

CONTEXT

This work analyzes the isomerization effects and solvent contributions to the stability, electronic excitations, reactivity, and non-linear optical properties (NLO) of resveratrol molecules within the formalism of the Density Functional Theory. The findings suggest that resveratrol solvatochromism is significantly influenced by solvent polarization. The electronic and free energies (E and G) indicate that trans is the most stable conformer. The system is classified as a strong nucleophile. However, the analysis of the Fukui functions and the Mulliken charges indicate that cis-trans isomerization jointly affects the reactive indices of the carbon and hydrogen atoms. The results also suggest that solvent is relevant to solvatochromism and the NLO response. Both cis and trans conformers present strong π - π ∗ excitations that undergo a visible hypsochromic change when the polarity of the solvent increases. Once the absorption spectra are connected to the first hyperpolarization ( β ) by the Oudar and Chemla relation, the hypsochromism of resveratrol is the reason for the drop in the generation of the second harmonic when the ambient polarity decreases. The CAM-B3LYP DFT results suggest that resveratrol is interesting for NLO applications. Depending on the choice of solvent, values ∼ 50 times those observed for urea ( β = 0.34 × 10 - 34 esu), which is a standard NLO material.

METHODS

The optimized geometries of cis and trans isomers of resveratrol in vacuum were obtained using Density Functional Theory (DFT) with the hybrid exchange-correlation function (CAM-B3LYP) and Pople basis set functions, specifically 6-311++G(d,p). The solvent effect on the geometries of both isomers was included using the polarizable continuum model (PCM) with the same level of QM calculation. Vibrational analysis was conducted to confirm that all optimized geometries correspond to the minimum energy. Various electronic properties, including dipole moments, molecular orbitals, transition energy, dipole polarizabilities, and global reactivity parameters, were calculated using both continuum and discrete solvation models based on the sequential QM/MM methodology. All QM calculations were performed with the Gaussian 09 program and the MC simulations with the DICE program. All NLO analysis was carried out using the Multiwfn code.

Reactivity properties of the HOSO and HSO2 isomers in liquid medium: a sequential Monte Carlo/quantum mechanics study

CONTEXT

The rationalization of acid rain formation steps is fundamental for mitigating its effects. It is believed the hydroxysulfinyl radical is an intermediate species for the production of atmospheric sulfuric acid. Two stable configurations HOSO and HSO₂ have been reported for such a radical in the gas phase. This work aims at studying these isomers in the aqueous medium. The electrical and reactivity quantities - electronic chemical potential ([Formula: see text]), chemical hardness ([Formula: see text]), and electrophilicity ([Formula: see text]) - are here calculated and compared. Considering first solvation shells (15 H₂O for HSO₂ and 23 H₂O for HOSO), an increase of 25% in the dipole moment of HSO₂ was obtained, while the dipole moment of HOSO decreases in 11%. Both solvated isomers grow softer ([Formula: see text] decreases) contrasted to the gas phase.

METHODS

HOSO and HSO₂ are studied through a sequential Monte Carlo/quantum mechanics approach. Lennard-Jones plus the Coulomb potentials were used to represent intermolecular potential interaction in the frame of the DICE package. Molecular structure calculations were performed at CASPT2/aug - cc - pV(T + d)Z level of theory using the MOLPRO suite of programs.

Structure-Dynamics Interrelation Governing Charge Transport in Cosolvated Acetonitrile/LiTFSI Solutions

Concentrated ionic solutions present a potential improvement for liquid electrolytes. However, their conductivity is limited by high viscosities, which can be attenuated via cosolvation. This study employs a series of experiments and molecular dynamics simulations to investigate how different cosolvents influence the local structure and charge transport in concentrated lithium bis(trifluoromethane-sulfonyl)imide (LiTFSI)/acetonitrile solutions. Regardless of whether the cosolvent's dielectric constant is low (for toluene and dichloromethane), moderate (acetone), or high (methanol and water), they preserve the structural and dynamical features of the cosolvent-free precursor. However, the dissimilar effects of each case must be individually interpreted. Toluene and dichloromethane reduce the conductivity by narrowing the distribution of Li⁺-TFSI^- interactions and increasing the activation energies for ionic motions. Methanol and water broaden the distributions of Li⁺-TFSI^- interactions, replace acetonitrile in the Li⁺ solvation, and favor short-range Li⁺-Li⁺ interactions. Still, these cosolvents strongly interact with TFSI^-, leading to conductivities lower than that predicted by the Nernst-Einstein relation. Finally, acetone preserves the ion-ion interactions from the cosolvent-free solution but forms large solvation complexes by joining acetonitrile in the Li⁺ solvation. We demonstrate that cosolvation affects conductivity beyond simply changing viscosity and provide fairly unexplored molecular-scale perspectives regarding structure/transport phenomena relation in concentrated ionic solutions.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.

Structures, intermolecular interactions, and chemical hardness of binary water-organic solvents: a molecular dynamics study

The evolution of structural properties, thermodynamics and averaged (dynamic) total hardness values as a function of the composition of binary water-organic solvents, was rationalized in view of the intermolecular interactions. The organic solvents considered were ethanol, acetonitrile, and isopropanol at 0.25, 0.5, 0.75, and 1 mass fractions, and the results were obtained using molecular dynamics simulations. The site-to-site radial distribution functions reveal a well-defined peak for the first coordination shell in all solvents. A characteristic peak of the second coordination shell exists in aqueous mixtures of acetonitrile, whereas in the water-alcohol solvents, a second peak develops with the increase in alcohol content. From the computed coordination numbers, averaged hydrogen bonds and their lifetimes, we found that water mixed with acetonitrile largely preserves its structural features and promotes the acetonitrile structuring. Both the water and alcohol structures in their mixtures are disturbed and form hydrogen bonds between molecules of different kinds. The dynamic hardness values are obtained as the average over the total hardness values of 1200 snapshots per solvent type, extracted from the equilibrium dynamics. The dynamic hardness profile has a non-linear evolution with the liquid compositions, similarly to the thermodynamic properties of these non-ideal solvents. Graphical abstract Computed dynamic total hardness, as a function of the cosolvent mass fraction for water-ethanol (EtOH), water-isopropanol (2PrOH) and water-acetonitrile (AN).

What is measured by hyper-Rayleigh scattering from a liquid?

Polarization and angle dependence of hyper-Rayleigh scattering (HRS) measured for liquid acetonitrile and dimethyl sulfoxide (DMSO) is analyzed in terms of contributions from randomly oriented molecules and additional contributions produced during intermolecular collisions and induced by the electric field of dissolved ions. All three contributions show the effect of long-range correlation, and the correlation functions are determined using the HRS observations combined with the results of molecular dynamics simulations. HRS from acetonitrile is polarized transverse to the scattering vector. This is due to long-range molecular orientation correlation produced by the dipole-dipole interaction, and correlation at distances r > 100 nm must be included to account for the HRS observations. Analysis of the HRS measurements for acetonitrile determines the length scale a = 0.185 nm for the long-range longitudinal and transverse orientation correlation functions B_L=-2B_T=a³/r³. Transverse polarized collision-induced HRS is also observed for acetonitrile, indicating long-range correlation of intermolecular modes. Strong longitudinal HRS is induced by the radial electric field of dissolved ions in acetonitrile. For DMSO, the angle between the molecular dipole and the vector part of the first hyperpolarizability tensor is about 100°. As a result, HRS from the randomly oriented molecules in DMSO is nearly unaffected by dipole correlation, and ion-induced HRS is weak. The strong longitudinal polarized HRS observed for DMSO is due to the collision-induced contribution, indicating long-range correlation of intermolecular modes. The HRS observations require correlation that has r^-3 long-range asymptotic form, for molecular orientation and for intermolecular vibration and libration, for both acetonitrile and DMSO.

Solvent Effects on the Dynamic Polarizability and Raman Response of Molecule-Metal Oxide Hybrid Clusters

Currently, there is considerable interest in the properties of semiconducting metal oxide nanoparticle substrates because of their utility in surface-enhanced Raman scattering, dye-sensitized solar cells, and photocatalysis. While the enhancement of Raman activities of molecules adsorbed on these nanoparticles is due to a large increase in the polarizability, because of charge transfer from the molecule to the semiconducting nanoparticle, little is known about the factors responsible for modulating the polarizability, particularly the influence of the solvent. Consequently, we have carried out Monte Carlo simulations of several hybrids to study the solvent effect on the dynamic polarizabilities and electronic spectra. Our results indicate that the presence of the solvent induces a shift and an increase in the polarization response that is dependent on the identity of the hybrid. The observed enhancement can be attributed to both the resonant character of the excitation and the participation of the solvent in the charge redistribution. The methodology employed in this work could be very valuable in both identifying and developing metal oxides as novel molecular sensors.

Theoretical study of the solvent effect on the aromaticity of benzene: a NICS analysis

Nucleus-independent chemical shift (NICS) quantities for benzene-benzene and benzene-water species were obtained and are discussed in gas phase and in solution. Besides standard polarizable continuum model (PCM) calculations, sequential Monte Carlo/quantum mechanics (S-MC/QM) were also performed. Benzene was shown to be slightly more aromatic in condensate phase when we considered the average solvent configuration (ASEC) approach with explicit molecules.

The UV-vis absorption spectrum of the flavonol quercetin in methanolic solution: A theoretical investigation

The UV-vis absorption spectrum of the solvated quercetin molecule in methanol was investigated theoretically by means of an elegant type of QM/MM scheme better known as sequential Monte Carlo/quantum mechanics (S-MC/QM) methodology. A set of 125 uncorrelated Monte Carlo molecular liquid structures were properly selected through the autocorrelation function of the energy in order to be used in the quantum mechanical calculations. These molecular liquid structures were obtained by means of the radial and minimum distance distribution functions. A detailed account of the pattern of hydrogen bond structures obtained in this study is also available. The computed results obtained here were directly compared with the available experimental data in order to validate our theoretical model and through this comparison a very good conformity between theoretical and available experimental results was found.

Structure and UV−Vis Spectrum of C60 Fullerene in Ethanol: A Sequential Molecular Dynamics/Quantum Mechanics Study

A molecular dynamics simulation combined with semiempirical quantum mechanics calculations has been performed to investigate the structure, dynamical, and electronic properties of pure C60 in liquid ethanol. The behavior of the fullerene alcoholic solution was obtained by using the NPT ensemble under ambient conditions, including one C60 fullerene immersed in 1000 ethanol molecules. Our analyzed center-of-mass pairwise radial distribution function indicated that, on average, there are 32, 72, 132, and 187 ethanol molecules around, respectively, the first, second, third, and fourth solvation shells of the C60 molecule. To investigate the UV-vis transition energies of C60 in the presence of ethanol, we have considered constituents of the time uncorrelated supramolecular structures of the first solvation shell, i.e., clusters of C60@{EtOH}32 types. The semiempirical calculations were performed at the intermediate neglect of differential overlap level with configuration interaction singles (INDO/CIS). Our results have pointed out that the characteristic C60 UV-vis absorbance peaks are slightly shifted to longer wavelengths, as compared to the isolated molecule. These findings are in connection with the weak donor-acceptor character of the interactions involving electron lone pairs of oxygen atoms on the solvent and the fullerene surface.

Band gap and density of states of the hydrated C60 fullerene system at finite temperature

We examine the electronic properties of the hydrated C60 fullerene under ambient conditions using a sequential Monte Carlo/density functional theory scheme. In this procedure, the average electronic properties of the first hydration shell of C60 equilibrate for ca. 40 uncorrelated configurations of the fullerene aqueous solution. We obtain a systematic red-shift of 0.8 eV in the band gap of the hydrated system, which is mainly attributed to the thermal fluctuations of the aqueous environment.

Solvent effects on the UV-visible absorption spectrum of benzophenone in water: a combined Monte Carlo quantum mechanics study including solute polarization

The entire ultraviolet-visible absorption spectrum of benzophenone in water is studied and compared with the same spectrum in gas phase. Five transitions are considered, and the corresponding solvatochromic shifts are obtained and compared to experiment. Using a sequential procedure of Monte Carlo simulations and quantum mechanical calculations, liquid configurations were generated and an averaged spectrum of the solution was calculated. The solute polarization was included by an iterative procedure where the atomic charges of the solute were obtained as an average with the solvent distribution. The calculated average dipole moment of benzophenone in water, with MP26-31++G(d,p), converges to the value of 5.84+/-0.05 D, 88% larger than the gas-phase value of 3.11 D. Using 100 statistically uncorrelated configurations and solvation shells with 235 explicit water molecules selected by a minimum-distance distribution of solvent shells, instead of the usual radial distribution, the excitation energies were obtained from solute-solvent all-valence-electron INDO/CIS calculations. The shift of the weak n-pi(*) transition is obtained as 2045+/-40 cm(-1) and the strong and broad pi-pi(*) shift as -1790+/-30 cm(-1). These results are in good agreement with the experimental values of 2200 and -1600 cm(-1), respectively. Standard procedure used by common force fields to generate atomic charges to describe the electrostatic moments of the solute, with HF6-31G(d), gives a dipole moment of 3.64 D. Using these standard charges in the simulation, the average shifts are calculated as 1395+/-35 and -1220+/-25 cm(-1), both about 600 cm(-1) smaller in magnitude than those obtained with the average converged fully polarized solute. The influence of the solute polarization in the solute-solvent interaction and, in particular, in solute-solvent hydrogen bonds is analyzed.

Calculation of the Microscopic and Macroscopic Linear and Nonlinear Optical Properties of Liquid Acetonitrile. II. Local Fields and Linear and Nonlinear Susceptibilities in Quadrupolar Approximation

A discrete model based on the multipolar expansion including terms up to hexadecapoles was employed to describe the electrostatic interactions in liquid acetonitrile. Liquid structures obtained form molecular dynamics simulations with different classical, nonpolarizable potentials were used to analyze the electrostatic interactions. The computed average local field was employed for the determination of the environmental effects on the linear and nonlinear electrical molecular properties. Dipole-dipole interactions yield the dominant contribution to the local field, whereas higher multipolar contributions are small but not negligible. Using the effective in-phase properties, macroscopic linear and nonlinear susceptibilities of the liquid were computed. Depending on the partial charges describing the Coulomb interactions of the force field employed, either the linear properties (refractive index and dielectric constant) were reproduced in good agreement with experiment or the nonlinear properties [third-harmonic generation (THG) and electric field induced second-harmonic (EFISH) generation] and the bulk density but never both sets of properties together. It is concluded that the partial charges of the force fields investigated are not suitable for reliable dielectric properties. New methods are probably necessary for the determination of partial charges, which should take into account the collective and long-range nature of electrostatic interactions more precisely.

Monte Carlo Simulation of Cisplatin Molecule in Aqueous Solution

The Lennard-Jones (12-6) parameters were obtained for all atoms of cisplatin molecule using the ab initio quantum mechanical potential energy surface for the water-cisplatin interaction as reference data. The parameters found were (epsilon/kcal.mol(-1) and sigma/angstroms) 1.0550, 3.6590 (Pt); 0.0381, 4.6272 (Cl); 0.0455, 3.3783 (N); and 0.0185, 0.0936 (H) and provided very good results for the description of the aqueous solution of cisplatin through Monte Carlo simulation. From statistical analysis of solute-solvent interactions, we observed that the NH3 groups are involved in 53% of the calculated hydrogen bonds with a significant contribution from chlorides (41%) and only 6% involving the Pt center. This is in agreement with the expected behavior for such molecules. Two hydration shells with 22 and 81 water molecules, respectively, centered around 4.6 and 7.3 angstroms were found from the center of mass pair correlation function analysis. The cisplatin atomic Lennard-Jones parameters are reported for the first time, and they might be useful for studying the structure, properties, and processes of cisplatin-like molecules in aqueous solution, including explicitly the solvent effect.

Are dipolar liquids ferroelectric?

VH and HV depolarized hyper-Rayleigh scattering spectra were measured for liquid solutions of dipolar CH3CN in nondipolar C2Cl4 at T=300 K. The VH spectrum contains a strong narrow peak due to a slowly relaxing longitudinal orientation mode. This peak is absent in the HV spectrum, and it disappears from the VH spectrum when the CH3CN concentration is reduced to 8%. This observation is consistent with a ferroelectric phase transition predicted to occur when rho mu0(2)=9epsilon0kT=49 D2 M.