A variational reformulation of molecular properties in electronic-structure theory

Conventional quantum-mechanical calculations of molecular properties, such as dipole moments and electronic excitation energies, give errors that depend linearly on the error in the wave function. An exception is the electronic energy, whose error depends quadratically on the error in wave function. We here describe how all properties may be calculated with a quadratic error, by setting up a variational Lagrangian for the property of interest. Because the construction of the Lagrangian is less expensive than the calculation of the wave function, this approach substantially improves the accuracy of quantum-chemical calculations without increasing cost. As illustrated for excitation energies, this approach enables the accurate calculation of molecular properties for larger systems, with a short time-to-solution and in a manner well suited for modern computer architectures.

Dynamical Effects of Solvation on Norbornadiene/Quadricyclane Systems

Molecules that can undergo reversible chemical transformations following the absorption of light, the so-called molecular photoswitches, have attracted increasing attention in technologies, such as solar energy storage. Here, the optical and thermochemical properties of the photoswitch are central to its applicability, and these properties are influenced significantly by solvation. We investigate the effects of solvation on two norbornadiene/quadricyclane photoswitches. Emphasis is put on the energy difference between the two isomers and the optical absorption as these are central to the application of the systems in solar energy storage. Using a combined classical molecular dynamics and quantum mechanical/molecular mechanical computational scheme, we showcase that the dynamic effects of solvation are important. In particular, it is found that standard implicit solvation models generally underestimate the energy difference between the two isomers and overestimate the strength of the absorption, while the explicit solvation spectra are also less red-shifted than those obtained using implicit solvation models. We also find that the absorption spectra of the two systems are strongly correlated with specific dihedral angles. Altogether, this highlights the importance of including the dynamic effects of solvation.

Cluster perturbation theory IX: Perturbation series for the coupled cluster singles and doubles ground state energy

In this paper, we develop and analyze a number of perturbation series that target the coupled cluster singles and doubles (CCSD) ground state energy. We show how classical Møller-Plesset perturbation theory series can be restructured to target the CCSD energy based on a reference CCS calculation and how the corresponding cluster perturbation series differs from the classical Møller-Plesset perturbation series. Subsequently, we reformulate these series using the coupled cluster Lagrangian framework to obtain series, where fourth and fifth order energies are determined only using parameters through second order. To test the methods, we perform a series of test calculations on molecular photoswitches of both total energies and reaction energies. We find that the fifth order reaction energies are of CCSD quality and that they are of comparable accuracy to state-of-the-art approximations to the CCSD energy based on local pair natural orbitals. The advantage of the present approach over local correlation methods is the absence of user defined threshold parameters for neglecting or approximating contributions to the correlation energy. Fixed threshold parameters lead to discontinuous energy surfaces, although this effect is often small enough to be ignored, but the present approach has a differentiable energy that will facilitate derivation and implementation of gradients and higher derivatives. A further advantage is that the calculation of the perturbation correction is non-iterative and can, therefore, be calculated in parallel, leading to a short time-to-solution.

Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies

Even though cluster perturbation theory has been shown to be a robust noniterative alternative to coupled cluster theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. We present a proof-of-concept strategy for implementing a cluster perturbation theory ground-state energy series for the coupled cluster singles and doubles energy with N4 computational scaling using tensor hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected double intermediates to the THC format. Using the outlined strategy, we showcase that the THC pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations, and we empirically verify the N4 scaling.

A Combined Quantum Mechanics and Molecular Mechanics Approach for Simulating the Optical Properties of DNA-Stabilized Silver Nanoclusters

DNA-stabilized silver nanoclusters have emerged as an intriguing type of nanomaterial due to their unique optical and electronic properties, with potential applications in areas such as biosensing and imaging. The development of efficient methods for modeling these properties is paramount for furthering the understanding and utilization of these clusters. In this study, a hybrid quantum mechanical and molecular mechanical approach for modeling the optical properties of a DNA-templated silver nanocluster is evaluated. The influence of different parameters, including ligand fragmentation, damping, embedding potential, basis set, and density functional, is investigated. The results demonstrate that the most important parameter is the type of atomic properties used to represent the ligands, with isotropic dipole-dipole polarizabilities outperforming the rest. This underscores the importance of an appropriate representation of the ligands, particularly through the selection of the properties used to represent them. Moreover, the results are compared to experimental data, showing that the applied methodology is reliable and effective for the modeling of DNA-stabilized silver nanoclusters. These findings offer valuable insights that may guide future computational efforts to explore and harness the potential of these novel systems.

Extended Tetrathiafulvalenes with Fluoreno[3,2-b]fluorene and Diindeno[1,2-b : 1',2'-i]anthracene Cores

In one-dimensional polycyclic aromatic hydrocarbons (PAHs) containing five- and six-membered rings fused together, one key question is whether the structures possess a quinoidal or aromatic diradical character. Here, we generate such PAHs by reversible oxidation of PAH-extended tetrathiafulvalenes (TTFs). Extended TTFs were thus prepared and studied for their geometrical properties (crystallography), redox properties, and UV/Vis/NIR/EPR characteristics as a function of charge state. The EPR measurements of radical cations showed unique features for each PAH-TTF. The dications, formally composed of fluoreno[3,2-b]fluorene and diindeno[1,2-b:1',2'-i]anthracene cores, were experimentally found to exhibit singlet ground states. For the latter, calculations reveal the closed shell, quinoid singlet state to be isoenergetic with the open shell singlet diradical. Each charge state exhibited unique optical properties with radical cations absorbing strongly in the NIR region with signatures from π-dimers for the large core. The experimental results were paralleled and supported by detailed computations, including spin density distribution calculations, EPR simulations, and nucleus independent chemical shift (NICS) xy scans.

