2026-05-28 | | Total: 9
Liquid water is fundamentally important, and its accurate computer simulation has been the driving force for myriad methodological developments. Ab initio molecular dynamics with forces obtained from density functional theory (DFT) is now a standard tool widely used by researchers. However, we reveal that previous studies of liquid water using the same widely-used density functional (revPBE-D3) exhibit significant discrepancies with one another, varying by over 20% in the diffusion coefficient and 10% in the density, raising fundamental questions about reproducibility. By combining modern long-range machine-learning interatomic potentials that enable robust statistical sampling with carefully converged DFT training data, we resolve these discrepancies, achieving consensus across six diverse community codes. Our predictions differ markedly from previous literature: we show that most previous results overestimate the density and underestimate the diffusion coefficient of revPBE-D3 water due to basis set incompleteness and pseudopotential inconsistencies, coupled with limitations in statistical sampling (in some cases). These benchmark values provide a reliable reference for validating current and future implementations of DFT-based ab initio molecular dynamics. Reaching agreement establishes confidence and credibility and serves as a prerequisite for the systematic assessment of new density functionals and numerical approximations.
Block tensor decomposition (BTD) and canonical polyadic decomposition (CPD) are combined into a unified $O(N^3)$-scaling framework for second-order perturbation theory (PT2), demonstrated on MP2 and renormalized PT2 (rPT2). BTD constructs the tensor hyper-contraction kernel at $O(N^3)$ via a dual-grid scheme; CPD factorizes the exchange channel through a block-based two-stage ALS. An asymmetric half-kernel design applies bare Coulomb to one vertex and coupling-constant-averaged screening to the other, capturing the SOSEX component of rPT2 without a frequency-dependent CPD. For MP2, BTD-CPD reproduces canonical RI-MP2 to 0.058~kcal/mol per heavy atom. For rPT2@PBE0 on the S66x8 benchmark, the mean absolute error is 0.36~kcal/mol (ME $-$0.19, RMSE 0.46) over 528 data points. The CPD-compressed intermediates yield $O(N^2)$ storage alongside $O(N^3)$ scaling.
Variational excited-state density functional theory (DFT) enables the calculation of excited states at a cost comparable to ground-state calculations, but single-configuration approaches often suffer from spin contamination. We implement restricted open-shell Kohn-Sham (ROKS) DFT, which recovers spin-pure singlet excitation energies via the variational minimization of a weighted combination of mixed-spin and triplet configurations, within the plane-wave projector augmented-wave framework of VASP. The energy functional is optimized using a preconditioned conjugate-gradient or a direct inversion in the iterative subspace algorithm, and analytical atomic forces are derived. The implementation is validated for eight organic molecules by comparison to the Q-Chem quantum chemistry code, yielding mean deviations of approximately $30$ eV. As a solid-state application, we investigate the three lowest lying excitations of MgO with a neutral oxygen vacancy. For a dielectric-dependent hybrid functional, vertical excitation energies from ROKS and time-dependent density functional theory (TDDFT) differ on average by about $0.21$ eV. The Franck-Condon shifts deviate on average by $0.14$ eV between the two methods and mass-weighted displacements between the excited states and the ground state by $0.12$ amu$^{1/2}$ Ang. Additional calculations at the PBE level reveal that these properties depend less strongly on the DFT functional for ROKS than for TDDFT. These results demonstrate that ROKS provides excitation energies and excited-state forces with an accuracy similar to TDDFT while retaining the favorable scaling of ground-state DFT, making it a promising approach for affordable excited-state simulations in extended systems.
Enforcing universal symmetries in machine learning (ML) models is a common strategy to mitigate data scarcity. We show that exploiting exact, as well as approximate, label symmetries can benefit scaling laws. We illustrate the idea for the s, p, d orbital densities of the electron in the hydrogen atom, for the three vibrational normal modes of the water molecule, as well as its full 3D potential energy hypersurface. Resulting ML models of electron density and potential energies exhibit superior learning curves, demonstrating improved generalization efficiency. When label symmetries are not exact, the same principles govern the observed learning behavior -- up to the convergence floors set by the degree to which the symmetry is approximate. For convex wells in the molecular potential energy surface, a Hessian-based correction suppresses the leading symmetry-breaking error in augmented labels.
Hot-exciton relaxation in semiconductor nanocrystals (NCs) is often described using perturbative theories, but their accuracy is difficult to assess for realistic exciton--phonon Hamiltonians. Here, we benchmark the perturbative quantum master equation (QME) and several mixed quantum--classical (MQC) methods against fully quantum mechanical dynamics. Using atomistically parameterized models for CdSe core and CdSe/CdS core--shell NCs, we find that bare CdSe exhibits an ultrafast initial decay followed by slower cooling, whereas the core--shell system is dominated by the slower component. Analysis of reduced models shows that the ultrafast component arises from rapid diabatic state mixing driven by thermal fluctuations of low-frequency phonons, rather than from nuclear-assisted energy relaxation. The QME captures the initial fast decay but can fail for the slower relaxation in the diabatic representation, while the mapping approach to surface hopping (MASH) gives the most consistent agreement with both benchmark dynamics and equilibrium populations. These results establish a benchmark for exciton-cooling dynamics in NCs and clarify the physical regimes in which widely used approximate methods are reliable.
Nonlinear optical light sources enable the generation of photons with quantum states that are intrinsically linked to underlying material dynamics, rather than imposed through external modulation. Here we investigate fundamental quantum communication metrics of four-wave-mixing signal fields generated by the two-dimensional perovskite (PEA)2PbI4. Using polarization-resolved interferometric measurements together with a microscopic nonlinear response model for the Bloch vector, we reconstruct effective single-photon polarization density matrices inferred from the experimental signal fields and evaluate the corresponding Holevo bound and effective secret-bit rates as a function of the coherence time, population time, and detection wavelength. We find that incorporating the coherence-time degree of freedom systematically lowers the Holevo bound by approximately 2.6-5.8% across the various excitonic resonances, indicating reduced distinguishability of the polarization states when the full multidimensional parameter space is sampled. To connect the polarization-state indistinguishability with experimentally achievable throughput, we further introduce an effective secret-bits-per-pulse metric that enables rapid evaluation of secure information throughput for candidate materials without requiring photon-number-resolved detection. For the present system, control of the population time via spin-dependent evolution yields substantially higher secret-bit rates than manipulation of the coherence time, while spectral regions associated with single-exciton and biexciton resonances define complementary operating regimes for secure communication. More broadly, this work positions nonlinear spectroscopy as a framework for exploring how emergent optical materials can generate and structure quantum states in ways that are advantageous for established quantum communication schemes.
Designing compact and accurate circuits for the variational quantum eigensolver (VQE) is a central challenge in near-term quantum chemistry. Existing adaptive methods such as ADAPT-VQE design circuits by iteratively selecting operators from a predefined pool guided by gradient information and greedy heuristics. In this work, we adopt differentiable quantum architecture search (DQAS) as a circuit design framework based on the UCCSD operator pool, and introduce two complementary strategies: a global mode that simultaneously optimizes all operator selections, and a layerwise mode that constructs circuits incrementally while preserving previously learned structure. By relaxing discrete operator selection into a continuous differentiable optimization, DQAS enables gradient-based exploration over the combinatorial space of UCC circuit architectures. Benchmarks on BeH2, H4, LiH, H6, and H2O (8-14 qubits) show that both strategies achieve higher accuracy and fewer CNOT gates than ADAPT-VQE in the compact circuit regime, with up to 2.7-fold accuracy improvement for H2O and CNOT reductions of 13-17% at equivalent circuit depths. Benchmarks on the qubit-excitation-based (QEB) operator pool confirm that both advantages generalize beyond UCCSD. These results demonstrate that differentiable architecture search provides an effective and generalizable framework for designing accurate and compact VQE circuits in near-term quantum chemistry.
Strong-field enhanced ionization (EI) is a phenomenon in which stretching of interatomic bonds into a distorted molecular geometry leads to an increase in the tunneling ionization rate driven by a strong field. Isolating the momentum distribution of the electrons involved in EI is critical to fully characterizing the phenomenon. We have measured this EI distribution in triple ionization of D$_2$O using 6-fs pulse pairs together with full fragment momentum imaging and electron-ion correlation methods. We find that the EI electron momentum distribution differs substantially from that of standard strong-field tunneling from molecules, exhibiting an increased yield of electrons with large momentum in the direction of the laser polarization, and a change from the expected Gaussian distribution. These observations indicate that the instantaneous EI tunneling rate is maximized at a critical value of the laser electric field, rather than at the peak of an optical cycle. This finding distinguishes EI from Keldysh tunneling rate predictions, where tunneling rate increases monotonically with field strength. These pronounced differences between EI and non-EI electron spectra are critical tests of models of enhanced ionization and suggest a route towards control of the sub-cycle timing of electron emission.
We investigate how proximity to a metal nanostructure, particularly to a flat mirror or a cavity confined between two mirrors, affects the vibronic structure of Raman scattering signals. We find that such proximity, particularly for the strong-coupling situation encountered in cavity environments, plays multiple roles in shaping Raman signals beyond the now-familiar signal enhancement known as surface-enhanced Raman scattering (SERS). First, in analogy to the electromagnetic SERS mechanism, near or between mirrors, the local field experienced by a molecule differ from that in vacuum. In particular, between mirrors, the cavity enhances the effective excited state population by trapping the EM field inside it. Second, the nearby metal surface provides a relaxation channel and a lineshape broadening mechanism, and inside a cavity this lineshape is inherited by the cavity polaritons. This relaxation results in a loss of yield but the associated broadening also leads to significant absorption over a larger frequency range. Third, near metallic interfaces interference between incident and reflected light can lead to a richly structured Raman spectrum. For instance, we find that the Rabi contraction (that results from depopulating the ground state) can interfere with Raman signals (and the effect appears to be the same order as Raman itself). These cavity effects are calculated by a full-scale FDTD simulation and highlight the convoluted but fascinating roles of photonic materials on optical signals.