Search

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn–Sham method. The light–matter interaction is included in the Jaynes–Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light–matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined. [ABSTRACT FROM AUTHOR]

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy. [ABSTRACT FROM AUTHOR]

The optimization of the quantum efficiency of single-molecule light-driven rotary motors typically relies on chemical modifications. While, in isolated conditions, computational methods have been frequently used to design more efficient motors, the role played by the solvent environment has not been satisfactorily investigated. In this study, we used multiscale nonadiabatic molecular dynamics simulations of the working cycle of a 2-stroke photon-only molecular rotary motor. The results, which display dynamics consistent with the available transient spectroscopy measurements, predict a considerable decrease in the isomerisation quantum efficiency in methanol solution with respect to the gas phase. The origin of such a decrease is traced back to the ability of the motor to establish hydrogen bonds with solvent molecules. The analysis suggests that a modified motor with a reduced ability to form hydrogen bonds will display increased quantum efficiency, therefore extending the set of engineering rules available for designing light-driven rotary motors. Increasing the rotational efficiency of single-molecule light-driven rotary motors often relies on chemical modifications aimed at eliminating the factors that hinder rotation. Using multiscale nonadiabatic simulations, the authors investigate the transient conformations assumed by the motor molecule during its operation in a solvent and examine possibilities for enhancing the motor's efficiency by blocking certain solvent-solute interactions that restrain successful completion of the rotational movement. [ABSTRACT FROM AUTHOR]

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm −1) originates from the effect of electron donating and electron withdrawing substituents rather than π‐delocalization and related (anti‐)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S 1/S 0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π‐delocalization effects in π‐conjugated compounds. [ABSTRACT FROM AUTHOR]

The formulations and performances of spin‐restricted ensemble‐referenced KS and mixed‐reference spin‐flip time‐dependent density functional theory incorporating strong correlations into density functional theories are documented. As a result of balanced dynamic and nondynamic correlation, they are capable of describing strongly correlated challenging systems including diradicals, bond dissociation, conical intersections, doubly excited states, and so on, overcoming the limitations of current density functional theories. [ABSTRACT FROM AUTHOR]

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B u + and 2 1 A g - , of s- trans -butadiene and s- trans -hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){ N o }, and DEA-EOMCC(4p2h){ N u } or DEA-EOMCC(3p1h,4p2h){ N u }, consistently predict that the 1 1 B u + /2 1 A g - crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B u + and 2 1 A g - potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.