Search

Quantum Monte Carlo Algebraic Diagrammatic Construction (QMCADC) has been proposed as a reformulation of the second-order ADC scheme for the polarization propagator within the projection quantum Monte Carlo formalism. Dense-sparse partitioning and importance ranking filtering strategies are now exploited to accelerate its convergence and to alleviate the sign problem inherent in such calculations. By splitting the configuration space into dense and sparse subsets, the corresponding projection operator is decomposed into four distinct blocks. Deterministic calculations handle the dense-to-dense and sparse-to-dense blocks, while the remaining blocks, dense-to-sparse and sparse-to-sparse, are stochastically evaluated. The dense set is efficiently stored in a fixed-size array, and the sparse set is represented through conventional floating random Monte Carlo walks. The stochastic projection is further refined through importance ranking criteria, enabling a reduction in the required number of walkers with a controllable bias. Our results demonstrate the integration of dense-sparse partitioning with importance ranking filtering to significantly enhance the efficiency of QMCADC, enabling large-scale molecular excited-state calculations. Furthermore, this novel approach maximizes the utilization of the sparsity of ADC(2), transforming QMCADC into a tailored framework for ADC calculations. [ABSTRACT FROM AUTHOR]

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. Calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

The surprisingly low rupture force and remarkable mechanical anisotropy of rubredoxin have been known for several years. Exploiting the first combination of steered molecular dynamics and the quantum chemical Judgement of Energy DIstribution (JEDI) analysis, the common belief that hydrogen bonds between neighboring amino acid backbones and the sulfur atoms of the central FeS4 unit in rubredoxin determine the low mechanical resistance of the protein is invalidated. The distribution of strain energy in the central part of rubredoxin is elucidated in real-time with unprecedented detail, giving important insights into the mechanical unfolding pathway of rubredoxin. While structural anisotropy as well as the contribution of angle bendings in the FeS4 unit have a significant influence on the mechanical properties of rubredoxin, these factors are insufficient to explain the experimentally observed low rupture force. Instead, the rupture mechanism of rubredoxin is far more complex than previously thought and requires more than just a hydrogen bond network.

In this paper we present the Markovian Multiagent Monte-Carlo Second Order Self-Consistent Field Algorithm (M3-SOSCF). This algorithm provides a highly reliable methodology for converging SCF calculations in single-reference methods using a modified differential evolution approach. Additionally, M3 is embarrassingly parallel and modular in regards to Newton–Raphson subroutines. We show that M3 is able to surpass contemporary SOSCFs in reliability, which is illustrated by a benchmark employing poor initial guesses and a second benchmark with SCF calculations which face difficulties using standard SCF algorithms. Furthermore, we analyse inherent properties of M3 and show that in addition to its robustness and efficiency, it is more user-friendly than current SOSCFs. [ABSTRACT FROM AUTHOR]

The intermediate state representation (ISR) formalism allows for the straightforward calculation of excited state properties and state-to-state transition moments using the algebraic-diagrammatic construction (ADC) scheme for the polarization propagator. Here, the derivation and implementation of the ISR in third-order perturbation theory for the one-particle operator are presented, enabling, for the first time, the calculation of consistent third-order ADC [ADC(3)] properties. The accuracy of ADC(3) properties is evaluated with respect to high-level reference data and compared to the previously used ADC(2) and ADC(3/2) schemes. Oscillator strengths and excited state dipole moments are computed, and typical response properties are considered: dipole polarizabilities, first-order hyperpolarizabilities, and two-photon absorption strengths. The consistent third-order treatment of the ISR leads to an accuracy similar to that of the mixed-order ADC(3/2) method; the individual performance, however, depends on the property and molecule under investigation. ADC(3) produces slightly improved results in the case of oscillator strengths and two-photon absorption strengths, while excited state dipole moments, dipole polarizabilities, and first-order hyperpolarizabilities exhibit similar accuracy at ADC(3) and ADC(3/2) levels. Taking the significant increase of central processing unit time and memory requirements of the consistent ADC(3) approach into account, the mixed-order ADC(3/2) scheme offers a better compromise between accuracy and efficiency for the properties considered. [ABSTRACT FROM AUTHOR]

