Fock Matrix Research Articles

Accurate Thermochemistry with Multireference Methods: A Stress Test for Internally Contracted Multireference Coupled-Cluster Theory.

The internally contracted multireference coupled-cluster method with single, double and perturbative triple excitations, icMRCCSD(T), was tested for its performance in the context of computational high-accuracy thermochemistry. The results were gauged against the standard single-reference coupled-cluster hierarchy with up to 5-fold excitations. The test set comprised of a selection of first-row dinuclear compounds and the three 3d-transition metal compounds MnH, FeH, and CoH. The results revealed two problems with the current formulation of icMRCCSD(T). First, the choice of the Dyall Hamiltonian as the zeroth-order Hamiltonian, which leads to a biased description of the different orbital subspaces and particularly poor results for the atomic correlation energies, and second, the tendency to overestimate the perturbative correction for triply excited clusters, in particular in the presence of open shells and correspondingly low orbital-energy gaps. The two problems could be solved by resorting to the effective Fock operator as zeroth-order Hamiltonian and by adopting a modified amplitude equation that includes terms quadratic in the pair clusters. A similar modification was recently proposed by Masios et al. (Phys. Rev. Lett. 2023, 131, 186401) in the context of applying single-reference coupled-cluster theory to systems with small or vanishing band gaps and we chose the acronym '(cT*) correction' in analogy to that work. In contrast to the work of Masios et al., additional terms including single excitation clusters were omitted, as these again lead to an overestimation of correlation effects in more difficult cases. We also tested another alternative for the zeroth-order Hamiltonian and additional higher-order corrections for the correlation energy. These extensions did not significantly improve the results and were also computationally more demanding. The improved icMRCCSD(cT*)F method yields very accurate results with errors, relative to accurate benchmarks, better than 2 kJ/mol for total energies and atomization energies for the entire set of examples considered in this work.

Coupled Perturbed Approach to Dual Basis Sets for Molecules and Solids. II: Energy and Band Corrections for Periodic Systems.

When trying to reach convergence of quantum chemical calculations toward the complete basis set limit, crystalline solids generally prove to be more challenging than molecules. This is due both to the closer packing of atoms─hence, to linear dependencies─and to the problematic behavior of Ewald techniques used for dealing with the infinite character of Coulomb sums. Thus, a dual basis set approach is even more desirable for periodic systems than for molecules. In such an approach, the self-consistent procedure is implemented in a small basis set, and the effect of the enlargement of the basis set is estimated a posteriori. In this paper, we extend to crystalline solids our previous coupled perturbed dual basis set approach [J. Chem. Theory Comput. 2020, 16, 1, 340-353] in which the basis set enlargement is treated as a perturbation. Among the notable features of this approach are (i) the possibility of obtaining not only a correction to the energy but also to energy bands and electron density; (ii) the absence of a diagonalization step for the full Fock matrix in the large basis set; and (iii) the possibility of extrapolating low order perturbation energy corrections to infinite order. We also present here the first periodic implementation of the dual basis set method of Liang and Head-Gordon [J. Phys. Chem. A 2004, 108, 3206-3210]. The effectiveness of both approaches is, then, compared on a small, but representative, set of solids.

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian.

A novel doubly iterative self-consistent field (SCF) approach using a semiempirical model Hamiltonian (denoted as the SMH algorithm) is proposed to accelerate the Hartree-Fock (HF) and density functional theory (DFT) calculations. This method first constructs the Fock matrix exactly in each SCF macroiteration, followed by a few SCF microiterations, where the Fock matrix is incrementally updated using an inexpensive semiempirical approximation. This leads to an improved wave function in each SCF macroiteration with minimal additional cost, and therefore a reduced number of exact Fock builds is required for SCF convergence. The SMH algorithm can be combined with conventional SCF convergence techniques such as level shifting, direct inversion in the iterative subspace (DIIS), and energy-DIIS (EDIIS). When integrated with DIIS, SMH enhances the convergence of simple organic molecules by approximately 10% compared to plain DIIS, with speedups of up to 60% for the more challenging transition metal systems compared to EDIIS + DIIS. Our results show that SMH is a reliable SCF accelerator that seldom deteriorates convergence and is highly robust.

