2026-04-17 | | Total: 5
We present a relativistic time-dependent density functional theory (TDDFT) approach for the simulation of resonant inelastic X-ray scattering (RIXS) spectra, based on both a full four-component (4c) Dirac-Coulomb Hamiltonian and a modern atomic mean-field exact two-component (amfX2C) Hamiltonian model. The approach builds on the pseudo-wavefunction formalism and a core-valence separation scheme, enabling the efficient evaluation of couplings between two manifolds of excited states relative to a common ground state, as required for solving the Kramers-Heisenberg equation for RIXS. The relativistic formulation provides a variational description of scalar and spin-orbit relativistic effects, which are essential for accurately describing inner-shell excitations involved in RIXS processes. Its transformation to the 2c regime via the amfX2C Hamiltonian significantly reduces the computational cost while offering 4c-quality results by accounting for two-electron and exchange-correlation picture-change effects arising from the X2C transformation. In addition to two-dimensional RIXS maps, the methodology enables the direct evaluation of high-energy-resolution fluorescence detection (HERFD) and resonant X-ray emission spectra (RXES). Applications to 2p3d and 3d4f RIXS maps of selected ruthenium and uranium complexes demonstrate that the amfX2C approach reproduces reference 4c results and experimental spectra with high accuracy, capturing all key spectral features and providing reliable peak assignments.
The pair-coupled-cluster doubles (pCCD) method has emerged as a viable approach for quantum-chemical studies of strongly correlated systems. Despite its lower formal scaling (O(N$^4$)) compared to other versions of coupled cluster (CC) theory, applications to large chemical structures are still expensive. Fragmentation and embedding strategies offer a viable approach in such cases. In this work, we present a simple and efficient density-embedding scheme based on pCCD electron densities. The main computational benefit arises from the fact that pCCD response $Λ$-equations are much cheaper to compute than those of standard CC methods, providing easy access to one-electron properties. The pCCD densities of the individual subsystems are used to generate static embedding potentials that capture the environment's effect on the embedded system. The individual fragment energies are then iteratively converged in a self-consistent fashion. We demonstrate the reliable performance of this scheme with the estimation of dipole moments of the weakly bound CO2$\cdots$Rg (Rg = He, Ne, Ar, and Kr) complexes and with the modeling of vertical excitations of some microsolvated molecules.
The origin of enhanced reactivity in aqueous microdroplets remains debated, with interfacial electric fields (IEFs) often invoked as catalytic drivers. Here, we provide a quantum-mechanical, spatially resolved characterization of the electric field at air-water interfaces by combining deep-learning molecular dynamics with \emph{ab initio} re-sampling. Across planar interfaces and nanodroplets of varying curvature and charge state, we find an outward-oriented field of $\sim 1.0$--$1.2$ V/Å along the intrinsic surface normal. Crucially, its magnitude scales linearly with the average number of hydrogen bonds per interfacial molecule, directly tying the field to the local hydrogen-bond network. Despite its large magnitude and contrary to common expectations, we find that curvature and pH exert only a minor influence on the IEF, becoming negligible at experimentally relevant droplet sizes and pH. Consequently, the reactivity differences observed in $μ$m-sized droplets cannot be ascribed to variations in the IEF, which changes by a factor of only $\sim10^{-5}$ between $3$ and $40μ$m-sized droplets. Moreover, the IEF is localized inside the interfacial region and rapidly vanishes within a few Å. This strong spatial confinement renders the IEF strongly tied to the local electronic structure, identifying it as a local property of the air-water boundary rather than an independent physical driver of ``on-water'' catalysis.
Modelling micro- and meso-scopic scale thermodynamic and transport properties of soft condensed matter hinges upon its representation. This is especially relevant for polar solvents such as water, since these require effective representation of their dielectric nature as driven by molecular charge distributions and molecular network structuring. The dielectric nature of a medium leads to complex phenomena such as local polarisability response and restructuring near interfaces in reaction to changes in local charge distributions. Inclusion of such phenomena when using larger-than-atomistic techniques such as coarse-grained molecular dynamics (CG-MD) and dissipative particle dynamics (DPD) is still an open question, to which we provide a novel way to consider and justify the necessary and suitable coarse-graining level, enabling us to compare new polar CG models' performance against that of an underlying atomistic model. We polarise our previous non-polar nDPD water model to prepare it for use in simulations of liquid electrolytes as well as solvated organic membranes and measure its fitness to serve as a dielectric medium by comparing its properties to those of the TIP3P water model, while simultaneously observing changes to properties already represented well by the non-polar model.
We address a fundamental question: under which conditions do the dynamics and thermodynamics of open chemical reaction networks (CRNs), grounded on the notion of idealized chemostats that exchange selected species, emerge from underlying closed CRNs? While open CRNs provide the standard framework to describe out-of-equilibrium chemical systems, real systems are finite and ultimately relax to equilibrium, leaving the status of this description conceptually unresolved. Here we show that open-CRN behavior arises as an asymptotic regime of closed CRNs when two minimal and physically transparent conditions are met: a time-scale separation, whereby fast reactions effectively act as exchange mechanisms, and an abundance separation, whereby a subset of species behaves as chemostats with diverging chemical capacity. In this regime, both the stochastic dynamics and the thermodynamic structure \ -- including local detailed balance, entropy production, and free-energy balance \ -- emerge to leading order from the underlying closed CRN. Our results apply to arbitrary stoichiometries. They show that chemostats need not be introduced as external idealizations, but instead arise as emergent thermodynamic structures within closed systems, providing a unified and physically grounded foundation for the nonequilibrium thermodynamics of CRNs.