RESUMO

CONTEXT: Push-pull compounds are model systems and have numerous applications. By changing their substituents, properties are modified and new molecules for different applications can be designed. The work investigates the gas-phase electronic absorption spectra of 15 derivatives of push-pull para-nitroaniline (pNA). This molecule has applications in pharmaceuticals, azo dyes, corrosion inhibitors, and optoelectronics. Both electron-donor and electron-withdrawing groups were investigated. Employing machine learning-derived Hammett's constants σm, σm0, σR, and σI, correlations between substituents and electronic properties were obtained. Overall, the σm0 constants presented the best correlation with HOMO and LUMO energies, whereas the σR constants best agreed with the transition energy of the first band and HOMO-LUMO energy gap. Electron-donors, which have lower σR values, redshift the absorption spectrum and reduce the HOMO-LUMO energy gap. Conversely, electron-withdrawing groups (higher σR's) blueshift the spectrum and increase the energy gap. The second band maximum energies, studied here for the first time, showed no correlation with σ but tended to increase with σ. A comprehensive charge transfer (CT) analysis of the main transition of all systems was also carried out. We found that donors (lower σ's) slightly enhance the CT character of the unsubstituted pNA, whereas acceptors (higher σ's) decrease it, leading to increased local excitations within the aromatic ring. The overall CT variation is not large, except for pNA-SO2H, which considerably decreases the total CT value. We found that the strong electron donors pNA-OH, pNA-OCH3, and pNA-NH2, which have the smallest HOMO-LUMO energy gaps and lowest σ's, have potential for optoelectronic applications. The results show that none of the studied molecules is fluorescent in the gas phase. However, pNA-NH2 and pNA-COOH in cyclohexane and water reveal fluorescence upon solvation. METHODS: We investigated theoretically employing the second-order algebraic diagrammatic construction (ADC(2)) ab initio wave function and time-dependent density functional theory (TDDFT) the gas-phase electronic absorption spectra of 15 derivatives of p-nitroaniline (pNA). The investigated substituents include both electron-donor (C6H5, CCH, CH3, NH2, OCH3, and OH,) and electron-withdrawing (Br, CCl3, CF3, Cl, CN, COOH, F, NO2, and SO2H) substituents.

RESUMO

CONTEXT: An accurate description of the molecular charge density is crucial for investigating intra- and inter-molecular properties. Among the different ways of describing and analyzing it, the widely used distributed multipole analysis (DMA) is an accurate method for decomposing the molecular charge density into atom-centered electric multipoles (monopole, dipole, quadrupole, and so on) that have a direct chemical interpretation. In this work, DMA was employed to decompose the molecular charge density of six chemically distinct molecules, namely, (2R)-2-amino-3-[(S)-prop-2-enylsulfinyl] propanoic acid (AAP), 4-amine-2-nitro-1,3,5 triazole (ANTA), (RS)-Propan-2-yl methylphosphonofluoridate (SARIN), chloromethane (CLMET), and 2-aminoacetic acid (GLY) into monopole, dipole, and quadrupole values. A hypothetical variation of ANTA built by exchanging all the nitrogen atoms with phosphorus that we named 4-phosphine-2-phosphite-1,3,5-phosphorine (ANTAP) was also studied. These molecules have different chemical structures bearing distinct carbon skeletons, electronegative atoms, and electron-withdrawing/donating groups. We found that although DFT multipole values can depend considerably on the exchange-correlation functional for specific atomic sites, the associated root-mean-square errors (RMSEs) compared to benchmark MP4 mainly were about [Formula: see text] The most significant variations were for monopoles and dipoles of sites highly polarized by adjacent atoms, and to a lesser degree, for the quadrupoles. The double hybrid B2PLYP and the hybrid meta M06-2X functionals, as expected in the framework of Jacob's ladder, overall give the most accurate results among the DFT methods. The MP2 DMA multipole values have an RMSE in relation to the MP4 benchmark mainly in the range [Formula: see text], thus representing a lower computational cost to obtain results with similar good accuracy without the ambiguity of choosing a DFT functional. The deviations of the HF multipoles from the benchmark in most cases were less than 20%, in agreement with the well-known fact that non-correlated charge densities have a slight dependence on the electronic correlation. We also confirmed that DMA values have a small dependence on the size of the basis set: deviations did not exceed 5% in most cases. However, the dependence of the DMA values on the size of the basis set increases with the rank of the electric multipole. To compute accurate values of DMA multipoles of an atom bonded to very electronegative atoms, especially dipoles and quadrupoles, a large basis set including diffuse functions is necessary. Despite that, for a given polarized basis set, the choice of the basis set to compute accurate DMA multipole values is not critical. METHODS: The molecular charge densities were computed using the electronic structure methods Hartree-Fock (HF), MP2, MP4, DFT/PBE, DFT/B3LYP, DFT/B3PW91, DFT/M06-2X, and DFT/B2PLYP implemented in the Gaussian 09 package. MP4 was the benchmark method. The DMA multipoles were obtained with the GDMA program of Stone. The 6-311G + + (d,p) basis set was used for the production calculations, and the augmented correlation-consistent Dunning's hierarchy of basis sets was employed to evaluate the dependence of the DMA multipoles on the basis set size.