Investigation of Solvent Effects on the Molecular Energy Storage and Optical Properties of Bicyclooctadiene/Tetracyclooctane Photoswitches

In this paper, we investigate the effects of solvation on the solar energy storage properties of bicyclooctadiene/tetracyclooctane (BOD/TCO) photoswitches. The solvent effects on the thermochemical and optical properties are studied in cyclohexane, toluene, dichloromethane, ethanol, acetonitrile, and a vacuum using density functional theory and coupled cluster theory. Our results show that the energy storage capacity of the BOD/TCO system increases as the solvent polarity increases, and the change is more significant with an unsubstituted system. The energy storage capacity of the substituted system is not dependent on the polarity of the solvent. As the solvent polarity increases, the absorption peaks shift away from each other and the absorption intensities increase. Overall, the solvents improve the performance of the optical properties and the energy storage capacities of the BOD/TCO molecular solar thermal systems.

Solvent-Induced Symmetry Breaking of the Photoinduced Charge Transfer Dynamics in the Bridged Perylene Dimer

The charge transfer dynamics of the bridged perylene dimer were investigated with the recently developed solvent coupling Redfield time propagation model. The results are compared with previous experimental findings to showcase the significance and applicability of the model. The charge transfer dynamics in vacuum showed that no breaking of the charge transfer direction symmetry occurred upon optical excitation, in perfect agreement with the experiment. Meanwhile, attractive solute-solvent interactions facilitated by the dipole moments of the polar solvents were observed to break the charge transfer direction symmetry. The conformational isomerism effect on the transfer dynamics manifested itself by promoting different transport channels upon solvation. Consequently, the solvent coupling Redfield time propagation model was indeed found to be able to quantitatively describe the charge transfer dynamics including exotic phenomena such as symmetry breaking of charge transfer direction.

Heterocyclic [9]Helicenes Exhibiting Bright Circularly Polarized Luminescence

We describe a concise synthetic strategy for the preparation of heterocyclic [9]helicenes and a simple preparative-scale protocol for the optical resolution of the resulting M- and P-enantiomers. The helicenes were characterized by single-crystal X-ray diffraction along with a range of spectroscopic and computational techniques. A fluorescence quantum yield of up to 65 % was observed, and the chiroptical properties of both M- and P-helicenes revealed large dissymmetry factors. The circularly polarized luminescence brightness reaches up to 17 M-1  cm-1 , as measured experimentally and verified computationally, which makes this the highest circularly polarized luminescence brightness among heterocyclic helicenes. We describe how chiroptical properties (both circular dichroism and circularly polarized luminescence) can be described and predicted using quantum chemical calculations. The synthetic approach also reveals by-products that originate from internal oxidation reactions, presumably mediated by the close proximity of the π-surfaces in the helicene structure.

Searching the Chemical Space of Bicyclic Dienes for Molecular Solar Thermal Energy Storage Candidates

Photoswitches are molecular systems that are chemically transformed subsequent to interaction with light and they find potential application in many new technologies. The design and discovery of photoswitch candidates require intricate molecular engineering of a range of properties to optimize a candidate to a specific applications, a task which can be tackled efficiently using quantum chemical screening procedures. In this paper, we perform a large scale screening of approximately half a million bicyclic diene photoswitches in the context of molecular solar thermal energy storage using ab initio quantum chemical methods. We further device an efficient strategy for scoring the systems based on their predicted solar energy conversion efficiency and elucidate potential pitfalls of this approach. Our search through the chemical space of bicyclic dienes reveals systems with unprecedented solar energy conversion efficiencies and storage densities that show promising design guidelines for next generation molecular solar thermal energy storage systems.

The Effects of Solvent Dynamics on the Back Reaction of Solar-Thermal Energy Storage Systems

We have investigated dynamic solvent effects on molecular solar-thermal energy storage systems using models describing the effects of frequency dependent viscosities and dielectric constants on chemical reaction rates. We have utilized the generalized Langevin model for understanding how the reactions are affected by the frequency dependent viscosities and dielectric constants. Our results show that the rate constants of the molecular solar-thermal energy storage systems depend strongly on the dielectric electric solvent properties and the frequency dependent viscosities of the solvents.

Computational investigation of photoswitch conjugates for molecular solar energy storage

Solar energy conversion and storage are vital for combating climate change. Molecular solar thermal systems offer a promising solution, where energy is stored in molecular compounds. This study investigates dyad molecular photoswitches by combining bicyclooctadiene/tetracyclooctane and dihydroazulene/vinylheptafulvene systems with phenyl and cyano groups. Density functional theory calculations were employed to determine molecular properties and consider solvation effects in toluene and dichloromethane. The results show that the combined systems have a predicted storage energy of up to 206.14 kJ mol-1 and an absorption peak at 390.26 nm with appreciable intensity. These dyad photoswitches exhibit favorable properties for molecular solar thermal storage and other applications. A comparison with individual photoswitches reveals advantages and disadvantages. The most effective conjugate has a slightly lower storage density than an equal mixture of individual systems, but it demonstrates better absorption characteristics, with improved overlap with the solar spectrum and higher absorption intensity. These findings contribute to the understanding of dyad molecular photoswitches, showcasing their potential for advanced energy storage and conversion technologies.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Elucidating the Mystery of DNA-Templating Effects on a Silver Nanocluster

