Structure of [18]Annulene Revisited: Challenges for Computing Benzenoid Systems.

For cyclic conjugated structures, erratic computational results have been obtained with Hartree-Fock (HF) molecular orbital (MO) methods as well as density functional theory (DFT) with large HF-exchange contributions. In this work, the reasons for this unreliability are explored. Extensive computations on [18]annulene and related compounds highlight the pitfalls to be avoided and the due diligence required for such computational investigations. In particular, a careful examination of the MO singlet-stability eigenvalues is recommended. The appearance of negative eigenvalues is not (necessarily) problematic, but near-zero (positive or negative) eigenvalues can lead to dramatic errors in vibrational frequencies and related properties. DFT approaches with a lower HF admixture generally appear more robust in this regard for the description of benzenoid structures, although they may exaggerate the tendency toward planarity and C-C bond-equalization. For the iconic [18]annulene, the results support a nonplanar equilibrium structure. The density-fitted frozen natural orbital coupled-cluster singles and doubles with perturbative triples [DF-FNO CCSD(T)] method of electron correlation with an aug-pVQZ/aug-pVTZ basis set places the C2 global minimum 1.1 kcal mol-1 below the D6h stationary point.

Reduced Scaling Real-Time Coupled Cluster Theory.

Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.

Challenges in the Use of Quantum Computing Hardware-Efficient Ansätze in Electronic Structure Theory.

Advances in quantum computation for electronic structure, and particularly heuristic quantum algorithms, create an ongoing need to characterize the performance and limitations of these methods. Here we discuss some potential pitfalls connected with the use of hardware-efficient Ansätze in variational quantum simulations of electronic structure. We illustrate that hardware-efficient Ansätze may break Hamiltonian symmetries and yield nondifferentiable potential energy curves, in addition to the well-known difficulty of optimizing variational parameters. We discuss the interplay between these limitations by carrying out a comparative analysis of hardware-efficient Ansätze versus unitary coupled cluster and full configuration interaction, and of second- and first-quantization strategies to encode Fermionic degrees of freedom to qubits. Our analysis should be useful in understanding potential limitations and in identifying possible areas of improvement in hardware-efficient Ansätze.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Modeling Complex Solvent Effects on the Optical Rotation of Chiral Molecules: A Combined Molecular Dynamics and Density Functional Theory Study.

The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.

Structure Elucidation and Confirmation of Phloroglucinols from the Roots of Garcinia dauphinensis by Comparison of Experimental and Calculated ECD Spectra and Specific Rotations.

Eight phloroglucinols from Garcinia dauphinensis were recently reported to have good to moderate antiplasmodial and anticancer activities, consistent with other phloroglucinol derivatives isolated from natural sources. Chiroptical properties were previously calculated and compared to experimental data for compound 2 as a means to deduce its absolute configuration. Tentative assignments for the remaining compounds were also reported based on these data. In order to arrive at stereochemical assignments for phloroglucinols 1 and 3-8, ECD spectra and specific rotations were computed for all stereoisomers of each compound. Molecular orbital analyses were also carried out for the most energetically favorable conformers of each compound. Absolute configurations are reported for all eight phloroglucinols for the first time.

Machine-Learning Coupled Cluster Properties through a Density Tensor Representation.

The introduction of machine-learning (ML) algorithms to quantum mechanics enables rapid evaluation of otherwise intractable expressions at the cost of prior training on appropriate benchmarks. Many computational bottlenecks in the evaluation of accurate electronic structure theory could potentially benefit from the application of such models, from reducing the complexity of the underlying wave function parameter space to circumventing the complications of solving the electronic Schrödinger equation entirely. Applications of ML to electronic structure have thus far been focused on learning molecular properties (mainly the energy) from geometric representations. While this line of study has been quite successful, highly accurate models typically require a "big data" approach with thousands of training data points. Herein, we propose a general, systematically improvable scheme for wave function-based ML of arbitrary molecular properties, inspired by the underlying equations that govern the canonical approach to computing the properties. To this end, we combine the established ML machinery of the t-amplitude tensor representation with a new reduced density matrix representation. The resulting model provides quantitative accuracy in both the electronic energy and dipoles of small molecules using only a few dozen training points per system.

Theory and implementation of a novel stochastic approach to coupled cluster.

We present a detailed discussion of our novel diagrammatic coupled cluster Monte Carlo (diagCCMC) [Scott et al. J. Phys. Chem. Lett. 10, 925 (2019)]. The diagCCMC algorithm performs an imaginary-time propagation of the similarity-transformed coupled cluster Schrödinger equation. Imaginary-time updates are computed by the stochastic sampling of the coupled cluster vector function: each term is evaluated as a randomly realized diagram in the connected expansion of the similarity-transformed Hamiltonian. We highlight similarities and differences between deterministic and stochastic linked coupled cluster theory when the latter is re-expressed as a sampling of the diagrammatic expansion and discuss details of our implementation that allow for a walker-less realization of the stochastic sampling. Finally, we demonstrate that in the presence of locality, our algorithm can obtain a fixed errorbar per electron while only requiring an asymptotic computational effort that scales quartically with system size, independent of the truncation level in coupled cluster theory. The algorithm only requires an asymptotic memory cost scaling linearly, as demonstrated previously. These scaling reductions require no ad hoc modifications to the approach.

