No Arabic abstract
A method for study of charge-transfer interactions between solute molecules and solvent based on the comparison of the ratios of spectral shifts of different electronic transitions in solute molecules in chemically inert solvent is proposed. The method is applicable to molecules that do not change their dipole moment on excitation. As an example, a presence of charge transfer interactions in higher electronic states of aromatic hydrocarbons (benzene, phenanthrene, and naphthalene) dissolved in water and alcohols was demonstrated.
The electronic and optical properties of polycyclic aromatic hydrocarbons (PAHs) present a strong dependence on their size and geometry. We tackle this issue by analyzing the spectral features of two prototypical classes of PAHs, belonging to D6h and D2h symmetry point groups and related to coronene as multifunctional seed. While the size variation induces an overall red shift of the spectra and a redistribution of the oscillator strength between the main peaks, a lower molecular symmetry is responsible for the appearance of new optical features. Along with broken molecular orbital degeneracies, optical peaks split and dark states are activated in the low-energy part of the spectrum. Supported by a systematic analysis of the composition and the character of the optical transitions, our results contribute in shedding light to the mechanisms responsible for spectral modifications in the visible and near UV absorption bands of medium-size PAHs.
We present the extension of the quantum/classical polarizable fluctuating charge model to the calculation of single residues of quadratic response functions, as required for the computational modeling of two-photon absorption cross-sections. By virtue of a variational formulation of the quantum/classical polarizable coupling, we are able to exploit an atomic orbital-based quasienergy formalism to derive the additional coupling terms in the response equations. Our formalism can be extended to the calculation of arbitrary order response functions and their residues. The approach has been applied to the challenging problem of one- and two-photon spectra of rhodamine 6G (R6G) in aqueous solution. Solvent effects on one- and two-photon spectra of R6G in aqueous solution have been analyzed by considering three different approaches, from a continuum (QM/PCM) to two QM/MM models (non-polarizable QM/TIP3P and polarizable QM/FQ). Both QM/TIP3P and QM/FQ simulated OPA and TPA spectra show that the inclusion of discrete water solvent molecules is essential to increase the agreement between theory and experiment. QM/FQ has been shown to give the best agreement with experiments.
Functionalized polcyclic aromatic hydrocarbons (PAHs) are an interesting class of molecules in which the electronic state of the graphene-like hydrocarbon part is tuned by the functional group. Searching for new types of donor and acceptor molecules, a set of new PAHs has recently been investigated experimentally using ultraviolet photoelectron spectroscopy (UPS). In this work, the electronic structure of the PAHs is studied numerically with the help of B3LYP hybrid density functionals. Using the DELTA-SCF method, electron binding energies have been determined which affirm, specify and complement the UPS data. Symmetry properties of molecular orbitals are analyzed for a categorization and an estimate of the related signal strength. While SIGMA-like orbitals are difficult to detect in UPS spectra of condensed film, calculation provides a detailed insight into the hidden parts of the electronic structure of donor and acceptor molecules. In addition, a diffuse basis set (6-311++G**) was used to calculate electron affinity and LUMO eigenvalues. The calculated electron affinity (EA) provides a classification of the donor/acceptor properties of the studied molecules. Coronene-hexaone shows a high EA, comparable to TCNQ, which is a well-known classical acceptor. Calculated HOMO-LUMO gaps using the related eigenvalues have a good agreement with the experimental lowest excitation energies. TD-DFT also accurately predicts the measured optical gap.
We present a new method to accurately describe the ionization fraction and the size distribution of polycyclic aromatic hydrocarbons (PAHs) within astrophysical sources. To this purpose, we have computed the mid-infrared emission spectra of 308 PAH molecules of varying sizes, symmetries, and compactness, generated in a range of radiation fields. We show that the intensity ratio of the solo CH out-of-plane bending mode in PAH cations and anions (referred to as the 11.0 $mu$m band, falling in the 11.0-11.3 $mu$m region for cations and anions) to their 3.3 $mu$m emission, scales with PAH size, similarly to the scaling of the 11.2/3.3 ratio with the number of carbon atoms (N$_{mathrm{C}}$) for neutral molecules. Among the different PAH emission bands, it is the 3.3 $mu$m band intensity which has the strongest correlation with N$_{mathrm{C}}$, and drives the reported PAH intensity ratio correlations with N$_{mathrm{C}}$ for both neutral and ionized PAHs. The 6.2/7.7 intensity ratio, previously adopted to track PAH size, shows no evident scaling with N$_{mathrm{C}}$ in our large sample. We define a new diagnostic grid space to probe PAH charge and size, using the (11.2+11.0)/7.7 and (11.2+11.0)/3.3 PAH intensity ratios respectively. We demonstrate the application of the (11.2+11.0)/7.7 - (11.2+11.0)/3.3 diagnostic grid for galaxies M82 and NGC 253, for the planetary nebula NGC 7027, and the reflection nebulae NGC 2023 and NGC 7023. Finally, we provide quantitative relations for PAH size determination depending on the ionization fraction of the PAHs and the radiation field they are exposed to.
Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions - described by a truncated Lennard-Jones (LJ) potential - and electrostatic interactions - described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges (electronic continuum correction, ECC), possibly complemented with rescaling of LJ parameters (electronic continuum correction rescaled, ECCR), have shown promising results in bulk, but their performance at interfaces has been less explored. Here we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk also describes accurately image charge repulsion at the liquid-vapor interface, and that the rescaling of LJ parameters in ECCR models - aimed at capturing local ion-ion and ion-water interactions in bulk - also describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces, and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.