This presentation considers the effects that DNA-templating has on the optical properties of a 16-atom silver cluster. To accomplish this, hybrid quantum mechanical and molecular mechanical simulations of a Ag₁₆-DNA complex have been carried out and compared with pure time-dependent density functional theory calculations of two Ag₁₆ clusters in vacuum. The presented results show that the templating DNA polymers both red-shift the one-photon absorption of the silver cluster and increase its intensity. This occurs through a change in cluster shape prompted by the structural constraints of the DNA ligands combined with silver-DNA interactions. The overall charge of the cluster also contributes to the observed optical response, as oxidation of the cluster results in a simultaneous blue-shift of the one-photon absorption and a decrease in intensity. Additionally, the changes in shape and environment also lead to a blue-shift and enhancement of the two-photon absorption.

Coupled cluster theory on modern heterogeneous supercomputers

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

Calculation of electric field gradients in Cd(II) model complexes of the CueR protein metal site

With this work we first test various DFT functionals against CCSD(T) for calculation of EFGs at the position of Cd(II) in a very small model system, Cd(SCH₃)₂. Moreover, the available basis sets in ADF are tested in terms of basis set convergence, and the effect of including relativistic effects using the scalar relativistic and spin orbit ZORA Hamiltonians is explored. The results indicate that an error of up to around 10% on the calculated EFG may be expected using spin-orbit ZORA and the BHandHLYP functional with a locally dense basis set. Next, this method was applied to model systems of the CueR protein, aiming to interpret ¹¹¹Ag-PAC spectroscopic data. Note that ¹¹¹Ag decays to ¹¹¹Cd on which the PAC data are recorded. Surprisingly, model systems truncated - as is often done - at the first C-C bond from the central Cd(II) are inadequate in size, and larger model systems must be employed to achieve reliable EFG calculations. The calculated EFGs agree well with experimental PAC data, and indicate that shortly after the nuclear decay the structure relaxes from linear two-coordinate AgS₂ in the native protein, to a structure (or structures) where Cd(II) recruits additional ligands such as backbone carbonyl oxygens to achieve higher coordination number(s).

Massively Parallel GPU Enabled Third Order Cluster Perturbation Excitation Energies for Cost-Effective Large Scale Excitation Energy Calculations

We present here a massively parallel implementation of the recently developed CPS(D-3) <p>excitation energy model which is based on Cluster Perturbation Theory. The new algorithm</p> <p>extends the one developed in [P. Baudin, et al, J. Chem. Phys., 150.13, 134110 (2019)] to</p> <p>leverage multiple nodes and to utilize graphical processing units for acceleration of heavy</p> <p>tensor contractions. Furthermore, we show that the extended algorithm scales efficiently</p> <p>with increasing amounts of computational resources and that the developed code enables</p> <p>CPS(D-3) excitation energy calculations on large molecular systems with a low time-tosolution.</p> <p>More specifically, calculations on systems with over 100 atoms and 1000 basis</p> <p>functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation</p> <p>energies as a computationally efficient alternative to those obtained from the Coupled</p> <p>Cluster Singles and Doubles model.

Indirect nuclear spin-spin couplings with third order contributions added to the SOPPA method

In this article, a modification of the second order polarization propagator approximation (SOPPA) method is introduced and illustrated for the calculation of the indirect nuclear spin-spin couplings. The standard SOPPA method, although cheaper in terms of computational cost, offers less accurate results than the ones obtained with coupled cluster methods. A new method, named SOPPA+A3-3, was therefore developed by adding the terms of the third order A matrix that rely on the second order double amplitudes. The performance of this third order contribution was studied using the CCSD method as a reference, calculating the spin-spin couplings of molecules of diverse sizes and compositions, and comparing them to the SOPPA method. The results show that inclusion of this third order contribution gives more accurate results than the standard SOPPA method with a level of accuracy close to that of the coupled cluster method with only small increase of computational cost of the response calculation that dominates the computational cost for small to medium sized molecules. The implementation of the first contributions to the third order polarization propagator approximation (TOPPA) in the Dalton program thus already shows a significant change in these molecular properties over those obtained with the standard SOPPA method.

Imparting Stability to Organic Photovoltaic Components through Molecular Engineering: Mitigating Reactions with Singlet Oxygen

One key challenge in the development of viable organic photovoltaic devices is to design component molecules that do not degrade during combined exposure to oxygen and light. Such molecules should thus remain comparatively unreactive towards singlet molecular oxygen and not act as photosensitizers for the generation of this undesirable species. Here, novel redox-active chromophores that combine these two properties are presented. By functionalizing indenofluorene-extended tetrathiafulvalenes (IF-TTFs) with cyano groups at the indenofluorene core using Pd-catalyzed cyanation reactions, we find that the reactivity of the exocyclic fulvene carbon-carbon double bonds towards singlet oxygen is considerably reduced. The new cyano-functionalized IF-TTFs were tested in non-fullerene acceptor based organic photovoltaic proof-of-principle devices, revealing enhanced device stability.

A Benchmark Study of Aromaticity Indexes for Benzene, Pyridine, and the Diazines - II. Excited State Aromaticity