Use of Multigrids to Reduce the Cost of Performing Interpolative Separable Density Fitting.

In this article, we present an interpolative separable density fitting (ISDF)-based algorithm to calculate the exact exchange in periodic mean field calculations. In the past, decomposing the two-electron integrals into the tensor hypercontraction (THC) form using ISDF was the most expensive step of the entire mean field calculation. Here, we show that by using a multigrid-ISDF algorithm, both the memory and the CPU cost of this step can be reduced. The CPU cost is brought down from cubic scaling to quadratic scaling with a low computational prefactor which reduces the cost by almost 2 orders of magnitude. Thus, in the new algorithm, the cost of performing ISDF is largely negligible compared to other steps. Along with the CPU cost, the memory cost of storing the factorized two-electron integrals is also reduced by a factor of up to 35. With the current algorithm, we can perform Hartree-Fock calculations on a diamond supercell containing more than 17,000 basis functions and more than 1500 electrons on a single node with no disk usage. For this calculation, the cost of constructing the exchange matrix is only a factor of 4 slower than the cost of diagonalizing the Fock matrix. Augmenting our approach with linear scaling algorithms can further speed up the calculations.

A modular, composite framework for the utilization of reduced-scaling Coulomb and exchange construction algorithms: Design and implementation.

Multiple algorithms exist for calculating Coulomb (J) or exchange (K) contributions to Fock-like matrices, and it is beneficial to develop a framework that allows the seamless integration and combination of different J and K construction algorithms. In Psi4, we have implemented the "CompositeJK" formalism for this purpose. CompositeJK allows for the combination of any J and K construction algorithms for any quantum chemistry method formulated in terms of J-like or K-like matrices (including, but not limited to, Hartree-Fock and density functional theory) in a highly modular and intuitive fashion, which is simple to utilize for both developers and users. Using the CompositeJK framework, Psi4 was interfaced to the sn-LinK implementation in the GauXC library, adding the first instance of noncommercial graphics processing unit (GPU) support for the construction of Fock matrix elements to Psi4. On systems with hundreds of atoms, the interface to the CPU sn-LinK implementation displays a higher performance than all the alternative JK construction methods available in Psi4, with up to x2.8 speedups compared to existing Psi4JK implementations. The GPU sn-LinK implementation, harnessing the power of GPUs, improves the observed performance gains to up to x7.0.

Revisiting Artifacts of Kohn-Sham Density Functionals for Biosimulation.

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn-Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO-LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver [ Phys. Chem. Chem. Phys. 2024, 26, 6557] to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree-Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping "frontier" orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO-LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

An improved DIIS method using a versatile residual matrix to accelerate SCF starting from a crude guess.

The minimization of the commutator of the Fock and density matrices as the error matrix in the direct inversion of the iterative subspace (CDIIS) developed by Pulay is a powerful self-consistent field (SCF) acceleration technique for the construction of optimum Fock matrix, if initiated with a fair initial guess. In this work, we present an alternative minimized error matrix to the commutator in the CDIIS, namely the residual or the gradient of the energy-functional for a Slater determinant subject to the orthonormality constraints among orbitals, representing the search for a newly improved Fock matrix in the direction of the residual in the direct inversion of the iterative subspace (RDIIS). Implemented in the computational chemistry package GAMESS, the RDIIS is compared with the standard CDIIS and the second order SCF orbital optimization (SOSCF) for tested molecules started with a crude guess. As a result, the RDIIS stably and efficiently performs the SCF convergence acceleration. Furthermore, the RDIIS is considerably independent on the subspace size with the concentrated linear coefficients accounting proportionally for the Fock matrices close to the current iteration.

Synthesis, structural, spectral, Anticancer activity, and density functional theory investigations of 2- [hydrazinylidene(phenyl)methyl]pyridine