We present an implementation for the calculation of molecular response properties using the algebraic-diagrammatic construction (ADC)/intermediate state representation approach. For the second-order ADC model [ADC(2)], a memory-efficient ansatz avoiding the storage of double excitation amplitudes is investigated. We compare the performance of different numerical algorithms for the solution of the underlying response equations for ADC(2) and show that our approach also strongly improves the convergence behavior for the investigated algorithms compared with the standard implementation. All routines are implemented in an open-source Python library. [ABSTRACT FROM AUTHOR]

Molecular photoswitches are widely used in material sciences, physics, chemistry, and biology. As needs grow more complex, materials have to react more than one‐dimensionally. The use of multiple photoswitches at once opens manifold opportunities for further improved and more complicated systems. However, this requires independent addressability, i.e., orthogonality, and reversible processes. Herein, the first study on ultrafast excited state dynamics of two orthogonal photoswitches, a push‐pull azobenzene and a donor‐acceptor Stenhouse adduct is reported. In order to gain detailed insight in their interactions and mutual influences on their photoswitching behavior, they are addressed individually and simultaneously via transient absorption spectroscopy supported by quantum chemical calculations. They show reversible photoswitchability and in addition, can be used in 4D printing to provide easy access to a plethora of functional devices. Furthermore, environmental influences on the excited state dynamics are examined using different solvents and thin films. Both compounds photoisomerize independently when addressed individually or simultaneously and only little impacts on the excited state dynamics are found. Especially the vibrational relaxation is affected by different surroundings changing the energy dissipation while hardly affecting the electronic states involved. The orthogonal and simultaneous addressability is thereby crucial for their usage in 4D printed microactuators. [ABSTRACT FROM AUTHOR]

Compounds with multiple photoswitching units are appealing for complex photochemical control of molecular materials and nanostructures. Herein, we synthesized novel meta‐ and para‐ connected (related to the nitrogen of the indoline) azobenzene‐spiropyran dyads, in which the central benzene unit is shared by both switches. We investigated their photochemistry using static and time‐resolved transient absorption spectroscopy as well as quantum chemical calculations. In the meta‐compound, the individual components are photochemically decoupled due to the meta‐pattern. In the para‐compound the spiro‐connectivity leads to a bifunctional photoswitchable system with a red‐shifted absorption. The azobenzene and the spiropyran can thus be addressed and switched independently by light of appropriate wavelength. Through the different connectivity patterns two different orthogonally photoswitchable systems have been obtained which are promising candidates for complex applications of light control. [ABSTRACT FROM AUTHOR]

A π‐expanded X‐type double [5]helicene comprising dihydropyracylene moieties was synthesized from commercially available acenaphthene. X‐ray crystallographic analysis revealed the unique highly twisted structure of the compound resulting in the occurrence of two enantiomers which were separated by chiral HPLC, owing to their high conformational stability. The compound shows strongly bathochromically shifted UV/vis absorption and emission bands with small Stokes shift and considerable photoluminescence quantum yield and circular polarized luminescence response. The electrochemical studies revealed five facilitated reversible redox events, including three reductions and two oxidations, thus qualifying the compound as chiral multistage redox amphoter. The experimental findings are in line with the computational studies based on density functional theory pointing towards increased spatial extension of the frontier molecular orbitals over the polycyclic framework and a considerably narrowed HOMO–LUMO gap. [ABSTRACT FROM AUTHOR]

