Time-Dependent Multilevel Density Functional Theory

We present a novel three-layer approach based on multilevel density functional theory (MLDFT) and polarizable molecular mechanics to simulate the electronic excitations of chemical systems embedded in an external environment within the time-dependent DFT formalism. In our method, the electronic structure of a target system, the chromophore, is determined in the field of an embedded inactive layer, which is treated as frozen. Long-range interactions are described by employing the polarizable fluctuating charge (FQ) force field. The resulting MLDFT/FQ thus accurately describes both electrostatics (and polarization) and non-electrostatic target-environment interactions. The robustness and reliability of the approach are demonstrated by comparing our results with experimental data reported for various organic molecules in solution.

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

Antitumor efficacy of synthesized Ag-Au nanocomposite loaded with PEG and ascorbic acid in human lung cancer stem cells

Lung cancer is the leading cause of mortality worldwide of which non-small cell lung carcinoma constitutes majority of the cases. High mortality is attributed to early metastasis, late diagnosis, ineffective treatment and tumor relapse. Chemotherapy and radiotherapy form the mainstay of its treatment. However, their associated side effects involving kidneys, nervous system, gastrointestinal tract, and liver further adds to dismal outcome. These disadvantages of conventional treatment can be circumvented by use of engineered nanoparticles for improved effectiveness with minimal side effects. In this study we have synthesized silver gold nanocomposite (Ag-Au NC) using polyethylene glycol and l-ascorbic acid as surfactant and reducing agent respectively. Synthesized nanocomposite was characterized by ultraviolet-visible absorption, dynamic light scattering, scanning and transmission electron microscopy. Compositional analysis was carried out by energy dispersive X-ray analysis and average pore diameter was estimated using Barrett-Joyner-Halenda method. In-silico molecular docking analysis of the synthesized NC against active regions of epidermal growth factor receptor revealed good binding energy. Subsequently, we investigated the effect of NC on growth and stem cell attributes of A549 lung cancer cells. Results showed that NC was effective in inhibiting A549 cell proliferation, induced DNA damage, G2/M phase arrest and apoptosis. Further, tumor cell migration and spheroid formation were also negatively affected. NC also enhanced reactive oxygen species generation and mitochondrial depolarization. In addition, the effect of NC on putative cancer stem cells in A549 cells was evaluated. We found that Ag-Au NC at IC50 targeted CD44, CD24, CD166, CD133 and CD326 positive cancer stem cells and induced apoptosis. CD166 positive cells were relatively resistance to apoptosis. Together our results demonstrate the anticancer efficacy of Ag-Au NC mediated by a mechanism involving cell cycle arrest and mitochondrial derangement.

Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment

In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.

Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.

Solute - solvent repulsion effects on the absorption spectra of anthracene in n-hexane investigated under high pressure

The band positions in the UV-VIS absorption spectra of compressed solution of anthracene in n-hexane significantly depend not only on the dispersive but also on the repulsive solute-solvent interactions, what has so far been omitted. Their strength is determined not only by the solvent polarity but also by Onsager cavity radius changing with pressure. The results obtained for anthracene show that repulsive interactions should be included in the interpretation of barochromic and solvatochromic results of aromatic compounds. We show that the barochromic studies in the liquid solvent can be an alternative to solvatochromic studies, e.g. to determine the polarizability of organic molecules in the electronic excited state. The pressure-induced polarity change in n-hexane exceeds that induced by the exchange of n-alkane solvents between n-pentane and n-hexadecane.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening

The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ^★ S₀ → S₁ transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm^-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.

Charge-Transfer Luminescence in a Molecular Donor-Acceptor Complex: Computational Insights

Donor-acceptor molecular complexes are a popular class of materials utilizing charge-transfer states for practical applications. A recent class of donor-acceptor dyads based on the fluorescent BODIPY functionalized with triphenylamine (TPA) shows the peculiar property of dual fluorescence. It is hypothesized that instead of the sensitized charge-transfer state being optically dark, it provides an additional bright radiative pathway. Here we use time-dependent density functional theory to characterize the energetic alignment of excitonic and charge-transfer states in a BODIPY-TPA molecular complex. We observe that using a long-range exchange corrected functional in combination with state-specific solvation scheme gives a qualitatively correct alignment of the exciton and charge-transfer states and an enhancement in oscillator strength for the equilibrium solvated charge-transfer state, in agreement with experiment. This work provides rationalization of charge-transfer state emission and provides a foundation to explore charge-transfer using ab initio excited-state nonadiabatic dynamics.