RESUMO

Two-dimensional (2D) materials formed by thin-films of metal oxides that grow on metal supports are commonly used in heterogeneous catalysis and multilayer electronic devices. Despite extensive research on these systems, the effects of charged defects at supported oxides on surface processes are still not clear. In this work, we perform spin-polarized density-functional theory (DFT) calculations to investigate formation and interaction of charged magnesium and oxygen vacancies, and Al dopants on MgO(001)/Ag(001) surface. The results show a sizable interface compressive effect that decreases the metal work function as electrons are added on the MgO surface with a magnesium vacancy. This surface displays a larger formation energy in a water environment (O-rich condition) even with additional Al-doping. Under these conditions, we found that a polar molecule such as CO is more strongly adsorbed on the low-coordination oxygen sites due to a larger contribution of the channeled electronic transport with the silver interface regardless of the surface charge. Therefore, these findings elucidate how surface intrinsic vacancies can influence or contribute to charge transfer, which allows one to explore more specific reactions at different surface topologies for more efficient catalysts for CO2 conversion.

RESUMO

The alternant polycyclic aromatic hydrocarbon pyrene has photophysical properties that can be tuned with different donor and acceptor substituents. Recently, a D (donor)-Pyrene (bridge)-A (acceptor) system, DPA, with the electron donor N,N-dimethylaniline (DMA), and the electron acceptor trifluoromethylphenyl (TFM), was investigated by means of time-resolved spectroscopic measurements (J. Phys. Chem. Lett. 2021, 12, 2226-2231). DPA shows great promise for potential applications in organic electronic devices. In this work, we used the ab initio second-order algebraic diagrammatic construction method ADC(2) to investigate the excited-state properties of a series of analogous DPA systems, including the originally synthesized DPAs. The additionally investigated substituents were amino, fluorine, and methoxy as donors and nitrile and nitro groups as acceptors. The focus of this work was on characterizing the lowest excited singlet states regarding charge transfer (CT) and local excitation (LE) characters. For the DMA-pyrene-TFM system, the ADC(2) calculations show two initial electronic states relevant for interpreting the photodynamics. The bright S1 state is locally excited within the pyrene moiety, and an S2 state is localized ~0.5 eV above S1 and characterized as a donor to pyrene CT state. HOMO and LUMO energies were employed to assess the efficiency of the DPA compounds for organic photovoltaics (OPVs). HOMO-LUMO and optical gaps were used to estimate power conversion and light-harvesting efficiencies for practical applications in organic solar cells. Considering the systems using smaller D/A substituents, compounds with the strong acceptor NO2 substituent group show enhanced CT and promising properties for use in OPVs. Some of the other compounds with small substituents are also found to be competitive in this regard.

RESUMO