In this work, one geometrical aromaticity index and four electron sharing indexes are benchmarked for their application in excited state aromaticity calculations. Two computationally feasible and reliable procedures are identified, namely, CAM-B3LYP/cc-pVTZ and ωB97X-D/cc-pVTZ. Topological effects on the first excited singlet and triplet electronic manifold were investigated, and the latter was in general found to display more aromatic character compared to the S₁ surface. Besides, geometrical relaxation on each of the manifolds was observed to hamper the aromaticity, thereby resulting in more antiaromatic character. The relative order of excited state aromaticity within the studied molecules was noted to resemble the reversed version of the relative order of ground state aromaticity. Thereby, the following generalization was postulated: The more aromatic a molecule is in its ground state, the more antiaromatic it will be in its electronic first excited manifolds.

Redfield Propagation of Photoinduced Electron Transfer Reactions in Vacuum and Solution

Dynamical simulations of ultrafast electron transfer reactions are of utmost interest. To allow for energy dissipation directly into an external surrounding environment, a solvent coupling model has been deduced, implemented, and utilized to describe the photoinduced electron transfer dynamics within a model triad system herein. The model is based on Redfield theory, and the environment is represented by harmonic oscillators filled with bosonic quanta. To imitate real solvents, the oscillators have been equipped with frequencies and polarization lifetimes characteristic of the corresponding solvent. The population was found to transfer through the energetically lowest electron transfer route regardless of the medium. The condensed population transfer dynamics were observed to be highly dependent on the solvent parameters. In particular, an increase in the solvent coupling entailed a detainment in the population transfer from the initially prepared diabatic state and a promotion in the population transfer through the other electron transfer route. Two explanations based on the diagonal and off-diagonal matrix elements of the Kohn-Sham Fock matrix, respectively, have been provided. The lifetime of the populated partially charge-separated state was prolonged with increasing solvent polarity, and it was explained in terms of attractive interactions between the solvent's dipole moments and the fragments' charges. The high-frequency vibrational fine-structure in the correlation function was demonstrated to be important for the transfer dynamics, and the importance of dephasing effects in polar solvents was verified and precised to concern the optical polarization of the solvents.

Machine learning the frontier orbital energies of SubPc based triads

Organic photovoltaic devices are promising candidates for efficient energy harvesting from sunlight. Designing new dye molecules suitable for such devices is a challenging task restricted by the rapid increase of computational cost with system size. Solar cell material properties are closely related to the electronic structure of the dye, and an effective molecular orbital energy screening method for a family of dyes is therefore desired. In this work, a machine learning approach is used to sort through the chemical space of peripheral double-substituted boron-Subphthalocyanine dyes. A database of 12,102 PM6 optimized structures was built and for each of the structures time-dependent density functional theory (LC-[Formula: see text]HPBE/6-31+G(d)) calculations were performed. We investigated the changes of the molecular orbital energies of the molecular orbitals related to reduction and oxidation of the compounds. With the Electrotopological-state index moleculear representation all the tested algorithms, Support Vector Machine, Random Forest Regression, Neural Network, and Simple Linear Regression, captured the calculated frontier orbital energies with a prediction root-mean-square-error in the order of 0.05 eV. Finally, frontier orbital energies were predicted for more than 40,000 new structures by the trained Support Vector Machine algorithm. Compared to the parent boron-Subphthalocyanine structure, 237 and 132 functionalized dyes were predicted to have upshifted molecular orbital energies using the Electrotopological-state index and OneHot encoding feature vector, respectively. Out of 27 investigated donor and acceptor ligands, the acetamide and hydroxyl ligands gave rise to the desired increase in frontier molecular orbital energy.

Modeling Absorption and Emission Spectroscopies of Symmetric and Asymmetric Azaoxahelicenes in Vacuum and Solution

Helicenes are of general interest due to the significant chiral signals in both absorption- and emission-based spectroscopy. Herein, the spectroscopic properties of four recently synthesized azaoxahelicenes are studied using density functional theory methods. The azaoxahelicenes have 7, 9, 10, and 13 units and one to two complete turns of the structure. UV-vis absorption and electronic circular dichroism spectra are determined both in vacuum and in solution using explicit solvation through a combined molecular dynamics/polarizable embedding framework. Additionally, emission and circularly polarized luminescence spectra are determined based on vibronic calculations. The resulting spectra are in good agreement with the experimentally available data, highlighting that both absorption- and emission-based spectra of the systems can be modeled computationally such that reliable predictions can be made for systems that are yet to be synthesized.

Cluster perturbation theory. VIII. First order properties for a coupled cluster state

We have extended cluster perturbation (CP) theory to comprehend the calculation of first order properties (FOPs). We have determined CP FOP series where FOPs are determined as a first energy derivative and also where the FOPs are determined as a generalized expectation value of the external perturbation operator over the coupled cluster state and its biorthonormal multiplier state. For S(D) orbital excitation spaces, we find that the CP series for FOPs that are determined as a first derivative in general in second order have errors of a few per cent in the singles and doubles correlation contribution relative to the targeted coupled cluster (CC) results. For a SD(T) orbital excitation space, we find that the CP series for FOPs determined as a generalized expectation value in second order have errors of about ten percent in the triples correlation contribution relative to the targeted CC results. These second order models therefore constitute viable alternatives for determining high quality FOPs.

Electric Properties of Photochromic Molecules Physisorbed on Silver and Copper Nanoparticles

