The use of Coulomb-attenuated methods for the calculation of electronic circular dichroism spectra

Origin invariant electronic circular dichroism in the length dipole gauge without London atomic orbitals

We present a method for obtaining origin-independent electronic circular dichroism (ECD) in the length-gauge representation LG(OI) without the usage of London atomic orbitals. This approach builds upon the work by Caricato [J. Chem. Phys. 153, 151101 (2020)] and is applied to rotatory strengths and ECD spectra from damped response theory. Numerical results are presented for time-dependent Hartree-Fock and density-functional theory, the second-order algebraic diagrammatic construction method, and linear-response coupled-cluster theory with singles and approximate doubles. We can support the finding that the common choice of placing the gauge origin in the center of mass of a molecule in conventional length-gauge calculations involving chiroptical properties might not be optimal and show that LG(OI) is a valuable alternative for the origin-independent calculation of ECD spectra. We show that, for a limited test set, the convergence of the rotatory strengths calculated with the LG(OI) approach toward the basis-set limit tends to be faster than for the established velocity gauge representation. Relationships between the sum-over-states expression of the optical rotation in the LG(OI) framework and its representation in terms of response functions are analyzed.

TD-DFT calculations of one- and two-photon absorption in Coumarin C153 and Prodan: attuning theory to experiment

We use TD-DFT to calculate the one-photon absorption (1PA) and two-photon absorption (2PA) properties of C153 and Prodan in toluene and DMSO, and benchmark different methods relative to accurate experimental data available from the literature on these particular systems. As the first step, we modify the range-separated TD-DFT to provide the best prediction for the peak 1PA wavelength, and then apply the optimized functionals to achieve quantitative predictions of the corresponding two-photon absorption cross section, σ_2PA, with an accuracy ∼10-20% in C153 and ∼20-30% in Prodan. To elucidate the origin of residual discrepancies between the theory and experimental observations, we invoked the two essential states model for σ_2PA, which allows us to verify not only the transition wavelength and the σ_2PA value, but also to quantitatively benchmark the calculation of key molecular parameters such as the transition dipole moment and the change of the permanent dipole moment. Such comprehensive cross-checking indicates that a larger discrepancy in Prodan is most likely caused by a noted failure of DFT to predict the relative intensity and relative ordering of closely lying excited states with different degrees of intramolecular charge transfer, which we further support by analyzing the predictions obtained by high-level coupled-cluster calculations in the gas phase. Our results highlight the utility of benchmarking the calculations not only relative to other theoretical methods, but also in comparison to the experimental measurements, wherever such data are available.

Systematic study on the TD-DFT calculated electronic circular dichroism spectra of chiral aromatic nitro compounds: A comparison of B3LYP and CAM-B3LYP

B3LYP is one of the most widely used functional for the prediction of electronic circular dichroism spectra, however if the studied molecule contains aromatic nitro group computations may fail to produce reliable results. A test set of molecules of known stereochemistry were synthesized to study this phenomenon in detail. Spectra were computed by B3LYP and CAM-B3LYP functionals with 6-311++G(2d,2p) basis set. It was found that the range separated CAM-B3LYP gives better predictions than B3LYP for all test molecules. Fragment population analysis revealed that the nitro groups form highly localized molecule orbitals but the exact composition depends on the functional. CAM-B3LYP allows sufficient spatial overlap between the nitro group and distant parts of the molecule, which is necessary for the accurate description of excited states especially for charge transfer states. This phenomenon and the synthesized test molecules can be used to benchmark theoretical methods as well as to help the development of new functionals intended for spectroscopical studies.

Gas Phase Detection of the NH-P Hydrogen Bond and Importance of Secondary Interactions

