1
|
A variational reformulation of molecular properties in electronic-structure theory. SCIENCE ADVANCES 2024; 10:eadn3454. [PMID: 38657075 PMCID: PMC11042728 DOI: 10.1126/sciadv.adn3454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Accepted: 03/19/2024] [Indexed: 04/26/2024]
Abstract
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.
Collapse
|
2
|
Pyrene Functionalized Norbornadiene-Quadricyclane Fluorescent Photoswitches: Characterization of their Spectral Properties and Application in Imaging of Amyloid Beta Plaques. Chemistry 2024:e202400322. [PMID: 38629212 DOI: 10.1002/chem.202400322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Indexed: 05/23/2024]
Abstract
This study presents the synthesis and characterization of two fluorescent norbornadiene (NBD) photoswitches, each incorporating two conjugated pyrene units. Expanding on the limited repertoire of reported photoswitchable fluorescent NBDs, we explore their properties with a focus on applications in bioimaging of amyloid beta (Aβ) plaques. While the fluorescence emission of the NBD decreases upon photoisomerization, aligning with what has been previously reported, for the first time we observed luminescence after irradiation of the quadricyclane (QC) isomer. We deduce how the observed emission is induced by photoisomerization to the excited state of the parent isomer (NBD) which is then the emitting species. Thorough characterizations including NMR, UV-Vis, fluorescence, X-ray structural analysis and density functional theory (DFT) calculations provide a comprehensive understanding of these systems. Notably, one NBD-QC system exhibits exceptional durability. Additionally, these molecules serve as effective fluorescent stains targeting Aβ plaques in situ, with observed NBD/QC switching within the plaques. Molecular docking simulations explore NBD interactions with amyloid, unveiling novel binding modes. These insights mark a crucial advancement in the comprehension and design of future photochromic NBDs for bioimaging applications and beyond, emphasizing their potential in studying and addressing protein aggregates.
Collapse
|
3
|
Dynamical Effects of Solvation on Norbornadiene/Quadricyclane Systems. J Phys Chem A 2024; 128:2602-2610. [PMID: 38511966 DOI: 10.1021/acs.jpca.4c00045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
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.
Collapse
|
4
|
Cluster perturbation theory IX: Perturbation series for the coupled cluster singles and doubles ground state energy. J Chem Phys 2024; 160:104108. [PMID: 38477336 DOI: 10.1063/5.0192388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 02/17/2024] [Indexed: 03/14/2024] Open
Abstract
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.
Collapse
|
5
|
Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies. J Chem Theory Comput 2024; 20:1932-1943. [PMID: 38380846 DOI: 10.1021/acs.jctc.3c01038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
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.
Collapse
|
6
|
Extended Tetrathiafulvalenes with Fluoreno[3,2-b]fluorene and Diindeno[1,2-b : 1',2'-i]anthracene Cores. Chemistry 2024; 30:e202302688. [PMID: 37930277 DOI: 10.1002/chem.202302688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 11/03/2023] [Accepted: 11/06/2023] [Indexed: 11/07/2023]
Abstract
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.
Collapse
|
7
|
Investigation of Solvent Effects on the Molecular Energy Storage and Optical Properties of Bicyclooctadiene/Tetracyclooctane Photoswitches. J Phys Chem A 2024; 128:41-50. [PMID: 38152898 DOI: 10.1021/acs.jpca.3c04760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
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.
Collapse
|
8
|
Searching the Chemical Space of Bicyclic Dienes for Molecular Solar Thermal Energy Storage Candidates. Angew Chem Int Ed Engl 2023; 62:e202309543. [PMID: 37489860 DOI: 10.1002/anie.202309543] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Revised: 07/21/2023] [Accepted: 07/25/2023] [Indexed: 07/26/2023]
Abstract
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.
Collapse
|
9
|
The Effects of Solvent Dynamics on the Back Reaction of Solar-Thermal Energy Storage Systems. J Phys Chem A 2023; 127:7058-7069. [PMID: 37607346 DOI: 10.1021/acs.jpca.3c02294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/24/2023]
Abstract
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.
Collapse
|
10
|
Computational investigation of photoswitch conjugates for molecular solar energy storage. Phys Chem Chem Phys 2023; 25:21964-21969. [PMID: 37554092 DOI: 10.1039/d3cp02555a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
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.
Collapse
|
11
|
Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1256510. [PMID: 37654900 PMCID: PMC10466216 DOI: 10.3389/fchem.2023.1256510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 07/11/2023] [Indexed: 09/02/2023] Open
Abstract
[This corrects the article DOI: 10.3389/fchem.2023.1154526.].
Collapse
|
12
|
Abstract
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 Ag16-DNA complex have been carried out and compared with pure time-dependent density functional theory calculations of two Ag16 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.
Collapse
|
13
|
Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1154526. [PMID: 37388945 PMCID: PMC10303140 DOI: 10.3389/fchem.2023.1154526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 05/11/2023] [Indexed: 07/01/2023] Open
Abstract
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.
Collapse
|
14
|
Massively Parallel GPU Enabled Third Order Cluster Perturbation Excitation Energies for Cost-Effective Large Scale Excitation Energy Calculations. J Chem Phys 2023; 158:144111. [PMID: 37061462 DOI: 10.1063/5.0142780] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2023] Open
Abstract
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.
Collapse
|
15
|
Indirect nuclear spin-spin couplings with third order contributions added to the SOPPA method. J Chem Phys 2023; 158:124118. [PMID: 37003784 DOI: 10.1063/5.0140117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023] Open
Abstract
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.
Collapse
|
16
|
Modeling Absorption and Emission Spectroscopies of Symmetric and Asymmetric Azaoxahelicenes in Vacuum and Solution. J Phys Chem A 2022; 126:6467-6472. [PMID: 36095036 DOI: 10.1021/acs.jpca.2c05721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
|
17
|
Cluster perturbation theory. VIII. First order properties for a coupled cluster state. J Chem Phys 2022; 157:024108. [DOI: 10.1063/5.0082585] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
|
18
|
Electric Properties of Photochromic Molecules Physisorbed on Silver and Copper Nanoparticles. J Phys Chem A 2022; 126:3145-3156. [PMID: 35583037 DOI: 10.1021/acs.jpca.2c01003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
|
19
|
Prospects of Improving Molecular Solar Energy Storage of the Norbornadiene/Quadricyclane System through Bridgehead Modifications. J Phys Chem A 2022; 126:2670-2676. [PMID: 35467862 DOI: 10.1021/acs.jpca.2c00950] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
|
20
|
Cluster Perturbation Theory. VII. The convergence of Cluster Perturbation Expansions. J Chem Phys 2022; 157:024107. [DOI: 10.1063/5.0082584] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
|
21
|
Cluster perturbation theory. VI. Ground-state energy series using the Lagrangian. J Chem Phys 2022; 157:024106. [DOI: 10.1063/5.0082583] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
|
22
|
Optimization of the thermochemical properties of the norbornadiene/quadricyclane photochromic couple for solar energy storage using nanoparticles. Phys Chem Chem Phys 2022; 24:5506-5521. [PMID: 35171973 DOI: 10.1039/d2cp00226d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
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.
Collapse
|
23
|
The effects of solvation on the back reaction and storage capabilities of solar thermal energy storage systems. Phys Chem Chem Phys 2022; 24:5564-5577. [PMID: 35174838 DOI: 10.1039/d2cp00401a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
|
24
|
Synthesis, characterization and computational evaluation of bicyclooctadienes towards molecular solar thermal energy storage. Chem Sci 2022; 13:834-841. [PMID: 35173948 PMCID: PMC8768882 DOI: 10.1039/d1sc05791j] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Accepted: 12/20/2021] [Indexed: 01/17/2023] Open
Abstract
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.
Collapse
|
25
|
High Throughput Screening of Norbornadiene/Quadricyclane Derivates for Molecular Solar Thermal Energy Storage. Phys Chem Chem Phys 2022; 24:28956-28964. [DOI: 10.1039/d2cp03032b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
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...
Collapse
|
26
|
Excited State Dynamics and Conjugation Effects of the Photoisomerization Reactions of Dihydroazulene. Phys Chem Chem Phys 2022; 24:28934-28943. [DOI: 10.1039/d2cp02706b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
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...
Collapse
|
27
|
Investigation of the Structural and Thermochemical Properties of [2.2.2]-Bicyclooctadiene Photoswitches. J Phys Chem A 2021; 125:10330-10339. [PMID: 34809434 DOI: 10.1021/acs.jpca.1c07737] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
|
28
|
Abstract
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).
Collapse
|
29
|
Fulvalene-Based Polycyclic Aromatic Hydrocarbon Ladder-Type Structures: Synthesis and Properties. Chemistry 2021; 27:8315-8324. [PMID: 33856724 DOI: 10.1002/chem.202100984] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Indexed: 12/11/2022]
Abstract
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.
Collapse
|
30
|
Promoting the thermal back reaction of vinylheptafulvene to dihydroazulene by physisorbtion on nanoparticles. Phys Chem Chem Phys 2021; 23:12889-12899. [PMID: 34075905 DOI: 10.1039/d0cp02893b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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 TiO2 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 106. This illustrates the prospects of utilizing nanoparticles for controlling the release of the stored thermal energy.
Collapse
|
31
|
Molecular solar thermal energy storage properties of photochromic molecules physisorbed onto nanoparticles. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136661] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
32
|
The influence of nanoparticles on the excitation energies of the photochromic dihydroazulene/vinylheptafulvene system. Phys Chem Chem Phys 2019; 21:6689-6698. [DOI: 10.1039/c8cp06539j] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
This paper studies how nanoparticles affect photochromic systems, focusing on the influence of gold nanoparticles on the optical properties of the dihydroazulene/vinylheptafulvene (DHA/VHF) system.
Collapse
|
33
|
The influence of gold nanoparticles on the two photon absorption of photochromic molecular systems. Phys Chem Chem Phys 2019; 21:18577-18588. [DOI: 10.1039/c9cp03283e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
In this study, we investigate the influence of gold nanoparticles on the nonlinear optical properties of the dihydroazulene/vinylheptafulvene photo- and thermochromic system.
Collapse
|