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]

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems. [ABSTRACT FROM AUTHOR]

It is estimated that 72% of the worldwide primary energy consumption is lost as waste heat. Thermoelectric Generators (TEGs) are a possible solution to convert a part of this energy into electricity and heat for space heating. However, for their deployment a proven long-term operation is required. Therefore, this research investigates the long-term stability of TEGs on system level in air and argon atmosphere under thermal cycling up to 543 K. The layout of the examined test objects resembles a TEG in stack design. The results show that the maximal output power of the test object in air reaches a plateau at 57% of the initial power after 50 cycles caused by an increased electrical resistance of the system. Whereas the test object in argon atmosphere shows no significant degradation of electrical power or resistance. The findings represent a step towards the understanding of the long-term stability of TEGs and can be used as a guideline for design decisions. [ABSTRACT FROM AUTHOR]

Molecular emitters simultaneously generating light at different wavelengths have wide applications. With a small molecule, however, it is challenging to realize two independent radiative pathways. We invented the first examples of dual‐emissive single‐benzene fluorophores (SBFs). Two emissive tautomers are generated by synthetic modulation of the hydrogen bond acidity, which opens up pathways for excited‐state proton transfer. White light is produced by a delicate balance between the energy and intensity of the emission from each tautomer. We show that the excited‐state antiaromaticity of the benzene core itself dictates the proton movements driving the tautomer equilibrium. Using this simple benzene platform, a fluorinated SBF was synthesized with a record high solubility in perfluorocarbon solvents. White light‐emitting devices and multicolor imaging of perfluorocarbon nanodroplets in live cells demonstrate the practical utility of these molecules. [ABSTRACT FROM AUTHOR]

The use of the mixed reference (MR) reduced density matrix, which combines reduced density matrices of the MS = +1 and -1 triplet-ground states, is proposed in the context of the collinear spin-fliptime-dependent density functional theory (SF-TDDFT) methodology. The time-dependent Kohn-Sham equation with the mixed state is solved by the use of spinor-like open-shell orbitals within the linear response formalism, which enables to generate additional configurations in the realm of TD-DFT. The resulting MR-SF-TDDFT computational scheme has several advantages before the conventional collinear SF-TDDFT. The spin-contamination of the response states of SF-TDDFT is nearly removed. This considerably simplifies the identification of the excited states, especially in the "black-box" type applications, such as the automatic geometry optimization, reaction path following, or molecular dynamics simulations. With the new methodology, the accuracy of the description of the excited states is improved as compared to the collinear SF-TDDFT. Several test examples, which include systems typified by strong non-dynamic correlation, orbital (near) degeneracy, and conical intersections, are given to illustrate the performance of the new method. [ABSTRACT FROM AUTHOR]

The rational engineering of photoresponsive materials, e.g., light-driven molecular motors, is a challenging task. Here, we use structure-related design rules to prepare a prototype molecular rotary motor capable of completing an entire revolution using, exclusively, the sequential absorption of two photons; i.e., a photon-only two-stroke motor. The mechanism of rotation is then characterised using a combination of non-adiabatic dynamics simulations and transient absorption spectroscopy measurements. The results show that the rotor moiety rotates axially relative to the stator and produces, within a few picoseconds at ambient T, an intermediate with the same helicity as the starting structure. We discuss how such properties, that include a 0.25 quantum efficiency, can help overcome the operational limitations of the classical overcrowded alkene designs. Improving the efficiency of light-driven molecular rotary motors is a challenging task. Here, the authors combine theoretical modeling, synthesis and spectroscopy to prepare a prototype molecular motor capable of avoiding inefficient thermally activated motion; thus offering prospects to implement a 2-stroke photon-only molecular motor. [ABSTRACT FROM AUTHOR]

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the newmethod correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena. [ABSTRACT FROM AUTHOR]

The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham(REKS) method and its state interaction (SI) extension, SI-SA-REKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. In this work, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. The implementation and benchmarking of the derived formalism will be described in a subsequent article of this series. [ABSTRACT FROM AUTHOR]

The spin-restricted ensemble-referenced Kohn-Sham (REKS) method is based on an ensemble representation of the density and is capable of correctly describing the non-dynamic electron correlation stemming from (near-)degeneracy of several electronic configurations. The existing REKS methodology describes systems with two electrons in two fractionally occupied orbitals. In thiswork, the REKS methodology is extended to treat systems with four fractionally occupied orbitals accommodating four electrons and self-consistent implementation of the REKS(4,4) method with simultaneous optimization of the orbitals and their fractional occupation numbers is reported. The new method is applied to a number of molecular systems where simultaneous dissociation of several chemical bonds takes place, as well as to the singlet ground states of organic tetraradicals 2,4-didehydrometaxylylene and 1,4,6,9-spiro[4.4]nonatetrayl. [ABSTRACT FROM AUTHOR]

A Kohn--Sham-type computational scheme capable of treating systems with strong nondynamic correlation is examined. The experiment results of REKS calculation are compared with the available Kohn--Sham solutions for cases in which the exact density is known, as well as with results of conventional multireference ab initio methods.

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]

