Search

Drachmann's regularization approach is implemented for floating explicitly correlated Gaussians (fECGs) and molecular systems. Earlier applications of drachmannized relativistic corrections for molecular systems were hindered due to the unknown analytic matrix elements of 1/rix1/rjy-type operators with fECGs. In the present work, one of the 1/r factors is approximated by a linear combination of Gaussians, which results in calculable integrals. The numerical approach is found to be precise and robust over a range of molecular systems and nuclear configurations, and thus, it opens the route toward an automated evaluation of high-precision relativistic corrections over potential energy surfaces of polyatomic systems. Furthermore, the newly developed integration approach makes it possible to construct the matrix representation of the square of the electronic Hamiltonian relevant for energy lower-bound as well as time-dependent computations of molecular systems with a flexible and high-precision fECG basis representation., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

The aim of this study was to investigate the transition from non-covalent reversible over covalent reversible to covalent irreversible inhibition of cysteine proteases by making delicate structural changes to the warhead scaffold. To this end, dipeptidic rhodesain inhibitors with different N -terminal electrophilic arenes as warheads relying on the S N Ar mechanism were synthesized and investigated. Strong structure-activity relationships of the inhibition potency, the degree of covalency, and the reversibility of binding on the arene substitution pattern were found. The studies were complemented and substantiated by molecular docking and quantum-mechanical calculations of model systems. Furthermore, the improvement in the membrane permeability of peptide esters in comparison to their corresponding carboxylic acids was exemplified.

Interactions in atomic and molecular systems are dominated by electromagnetic forces and the theoretical framework must be in the quantum regime. The physical theory for the combination of quantum mechanics and electromagnetism, quantum electrodynamics has been "established" by the mid-twentieth century, primarily as a scattering theory. To describe atoms and molecules, it is important to consider bound states. In the nonrelativistic quantum mechanics framework, bound states can be efficiently computed using robust and general methodologies with systematic approximations developed for solving wave equations . With the sight of the development of a computational quantum electrodynamics framework for atomic and molecular matter, the field theoretic Bethe-Salpeter wave equation expressed in space-time coordinates, its exact equal-time variant, and emergence of a relativistic wave equation, is reviewed. A computational framework, with initial applications and future challenges in relation with precision spectroscopy, is also highlighted., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

A general computational scheme for the (nonrelativistic) Bethe logarithm is developed, opening the route to "routine" evaluation of the leading-order quantum electrodynamics correction (QED) relevant for spectroscopic applications for small polyatomic and polyelectronic molecular systems. The implementation relies on the Schwartz method and minimization of a Hylleraas functional. In relation with electronically excited states, a projection technique is considered, which ensures positive definiteness of the functional over the entire parameter (photon momentum) range. Using this implementation, the Bethe logarithm is converged to a relative precision better than 1:10 3 for selected electronic states of the two-electron H 2 and H 3 + , and the three-electron He 2 + and H+H 2 molecular systems. The present work focuses on nuclear configurations near the local minimum of the potential energy surface, but the computations can be repeated also for other structures.

Variational and perturbative relativistic energies are computed and compared for two-electron atoms and molecules with low nuclear charge numbers. In general, good agreement of the two approaches is observed. Remaining deviations can be attributed to higher-order relativistic, also called non-radiative quantum electrodynamics (QED), corrections of the perturbative approach that are automatically included in the variational solution of the no-pair Dirac-Coulomb-Breit (DCB) equation to all orders of the α fine-structure constant. The analysis of the polynomial α dependence of the DCB energy makes it possible to determine the leading-order relativistic correction to the non-relativistic energy to high precision without regularization. Contributions from the Breit-Pauli Hamiltonian, for which expectation values converge slowly due the singular terms, are implicitly included in the variational procedure. The α dependence of the no-pair DCB energy shows that the higher-order (α 4 E h ) non-radiative QED correction is 5% of the leading-order (α 3 E h ) non-radiative QED correction for Z = 2 (He), but it is 40% already for Z = 4 (Be 2+ ), which indicates that resummation provided by the variational procedure is important already for intermediate nuclear charge numbers.

The Breit interaction is implemented in the no-pair variational Dirac-Coulomb (DC) framework using an explicitly correlated Gaussian basis reported in the previous paper [P. Jeszenszki, D. Ferenc, and E. Mátyus, J. Chem. Phys. 156, 084111 (2022)]. Both a perturbative and a fully variational inclusion of the Breit term are considered. The no-pair DC plus perturbative Breit and the no-pair DC-Breit energies are compared with perturbation theory results including the Breit-Pauli Hamiltonian and leading-order non-radiative quantum electrodynamics corrections for low Z values. Possible reasons for the observed deviations are discussed.

The Dirac-Coulomb equation with positive-energy projection is solved using explicitly correlated Gaussian functions. The algorithm and computational procedure aims for a parts-per-billion convergence of the energy to provide a starting point for further comparison and further developments in relation with high-resolution atomic and molecular spectroscopy. Besides a detailed discussion of the implementation of the fundamental spinor structure, permutation, and point-group symmetries, various options for the positive-energy projection procedure are presented. The no-pair Dirac-Coulomb energy converged to a parts-per-billion precision is compared with perturbative results for atomic and molecular systems with small nuclear charge numbers. Paper II [D. Ferenc, P. Jeszenszki, and E. Mátyus, J. Chem. Phys. 156, 084110 (2022).] describes the implementation of the Breit interaction in this framework.

The quantum harmonic model and the two-phase thermodynamic method (2PT) are widely used to obtain quantum-corrected properties such as isobaric heat capacities or molar entropies. 2PT heat capacities were calculated inconsistently in the literature. For water, the classical heat capacity was also considered, but for organic liquids, it was omitted. We reanalyzed the performance of different quantum corrections on the heat capacities of common organic solvents against experimental data. We have pointed out serious flaws in previous 2PT studies. The vibrational density of states was calculated incorrectly causing a 39% relative error in diffusion coefficients and 45% error in the 2PT heat capacities. The wrong conversion of isobaric and isochoric heat capacities also caused about 40% error but in the other direction. We have introduced the concept of anharmonic correction (AC), which is simply the deviation of the classical heat capacity from that of the harmonic oscillator model. This anharmonic contribution is around +30 to 40 J/(mol K) for water depending on the water model and -8 to -10 J/(mol K) for hydrocarbons and halocarbons. AC is unrealistically large, +40 J/(K mol) for alcohols and amines, indicating some deficiency of the OPLS force field. The accuracy of the computations was also assessed with the determination of the self-diffusion coefficients.