Accurate predictions of the electronic excited states of BODIPY based dye sensitizers using spin-component-scaled double-hybrid functionals: a TD-DFT benchmark study

The vertical excitation energies of 13 BODIPY based dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals from different DFT rungs. Most TD-DFT results were found to overestimate the excitation energies, and show mean absolute error (MAE) values in the range 0.2–0.5 eV. The dispersion-corrected, spin-component-scaled, double-hybrid (DSD) functionals DSD-BLYP and DSD-PBEP86 were found to have the smallest MAE values of 0.083 eV and 0.106 eV, respectively, which is close to the range of average errors found in the more expensive coupled-cluster methods. Moreover, DSD-BLYP and DSD-PBEP86 functionals show excellent consistency and quality of results (standard deviation = 0.048 eV and 0.069 eV respectively). However, the range separated hybrid (RSH) and the range separated double hybrid (RSDH) functionals were found to provide the best predictability (linear determination coefficient R² > 0.97 eV).

The excitation energies of 13 BODIPY dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals. Spin-component-scaled double-hybrid (DSD) functionals are found to show the best performance.

Assessment of time-dependent density functionals for the electronic excitation energies of organic dyes used in DSSCs

The absorption spectra modeled as the vertical excitation energies of 13 dye sensitizers used in dye-sensitized solar cells (DSSCs) are benchmarked by means of time-dependent (TD)-DFT, using 36 functionals from different DFT rungs.

Single-ion solvation free energy: A new cluster-continuum approach based on the cluster expansion method

Accurate calculation of the solvation free energy of single ions remains an important goal, involving development in the dielectric continuum solvation models, and statistical mechanics with explicit solvent and hybrid discrete-continuum methods. In the last case, many of the research studies involve a quasi-chemical approach using the monomer cycle or the cluster cycle to calculate the solvation free energy of single ions. In this work, a different cluster-continuum approach based on the cluster expansion method was tested for solvation of 16 cations and 32 anions in aqueous solution. The SMD model was used for the dielectric continuum part and three explicit water molecules were introduced in the region of the solute with the highest interaction energy. Harmonic frequency calculations and molecular dynamics sampling of configurations are not required. An empirical γ_N parameter for cations and another for anions is introduced. The method produces a substantial improvement of the SMD model with a mean absolute deviation of 2.3 kcal mol^-1 for cations and 2.9 kcal mol^-1 for anions. The analysis of the correlation between theoretical and experimental data produces a linear regression line with a slope of 1.09 for cations and 1.01 for anions. The good results of this approximated cluster expansion approach suggest that the method could be further improved by including more solvent molecules and sampling the configurations.

Are all charge-transfer parameters created equally? A study of functional dependence and excited-state charge-transfer quantification across two dye families

Small molecule organic dyes have many potential uses in medicine, textiles, forensics, and light-harvesting technology. Being able to computationally predict the spectroscopic properties of these dyes could greatly expedite screening efforts, saving time and materials. Time-dependent density functional theory (TD-DFT) has been shown to be a good tool for this in many instances, but characterizing electronic excitations with charge-transfer (CT) character has historically been challenging and can be highly sensitive to the chosen exchange-correlation functional. Here we present a combined experimental and computational study of the excited-state electronic structure of twenty organic dyes obtained from the Max Weaver Dye Library at NCSU. Results of UV-vis spectra calculations on these dyes with six different exchange-correlation functionals, BP86, B3LYP, PBE0, M06, BH and HLYP, and CAM-B3LYP, were compared against their measured UV-vis spectra. It was found that hybrid functionals with modest amounts (20-30%) of included Hartree-Fock exchange are the most effective at matching the experimentally determined λ_max. The interplay between the observed error, the functional chosen, and the degree of CT was analyzed by quantifying the CT character of λ_max using four orbital and density-based metrics, Λ, Δr, S_C and D_CT, as well as the change in the dipole moment, Δμ. The results showed that the relationship between CT character and the functional dependence of error is not straightforward, with the observed behavior being dependent both on how CT was quantified and the functional groups present in the molecules themselves. It is concluded that this may be a result of the examined excitations having intermediate CT character. Ultimately it was found that the nature of the molecular "family" influenced how a given functional behaved as a function of CT character, with only two of the examined CT quantification methods, Δr and D_CT, showing consistent behavior between the different molecular families. This suggests that further work needs to be done to ensure that currently used CT quantification methods show the same general trends across large sets of multiple dye families.

Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview

