MiMiC: A high-performance framework for multiscale molecular dynamics simulations.

MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy [Kanungo, B. J. Phys. Chem. Lett. 2024, 15, 323-328], we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r2SCAN 50, a 50% global hybrid of exact exchange with r2SCAN, which may be less correct than r2SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r2SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree-Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r2SCAN@HF, as in the composite method HF-r2SCAN-DC4, makes its FE more negative than in r2SCAN@HF, permitting a near-perfect cancellation of FE and DE. r2SCAN by itself (and even more so, r2SCAN evaluated on the r2SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r2SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Finally, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn-Sham inversion of the CCSD(T) density.

Tucker Tensor Approach for Accelerating Fock Exchange Computations in a Real-Space Finite-Element Discretization of Generalized Kohn-Sham Density Functional Theory.

The evaluation of Fock exchange is often the computationally most expensive part of hybrid functional density functional theory calculations in a systematically improvable, complete basis. In this work, we employ a Tucker tensor based approach that substantially accelerates the evaluation of the action of Fock exchange by transforming three-dimensional convolutional integrals into a tensor product of one-dimensional convolution integrals. Our numerical implementation uses a parallelization strategy that balances the memory and communication bottlenecks, alongside overlapping compute and communication operations to enhance computational efficiency and parallel scalability. The accuracy and computational efficiency are demonstrated on various systems, including Pt clusters of various sizes and a TiO2 cluster with 3684 electrons.

Unconventional Error Cancellation Explains the Success of Hartree-Fock Density Functional Theory for Barrier Heights.

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, which has been known for ∼30 years. The conventional explanation is that the Hartree-Fock theory yields the more accurate density. This work presents a benchmark Kohn-Sham inversion of accurate coupled-cluster densities for the reaction H2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree-Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532-543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.

Exact and Model Exchange-Correlation Potentials for Open-Shell Systems.

The conventional approaches to the inverse density functional theory problem typically assume nondegeneracy of the Kohn-Sham (KS) eigenvalues, greatly hindering their use in open-shell systems. We present a generalization of the inverse density functional theory problem that can seamlessly admit degenerate KS eigenvalues. Additionally, we allow for fractional occupancy of the Kohn-Sham orbitals to also handle noninteracting ensemble-v-representable densities, as opposed to just noninteracting pure-v-representable densities. We present the exact exchange-correlation (XC) potentials for six open-shell systems─four atoms (Li, C, N, and O) and two molecules (CN and CH2)─using accurate ground-state densities from configuration interaction calculations. We compare these exact XC potentials with model XC potentials obtained using nonlocal (B3LYP, SCAN0) and local/semilocal (SCAN, PBE, PW92) XC functionals. Although the relative errors in the densities obtained from these DFT functionals are of O(10-3 to 10-2), the relative errors in the model XC potentials remain substantially large─O(10-1 to 100).

Exchange correlation potentials from full configuration interaction in a Slater orbital basis.

Ryabinkin-Kohut-Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. Examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

Efficient All-Electron Time-Dependent Density Functional Theory Calculations Using an Enriched Finite Element Basis.

We present an efficient and systematically convergent approach to all-electron real-time time-dependent density functional theory (TDDFT) calculations using a mixed basis, termed as enriched finite element (EFE) basis. The EFE basis augments the classical finite element basis (CFE) with a compactly supported numerical atom-centered basis, obtained from atomic ground-state DFT calculations. Particularly, we orthogonalize the enrichment functions with respect to the classical finite element basis to ensure good conditioning of the resultant basis. We employ the second-order Magnus propagator in conjunction with an adaptive Krylov subspace method for efficient time evolution of the Kohn-Sham orbitals. We rely on a priori error estimates to guide our choice of an adaptive finite element mesh as well as the time step to be used in the TDDFT calculations. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretizations. Notably, we attain a 50-100 times speedup for the EFE basis over the CFE basis. We also demonstrate the efficacy of the EFE basis for both linear and nonlinear responses by studying the absorption spectra in sodium clusters, the linear to nonlinear response transition in the green fluorescence protein chromophore, and the higher harmonic generation in the magnesium dimer. Lastly, we attain good parallel scalability of our numerical implementation of the EFE basis for up to ∼1000 processors, using a benchmark system of a 50-atom sodium nanocluster.

A Comparison of Exact and Model Exchange-Correlation Potentials for Molecules.

Accurate exchange-correlation (XC) potentials for three-dimensional systems─via solution of the inverse density functional theory (DFT) problem─are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules─dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne─are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10-3-10-2), the model XC potentials have much higher errors of O(10-1-100). Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

Backbone charge transport in double-stranded DNA.