Psi4 1.4: Open-source software for high-throughput quantum chemistry.

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

Basis Set Superposition Errors in the Many-Body Expansion of Molecular Properties.

The underlying reasons for the poor convergence of the venerated many-body expansion (MBE) for higher-order response properties are investigated, with a particular focus on the impact of basis set superposition errors. Interaction energies, dipole moments, dynamic polarizabilities, and specific rotations are computed for three chiral solutes in explicit water cages of varying sizes using the MBE including corrections based on the site-site function counterpoise (or "full-cluster" basis) approach. In addition, we consider other possible causes for the observed oscillatory behavior of the MBE, including numerical precision, basis set size, choice of density functional, and snapshot geometry. Our results indicate that counterpoise corrections are necessary for damping oscillations and achieving reasonable convergence of the MBE for higher order properties. However, oscillations in the expansion cannot be completely eliminated for chiroptical properties such as specific rotations due to their inherently nonadditive nature, thus limiting the efficacy of the MBE for studying solvated chiral compounds.

Phloroglucinols from the Roots of Garcinia dauphinensis and Their Antiproliferative and Antiplasmodial Activities.

Garcinia dauphinensis is a previously uninvestigated endemic plant species of Madagascar. The new phloroglucinols dauphinols A-F and 3'-methylhyperjovoinol B (1-7) and six known phloroglucinols (8-13) together with tocotrienol 14 and the three triterpenoids 15-17 were isolated from an ethanolic extract of G. dauphinensis roots using various chromatographic techniques. The structures of the isolated compounds were elucidated by NMR, MS, optical rotation, and ECD data. Theoretical ECD spectra and specific rotations for 2 were calculated and compared to experimental data in order to assign its absolute configuration. Among the compounds tested, 1 showed the most promising growth inhibitory activity against A2870 ovarian cancer cells, with IC50 = 4.5 ± 0.9 µM, while 2 had good antiplasmodial activity against the Dd2 drug-resistant strain of Plasmodium falciparum, with IC50 = 0.8 ± 0.1 µM.

Calculating Optical Rotatory Dispersion Spectra in Solution Using a Smooth Dielectric Model.

The calculation of specific rotation of molecules in solution is probed at the coupled cluster (CC) level utilizing a continuum dielectric model based on a definition of the dielectric permittivity as a smooth function of electron density. Solvation effects are captured through polarization of Hartree-Fock (HF) molecular orbitals before subsequent calculations with the coupled cluster singles and doubles (CCSD) method. For the challenging ( S)-methyloxirane molecule, CCSD specific rotations yield an incorrect sign for the rotation in water, and the continuum model is unable to predict the wide variations in the optical rotatory dispersion (ORD) curves seen for nonpolar solvents of similar dielectric constant. In two molecules, (1 S,4 S)-norbornenone and ( S)-2-chloropropionitrile, specific rotations computed with CCSD in conjunction with implicit solvent fail to provide solvent shifts of the correct order of magnitude, indicating that the solvent response is a major contribution to the overall solvation effect.

Performance of Property-Optimized Basis Sets for Optical Rotation with Coupled Cluster Theory.

The effectiveness of the optical rotation prediction (ORP) basis set for computing specific rotations at the coupled cluster (CC) level has been evaluated for a test set of 14 chiral compounds. For this purpose, the ORP basis set has been developed for the second-row atoms present in the investigated systems (that is, for sulfur, phosphorus, and chlorine). The quality of the resulting set was preliminarily evaluated for seven molecules using time-dependent density-functional theory (TD-DFT). Rotations were calculated with the coupled cluster singles and doubles method (CCSD) as well as the second-order approximate coupled cluster singles and doubles method (CC2) with the correlation-consistent aug-cc-pVDZ and aug-cc-pVTZ basis sets and extrapolated to estimate the complete basis-set (CBS) limit for comparison with the ORP basis set. In the compounds examined here, the ORP calculations on molecules containing only first-row atoms compare favorably with results from the larger aug-cc-pVTZ basis set, in some cases lying closer to the estimated CBS limit, while results for molecules containing second-row atoms indicate that larger correlation-consistent basis sets are necessary to obtain reliable estimates of the CBS limit.

Perspective: Computational chemistry software and its advancement as illustrated through three grand challenge cases for molecular science.

The field of computational molecular sciences (CMSs) has made innumerable contributions to the understanding of the molecular phenomena that underlie and control chemical processes, which is manifested in a large number of community software projects and codes. The CMS community is now poised to take the next transformative steps of better training in modern software design and engineering methods and tools, increasing interoperability through more systematic adoption of agreed upon standards and accepted best-practices, overcoming unnecessary redundancy in software effort along with greater reproducibility, and increasing the deployment of new software onto hardware platforms from in-house clusters to mid-range computing systems through to modern supercomputers. This in turn will have future impact on the software that will be created to address grand challenge science that we illustrate here: the formulation of diverse catalysts, descriptions of long-range charge and excitation transfer, and development of structural ensembles for intrinsically disordered proteins.