Absorption and Fluorescence Spectra of Uracil in the Gas Phase and in Aqueous Solution: A TD-DFT Quantum Mechanical Study

Chalcone-based Turn-Off Chemosensor for Selective and Susceptible Detection of Fe2+ Ions: Spectroscopic and DFT Investigations

Herein, in this report we are introducing newly synthesized chalcone derivative, "(E)-1-phenyl-3-(4-((5-(((Z)-thiophen-2-ylmethylene)amino)-1,3,4-thiadiazol-2-yl)thio)phenyl)prop-2-en-1-one" (5), as a chemosensor to detect Fe2+ metal ions in HEPES buffer solution of pH 7.5. Spectroscopic techniques were used to confirm the synthesized sensor. To determine the chemical reactivity and molecular stability of the probe, a frontier molecular orbitals investigation was carried out. A molecular electrostatic potential map was investigated to know the binding site of 5 for metal ion coordination. The theoretical absorption and fluorescence emission properties were estimated and correlated with the experimental observations. The sensor showed excellent selectivity for Fe2+ compared to all other studied metal ions. The fluorescence binding studies were carried out by adding different amounts of Fe2+ ions for a fixed concentration of probe 5. The inclusion of Fe2+ ions resulted in a decrease in fluorescence intensity with a bathochromic shift of emission wavelength of 5 due to the 5-Fe2+ complexation. The binding affinity value for the probe was found to be 576.2 M-1 with the help of the Stern-Volmer plot. The Job's plot and mass spectra supported the 2:1 (5: Fe2+) stoichiometry of complex formation. The detection limit and limit of quantification of 5 for Fe2+ were calculated to be 4.79 × 10-5 M and 14.54 × 10-5 M. Further, in addition to this, the photophysical parameters such as fluorescence lifetime of 5 and 5-Fe2+ complex measured to be 0.1439 and 0.1574 ns. The quantum yield of 5 and 5-Fe2+ was found to be 0.0398 and 0.0376. All these experimental findings revealed that probe 5 has excellent selectivity and sensitivity for Fe2+ ions.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain onπ π * state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained byn π * trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of theπ π * state leading to a wrongπ π * / n π * energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Is 1-methylcytosine a faithful model compound for ultrafast deactivation dynamics of cytosine nucleosides in solution?

1-Methylcytosine (1mCyt) is the base for nucleoside N1-methylpseudodeoxycytidine of Hachimoji nucleic acids and a frequently used model compound for theoretical studies on excited states of cytosine nucleosides. However, there is little experimental characterization of spectra and photo-dynamic properties of 1mCyt. Herein, we report a comprehensive investigation into excited state dynamics and effects of solvents on fluorescence dynamics of 1mCyt in both water and acetonitrile. The study employed femtosecond broadband time-resolved fluorescence, transient absorption, and steady-state spectroscopy, along with density functional theory and time-dependent density functional theory calculations. The results obtained provide the first experimental evidence for identifying a dark-natured ∼5.7 ps lifetime nπ* state in the ultrafast non-radiative deactivation with 1mCyt in aqueous solution. This study also demonstrates a significant effect of the solvent on 1mCyt's fluorescence emission, which highlights the crucial role of solute-solvent hydrogen bonding in altering structures and reshaping the radiative as well as nonradiative dynamics of the 1mCyt's ππ* state in the aprotic solvent compared to the protic solvent. The solvent effect exhibited by 1mCyt is distinctive from that known for deoxycytidine, indicating the need for caution in using 1mCyt for modelling the ultrafast dynamics of Cyt nucleosides in solvents with varying properties. Overall, our study unveils a deactivation mechanism that confers a high degree of photo-stability for 1mCyt in solution, shedding light on the molecular basis for solvent-induced effects on the excited state dynamics of nucleobases and derivatives.

Role of Inner Molecular Orbitals in High-Harmonic Generation Spectra of Aligned Uracil