Diketopyrrolopyrrole (DPP) systems have promising applications in different organic electronic devices. In this work, we investigated the effect of 20 different substituent groups on the optoelectronic properties of DPP-based derivatives as the donor ( D )-material in an organic photovoltaic (OPV) device. For this purpose, we employed Hammett's theory (HT), which quantifies the electron-donating or -withdrawing properties of a given substituent group. Machine learning (ML)-based σ m , σ p , σ m 0 , σ p 0 , σ p + , σ p - , σ I , and σ R Hammett's constants previously determined were used. Mono- (DPP-X1 ) and di-functionalized (DPP-X2 ) DPPs, where X is a substituent group, were investigated using density functional theory (DFT), time-dependent DFT (TDDFT), and ab initio methods. Several properties were computed using CAM-B3LYP and the second-order algebraic diagrammatic construction, ADC(2), an ab initio wave function method, including the adiabatic ionization potential ( I P A ), the electron affinity ( E A A ), the HOMO-LUMO gaps ( E g ), and the maximum absorption wavelengths ( λ max ), the first excited state transition 1 S0 → 1 S1 energies ( ∆ E ) (the optical gap), and exciton binding energies. From the optoelectronic properties and employing typical acceptor systems, the power conversion efficiency ( PCE ), open-circuit voltage ( V OC ), and fill factor ( FF ) were predicted for a DPP-based OPV device. These photovoltaic properties were also correlated with the machine learning (ML)-based Hammett's constants. Overall, good correlations between all properties and the different types of σ constants were obtained, except for the σ I constants, which are related to inductive effects. This scenario suggests that resonance is the main factor controlling electron donation and withdrawal effects. We found that substituent groups with large σ values can produce higher photovoltaic efficiencies. It was also found that electron-withdrawing groups (EWGs) reduced E g and ∆ E considerably compared to the unsubstituted DPP-H. Moreover, for every decrease (increase) in the values of a given optoelectronic property of DPP-X1 systems, a more significant decrease (increase) in the same values was observed for the DPP-X2 , thus showing that the addition of the second substituent results in a more extensive influence on all electronic properties. For the exciton binding energies, an unsupervised machine learning algorithm identified groups of substituents characterized by average values (centroids) of Hammett's constants that can drive the search for new DDP-derived materials. Our work presents a promising approach by applying HT on molecular engineering DPP-based molecules and other conjugated molecules for applications on organic optoelectronic devices.

RESUMO

The A-series is the most recent generation of chemical warfare nerve agents (CWA) which act directly on the inhibition of the human acetylcholinesterase (HssAChE) enzyme. These compounds lack accurate experimental data on their physicochemical properties, and there is no evidence that traditional antidotes effectively reactivate HssAChE inhibited by them. In the search for potential antidotes, we employed virtual screening, molecular docking, and molecular dynamics (MD) simulations for the theoretical assessment of the performance of a library of Mannich phenols as potential reactivators of HssAChE inhibited by the Novichok agents A-230, A-232, and A-234, in comparison with the commercial oximes pralidoxime (2-PAM), asoxime (HI-6), trimedoxime (TMB-4), and obidoxime. Following the near-attack conformation (NAC) approach, our results suggest that the compounds assessed would face difficulties in triggering the proposed nucleophilic in-line displacement mechanism. Despite this, it was observed that certain Mannich phenols presented similar or superior results to those obtained by reference oximes against A-232 and A-234 model, suggesting that these compounds can adopt more favourable conformations. Additional binding energy calculations confirmed the stability of the model/ligands complexes and the reactivating potential observed in the molecular docking and MD studies. Our findings indicate that the Mannich phenols could be alternative antidotes and that their efficacy should be evaluated experimentally against the A-series CWA.

Assuntos

Substâncias para a Guerra Química , Reativadores da Colinesterase , Agentes Neurotóxicos , Humanos , Antídotos/farmacologia , Reativadores da Colinesterase/farmacologia , Acetilcolinesterase/metabolismo , Simulação de Acoplamento Molecular , Inibidores da Colinesterase/farmacologia , Inibidores da Colinesterase/química , Oximas/farmacologia , Oximas/química , Trimedoxima/química , Trimedoxima/farmacologia , Substâncias para a Guerra Química/farmacologia , Compostos de Piridínio/farmacologia