The analytical gradient for the two-component Normalized Elimination of the Small Component (2c-NESC) method is presented. The 2c-NESC is a Dirac-exact method that employs the exact two-component one-electron Hamiltonian and thus leads to exact Dirac spin-orbit (SO) splittings for one-electron atoms. For many-electron atoms and molecules, the effect of the two-electron SO interaction is modeled by a screened nucleus potential using effective nuclear charges as proposed by Boettger [Phys. Rev. B 62, 7809 (2000)]. The effect of spin-orbit coupling (SOC) on molecular geometries is analyzed utilizing the properties of the frontier orbitals and calculated SO couplings. It is shown that bond lengths can either be lengthened or shortened under the impact of SOC where in the first case the influence of low lying excited states with occupied antibonding orbitals plays a role and in the second case the j j-coupling between occupied antibonding and unoccupied bonding orbitals dominates. In general, the effect of SOC on bond lengths is relatively small (=5% of the scalar relativistic changes in the bond length). However, large effects are found for van der Waals complexes Hg2 and Cn2, which are due to the admixture of more bonding character to the highest occupied spinors. [ABSTRACT FROM AUTHOR]

State-averaged (SA) variants of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, SA-REKS and state-interaction (SI)-SA-REKS, implement ensemble density functional theory for variationally obtaining excitation energies of molecular systems. In this work, the currently existing version of the SA-REKS method, which included only one excited state into the ensemble averaging, is extended by adding more excited states to the averaged energy functional. A general strategy for extension of the REKS-type methods to larger ensembles of ground and excited states is outlined and implemented in extended versions of the SA-REKS and SI-SA-REKS methods. The newly developed methods are tested in the calculation of several excited states of ground-state multi-reference systems, such as dissociating hydrogen molecule, and excited states of donor-acceptor molecular systems. For hydrogen molecule, the new method correctly reproduces the distance dependence of the lowest excited state energies and describes an avoided crossing between the doubly excited and singly excited states. For bithiophene-perylenediimide stacked complex, the SI-SA-REKS method correctly describes crossing between the locally excited state and the charge transfer excited state and yields vertical excitation energies in good agreement with the ab initio wavefunction methods. [ABSTRACT FROM AUTHOR]

Electron transfer in the ground and excited states of a model donor-acceptor (D-A) system is investigated using the single-reference and multi-reference density functional theory (DFT) methods. To analyze the results of the calculations, a simple two-site multi-reference model was derived that predicts a stepwise electron transfer in the S0 state and a wave-like dependence of the S1 electron transfer on the external stimulus. The standard single-reference Kohn-Sham (KS) DFT approach and the time-dependent DFT (TDDFT) method failed to describe the correct dependence of the S0 and S1 electron transfer on the external electric field applied along the donor-acceptor system. The multi-reference DFT approach, the spin-restricted ensemble-referenced KS (REKS) method, was able to successfully reproduce the correct behavior of the S0 and S1 electron transfer on the applied field. The REKS method was benchmarked against experimentally measured gas phase charge transfer excitations in a series of organic donor-acceptor complexes and displayed its ability to describe this type of electronic transitions with a very high accuracy, mean absolute error of 0.05 eV with the use of the standard range separated density functionals. On the basis of the calculations undertaken in this work, it is suggested that the non-adiabatic coupling between the S0 and S1 states may interfere with the electron transfer in a weakly coupled donor-acceptor system. It is also suggested that the electronic excitation of a D+-A- system may play a dual role by assisting the further electron transfer at certain magnitudes of the applied electric field and causing the backward transfer at lower electric field strengths. [ABSTRACT FROM AUTHOR]

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 À or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems. [ABSTRACT FROM AUTHOR]

Computational investigation of the longest wavelength excitations in a series of cyanines and linear n-acenes is undertaken with the use of standard spin-conserving linear response time-dependent density functional theory (TD-DFT) as well as its spin-flip variant and a ΔSCF method based on the ensemble DFT. The spin-conserving linear response TD-DFT fails to accurately reproduce the lowest excitation energy in these π-conjugated systems by strongly overestimating the excitation energies of cyanines and underestimating the excitation energies of n-acenes. The spin-flip TD-DFT is capable of correcting the underestimation of excitation energies of n-acenes by bringing in the non-dynamic electron correlation into the ground state; however, it does not fully correct for the overestimation of the excitation energies of cyanines, for which the non-dynamic correlation does not seem to play a role. The ensemble DFT method employed in this work is capable of correcting for the effect of missing non-dynamic correlation in the ground state of n-acenes and for the deficient description of differential correlation effects between the ground and excited states of cyanines and yields the excitation energies of both types of extended π-conjugated systems with the accuracy matching high-level ab initio multireference calculations. [ABSTRACT FROM AUTHOR]

The analytic derivative algorithm for the effective contact densities obtained by the exact two-component (X2C) relativistic Hamiltonian is extended to the local approximations to X2C to achieve a higher computational efficiency without losing accuracy. The new algorithm has been implemented in a standalone program, which can utilize the molecular orbitals from state-of-the-art ab initio or density functional theory (DFT) calculations by other quantum chemistry programs in connection with various relativistic Hamiltonians. With the help of the utility program, the effective contact densities as well as the related Mössbauer isomer shifts can be studied by various advanced single-reference and multi-reference ab initio methods as long as the canonical or natural orbitals are available. Using the developed algorithm, the effective contact densities and the Mössbauer isomer shifts in a series of iron compounds and in HgFn (n = 1, 2, 4, and 6) molecules were studied. The obtained results show that (1) adequate account of the static electron correlation significantly improves the agreement of the theoretical 57Fe effective contact densities with the experimental isomer shifts, and (2) the non-monotonous changes of the effective contact density in a series of HgFn compounds are caused by the increasing screening effect due to shrinking of the Hg 5d orbitals. [ABSTRACT FROM AUTHOR]