In this work, we decompose the high-harmonic generation (HHG) signal of aligned gas-phase uracil into single molecular-orbital (MO) contributions. We compute HHG spectra for a pulse linearly polarized perpendicular to the molecular plane, with an intensity of 0.6 and 0.85 × 1014 W/cm2 and a wavelength of 800 nm. We use the real-time time-dependent Configuration Interaction with singles method, coupled to a Gaussian-based representation of the time-dependent wavefunction. The strong-field dynamics is affected by the energy of the ionization/recombination channels and by the coupling between the orbital symmetry and laser polarization. In the configuration studied here, we expect that π-type MOs favorably couple with the incoming pulse and play a substantial role in generating the HHG spectrum. Indeed, we show that HOMO, HOMO - 1, and HOMO - 4, which all are π-like, determine the intensity of harmonic peaks at different energies, while HOMO - 2 and HOMO - 3 provide a smaller contribution. It is worth mentioning that HOMO - 4 produces a stronger signal than that from HOMO - 1, even though the corresponding ionization energy, in an one-electron picture, is around 2.5 eV larger and more than 4 eV larger than the HOMO one.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects

We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.

Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials

The direct energy functional minimization method using the orbital transformation (OT) scheme in the program package CP2K has been employed for Δ self-consistent field (ΔSCF) calculations. The OT method for non-uniform molecular orbitals occupations allows us to apply the ΔSCF method for various kinds of molecules and periodic systems. Vertical excitation energies of heteroaromatic molecules and condensed phase systems, such as solvated ethylene and solvated uracil obeying periodic boundary conditions, are reported using the ΔSCF method. In addition, a Re-phosphate molecule attached to the surface of anatase (TiO₂) has been investigated. Additionally, we have implemented a recently proposed state-targeted energy projection ΔSCF algorithm [K. Carter-Fenk and J. M. Herbert, J. Chem. Theory Comput. 16(8), 5067-5082 (2020)] for diagonalization based SCF in CP2K. It is found that the OT scheme provides a smooth and robust SCF convergence for all investigated excitation energies and (non-)periodic systems.

Designing small organic non-fullerene acceptor molecules with diflorobenzene or quinoline core and dithiophene donor moiety through density functional theory