The condensation reaction of equimolar amounts of 2-benzoylpyridine and hydrazine hydrate to produce hydrazone, 2- [hydrazinylidene(phenyl)methyl]pyridine (HPMP) with a mass spectral peak at m/z = 197. The compound showed 15N NMR peaks at 330.79, 306.64, and 110.45 ppm due to CN (azomethine), pyridine nitrogen and —NH2 nitrogen atoms. The crystal structure data, IR, and Raman studies were correlated with DFT/B3LYP/6–311++G(d,p) calculations. The compound exhibited strong fluorescence at 592 and 481 nm when excited at 290 and 210 nm, respectively. NBO Fock matrix analyses indicated the existence of different π and π* energy levels, and the molecule can undergo internal charge transfer and vibrational relaxation. The IC-50 value for HPMP was found to be 62.5 µg/ml against the MCF-7 breast cancer cell line as shown by the MTT assay. The IC-50 value for HPMP was higher than the standard doxorubicin but expected to have fewer side effects. Fluorescence microscopy analyses showed significant cell mortality compared to that of the control. The NH2 group in HPMP has been shown to interact with Asp73 (A) (2.89 Å) and ASM46 (A) (2.90 Å) through H-bonding with the DNA gyrase protein (Protein ID: 1KZN). The other selected protein residues were made hydrophobic by pyridine and phenyl rings to inactivate DNA gyrase protein. Therefore, HPMP can be used as an anticancer drug and for optical LED applications.

Unified construction of relativistic Hamiltonians.

It is shown that the four-component (4C), quasi-four-component (Q4C), and exact two-component (X2C) relativistic Hartree-Fock equations can be implemented in a unified manner by making use of the atomic nature of the small components of molecular 4-spinors. A model density matrix approximation can first be invoked for the small-component charge/current density functions, which gives rise to a static, pre-molecular mean field to be combined with the one-electron term. As a result, only the nonrelativistic-like two-electron term of the 4C/Q4C/X2C Fock matrix needs to be updated during the iterations. A "one-center small-component" approximation can then be invoked in the evaluation of relativistic integrals, that is, all atom-centered small-component basis functions are regarded as extremely localized near the position of the atom to which they belong such that they have vanishing overlaps with all small- or large-component functions centered at other nuclei. Under these approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians share precisely the same structure and accuracy. Beyond these is the effective quantum electrodynamics Hamiltonian that can be constructed in the same way. Such approximations lead to errors that are orders of magnitude smaller than other sources of errors (e.g., truncation errors in the one- and many-particle bases as well as uncertainties of experimental measurements) and are, hence, safe to use for whatever purposes. The quaternion forms of the 4C, Q4C, and X2C equations are also presented in the most general way, based on which the corresponding Kramers-restricted open-shell variants are formulated for "high-spin" open-shell systems.

Toward a Stochastic Complete Active Space Second-Order Perturbation Theory.

In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.

Quantum capacitances of transition metal-oxides (CoO, CuO, NiO, and ZnO) doped graphene oxide nanosheet: Insight from DFT computation

Density functional theory (DFT) computation has been utilized to explore the effects of the transition metal oxides: CoO, CuO, NiO, and ZnO doping on the electronic properties, structural, and quantum capacitances of graphene oxide nanosheet. From the magnetic moment analysis CoO@GO was observed to have higher magnetic moment of 11.688 μB compared to the studied the transition metal oxide doped systems. Investigation into the electronic properties revealed that NiO@GO attained higher energy gap with value of 0.144 eV. It was observed that the GO O/C affects the bandgaps of the modelled systems. Perturbation theory analysis of fock matrix showed that CoO@GO and CuO@GO possessed higher second order stabilization energy with values 238.56 kcal/mol and 208.94 kcal/mol respectively. From the quantum capacitance studies, it was observed that the value of CQ for graphene oxide (GO) increased slightly from 72.276 µF/cm2 to ZnO@GO (121.550 µF/cm2) > NiO@GO (93.870 µF/cm2) > CoO@GO (90.52 µF/cm2) > CuO@GO (89.375 µF/cm2). The results obtained herein can provide an effective and simple new idea for the design of graphene-based supercapacitors that possess high energy density.

