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.

Molecules with an inverted singlet–triplet gap (STG) between the first excited singlet and triplet states, for example, heptazine, have recently been reported and gained substantial attention since they violate the famous Hund’s rule. Utilizing state-of-the-art high-level ab initio methods, the singlet–triplet gap vanishes and approaches zero from below whatever is improved in the theoretical description of the molecules: the basis set or the level of electron correlation. Seemingly, the phenomenon of inverted singlet–triplet gaps tends to vanish the closer we observe.

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.

Orbital-free density functional theory (OF-DFT) holds promise to compute ground state molecular properties at minimal cost. However, it has been held back by our inability to compute the kinetic energy as a functional of electron density alone. Here, we set out to learn the kinetic energy functional from ground truth provided by the more expensive Kohn–Sham density functional theory. Such learning is confronted with two key challenges: Giving the model sufficient expressivity and spatial context while limiting the memory footprint to afford computations on a GPU and creating a sufficiently broad distribution of training data to enable iterative density optimization even when starting from a poor initial guess. In response, we introduce KineticNet, an equivariant deep neural network architecture based on point convolutions adapted to the prediction of quantities on molecular quadrature grids. Important contributions include convolution filters with sufficient spatial resolution in the vicinity of nuclear cusp, an atom-centric sparse but expressive architecture that relays information across multiple bond lengths, and a new strategy to generate varied training data by finding ground state densities in the face of perturbations by a random external potential. KineticNet achieves, for the first time, chemical accuracy of the learned functionals across input densities and geometries of tiny molecules. For two-electron systems, we additionally demonstrate OF-DFT density optimization with chemical accuracy. [ABSTRACT FROM AUTHOR]

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]

Abstract Background Contact tracing is one of the main public health tools in the control of coronavirus disease 2019 (COVID-19). A centralized contact tracing system was developed in Belgium in 2020. We aim to evaluate the performance and describe the results, between January 01, 2021, and September 30, 2021. The characteristics of COVID-19 cases and the impact of COVID-19 vaccination on testing and tracing are also described. Methods We combined laboratory diagnostic test data (molecular and antigen test), vaccination data, and contact tracing data. A descriptive analysis was done to evaluate the performance of contact tracing and describe insights into the epidemiology of COVID-19 by contact tracing. Results Between January and September 2021, 555.181 COVID-19 cases were reported to the central contact center and 91% were contacted. The average delay between symptom onset and contact tracing initiation was around 5 days, of which 4 days corresponded to pre-testing delay. High-Risk Contacts (HRC) were reported by 49% of the contacted index cases. The mean number of reported HRC was 2.7. In total, 666.869 HRC were reported of which 91% were successfully contacted and 89% of these were tested at least once following the interview. The estimated average secondary attack rate (SAR) among the contacts of the COVID-19 cases who reported at least one contact, was 27% and was significantly higher among household HRC. The proportion of COVID-19 cases who were previously identified as HRC within the central system was 24%. Conclusions The contact-tracing system contacted more than 90% of the reported COVID-19 cases and their HRC. This proportion remained stable between January 1 2021 and September 30 2021 despite an increase in cases in March–April 2021. We report high SAR, indicating that through contact tracing a large number of infections were prospectively detected. The system can be further improved by (1) reducing the delay between onset of illness and medical consultation (2) having more exhaustive reporting of HRC by the COVID-19 case.

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]

Abstract Triarylboranes containing linear or angular benzodithiophene moieties and bearing one or two dimesitylboron units were synthesized. The electrochemical and optical features of these compounds were investigated by cyclic voltammetry, UV/Vis and fluorescence spectroscopy while DFT calculations were run to analyze the energetic landscape of these systems. For both linear and angular benzodithiophenes, symmetrical disubstitution leads to the highest photoluminescence yields. The linear benzodithiophene disubstituted with two dimesitylboron units proved to be the most interesting and promising molecule as an electron‐transport material for organic electronics owing to its LUMO energy level of −2.84 eV which is close to those of commonly used electron transport materials like bathocuproine or bathophenantroline.

Until today, perturbation-theoretical consistent algebraic diagrammatic construction (ADC) schemes for the polarization propagator had been derived and implemented up to third order. They have turned out to be versatile and reliable ab initio single-reference methods for the quantum chemical investigation of electronic transitions as well as excited-state properties. Here we present, for the first time, the derivation of consistent fourth-order ADC(4) schemes exploiting novel techniques of automated equation and code generation. The accuracies of the resulting ADC(4) excitation energies have been benchmarked against recent high-level, near exact reference data. The mean absolute error for singly and doubly excited states turns out to be smaller than 0.1 and 0.5 eV, respectively. These developments open also new avenues toward highly accurate ADC methods for electron-detached and attached states. [ABSTRACT FROM AUTHOR]

An atomic-orbital reformulation of the Laplace-transformed scaled opposite-spin (SOS) coupled cluster singles and doubles (CC2) model within the resolution of the identity (RI) approximation (SOS-RI-CC2) is presented that extends its applicability to molecules with several hundreds of atoms and triple-zeta basis sets. We exploit sparse linear algebra and an attenuated Coulomb metric to decrease the disk space demands and the computational efforts. In this way, an effective sub-quadratic computational scaling is achieved with our ω-SOS-CDD-RI-CC2 model. Moreover, Cholesky decomposition of the ground-state one-electron density matrix reduces the prefactor, allowing for an early crossover with the molecular orbital formulation. The accuracy and performance of the presented method are investigated for various molecular systems. [ABSTRACT FROM AUTHOR]

Multidimensional spectroscopic tools are important to explore the details of molecular dynamics. Here the authors use shaped pulses to demonstrate a 3D fluorescence spectroscopy method to extract the fourth and higher-order nonlinear responses in light-molecule interaction.

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]