The non-fullerene acceptors A1-A5 with diflourobenzene or quinoline core (bridge) unit, donor cyclopenta[1,2-b:3,4-b']dithiophene unit and 2-(2-methylene-3-oxo-2,3-dihydro-1H-inden-1-ylidene)malononitrile as acceptor unit with additional phenyl, fulvene or thieno[3,2-d]pyrimidinyl 5-oxide groups have been designed through DFT calculations. The optimization of molecular geometries were performed with density functional theory (DFT) at B3LYP 6-31G (d,p) level of theory. The frontier molecular orbital (FMO) energies, band gap energies and dipole moments (ground and excited state) have been calculated to probe the photovoltaic properties. The band gap (1.42-2.01 eV) and dipole moment values (5.5-18. Debye) showed that these designed acceptors are good candidates for organic solar cells. Time-Dependent Density Functional Theory (TD-DFT) results showed λmax (wave length at maximum absorption) value (611-837 nm), oscillator strength (f) and excitation energies (1.50-2.02 eV) in gas phase and in CHCl3 solvent (1.48-1.89 eV) using integral equation formalism variant (IEFPCM) model. The λmax in CHCl3 showed marginal red shift for all designed acceptors compared with gas phase absorption. The partial density of states (PDOS) has been plotted by using multiwfn which showed that all the designed molecules have more electronic distribution at the donor moiety and lowest at the central bridge. The reorganization energies of electron (λe) (0.0007 eV to 0.017 eV), and the hole reorganization energy values (0.0003 eV to - 0.0403 eV) were smaller which suggested that higher charged motilities. The blends of acceptors A1-A5 with donor polymer D1 provided open circuit voltage (Voc) and ∆HOMO off-set of the HOMO of donor and acceptors. These blends showed 1.04 to 1.5 eV values of Voc and 0 to 0.38 eV ∆HOMO off set values of the donor-acceptor bends which indicate improved performance of the cell. Finally, the blend of D1-A4 was used for the study of distribution of HOMO and LUMO. The HOMO were found distributed on the donor polymer (D1) while the A4 acceptor was found with LUMO distribution. Based on λmax values, and band gap energies (Eg), excitation energies (Ex), reorganization energies; the A3 and A4 will prove good acceptor molecules for the development of organic solar cells.

A new oxygen-rich energetic salt dihydrazine tetranitroethide: a promising explosive alternative with high density and good performance

A novel high-energy salt with good oxygen balance, dihydrazine tetranitroethide (5), has been synthesized and characterized by FT-IR spectroscopy, NMR spectroscopy, elemental analysis, and X-ray single crystal diffraction. Compound 5 exhibits high crystal density (1.81 g cm-3) and impressive detonation velocity (9508 m s-1) and detonation pressure (37.9 GPa), showing potential applications as a high performance explosive and a promising additive of propellants.

Molecular designing of naphthalene diimide based fullerene-free small organic solar cell - Acceptors with high photovoltaic performance by density functional theory

With the help of computational chemistry tools, three non-fullerene acceptors, which are 2-methylene-malononitrile (M-1), 2-(3-methyl-5-methylene-2-thioxothiazolidin-4-ylidene) malononitrile (M-2) and 1-methyl-5-methylene-2,6-dioxo-1,2,5,6-tetrahydropyridine-3-carbonitrile (M-3), are designed with naphthalene diimide (NDI) central unit. Their different photovoltaic and optoelectronic properties like absorption spectrum, electrons density, solubility strength, reorganization energies, % ETC from donor to acceptor part, excitation energies, oscillating strength, morphology and crystallinity of device for constructing the thin film bulk hetro junction devices were computed at the WB97XD/6-31 G (d, p) level of density functional theory (DFT). Expected open circuit voltages of designed molecules are high as 4.05 eV to 4.49 eV, which are significantly larger than that of the previously reported 3-methyl-5-methylene-2-thioxothiazolidin-4-one (R) with the value of 3.60 eV at the zero current level. Charge carrier mobilities of designed molecules are high due to having low re-organization energies varying from 0.0163 eV to 0.0280 eV for electron and 0.0160 eV to 0.0190 eV for hole, strong absorption properties between the 420 nm to 550 nm in chloroform and 400 nm to 540 nm in gas phase conditions, respectively.

Solvent effect on the photophysical properties of thermally activated delayed fluorescence molecules

As the third-generation organic electroluminescent materials, thermally activated delayed fluorescence (TADF) molecules have become the research focus recently. Significant solvent effect on TADF molecules were found experimentally, while theoretical investigations are quite limited. In this work, the solvent effect on photophysical properties of DCBPy and DTCBPy are investigated with first-principles calculations. The solvent polarity has slight influence on the molecular geometries and orbitals, while it can decrease the energy gap between the first singlet excited state (S1) and first triplet excited state (T1) significantly. Both the oscillator strength and the radiation rates of S1 increase with larger solvent polarity. The large energy gap between S1 and T1 induce negligible intersystem crossing (ISC) and reverse ISC rates between them, which also indicates higher triplet excited states are involved in the up-conversion process. Our results provide valuable information about solvent influence on the light-emitting properties of TADF molecules, which could help one better understand the light-emitting mechanism of them and favor the design of TADF molecules.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Luminescent and DFT Study of Keto-Enol Tautomers of 5-Fluorouracil and Its Derivatives in Aqueous Solutions

Fluorescence (FL) spectra were recorded and identified for the following 5-fluorouracil (FU) tautomers: 5-fluoropyrimidine-2,4(1H,3H)-dione (A), 5-fluoro-2-hydroxypyrimidine-4(3H)-one (B), and 5-fluoro-4-hydroxypyrimidine-2(1H)-one (D), as well as the corresponding tautomers A and D of 1-(tetrahydrofuranyl-2)-5-fluorouracil and 1-methyl-5-fluorouracil, tautomers A and B of 3-methyl-5-fluorouracil, and the diketo tautomer of 1,3-dimethyl-5-fluorouracil. It was shown that the FL of rare tautomers of FU derivatives occurs by the excitation of uracil homoassociates followed by intramolecular proton transfer (IPT), formation of a pair of rare tautomers, and radiative deactivation of one of them. The FL quantum yields φ_i were estimated. In the case of rare tautomers, the φ_i value includes the probability of IPT in the homoassociate, and density functional theory calculations in the TPSSTPSS/6-311+G(d,p) approximation were used to determine the relative energies of tautomers, associated FU forms, and their spectral properties that agree with the observed luminescence pattern.

A "turn-on" fluorescent probe for the detection of permanganate in aqueous media

A novel carbazole-based probe was synthesized and found to exhibit fluorescence upon introduction to potassium permanganate in aqueous media. Guided by a series of experimental studies and DFT calculations, the underlying mechanism is proposed to proceed through a series of tandem proton and electron transfer processes.

A “turn-on” fluorescent probe for highly selective discrimination of hypochlorite (ClO−) from oxidants including dichromates (Cr2O72−) in aqueous media

A carbazole-based probe was synthesized and found to fluoresce (λ_em,max = 440 nm) in the presence of hypochlorite.

Solvent effects for vertical absorption and emission processes in solution using a self-consistent state specific method based on constrained equilibrium thermodynamics

A self-consistent state specific (SS) method in the framework of TDDFT is presented to account for solvent effects on absorption and emission processes for molecules in solution. In these processes, the initial state is an equilibrium state, while the polarization of the solvent is in nonequilibrium with the electron density of the solute in the final state. Nonequilibrium solvation free energy is calculated based on a novel nonequilibrium solvation model with constrained equilibrium manipulation. The bulk solvent effects are considered using the polarizable continuum method (PCM), where the solvent-solute interaction is described with a reaction field. Molecular orbitals and orbital energies in the presence of the reaction field corresponding to the excited state are employed and the response of the solvent is not included in the TDDFT calculations. A self-consistent procedure is designed to obtain the excited state reaction field. The equations based on this new nonequilibrium solvation model in the framework of the self-consistent SS-PCM/TDDFT method for calculation of vertical absorption and emission energies are presented and implemented in the Q-Chem package. Vertical absorption and emission energies for several small molecules in solution using the newly developed code are calculated and compared with available experimental data and the results of other theoretical studies. Solvent shifts of absorption and emission energies are reasonably reproduced with this approach. The new model is a promising approach to study nonequilibrium absorption and emission processes in solution.

Excitation-dependent fluorescence from atomic/molecular layer deposited sodium-uracil thin films

Atomic/molecular layer deposition (ALD/MLD) offers unique possibilities in the fabrication of inorganic-organic thin films with novel functionalities. Especially, incorporating nucleobases in the thin-film structures could open new avenues in the development of bio-electronic and photonic devices. Here we report an intense blue and widely excitation-dependent fluorescence in the visible region for ALD/MLD fabricated sodium-uracil thin films, where the crystalline network is formed from hydrogen-bonded uracil molecules linked via Na atoms. The excitation-dependent fluorescence is caused by the red-edge excitation shift (REES) effect taking place in the red-edge of the absorption spectrum, where the spectral relaxation occurs in continuous manner as demonstrated by the time-resolved measurements.

Characterization of the Photochemical Properties of 5-Benzyluracil via Time-Dependent Density Functional Theory

We present a detailed study of the excited state properties of 5-benzyluracil (5BU) in the gas phase and in implicit solvent using different electronic structure approaches ranging from time-dependent density functional theory in the linear response regime (LR-TDDFT) to a set of different wave-function-based methods for excited states, namely perturbed coupled cluster (CC2), algebraic diagrammatic construction method to second order (ADC(2)), and perturbed configuration interaction (CIS(D)). 5BU has been used to investigate DNA base-amino acid interactions. In particular, it served as a model of protein-DNA photoinduced cross-linking. While LR-TDDFT is computationally the most efficient first-principles approach for static and dynamic simulations of this bichromophoric system, its accuracy is difficult to assess due to the presence of excited states with charge transfer character. In this work, the performance of different exchange correlation functionals is compared against accurate benchmarks obtained either from high level wave-function-based methods or directly from experimental absorption spectra. Our investigation shows that accurate results for the excitation energies can be obtained using the hybrid meta-GGA functional M06. In view of dynamical studies of the relaxation of 5BU after photoexcitation, we also show that the PBE functional, while failing in the Franck-Condon region, provides qualitatively good results for the characterisation of a possible photocyclization path.

Excited-State Dynamics of the Thiopurine Prodrug 6-Thioguanine: Can N9-Glycosylation Affect Its Phototoxic Activity?

6-Thioguanine, an immunosuppressant and anticancer prodrug, has been shown to induce DNA damage and cell death following exposure to UVA radiation. Its metabolite, 6-thioguanosine, plays a major role in the prodrug's overall photoreactivity. However, 6-thioguanine itself has proven to be cytotoxic following UVA irradiation, warranting further investigation into its excited-state dynamics. In this contribution, the excited-state dynamics and photochemical properties of 6-thioguanine are studied in aqueous solution following UVA excitation at 345 nm in order to provide mechanistic insight regarding its photochemical reactivity and to scrutinize whether N9-glycosylation modulates its phototoxicity in solution. The experimental results are complemented with time-dependent density functional calculations that include solvent dielectric effects by means of a reaction-field solvation model. UVA excitation results in the initial population of the S₂(ππ*) state, which is followed by ultrafast internal conversion to the S₁(nπ*) state and then intersystem crossing to the triplet manifold within 560 ± 60 fs. A small fraction (ca. 25%) of the population that reaches the S₁(nπ*) state repopulates the ground state. The T₁(ππ*) state decays to the ground state in 1.4 ± 0.2 μs under N₂-purged conditions, using a 0.2 mM concentration of 6-thioguanine, or it can sensitize singlet oxygen in 0.21 ± 0.02 and 0.23 ± 0.02 yields in air- and O₂-saturated solution, respectively. This demonstrates the efficacy of 6-thioguanine to act as a Type II photosensitizer. N9-glycosylation increases the rate of intersystem crossing from the singlet to triplet manifold, as well as from the T₁(ππ*) state to the ground state, which lead to a ca. 40% decrease in the singlet oxygen yield under air-saturated conditions. Enhanced vibronic coupling between the singlet and triplet manifolds due to a higher density of vibrational states is proposed to be responsible for the observed increase in the rates of intersystem crossing in 6-thioguanine upon N9-glycosylation.

Multiconfigurational Effects in Theoretical Resonance Raman Spectra

We analyze resonance Raman spectra of the nucleobase uracil in the short‐time approximation calculated with multiconfigurational methods. We discuss the importance of static electron correlation by means of density‐matrix renormalization group self‐consistent field (DMRG‐SCF) calculations. Our DMRG‐SCF results reveal that a minimal active orbital space that leads to a qualitatively correct description of the resonance Raman spectrum of uracil should encompass parts of the σ/σ* bonding/anti‐bonding orbitals of the pyrimidine ring. We trace these findings back to the considerable entanglement between the σ/σ* bonding/anti‐bonding as well as valence π/π* orbitals in the excited‐state electronic structure of uracil, which indicates non‐negligible non‐dynamical correlation effects that are less pronounced in the electronic ground state.

Optical absorption and magnetic circular dichroism spectra of thiouracils: a quantum mechanical study in solution

The excited electronic states of thiouracils, the analogues of uracil where the carbonyl oxygens are substituted by sulphur atoms, have been investigated by computing the magnetic circular dichroism (MCD) and one-photon absorption (OPA) spectra at the TD-DFT level of theory.

Electrostatic Interaction-Induced Room-Temperature Phosphorescence in Pure Organic Molecules from QM/MM Calculations

Room temperature phosphorescence (RTP) from pure organic material is rare due to the low phosphorescence quantum efficiency. That is why the recent discovery of crystallization induced RTP for several organic molecules aroused strong interests. Through a combined quantum and molecular mechanics CASPT2/AMBER scheme taking terephthalic acid (TPA) as example, we found that electrostatic interaction not only can induce an enhanced radiative decay T1 → S0 through the dipole-allowed S1 intermediate state, but also can hinder the nonradiative decay process upon crystallization. From gas phase to crystal, the nature of S1 state is converted to (1)(π,π*) from (1)(n,π*) character, enhancing transition dipole moment and serving as an efficient intermediate radiative pathway for T1 → S0 transition, and eventually leading to a boosted RTP. The intermolecular packing also blocks the nonradiative decay channel of the high-frequency C═O stretching vibration with large vibronic coupling, rather than the conventional low-frequency aromatic rotation in crystal. This mechanism also holds for other organic compounds that contain both ketones and aromatic rings.

Toward a Theoretical Quantitative Estimation of the λmax of Anthraquinones-Based Dyes

We have computed the absorption spectra of a large series of anthraquinone dyes by using the time-dependent density functional theory (TD-DFT) for the excited-state calculations and the polarizable continuum model (PCM) for evaluating bulk solvent effects. On one hand, we compare the results obtained with the B3LYP and the PBE0 hybrid functionals, combined with different atomic basis sets. On the other hand, using multiple linear regression, we take advantage of the λmax predicted by these two functionals in order to reach the best agreement between theoretical estimates and experimental measurements. It turns out that 1. PBE0 provides more accurate results than B3LYP; in addition the average errors provided by the former are less basis set dependent. 2. Multiple linear regression provides excited state spectra in better agreement with experiment than any simple linear fit that could be performed. 3. Using our best fitting procedure, we obtained a mean absolute error of 6 nm for a set of 66 anthraquinones, with no deviations exceeding 25 nm. The related standard deviation, useful for predictions, is only 8 nm, i.e.,[Formula: see text] =[Formula: see text] ± 8 nm (or ±0.05 eV) for unknown anthraquinone compounds.

Polarizable Embedded RI-CC2 Method for Two-Photon Absorption Calculations

We present a novel polarizable embedded resolution-of-identity coupled cluster singles and approximate doubles (PERI-CC2) method for calculation of two-photon absorption (TPA) spectra of large molecular systems. The method was benchmarked for three types of systems: a water-solvated molecule of formamide, a uracil molecule in aqueous solution, and a set of mutants of the channelrhodopsin (ChR) protein. The first test case shows that the PERI-CC2 method is in excellent agreement with the PE-CC2 method and in good agreement with the PE-CCSD method. The uracil test case indicates that the effects of hydrogen bonding on the TPA of a chromophore with the nearest environment is well-described with the PERI-CC2 method. Finally, the ChR calculation shows that the PERI-CC2 method is well-suited and efficient for calculations on proteins with medium-sized chromophores.

Design of new aliphatic azido nitro compounds as plasticizer: an initial exploration on AFCTEE (1-azido-formic acid 1,1,1-trinitro-ethyl ester)

To explore new high-energy azido nitro compounds as plasticizers for propellants, AFCTEE (1-azido-formic acid 1,1,1-trinitro-ethyl ester) was designed and studied using density functional theory. The predicted density of AFCTEE, 1.90 gcm−3, is comparable to that of HMX (1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane) and much higher than that of general organic azido compounds. AFCTEE possesses higher energetic properties and chemical stability than the promising azido nitro plasticizer DAMNP (1,3-diazido-2-methyl-2-nitropropane) and the conventional plasticizer NG (nitroglycerine), and it has a moderate thermal stability. The pyrolysis of AFCTEE starts from the rupture of C–NO2 and then the breakage of N–N2 via Curtius rearrangement. This work is the initial exploration for AFCTEE, aiming at the energetics, spectra (IR, NMR, and UV), stability, and decomposition mechanism. Compared with DAMNP, the advantages of superior energetic properties and chemical stability suggest AFCTEE is a promising energetic azido nitro compound and is worth further investigation.

Mode-specific vibrational relaxation of photoexcited guanosine 5'-monophosphate and its acid form: a femtosecond broadband mid-IR transient absorption and theoretical study

UV-pump/broadband-mid-IR-probe transient absorption (TA) experiments and ab initio quantum mechanical (QM) calculations were used to investigate the photophysics in heavy water of the neutral and acid forms of guanosine 5'-monophosphate (GMP and GMPD(+), respectively). Excited GMP undergoes ultrafast internal conversion (IC) and returns to the electronic ground state in less than one picosecond with a large amount of excess vibrational energy. The spectroscopic signals are dominated by vibrational cooling - a process in which the solute dissipates vibrational energy to the solvent. For neutral GMP, cooling proceeds with a time constant of 3 ps. Following IC, at least some medium-frequency modes such as the carbonyl stretch and an in-plane ring vibration are excited, suggesting that the vibrational energy distribution is non-statistical. This is consistent with predicted structural changes upon passage through the S1/S0 conical intersection. GMPD(+) differs from GMP by a single deuteron at the N7 position, but has a dramatically longer lifetime of 200 ps. Vibrational cooling of the S1 state of GMPD(+) was monitored via several medium-frequency modes that were assigned using QM calculations. These medium-frequency modes are also vibrationally excited in a non-statistical fashion. Excitation of these modes is in line with the change in geometry at the S1 minimum of GMPD(+) predicted by QM calculations. Furthermore, these modes relax at different rates, fully consistent with QM calculations, which predict that excited vibrational states of the carbonyl stretch couple strongly to the D2O solvent and thus deactivate via intermolecular energy transfer (IET). In contrast, the ring stretch couples strongly to other ring modes of the guanine chromophore and appears to decay via intramolecular vibrational energy redistribution (IVR).

Theoretical study on excited states of bacteriochlorophyll a in solutions with density functional assessment

The excited-state properties of bacteriochlorophyll (BChl) a in triethylamine, 1-propanol, and methanol are investigated with the time-dependent density functional theory by using the quantum mechanical and molecular mechanical reweighting free energy self-consistant field method. It is found that no prevalent density functionals can reproduce the experimental excited-state properties, i.e., the absorption and reorganization energies, of BChl a in the solutions. The parameter μ in the range-separated hybrid functional is therefore optimized to reproduce the differences of the absorption energies in the solutions. We examine the origin of the differences of the absorption energies in the solutions and find that sensitive balance between contributions of structural changes and solute-solvent interactions determines the differences. The accurate description of the excitation with the density functional with the adjusted parameter is therefore essential to the understanding of the excited-state properties of BChl a in proteins and also the mechanism of the photosynthetic systems.

A theoretical study on the stability and intramolecular interaction in 5-nitrotetrazolates with the DFT and DFT-D methods

Two new salts 3,5-diazido-1,2,4-triazolium 5-nitrotetrazolate and 1-methyl-3,5-diazido-1,2,4-triazolium 5-nitrotetrazolate were designed based on the structures of experimentally synthesized 3-azido-1,2,4-triazolium 5-nitro-tetrazolate and 1-methyl-3-azido-1,2,4-triazolium 5-nitro-tetrazolate, to explore new promising candidates for energetic materials and to investigate the influences of the substituents (- CH 3 and - N 3) and solvent (water) on the intramolecular interactions and properties. The intramolecular hydrogen bonding interactions were investigated by the natural bond orbital (NBO) and the quantum theory of atoms in molecules (QTAIM) analyses using the density functional theory (DFT) and the dispersion correction DFT (DFT-D) methods. The low-lying singlet electronic transitions were estimated using the time-dependent DFT. All four examined salts exist as ionic structures in aqueous solution while acid–base molecular complexes form in gas phase. The hydrogen bond energy (E H ) obtained with the DFT-D method is larger than that obtained with the DFT method, but the trend is consistent, i.e. - N 3 increases while - CH 3 decreases E H . In addition, the position of the strongest electronic absorption peak has a little correlation with the number of - N 3 and - CH 3 groups. 3,5-diazido-1,2,4-triazolium 5-nitrotetrazolate is a valuable energetic salt with the highest nitrogen content, oxygen coefficient and density and the second highest heat of formation and chemical stability.