This paper investigates the electric properties of the photochromic dihydroazulene/vinylheptafulvene system as it is physisorbed onto silver and copper nanoparticles. Our focus is on how the polarizability and hyperpolarizability of the dihydroazulene, s-cis-vinylheptafulvene, and s-trans-vinylheptafulvene molecules depend on molecular orientation with respect to the nanoparticles, the molecule-cluster separation, and the type of nanoparticle. The computational approach utilizes a combined quantum mechanical/molecular mechanical method in which the molecules are treated quantum mechanically while the nanoparticles are treated with a simpler classical method. The molecules are described with density functional theory. The electric properties are calculated using response theory utilizing the long-range-corrected functional CAM-B3LYP and the correlation consistent basis set aug-cc-pVDZ. The atoms of the nanoparticles are represented using atomic polarizabilities. The interactions between the nanoparticles and the molecular systems are calculated using a polarizable embedding scheme after which the molecular properties are calculated with time-dependent density functional theory. The results show that the electric properties are indeed affected by the presence of the nanoparticles. It is also clear that it is the hyperpolarizabilities that change the most while the polarizabilities are less affected. Furthermore, the influence of the nanoparticles on the molecules depends heavily on the relative molecular orientation with respect to the nanoparticles and molecular conformation. Finally, it is observed that a copper nanoparticle has a larger influence on the molecular systems than a silver nanoparticle.

Perturbation of the UV transitions of formaldehyde by TiO2 photocatalysts and Aun nanoclusters

In the gas phase, formaldehyde has an electric-dipole forbidden transition that becomes allowed by vibronic coupling. In this paper we explore whether perturbation by surfaces could also enhance light absorption by CH₂O. We investigate the electronic transitions of formaldehyde in the gas phase and interacting with rutile (110) TiO₂, Au_n nanoclusters, and Au_n on (110)-TiO₂. These surfaces are chosen as being representative of metals and metal-oxide minerals, and also because of specific interest in photocatalysts and noble metal nanocluster catalysts. The oscillator strength of the forbidden n → π* transition of formaldehyde in vacuum is investigated by modelling vibrational coupling to the electronic transition with equation-of-motion coupled cluster theory. The excitation energies and oscillator strengths of formaldehyde are calculated for different orientations and distances to the surfaces using the coupled cluster singles and doubles linear response method within the Quantum Mechanical and Molecular Mechanical (QM/MM) model using the aug-cc-pVTZ basis set and compared with the values calculated in vacuo. The electronic transitions of formaldehyde vary very little when placed near a pure TiO₂-surface with only minor variations depending on the orientation of formaldehyde. Introducing a gold nanoparticle (by itself or supported by TiO₂) induces dramatic changes in the absorption properties. This is due to vibronic interactions and the effect of the broken symmetry on the n → π* transition. We see a large redshift in the transition of 90 nm and oscillator strengths larger than 1.0 × 10^-4 for CH₂O interacting with Au_n.

Prospects of Improving Molecular Solar Energy Storage of the Norbornadiene/Quadricyclane System through Bridgehead Modifications

We have investigated novel bicyclic diene molecular solar thermal energy storage systems that presently are the ones with the highest predicted energy density. Using a variety of different ab initio quantum chemical methods, we report storage energies, absorption spectra, and reaction barriers for the release of stored energy for a series of bicyclic dienes. The bicyclic dienes are all constructed by modifying the bridgehead of the well-known norbornadiene/quadricyclane (NBD/QC) system. In conclusion, we find it promising that it is possible to significantly amplify the storage energy of the NBD/QC system without seriously compromising other crucial properties by introducing simple modifications to the bridgehead.

Subphthalocyanine-triangulene dyads: Property tuning for light-harvesting device applications

Organic photovoltaics relies on the development of stable chromophores and redox-active organic molecules with tailor-made HOMO/LUMO energies. Here, we present the synthesis and properties of novel dyads composed of boron subphthalocyanine (SubPc) and triangulene units, connected either at the peripheral position of the subphthalocyanine or at the axial boron. The connectivity has strong implications for the absorption and fluorescence properties of the dyads, as well as their redox properties. While the SubPc unit has a bowl shape, triangulene is a planar structural unit that allows dyads to dimerize in the solid state on account of π-stacking interactions as shown by X-ray crystallography of one of the dyads. The electronic properties were also studied computationally by density functional theory methods. Excellent agreement between experimental and computed data were obtained, showing that our computational method is a strong tool in the rational design of optimum molecules to ultimately obtain finely tuned molecules for device applications.

Cluster Perturbation Theory. VII. The convergence of Cluster Perturbation Expansions

The convergence of the recently developed cluster perturbation CP expansions (Pawlowski et al, J. Chem. Phys. 150 134108(2019)) is analyzed with the double purpose of developing the mathematical tools and concepts needed to describe these expansions at general order and to identify the factors that define the rate of convergence of CP series. To this end, the CP energy, amplitude, and Lagrangian multiplier equations as functions of the perturbation strength are developed. By determining the critical points, defined as the perturbation strengths for which the Jacobian become singular, the rate of convergence as well as the intruder and critical states are determined for five simple molecules: BH, CO, H2O, NH3, and HF. To describe the patterns of convergence for these expansions at orders lower than the high-order asymptotic limit, a model is developed, where the perturbation corrections arise from two critical points. It is shown that this model allows rationalization of the behavior of the perturbation corrections at much lower order than required for the onset of the asymptotic convergence. For the H2O, CO, and HF molecules, the pattern and rate of convergence is defined by critical states where the Fock-operator underestimates the excitation energies, whereas the pattern and rate of convergence for BH is defined by critical states where the Fock-operator overestimates the excitation energy. For the NH3 molecule, both forms of critical points are required to describe the convergence behavior up to at least order 25.