RESUMO

Hammett's constants σ quantify the electron donor or electron acceptor power of a chemical group bonded to an aromatic ring. Their experimental values have been successfully used in many applications, but some are inconsistent or not measured. Therefore, developing an accurate and consistent set of Hammett's values is paramount. In this work,we employed different types of machine learning (ML) algorithms combined with quantum chemical calculations of atomic charges to predict theoretically new Hammett's constants σm, σp, σm0, σp0, σp+, σp-, σR, and σI for 90 chemical donor or acceptor groups. New σ values (219), including previously unknown 92, are proposed. The substituent groups were bonded to benzene and meta- and para-substituted benzoic acid derivatives. Among the charge methods (Mulliken, Löwdin, Hirshfeld, and ChelpG), Hirshfeld showed the best agreement for most kinds of σ values. For each type of Hammett constant, linear expressions depending on carbon charges were obtained. The ML approach overall showed very close predictions to the original experimental values, with meta- and para-substituted benzoic acid derivative values showing the most accurate values. A new consistent set of Hammett's constants is presented, as well as simple equations for predicting new values for groups not included in the original set of 90.

RESUMO

We decomposed density functional theory charge densities of 53 nitroaromatic molecules into atom-centered electric multipoles using the distributed multipole analysis that provides a detailed picture of the molecular electronic structure. Three electric multipoles, (the charge of the nitro groups), (the total dipole, i.e., polarization, of the nitro groups), (the total electron delocalization of the C ring atoms), and the number of explosophore groups (#NO2) were selected as features for a comprehensive machine learning (ML) investigation. The target property was the impact sensitivity h50 (cm) values quantified by drop-weight measurements, with a large h50 (e.g., 150 cm) indicating that an explosive is insensitive and vice versa. After a preliminary screening of 42 ML algorithms, four were selected based on the lowest root mean square errors: Extra Trees, Random Forests, Gradient Boosting, and AdaBoost. Compared to experimental data, the predicted h50 values of molecules having very different sensitivities for the four algorithms have differences in the range 19-28%. The most important properties for predicting h50 are the electron delocalization in the ring atoms and the polarization of the nitro groups with averaged weights of 39% and 35%, followed by the charge (16%) and number (10%) of nitro groups. A significant result is how the contribution of these properties to h50 depends on their actual sensitivities: for the most sensitive explosives (h50 up to ∼50 cm), the four properties contribute to reducing h50, and for intermediate ones (∼50 cm ≲ h50 ≲ 100 cm) #NO2 and contribute to increasing it and the other two properties to reducing it. For highly insensitive explosives (h50 ≳ 200 cm), all four properties essentially contribute to increasing it. These results furnish a consistent molecular basis of the sensitivities of known explosives that also can be used for developing safer new ones.

RESUMO

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S2) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S2 emission was not observed in some experiments.

RESUMO

Silole derivatives have been extensively employed for developing organic optoelectronics, but few studies focused on the photophysical properties of the silole molecule. In this work, we investigate these properties by computing the absorption spectra and performing nonadiabatic molecular dynamics of silole employing the algebraic diagrammatic construction [ADC(2)] and extended multi-state XMS-CASPT2 ab initio electronic structure methods. For vertical excitations and excited state optimizations, the equation of motion coupled-cluster singles and doubles (EOM-CCSD) was also used. The nuclear ensemble and the fewest-switches surface hopping molecular dynamics methods coupled with the first two high-level electronic structure methods were applied to probe the relaxation mechanisms of silole. We could reproduce the experimental first absorption maximum value and found an ultrafast relaxation process occurring exclusively through ring-puckering distortions without breaking ring bonds or hydrogen elimination. Minimum energy conical intersection optimizations were carried out and potential energy curves, including triplet states, were calculated to further elucidate the relaxation process of silole.