Methyl orange-doped thiourea barium chloride crystals: Natural bond orbital analysis of the fock matrix for laser applications

Pure thiourea barium chloride (TBC) crystals were synthesised using a slow cooling rate (0.05°C/day) and gradual evaporation from an aqueous solution. The growth process occurred over more than 19 days at a temperature of 39°C. Natural bond orbital (NBO) analysis has proven beneficial for studying charge transfer and conjugative interactions in molecular systems. Moreover, it sheds light on intermolecular bonding, and the interactions between these bonds have implications for optical applications. Additionally, the charge transfer interactions within the molecule that influence its catalytic activity can be elucidated by a lower HOMO and LUMO energy gap. This energy difference can be an indicator of a compound’s excitability; the more readily a compound can be excited, the narrower the energy gap. The energy gap might be significantly affected by the influence of thermal energy. The non-linear optical (NLO) efficiency of methyl orange-doped (MOD) TBC crystals was contrasted with that of the pure form, and the doped version exhibited greater efficiency, particularly for high-frequency applications.

PyBEST: Improved functionality and enhanced performance

PyBEST: Improved functionality and enhanced performance

A hybrid CPU/GPU method for Hartree-Fock self-consistent-field calculation.

The calculation of two-electron repulsion integrals (ERIs) is a crucial aspect of Hartree-Fock calculations. In computing the ERIs of varying angular momentum, both the central processing unit (CPU) and the graphics processing unit (GPU) have their respective advantages. To accelerate the ERI evaluation and Fock matrix generation, a hybrid CPU/GPU method has been proposed to maximize the computational power of both CPU and GPU while overlapping the CPU and GPU computations. This method employs a task queue where each task corresponds to ERIs with the same angular momentum. The queue begins with ERIs of low angular momentum, which are computationally efficient on GPUs, and ends with ERIs of high angular momentum, which are better suited for CPU computation. CPUs and GPUs dynamically grab and complete tasks from the start and end of the queue using OpenMP dynamic scheduling until all tasks are finished. The hybrid CPU/GPU computation offers the advantage of enabling calculations with arbitrary angular momentum. Test calculations showed that the hybrid CPU/GPU algorithm is more efficient than "GPU-only" when using a single GPU. However, as more GPUs are involved, the advantage diminishes or disappears. The scaling exponents of the hybrid method were slightly higher than "GPU-only," but the pre-exponent factor was significantly lower, making the hybrid method more effective overall.

Anti Microbial Activity, Vibrational, Electronic and Topology Analysis of 3, 4, 5-Trimethoxybenzaldehyde

The density functional theory (DFT) approach has become one of the most gainful means to scrutinize the molecular structure and vibrational spectra are finding well-known use in the applications related to biological systems. The experimental and theoretical on the molecular structure, electronic, and characteristics of 3, 4, 5-Trimethoxybenzaldehyde (TMB) are reported. The structure of the molecule was optimized and structural characteristics were determined by DFT using the wb97xd method with 6-311 G++ (d, p) basis sets. The vibrational assignments were made on the basis of potential energy distribution. A good consistency between the observed and calculated spectra was achieved. The intramolecular charge transfer and second-order perturbation theory of Fock Matrix in Natural bond orbital analysis by means NBO analysis. It shows the electronic stability of the compound arising from hyper conjugative interactions and charge delocalization of that molecule. Besides, the HOMO-LUMO, MEP, and Fukui were performed. In this review, the present benzaldehyde derivatives for skirmishing against bacteria strain Escherichia coli and fungi strain Candida Sp. Spectral and quantum chemical calculations are also discussed.