Cluster perturbation theory. VI. Ground-state energy series using the Lagrangian

In this paper, we derive alternative cluster perturbation series to the energy cluster perturbation (ECP) series of Paper I. The ECP series were derived using the standard coupled cluster energy framework. Here, we use the coupled cluster Lagrangian framework to derive the Lagrangian cluster perturbation (LCP) series and show that a slightly modified order concept means that the ECP and the LCP series become identical. Using the Lagrangian, we also derive a perturbation series where total cluster amplitudes and multipliers are determined through the same orders as dictated by the 2n+1/2n+2 rule, the VCP series. The VCP energies have errors that are bilinear in the errors of the total cluster amplitudes and multipliers. Test calculations have been performed for S(D) and SD(T) orbital excitation spaces. The convergence of the test calculations can be divided into two groups. The first group contains molecules that exhibit a slow monotonic geometric convergence pattern in both orbital spaces. The second group, containing the rest of the molecules, exhibits a fast low order convergence. For the S(D) calculations, the deviations in fourth order from the coupled cluster singles and doubles energy are below 1.3 per cent for the LCP series and 0.3 per cent for the VCP series, respectively. For the SD(T) calculations, the second group hardly shows any difference between the low order convergence of the three series ECP, LCP, and LCP and at fourth order the deviations from the triples correlation energy are less than 1 per cent.

A Neural Network Approach for Property Determination of Molecular Solar Cell Candidates

The dihydroazulene/vinylheptafulvene (DHA/VHF) photocouple is a promising candidate for molecular solar heat batteries, storing and releasing energy in a closed cycle. Much work has been done on improving the energy storage capacity and the half-life of the high-energy isomer via substituent functionalization, but similarly important is keeping these improved properties in common polar solvents, along with being soluble in these, which is tied to the dipole properties. However, the number of possible derivatives makes an overview of this combinatorial space impossible both for experimental work and traditional computational chemistry. Due to the time-consuming nature of running many thousands of computations, we look to machine learning, which bears the advantage that once a model has been trained, it can be used to rapidly estimate approximate values for the given system. Applying a convolutional neural network, we show that it is possible to reach good agreement with traditional computations on a scale that allows us to rapidly screen tens of thousands of the DHA/VHF photocouple, eliminating bad candidates and allowing computational resources to be directed toward meaningful compounds.

Optimization of the thermochemical properties of the norbornadiene/quadricyclane photochromic couple for solar energy storage using nanoparticles

In this paper, we present an investigation concerning the prospects of using nanoparticles to improve solar energy storage properties of three different norbornadiene/quadricyclane derivatives. Computationally, we study how different nanoparticles influence the properties of the systems that relate to the storage of solar energy, namely, the storage energy and the back reaction barrier. Our approach employs hybrid quantum mechanical/molecular mechanical calculations in which the molecular systems are described using density functional theory while the nanoparticles are described using molecular mechanics. The interactions between the two subsystems are determined using polarization dynamics. The results show that the influence of the nanoparticles on the thermochemical properties largely depends on the type of nanoparticle used, the relative orientation with respect to the nanoparticle, and the distance between the the nanoparticle and the molecular system. Additionally, we find indications that copper and/or titanium dioxide nanoparticles can lower the energy barrier of the back reaction for all of the studied systems without significantly lowering the storage capability of the systems. Consequently, the study shows that nanoparticles can potentially be employed in the optimization of molecular photoswitches towards solar energy storage.

The effects of solvation on the back reaction and storage capabilities of solar thermal energy storage systems

Solvent effects on molecular solar thermal energy storage systems have been investigated using density functional theory combined with solvent models describing the effects of viscosities and dielectric constants on chemical reaction rates. We have addressed the following issues concerning how solvents influence both the thermochemical properties and the thermal relaxation kinetics of the studied systems, how the friction of the solvent influences the recrossing of the reactions along with the dynamics and force constants of the transition state. We observe that the rate constants for the chemical reactions of the molecular solar thermal energy storage systems depend strongly on the dielectric solvent properties and the viscosities of the solvents.

Synthesis, characterization and computational evaluation of bicyclooctadienes towards molecular solar thermal energy storage

Molecular solar-thermal energy storage (MOST) systems are based on photoswitches that reversibly convert solar energy into chemical energy. In this context, bicyclooctadienes (BODs) undergo a photoinduced transformation to the corresponding higher energy tetracyclooctanes (TCOs), but the photoswitch system has not until now been evaluated for MOST application, due to the short half-life of the TCO form and limited available synthetic methods. The BOD system degrades at higher temperature via a retro-Diels-Alder reaction, which complicates the synthesis of the compounds. We here report a cross-coupling reaction strategy that enables an efficient synthesis of a series of 4 new BOD compounds. We show that the BODs were able to switch to the corresponding tetracyclooctanes (TCOs) in a reversible way and can be cycled 645 times with only 0.01% degradation. Half-lives of the TCOs were measured, and we illustrate how the half-life could be engineered from seconds to minutes by molecular structure design. A density functional theory (DFT) based modelling framework was developed to access absorption spectra, thermal half-lives, and storage energies which were calculated to be 143-153 kJ mol-1 (0.47-0.51 MJ kg-1), up to 76% higher than for the corresponding norbornadiene. The combined computational and experimental findings provide a reliable way of designing future BOD/TCO systems with tailored properties.