The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.

Simple Protocol for Capturing Both Linear-Response and State-Specific Effects in Excited-State Calculations with Continuum Solvation Models

We present an effective computational protocol (cLR²) to describe both solvatochromism and fluorosolvatochromism. This protocol, which couples the polarizable continuum model to time-dependent density functional theory, simultaneously accounts for both linear-response and state-specific solvation effects. A series of test cases, including solvatochromic and fluorosolvatochromic compounds and excited-state intramolecular proton transfers, are used to highlight that cLR² is especially beneficial for modeling bright excitations possessing a significant charge-transfer character, as well as cases in which an accurate balance between states of various polarities should be restored.

A polarizable three-layer frozen density embedding/molecular mechanics approach

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Thermally activated delayed fluorescence: A critical assessment of environmental effects on the singlet-triplet energy gap

The effective design of dyes optimized for thermally activated delayed fluorescence (TADF) requires the precise control of two tiny energies: the singlet-triplet gap, which has to be maintained within thermal energy, and the strength of spin-orbit coupling. A subtle interplay among low-energy excited states having dominant charge-transfer and local character then governs TADF efficiency, making models for environmental effects both crucial and challenging. The main message of this paper is a warning to the community of chemists, physicists, and material scientists working in the field: the adiabatic approximation implicitly imposed to the treatment of fast environmental degrees of freedom in quantum-classical and continuum solvation models leads to uncontrolled results. Several approximation schemes were proposed to mitigate the issue, but we underline that the adiabatic approximation to fast solvation is inadequate and cannot be improved; rather, it must be abandoned in favor of an antiadiabatic approach.

The vibronic absorption spectra and electronic states of acridine orange in aqueous solution

The time-dependent density functional theory (TD-DFT) was used to obtain vibronic absorption spectra of acridine orange dye (AO) in an aqueous solution that were in good agreement with the experiment. The protonated and neutral forms of the dye have been investigated. The results of calculations using various functionals and basis sets have been analyzed. The best agreement with experiment was given by the level of theory X3LYP/6-31G(d,p). AO molecular orbitals involved in electronic transitions due light absorption in the visible region of the spectrum have been obtained. The dipole moments and atomic charges of the ground and excited states of the AO molecule have been calculated. Maps of the electrostatic potential have been drawn. An insignificant photoinduced electron transfer was found in the central ring of the chromophore of the dye molecule. According to our calculations, the vibronic coupling and the Boltzmann distribution play a significant role in the absorption spectra of the AO.

Temperature-dependent oxidation of BSCAPE molecule in methanol medium

Temperature-dependent solvation free energy and oxidation by free energy of ionization of 2-Phenylethyl (2E)-3-(1-benzenesulfonyl-4,5-dihydroxyphenyl) acrylate (BSCAPE) in methanol medium are the concerns of the present work. This molecule is a relevant phenolic acid enclosing multiple bioactivities. The explicit, implicit and discrete-continuum models of solvation were used. The methanol molecules were coordinated to this acid to form cluster complexes. The dual method M06-2X/6-31++G(d,p)//B3LYP/6-31G(d) was employed along with basis set superposition error correction. The results show that, the free energy of coordination and solvation are distant. Both quantities increase with temperature. From discrete-continuum treatment, there is non-spontaneity of solvation process, while coordination yielded spontaneity and non-spontaneity at cold and hot room temperatures, respectively. The ionization potential in gas phase, decreases with temperature. All the solvation models yielded lower ionization potential than that of gas phase. Thus, it follows that, the increase of temperature and methanol medium favours the oxidation of BSCAPE. Consequently, this favours its metabolism processes.

Designing of benzodithiophene core-based small molecular acceptors for efficient non-fullerene organic solar cells