Understanding charge transport in DNA molecules is a long-standing problem of fundamental importance across disciplines1,2. It is also of great technological interest due to DNA's ability to form versatile and complex programmable structures. Charge transport in DNA-based junctions has been reported using a wide variety of set-ups2-4, but experiments so far have yielded seemingly contradictory results that range from insulating5-8 or semiconducting9,10 to metallic-like behaviour11. As a result, the intrinsic charge transport mechanism in molecular junction set-ups is not well understood, which is mainly due to the lack of techniques to form reproducible and stable contacts with individual long DNA molecules. Here we report charge-transport measurements through single 30-nm-long double-stranded DNA (dsDNA) molecules with an experimental set-up that enables us to address individual molecules repeatedly and to measure the current-voltage characteristics from 5 K up to room temperature. Strikingly, we observed very high currents of tens of nanoamperes, which flowed through both homogeneous and non-homogeneous base-pair sequences. The currents are fairly temperature independent in the range 5-60 K and show a power-law decrease with temperature above 60 K, which is reminiscent of charge transport in organic crystals. Moreover, we show that the presence of even a single discontinuity ('nick') in both strands that compose the dsDNA leads to complete suppression of the current, which suggests that the backbones mediate the long-distance conduction in dsDNA, contrary to the common wisdom in DNA electronics2-4.

Determining plasmonic hot-carrier energy distributions via single-molecule transport measurements.

Hot carriers in plasmonic nanostructures, generated via plasmon decay, play key roles in applications such as photocatalysis and in photodetectors that circumvent bandgap limitations. However, direct experimental quantification of steady-state energy distributions of hot carriers in nanostructures has so far been lacking. We present transport measurements from single-molecule junctions, created by trapping suitably chosen single molecules between an ultrathin gold film supporting surface plasmon polaritons and a scanning probe tip, that can provide quantification of plasmonic hot-carrier distributions. Our results show that Landau damping is the dominant physical mechanism of hot-carrier generation in nanoscale systems with strong confinement. The technique developed in this work will enable quantification of plasmonic hot-carrier distributions in nanophotonic and plasmonic devices.

Exact exchange-correlation potentials from ground-state electron densities.

The quest for accurate exchange-correlation functionals has long remained a grand challenge in density functional theory (DFT), as it describes the many-electron quantum mechanical behavior through a computationally tractable quantity-the electron density-without resorting to multi-electron wave functions. The inverse DFT problem of mapping the ground-state density to its exchange-correlation potential is instrumental in aiding functional development in DFT. However, the lack of an accurate and systematically convergent approach has left the problem unresolved, heretofore. This work presents a numerically robust and accurate scheme to evaluate the exact exchange-correlation potentials from correlated ab-initio densities. We cast the inverse DFT problem as a constrained optimization problem and employ a finite-element basis-a systematically convergent and complete basis-to discretize the problem. We demonstrate the accuracy and efficacy of our approach for both weakly and strongly correlated molecular systems, including up to 58 electrons, showing relevance to realistic polyatomic molecules.

Effect of length and contact chemistry on the electronic structure and thermoelectric properties of molecular junctions.

We present a combined experimental and computational study that probes the thermoelectric and electrical transport properties of molecular junctions. Experiments were performed on junctions created by trapping aromatic molecules between gold electrodes. The end groups (-SH, -NC) of the aromatic molecules were systematically varied to study the effect of contact coupling strength and contact chemistry. When the coupling of the molecule with one of the electrodes was reduced by switching the terminal chemistry from -SH to -H, the electrical conductance of molecular junctions decreased by an order of magnitude, whereas the thermopower varied by only a few percent. This has been predicted computationally in the past and is experimentally demonstrated for the first time. Further, our experiments and computational modeling indicate the prospect of tuning thermoelectric properties at the molecular scale. In particular, the thiol-terminated aromatic molecular junctions revealed a positive thermopower that increased linearly with length. This positive thermopower is associated with charge transport primarily through the highest occupied molecular orbital, as shown by our computational results. In contrast, a negative thermopower was observed for a corresponding molecular junction terminated by an isocyanide group due to charge transport primarily through the lowest unoccupied molecular orbital.

Role of macroscopic deformations in energetics of vacancies in aluminum.

Electronic structure calculations on million-atom samples are employed to investigate the effect of macroscopic deformations on energetics of vacancies in aluminum. We find that volumetric strain associated with a deformation largely governs the formation energies of monovacancies and divacancies. The calculations suggest that nucleation of these defects is increasingly favorable under volumetric expansion, so much to the point that they become thermodynamically favorable under large positive volumetric strains. On the contrary, on an average, existing vacancies are found to bind preferentially under compressive volumetric strains. Shear deformations did not affect the formation energies of vacancies, but strongly influenced the 110 divacancy binding energies, causing them to orient under energetically preferential directions.