Insights into structural, spectroscopic, and hydrogen bonding interaction patterns of nicotinamide-oxalic acid (form I) salt by using experimental and theoretical approaches.

In the present work, nicotinamide-oxalic acid (NIC-OXA, form I) salt was crystallized by slow evaporation of an aqueous solution. To understand the molecular structure and spectroscopic properties of NIC after co-crystallization with OXA, experimental infrared (IR), Raman spectroscopic signatures, X-ray powder diffraction (XRPD), and differential scanning calorimetry (DSC) techniques were used to characterize and validate the salt. The density functional theory (DFT) methodology was adopted to perform all theoretical calculations by using the B3LYP/6-311++G (d, p) functional/basis set. The experimental geometrical parameters were matched in good correlation with the theoretical parameters of the dimer than the monomer, due to the fact of covering the nearest hydrogen bonding interactions present in the crystal structure of the salt. The IR and Raman spectra of the dimer showed the red (downward) shifting and broadening of bands among (N15-H16), (N38-H39), and (C13=O14) bonds of NIC and (C26=O24), (C3=O1), and (C26=O25) groups of OXA, hence involved in the formation of NIC-OXA salt. The atoms in molecules (AIM) analysis revealed that (N8-H9···O24) is the strongest (conventional) intermolecular hydrogen bonding interaction in the dimer model of salt with the maximum value of interaction energy -12.1kcal mol-1. Furthermore, the natural bond orbital (NBO) analysis of the Fock matrix showed that in the dimer model, the (N8-H9···O24) bond is responsible for the stabilization of the salt with an energy value of 13.44kcal mol-1. The frontier molecular orbitals (FMOs) analysis showed that NIC-OXA (form I) salt is more reactive and less stable than NIC, as the energy gap of NIC-OXA (form I) salt is less than that of NIC. The global and local reactivity descriptor parameters were calculated for the monomer and dimer models of the salt. The electrophilic, nucleophilic, and neutral reactive sites of NIC, OXA, monomer, and dimer models of salt were visualized by plotting the molecular electrostatic potential (MESP) surface. The study provides valuable insights into combining both experimental and theoretical results that could define the physicochemical properties of molecules.

Analytical excited state gradients for time-dependent density functional theory plus tight binding (TDDFT + TB).

Understanding photoluminescent mechanisms has become essential for photocatalytic, biological, and electronic applications. Unfortunately, analyzing excited state potential energy surfaces (PESs) in large systems is computationally expensive, and hence limited with electronic structure methods such as time-dependent density functional theory (TDDFT). Inspired by the sTDDFT and sTDA methods, time-dependent density functional theory plus tight binding (TDDFT + TB) has been shown to reproduce linear response TDDFT results much faster than TDDFT, particularly in large nanoparticles. For photochemical processes, however, methods must go beyond the calculation of excitation energies. Herein, this work outlines an analytical approach to obtain the derivative of the vertical excitation energy in TDDFT + TB for more efficient excited state PES exploration. The gradient derivation is based on the Z vector method, which utilizes an auxiliary Lagrangian to characterize the excitation energy. The gradient is obtained when the derivatives of the Fock matrix, the coupling matrix, and the overlap matrix are all plugged into the auxiliary Lagrangian, and the Lagrange multipliers are solved. This article outlines the derivation of the analytical gradient, discusses the implementation in Amsterdam Modeling Suite, and provides proof of concept by analyzing the emission energy and optimized excited state geometry calculated by TDDFT and TDDFT + TB for small organic molecules and noble metal nanoclusters.

Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy.