Nowadays, organic solar cells (OSCs) with non-fullerene electron acceptors provide the highest efficiencies among all studied OSCs. To further improve the efficiencies of fullerene-free organic solar cells, end-capped acceptor modification is made with strong electron withdrawing groups. In this report, we have theoretically designed five new novel Benzodithiophene core-based acceptor molecules (H1-H5) with the aim to study the possible enhancement in photophysical, optoelectronic, and photovoltaic properties of newly designed molecules. The end-capped acceptor modification of famous and recently synthesized FBDIC molecule has been made with strong electron withdrawing groups. Density functional theory and time-dependent-density functional theory are extensively used to study the structural-property relationship, optical properties and various geometrical parameters like frontier molecular orbitals alignment, excitation and binding energy, transition density matrix along with open circuit voltage, density of states and dipole moment. Commonly, low reorganization energies (hole and electron) afford high charge mobility and our all designed systems are enriched in aspect (λ_e = 0.0044-0.0104 eV and λ_h = 0.0060-0.0090 eV). Moreover, H1-H5 molecules demonstrate red-shifting in absorption spectrum (λ_max = 741-812 nm) as compare to R (λ_max = 728 nm). Low excitation and binding energies with low HOMO (highest occupied molecular orbital)-LUMO (lowest unoccupied molecular orbital) energy gap of H1-H5 suggested that designed molecules are better and suitable candidates for high performance organic solar cell. Results of all analysis indicate that this theoretical framework demonstrates that end-capped acceptors modification is a simple and effective alternative strategy to achieve the desirable optoelectronic properties. Therefore, H1-H5 are recommended to experimentalist for out-looking future developments of highly efficient solar cells.

Understanding TADF: a joint experimental and theoretical study of DMAC-TRZ

TADF offers a promising way to harvest triplets in OLED for improved efficiency. To concurrently optimize the dye inside the matrix, a thorough experimental and theoretical study is presented of a the TADF dye addressing environmental effects.

Combining multilevel Hartree–Fock and multilevel coupled cluster approaches with molecular mechanics: a study of electronic excitations in solutions

We present the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable.

Simulating Absorption Spectra of Flavonoids in Aqueous Solution: A Polarizable QM/MM Study

We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.

Designing of near-infrared sensitive asymmetric small molecular donors for high-efficiency organic solar cells

Herein, we have designed four small molecular donors (SMDs) with Donor–Acceptor–Acceptor (D–Á–A) backbone having different acceptor units for highly efficient organic solar cells (OSCs). The specific molecular modeling has been made by replacing the additional acceptor unit (A) of recently synthesized TPA-DAA-MDN molecule (R) by employing different highly efficient acceptor units in order to improve the photovoltaic performances of the molecules. A theoretical approach (DFT and TD-DFT) has been applied to investigate the photophysical, opto-electronic and photovoltaic parameters of the designed molecules (DAA1–DAA4) and compared with the reference molecule (R). The red-shifting absorption of SMDs is the most important factor for highly efficient OSCs. Our all formulated molecules showed a red shifted absorption spectrum and also exhibit near IR sensitivity. Acceptor unit modification of R molecule causes reduction in HOMO-LUMO energy gap; therefore, all designed molecules offer better opto-electronic properties as compared to R molecule. A variety of certain critical factors essential for efficient SMDs like frontier molecular orbitals (FMOs), absorption maxima, dipole moment, exciton binding energy along with transition density matrix, excitation energy, open circuit voltages and charge mobilities of (DAA1–DAA4) and R have also been investigated. Generally, low values of reorganizational energy (hole and electron) offer high charge mobility and our all designed molecules are enriched in this aspect. High open circuit voltage values, low excitation energies, large dipole moment values indicate that our designed SMDs are suitable candidates for high-efficiency OSCs. Furthermore, conceptualized molecules are superior and thus are suggested to experimentalist for out-looking future progresses of highly efficient OSCs devices.

Optical Absorption of Atomically-Precise Sn14 Nanoclusters: The Antagonistic Interplay of Ligand Stabilization, Molecular Symmetry, and Solvatochromism

The synthesis of atomically precise clusters is nowadays well established. The study of isolated clusters in the gas phase has also become an approved field of research. Although both approaches examine the same research objects, namely nanoclusters, little is known about to what extent results from gas phase studies can be transferred to colloidal systems and vice versa. In particular, it is not yet sufficiently understood how ligands influence the geometric and electronic structure of clusters from an experimental point of view. By comparing a ligand-stabilized tin nanocluster in solution with an isolated species in the gas phase and considering different geometric arrangements with the same number of tin atoms, the impacts of ligand stabilization, molecular symmetry, and solvatochromism on the optical behavior are thoroughly worked out for the first time.

