Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers.

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

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.

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) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models.

The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue.

The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S1 and dark S2 surfaces near the S1 minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S1/S2 state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization.

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. In this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

New Photoelectron-Valence Electron Interactions Evident in the Photoelectron Spectrum of Gd2O.

Evidence of strong photoelectron-valence electron (PEVE) interactions has been observed in the anion photoelectron (PE) spectra of several lanthanide suboxide clusters, which are exceptionally complex from an electronic structure standpoint and are strongly correlated systems. The PE spectrum of Gd2O-, which should have relatively simple electronic structure because of its half-filled 4f subshell, exhibits numerous electronic transitions. The electron affinity determined from the spectrum is 0.26 eV. The intensities of transitions to excited states increase relative to the lower-energy states with lower photon energy, which is consistent with shakeup transitions driven by time-dependent electron-neutral interactions. A group of intense spectral features that lie between electron binding energies of 0.7 and 2.3 eV are assigned to transitions involving detachment of an electron from outer-valence σu and σg orbitals that have large Gd 6s contributions. The spectra show parallel transition manifolds in general, which is consistent with detachment from these orbitals. However, several distinct perpendicular transitions are observed adjacent to several of the vertical transitions. A possible explanation invoking interaction between the ejected electron and the high-spin neutral is proposed. Specifically, the angular momentum of electrons ejected from σu or σg orbitals, which is l = 1, can switch to l = 0, 2 with an associated change in the Ms of the remnant neutral, which is spin-orbit coupling between a free electron and the spin of a neutral.

A Macrocyclic 1,4-Diketone Enables the Synthesis of a p-Phenylene Ring That Is More Strained than a Monomer Unit of [4]Cycloparaphenylene.

The synthesis of a p-terphenyl-based macrocycle, containing a p-phenylene unit with 42.6 kcal/mol of strain energy (SE), is reported. The conversion of a macrocyclic 1,4-diketone to a highly strained arene system takes place over five synthetic steps, featuring iterative dehydrative reactions in the aromatization protocol. Spectroscopic data of the deformed benzenoid macrocycle are in excellent agreement with other homologues that have been reported, indicating that the central p-phenylene ring of 9 is aromatic.

Overcoming Strain-Induced Rearrangement Reactions: A Mild Dehydrative Aromatization Protocol for Synthesis of Highly Distorted p-Phenylenes.

A series of p-terphenyl-based macrocycles, containing highly distorted p-phenylene units, have been synthesized. Biaryl bonds of the nonplanar p-terphenyl nuclei were constructed in the absence of Pd-catalyzed or Ni-mediated cross-coupling reactions, using 1,4-diketones as surrogates to strained arene units. A streamlined synthetic protocol for the synthesis of 1,4-diketo macrocycles has been developed, using only 2.5 mol % of the Hoveyda-Grubbs second-generation catalyst in both metathesis and transfer hydrogenation reactions. Under protic acid-mediated dehydrative aromatization conditions, the central and most strained benzene ring of the p-terphenyl systems was susceptible to rearrangement reactions. To overcome this, a dehydrative aromatization protocol using the Burgess reagent was developed. Under these conditions, no strain-induced rearrangement reactions occur, delivering p-phenylene units with up to 28.4 kcal/mol strain energy and deformation angles that sum up to 40°.