Pair atomic density fitting (PADF) has been identified as a promising strategy to reduce the scaling with system size of quantum chemical methods for the calculation of the correlation energy like the direct random-phase approximation (RPA) or second-order Møller-Plesset perturbation theory (MP2). PADF can however introduce large errors in correlation energies as the two-electron interaction energy is not guaranteed to be bounded from below. This issue can be partially alleviated by using very large fit sets, but this comes at the price of reduced efficiency and having to deal with near-linear dependencies in the fit set. One posibility is to use global density fitting (DF), but in this work, we introduce an alternative methodology to overcome this problem that preserves the intrinsically favorable scaling of PADF. We first regularize the Fock matrix by projecting out parts of the basis set which gives rise to orbital products that are hard to describe by PADF. After having thus obtained a reliable self-consistent field solution, we then also apply this projector to the orbital coefficient matrix to improve the precision of PADF-MP2 and PADF-RPA. We systematically assess the accuracy of this new approach in a numerical atomic orbital framework using Slater type orbitals (STO) and correlation consistent Gaussian type basis sets up to quintuple-ζ quality for systems with more than 200 atoms. For the small and medium systems in the S66 database we show the maximum deviation of PADF-MP2 and PADF-RPA relative correlation energies to DF-MP2 and DF-RPA reference results to be 0.07 and 0.14 kcal/mol, respectively. When the new projector method is used, the errors only slightly increase for large molecules and also when moderately sized fit sets are used the resulting errors are well under control. Finally, we demonstrate the computational efficiency of our algorithm by calculating the interaction energies of large, non-covalently bound complexes with more than 1000 atoms and 20000 atomic orbitals at the RPA@PBE/CC-pVTZ level of theory.

Solvent potential effects (topological aspects, electron excitation), spectral characterization and biological attributes of NLO active 1-(2,4-dinitrophenyl)-2-((E)-3-phenylallylidene) hydrazine: Multiple anti tuberculosis agent

Density functional theory (DFT) method has many applications in relation to biological as well as NLO behaviour. It emerged as the most economically viable method for studying molecular structure. DNP2PAH is a synthesized compound with hydroxy derivatives that has undergone experimental and computational evaluation by DFT while being employed in a solvation research inquiry with a variety of solvents (protic, aprotic and non-polar solvents). DNP2PAH has been optimized and parameters were correlated for gas and solvents solvents [C2H3N, DMSO(aprotic), C2H6O, CH3OH(protic) and non-polar solvent (C6H6)]. FT-IR, FT-Raman, and UV–Vis spectrum has been computed and described in relation to practical analyses. It is revealed that observed and computational spectra are very similar. NBO technique utilized second order Fock matrix to assess the donor–acceptor correlations. The prominent interactions π*(C11-C12) - π*(C13-C14) with the highest stabilization energy of 208.65 kcal/ mol. Among solvents, benzene has the least softness (0.35496 (eV)) and the greatest energy gap (2.81721 (eV)). Global descriptors and molecular reactivity of the molecule were significantly changed by the liquids. Substance being studied possesses exceptional NLO (non-linear optical) properties. ELF, LOL and RDG enacted for both the gas phases and solvents. Whereas to convince the biological properties drug-likeness and biological score prediction were utilized. Additionally, this substance follows Lipinski's rule of five, which suggests that, in theory, taking an oral dose of such substances shouldn't be problematic. Results of molecular docking exemplified that the ligand under study is potent drug for anti-tuberculosis and anti-mycobacterium.

Quantum Mechanics/Molecular Mechanics Simulations on NVIDIA and AMD Graphics Processing Units.

We have ported and optimized the graphics processing unit (GPU)-accelerated QUICK and AMBER-based ab initio quantum mechanics/molecular mechanics (QM/MM) implementation on AMD GPUs. This encompasses the entire Fock matrix build and force calculation in QUICK including one-electron integrals, two-electron repulsion integrals, exchange-correlation quadrature, and linear algebra operations. General performance improvements to the QUICK GPU code are also presented. Benchmarks carried out on NVIDIA V100 and AMD MI100 cards display similar performance on both hardware for standalone HF/DFT calculations with QUICK and QM/MM molecular dynamics simulations with QUICK/AMBER. Furthermore, with respect to the QUICK/AMBER release version 21, significant speedups are observed for QM/MM molecular dynamics simulations. This significantly increases the range of scientific problems that can be addressed with open-source QM/MM software on state-of-the-art computer hardware.