A benchmark study of aromaticity indexes for benzene, pyridine and the diazines - I. Ground state aromaticity

Five different aromaticity indexes are benchmarked for benzene, pyridine and the diazines in their ground states. A basis set study was performed using the Pople style, Karlsruhe and Dunning's correlation consistent basis sets. Ten different DFT functionals, including LSDA, PBE, PBE0, B3LYP, CAM-B3LYP, wB97XD, M06-2X, SOGGA11X, M11 and MN15 were benchmarked by comparison with CCSD, CASSCF and MP2. Large out-of-plane imaginary frequencies were observed for some of the optimized structures at the correlated wavefunction level of theory. It was found that the DFT functionals in general predict the para-delocalization index, multicenter index and aromatic fluctuation index to be approximately 70%, 50% and 45% larger, respectively, compared to the CCSD method. Comparisons of the DFT functionals showed that the wB97XD, CAM-B3LYP and M06-2X functionals performed the best. Furthermore, the basis set dependence of the DFT functionals was found to be large for the electron sharing indexes. Based on these findings, it is recommended to perform ground state calculations of aromaticity indexes at the wB97XD, CAM-B3LYP or M06-2X level of theory utilizing a simple basis set of triple-ζ quality.

Five different aromaticity indexes are benchmarked for benzene, pyridine and the diazines in their ground states.

High Throughput Screening of Norbornadiene/Quadricyclane Derivates for Molecular Solar Thermal Energy Storage

We present a procedure for performing high throughput screening of molecular compounds for molecular solar thermal energy storage devices using extended tight binding (xTB) methods. In order to validate our...

Excited State Dynamics and Conjugation Effects of the Photoisomerization Reactions of Dihydroazulene

Herein, we present an investigation of the excited state dynamics of the Dihydroazulene photoswitch and its photoinduced reaction to Vinylheptafulvene. The focus is on how the introduction of a benzannulated...

Investigation of the Structural and Thermochemical Properties of [2.2.2]-Bicyclooctadiene Photoswitches

Molecular photoswitches can under certain conditions be used to store solar energy in the so-called molecular solar thermal storage systems, which is an interesting technology for renewable energy solutions. The current investigations focus on the performance of seven different density functional theory (DFT) methods (B3LYP, CAM-B3LYP, PBE0, M06-2X, ωB97X-D, B2PLYP, and PBE0DH) when predicting geometries and thermochemical properties of the [2.2.2]-bicyclooctadiene (BOD) photoswitch. We find that all of the investigated DFT methods provide geometries that are in good agreement with those obtained using coupled cluster singles and doubles (CCSD) calculations. The dependence on the employed basis set is not large when predicting geometries. With respect to the thermochemical properties, we find that the M06-2X, CAM-B3LYP, PBE0, and ωB97X-D functionals all predict thermochemical properties that are in good agreement with the results of the CCSD, the CCSD including perturbative triples (CCSD(T)), and the explicitly correlated CCSD-F12 and CCSD(T)-F12 models. Lastly, for energy calculations, we tested the newly developed fourth-order cluster perturbation theory singles and doubles CPS(D-4) model, which in this study provides energy differences that are of CCSD and sometimes also CCSD(T) quality at a relatively low cost. We find that the CPS(D-4) model is an excellent choice for further investigation of BOD derivatives because accurate energies can be obtained routinely using this methodology. From the results, we also note that the predicted storage energies and storage energy densities for the BOD photoswitch are very large compared to other molecular solar thermal storage systems and that these systems could be candidates for such applications.

Dynamics of nuclear recoil: QM-BOMD simulations of model systems following β-decay

The kinetic recoil energy received by the daughter nucleus in a nuclear decay is often large enough to affect the structure around the nucleus in chemical systems. The coinciding element change which typically occurs in a nuclear decay may additionally incur a structural reorganization. The effects of these phenomena on chemical systems where radio-isotopes are used are often little-known or neglected because the dynamics of nuclear decay is difficult to observe experimentally. In this work, QM-MD simulations are used to investigate local fs to ps dynamics following the β-decay of ¹¹¹Ag to ¹¹¹Cd in systems modelled on the metal-sensing CueR protein. An adiabatic approximation is applied, assuming that the electronic structure relaxes rapidly after the decay. PM7-MD simulations of recoil dynamics of the model systems show significant structural changes and bonding interactions that depend on the magnitude and direction of the recoil. We find that, in general, the kinetic recoil energy is rapidly distributed (<5 ps) uniformly throughout the systems in the studied scenarios.

Virtual screening of norbornadiene-based molecular solar thermal energy storage systems using a genetic algorithm

We present a computational methodology for the screening of a chemical space of 10²⁵ substituted norbornadiene molecules for promising kinetically stable molecular solar thermal (MOST) energy storage systems with high energy densities that absorb in the visible part of the solar spectrum. We use semiempirical tight-binding methods to construct a dataset of nearly 34 000 molecules and train graph convolutional networks to predict energy densities, kinetic stability, and absorption spectra and then use the models together with a genetic algorithm to search the chemical space for promising MOST energy storage systems. We identify 15 kinetically stable molecules, five of which have energy densities greater than 0.45 MJ/kg, and the main conclusion of this study is that the largest energy density that can be obtained for a single norbornadiene moiety with the substituents considered here, while maintaining a long half-life and absorption in the visible spectrum, is around 0.55 MJ/kg.

Dihydroazulene-Azobenzene-Dihydroazulene Triad Photoswitches