We have observed the NH···P hydrogen bond in a gas phase complex. The bond is identified in the dimethylamine-trimethylphosphine complex by a red shift of the fundamental NH-stretching frequency observed using Fourier transform infrared spectroscopy (FT-IR). On the basis of the measured NH-stretching frequency red shifts, we find that P is a hydrogen bond acceptor atom similar in strength to S. Both are stronger acceptors than O and significantly weaker acceptors than N. The hydrogen bond angle, ∠NHP, is found to be very sensitive to the functional employed in density functional theory (DFT) optimizations of the complex and is a possible parameter to assess the quality of DFT functionals. Natural bonding orbital (NBO) energies and results from the topological methods atoms in molecules (AIM) and noncovalent interactions (NCI) indicate that the sensitivity is caused by the weakness of the hydrogen bond compared to secondary interactions. We find that B3LYP favors the hydrogen bond and M06-2X favors the secondary interactions leading to under- and overestimation, respectively, of the hydrogen bond angle relative to a DF-LCCSD(T)-F12a calculated angle. The remaining functionals tested, B3LYP-D3, B3LYP-D3BJ, CAM-B3LYP, and ωB97X-D, as well as MP2, show comparable contributions from the hydrogen bond and the secondary interactions and are close to DF-LCCSD(T)-F12a results.

Electronic circular dichroism of chiral alkenes: B3LYP and CAM-B3LYP calculations

Time-dependent density functional theory (TDDFT) calculations of electronic circular dichroism (ECD) are widely used to determine absolute configurations (ACs) of chiral molecules. Two very popular DFT exchange-correlation functionals, one hybrid (B3LYP) and one long-range corrected (CAM-B3LYP), along with a hierarchical sequence of basis sets were investigated, and the ECD spectra predicted for eight alkenes and compared to gas-phase experimental spectra. Little variation in predicted ECD spectra was found with the basis set size enlargement, but the sensitivity to the functional is greater. Good agreement was obtained only with the CAM-B3LYP functional, leading to the conclusion that TDDFT calculations of ECD spectra can routinely provide reliable ACs if and only if an appropriate functional is used. For camphene, twistene, syn-(E)-bisfenchylidene, and phyllocladene, solvent effects were estimated.

On the origin of the very strong two-photon activity of squaraine dyes – a standard/damped response theory study

In the present work, we report the mechanism of a very large increase in the two-photon (TP) activity of squaraine based molecules upon changing the substituents.

Time-dependent density functional response theory for electronic chiroptical properties of chiral molecules

Methodology to calculate electronic chiroptical properties from time-dependent density functional theory (TDDFT) is outlined. Applications of TDDFT to computations of electronic circular dichroism, optical rotation, and optical rotatory dispersion are reviewed. Emphasis is put on publications from 2005 to 2010, but much of the older literature is also cited and discussed. The determination of the absolute configuration of chiral molecules by combined measurements and computations is an important application of TDDFT chiroptical methods and discussed in some detail. Raman optical activity (ROA) spectra are obtained from normal-mode derivatives of the optical rotation tensor and other linear response tensors. A few selected (ROA) benchmarks are reviewed.

The 559-to-600 nm shift observed in red fluorescent protein eqFP611 is attributed to cis-trans isomerization of the chromophore in an anionic protein pocket

Fluorescent proteins are commonly used as molecular labels, noninvasive markers of gene expression, and reporters of environmental conditions in live cells. We investigate the structural and spectroscopic properties of the chromophore of a far-red fluorescent protein eqFP611. Both the cis and trans isomers of the chromophore are examined within the protein for which both anionic and neutral states of protonation are considered. Spectroscopic properties are examined using time-dependent density functional theory (TDDFT), employing the B3LYP, PBE and B3PW91 density functionals. Intermolecular and long-range contributions to the structure and spectroscopy were treated using the own n-layered integrated molecular orbital and molecular mechanics (ONIOM) approach. The results indicated that the chromophore before excitation is in an anionic, protonated state, with the long-range contributions inducing a blue shift in the absorption and fluorescence maxima of the chromophore. Moreover, the calculated changes of the lowest pi-pi* excitation energy upon isomerization match the observed shift from 559 to 600 nm in the absorption maximum of the system following prolonged irradiation. Furthermore, decomposition analysis of the electrostatic contributions from individual residues indicated that the interactions from four residues Arg92, Lys67, Glu145, and His197 to the chromophore play a key role in the absorption and fluorescence spectra of eqFP611, suggesting that mutations at these sites should provide very useful mechanistic information.