A quantum study on novel azo-dyes containing a fullerene C60 unit as a smart material for optoelectronic applications

Quantum chemical calculations of some novel azo-dyes containing a fullerene C60 unit as a smart material have been carried out with the aims to determine their cis and trans electronic properties and to describe the change of their quantum parameters as a result of the trans/cis isomerization of these molecules. The effects of electron-withdrawing or electron-releasing groups on the R-position of these molecules on electronic, optical, spectroscopic, and other properties of these molecules have been considered with DFT and TDDFT calculations. The obtained results of the calculations show that compounds "b" and "c" with the strongest electron-releasing groups in the R-position of these molecules, particularly the trans isomers of these compounds, with higher chemical softness, higher electrophilicity index, higher thermodynamic properties, and higher charge transfer values, have the better electronic and optical properties and therefore the better chemical reactivity compared to the other compounds. Graphical abstract.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Antiadiabatic View of Fast Environmental Effects on Optical Spectra

An antiadiabatic approach is proposed to model how the refractive index of the surrounding medium affects optical spectra of molecular systems in condensed phases. The approach solves some of the issues affecting current implementations of continuum solvation models and more generally of effective models where a classical description is adopted for the molecular environment.

Zinc-Doped Boron Phosphide Nanocluster as Efficient Sensor for SO2

Adsorption of SO2 on pure B12P12 and Zn-doped B12P12 is investigated through density functional theory methods. Zn adsorption on BP delivers four optimized geometries: B-Top, P-top, b64, and ring-enlarged geometry with adsorption energies of −57.12 kJ/mol, −14.50 kJ/mol, −22.94 kJ/mol, and −14.83 kJ/mol, respectively. The adsorption energy of SO2 on pristine boron phosphide is −14.92 kJ/mol. Interaction of SO2 with Zn-doped boron phosphide gives four different geometries with adsorption energies of −69.76 kJ/mol, −9.82 kJ/mol, −104.92 kJ/mol, and −41.87 kJ/mol. Geometric parameters such as dipole moment, QNBO, frontier molecular orbital analysis, PDOS, and global indices of reactivity are performed to visualize the changes in electronic properties of B12P12 after Zn and SO2 adsorption.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Development of a 1,3a,6a-triazapentalene derivative as a compact and thiol-specific fluorescent labeling reagent

For the fluorescence imaging of biologically active small compounds, the development of compact fluorophores that do not perturb bioactivity is required. Here we report a compact derivative of fluorescent 1,3a,6a-triazapentalenes, 2-isobutenylcarbonyl-1,3a,6a-triazapentalene (TAP-VK1), as a fluorescent labeling reagent. The reaction of TAP-VK1 with various aliphatic thiols proceeds smoothly to afford the corresponding 1,4-adducts in high yields, and nucleophiles other than thiols do not react. After the addition of thiol groups in dichloromethane, the emission maximum of TAP-VK1 shifts to a shorter wavelength and the fluorescence intensity is substantially increased. The utility of TAP-VK1 as a compact fluorescent labeling reagent is clearly demonstrated by the labeling of Captopril, which is a small molecular drug for hypertension. The successful imaging of Captopril, one of the most compact drugs, in this study demonstrates the usefulness of compact fluorophores for mechanistic studies.

Optical spectra of organic dyes in condensed phases: the role of the medium polarizability

An effective model is presented to account for the effects of the medium electronic polarizability on spectral properties and on symmetry-breaking phenomena in charge-transfer dyes.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes

The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.

Quantum Confinement Effects on Solvatochromic Shifts of Molecular Solutes

We demonstrate the pivotal role of quantum mechanics density confinement effects on solvatochromic shifts. In particular, by resorting to a quantum mechanics/molecular mechanics (QM/MM) approach capable of accounting for confinement effects we successfully reproduce vacuo-to-water solvatochromic shifts for dark n → π* and bright π → π* transitions of acrolein and dark n → π* transitions of pyridine and pyrimidine without the need of including explicit water molecules in the QM portion. Remarkably, our approach is also able to dissect the effects of the single forces acting on the solute-solvent couple and allows for a rationalization of the experimental findings in terms of physicochemical quantities.