Different approaches to mixed-order algebraic-diagrammatic construction (ADC) schemes are investigated. The performance of two different strategies for scaling third-order contributions to the ADC secular matrix is evaluated. Both considered schemes employ a single tuning parameter and conserve general properties inherent to all ADC methods, such as hermiticity and size-consistency. The first approach, scaled-matrix ADC[(2) + x(3)], scales all contributions first occurring in ADC(3) equally and leads to an improvement of the accuracy of excitation energies compared to ADC(3) for x = 0.4–0.5. However, with respect to excited state dipole moments, this method provides lower accuracy than ADC(3). The second scaling approach, MP[(1) + x(2)] − ISR(3), scales the second order contributions of the ground-state wavefunction and derives a rigorous ADC scheme via the intermediate state representation formalism. Although the error in excitation energies is not improved, this method provides insight into the relevance of the individual terms of the ADC(3) matrix and indicates that the MP(2) wavefunction is, indeed, the optimal reference wavefunction for deriving a third-order single-reference ADC scheme. [ABSTRACT FROM AUTHOR]

This article describes a novel approach for the calculation of ionization potentials (IPs), or, more generally, electron-detachment energies, based on a unitary coupled-cluster (UCC) parameterization of the ground-state wave function. Explicit working equations for a scheme referred to as IP-UCC3 are given, providing electron-detachment energies and spectroscopic amplitudes of electron-detached states dominated by one-hole excitations correct through third order. In the derivation, an expansion of the UCC transformed Hamiltonian involving Bernoulli numbers as expansion coefficients is employed. Both the secular matrix and the effective transition moments are shown to be essentially equivalent to the strict third-order algebraic-diagrammatic construction scheme for the electron propagator (IP-ADC). Interestingly, due to the Bernoulli expansion, neglecting triple substitutions in the UCC expansion manifold does not affect the third-order consistency of the IP-UCC effective transition moments. Finally, the equivalence between ADC and UCC excited-state schemes is shown to not hold in fourth or higher order due to a different treatment of the correlated excited-state basis. [ABSTRACT FROM AUTHOR]

The performance of several methods for the calculation of vertical ionization potentials (IPs) or, more generally, electron-detachment energies based on unitary coupled-cluster (UCC) theory and the algebraic-diagrammatic construction (ADC) scheme is evaluated with respect to benchmark data computed at the level of equation-of-motion coupled-cluster theory, including single, double, and triple excitations (IP-EOM-CCSDT). Based on a statistical evaluation of about 200 electron-detached states of 41 molecules, the second-order methods IP-ADC(2) and IP-UCC2 show modest accuracies with IP-EOM-CCSDT as reference, exposing a mean signed error and a standard deviation of the error of −0.54 ± 0.50 and −0.49 ± 0.54 eV, respectively, accompanied by a mean absolute error (MAE) of 0.61 and 0.58 eV, respectively. The strict third-order IP-ADC method demonstrates an accuracy of 0.26 ± 0.35 eV (MAE = 0.35 eV), while the IP-UCC3 method is slightly more accurate with 0.24 ± 0.26 eV (MAE = 0.29 eV). Employing the static self-energy computed using the Dyson expansion method (DEM) improves the IP-ADC(3) performance to 0.27 ± 0.28 eV, with the mean absolute error of this method being 0.32 eV. However, employing the simpler improved fourth-order scheme Σ(4+) for the static self-energy provides almost identical results as the DEM. Based on the quality of the present benchmark results, it therefore appears not necessary to use the computationally more demanding DEM. [ABSTRACT FROM AUTHOR]

The quantum Monte Carlo (QMC) algebraic diagrammatic construction (ADC) method is introduced, which solves the eigenvalue problem of the second-order ADC scheme for the polarization propagator stochastically within the framework of QMC methodology allowing for massively parallel computations. As common virtue of the Monte Carlo integration techniques, quantum Monte Carlo algebraic diagrammatic construction (QMCADC) enables exploitation of the sparsity of the effective ADC matrix, and it reduces the memory requirements by storing only a portion of configurations at each iteration. Furthermore, distributing memory and processing loads to different computing nodes enables the use of fast developing parallel computing resources. Here, the theory and implementation of QMCADC is reported and its viability is demonstrated by the first proof-of-principle calculations. The focus lies on the first excited state and the reproduction of the corresponding lowest vertical excitation energy of various molecular systems. QMCADC is shown to be a genuine stochastic solution of the ADC eigenvalue problem, and exact ADC values can be obtained with a marginal controllable error. [ABSTRACT FROM AUTHOR]