Photoswitch triads comprising two dihydroazulene (DHA) units in conjugation with a central trans-azobenzene (AZB) unit were prepared in stepwise protocols starting from meta- and para-disubstituted azobenzenes. The para-connected triad had significantly altered optical properties and lacked the photoactivity of the separate photochromes. In contrast, for the meta-connected triad, all three photochromes could be photoisomerized to generate an isomer with two vinylheptafulvene (VHF) units and a cis-azobenzene unit. Ultrafast spectroscopy of the photoisomerizations revealed a fast DHA-to-VHF photoisomerization and a slower trans-to-cis AZB photoisomerization. This meta triad underwent thermal VHF-to-DHA back-conversion with a similar rate of all VHFs, independent of the identity of the neighboring units, and in parallel thermal cis-to-trans AZB conversion. The experimental observations were supported by computation (excitation spectra and orbital analysis of the transitions).

Fulvalene-Based Polycyclic Aromatic Hydrocarbon Ladder-Type Structures: Synthesis and Properties

Polycyclic aromatic hydrocarbons (PAHs) have found strong interest for their electronic properties and as model systems for graphene. While PAHs have been studied intensively as single units, here PAHs were constructed in ladder-type arrangements using cross-conjugated fulvalene and dithiafulvalene motifs as connecting units and moving forward a convenient synthetic approach for dimerizing (thio)ketones into olefins by the action of Lawesson's reagent. Some of the PAHs can also be regarded as "super-extended" tetrathiafulvalenes (TTFs) with some of the largest cores ever explored, being multi-redox systems that exhibit both reversible oxidations and reductions. Concomitant absorption redshifts were observed when expanding the ladder-type structures from one to two to three indenofluorene units, and optical and electrochemical HOMO-LUMO gaps were found to correlate linearly. Various conformations (and solid-state packing arrangements) were studied by X-ray crystallography and computations.

Promoting the thermal back reaction of vinylheptafulvene to dihydroazulene by physisorbtion on nanoparticles

We investigate the effects of nanoparticles on molecular solar thermal energy storage systems and how one can tune chemical reactivities of a molecular photo- and thermoswitch by changing the nanoparticles. We have selected the dihydroazulene/vinylheptafulvene system to illustrate the effects of the nanoparticles on the chemical reactivities of the molecular photo- and thermoswitch. We have utilized the following nanoparticles: a TiO₂ nanoparticle along with nanoparticles of gold, silver and copper. We calculate the rate constants for the release of the thermal energy utilizing a QM/MM method coupled to a transition state method. The molecular systems are described by density functional theory whereas the nanoparticles are given by molecular mechanics including electrostatic and polarization dynamics. In order to investigate whether the significant stabilization of the transitions state provided by the nanoparticles is general to the DHA/VHF system, we calculated the transition state rate constant of the parent- and 3-amino-substituted-DHA/VHF systems at 298.15 K in the four different orientations and at the three different separations. We observe that the transition state rate constant of the parent system is only increased as the cyano groups are oriented towards the nanoparticle while the presence of the nanoparticle actually impedes the reactions using the three other orientations. On the other hand, for the substituted system the nanoparticle generally leads to a significant increase in the rate of the reaction. We find that the nanoparticles can have a substantial effect on the calculated rate constants. We observe, depending on the nanoparticle and the molecular orientation, increases of the rate constants by a factor of 10⁶. This illustrates the prospects of utilizing nanoparticles for controlling the release of the stored thermal energy.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.

Simulating fullerene polyhedral formation from planar precursors

The synthesis path of the C₆₀-Buckyball fullerene from a planar precursor developed by Scott et al. [Science, 2002, 295, 5559] is investigated with density functional theory (DFT) methods. Various theoretically possible closing paths are analysed with respect to structural and energetic properties. The initial geometries were obtained by geometric interpolation of a cardboard-like model comprising rigid rings connected by hinges, which were then fully optimized with a selection of DFT-functionals. Analysis of the fully optimised geometries shows remarkable stability of face planarity, bond lengths and bond angles for all studied geometries, indicating soundness of the "cardboard with hinges"-model for approximating reaction paths for molecules of this type. This raises hope for development of a force field description of fullerene precursor molecules that can aid in discovery and analysis of good precursor candidates for rational synthesis of new fullerenes.

High throughput virtual screening of 230 billion molecular solar heat battery candidates

The dihydroazulene/vinylheptafulvene (DHA/VHF) thermocouple is a promising candidate for thermal heat batteries that absorb and store solar energy as chemical energy without the need for insulation. However, in order to be viable the energy storage capacity and lifetime of the high energy form (i.e., the free energy barrier to the back reaction) of the canonical parent compound must be increased significantly to be of practical use. We use semiempirical quantum chemical methods, machine learning, and density functional theory to virtually screen over 230 billion substituted DHA molecules to identify promising candidates. We identify a molecule with a predicted energy density of 0.38 kJ/g, which is significantly larger than the 0.14 kJ/g computed for the parent compound. The free energy barrier to the back reaction is 11 kJ/mol higher than the parent compound, which should correspond to a half-life of about 10 days—4 months. This is considerably longer than the 3–39 h (depending on solvent) observed for the parent compound and sufficiently long for many practical applications. Our paper makes two main important contributions: (1) a novel and generally applicable methodological approach that makes screening of huge libraries for properties involving chemical reactivity with modest computational resources, and (2) a clear demonstration that the storage capacity of the DHA/VHF thermocouple cannot be increased to >0.5 kJ/g by combining simple substituents.