2026-03-10 | | Total: 15
Hyperpolarisation using Signal Amplification By Reversible Exchange (SABRE) is a convenient method for high repeatability studies. The core of this technique is polarisation transfer to the target substrate during an on-going chemical exchange process. Typically, polarisation transfer is achieved as fast as possible. In this study we employ NMR sequences that on contrary slow down the polarisation transfer and yet demonstrate improved performance. Simulations confirm that such methods can lead to high polarisation yield in SABRE system that exhibit higher magnetic inequivalence and lower chemical exchange rate.
Excited atoms and molecules can utilize their excess energy to ionize a neighboring system by a process named interatomic and intermolecular Coulombic decay (ICD). ICD is ultrafast, in the femtosecond regime, and has many modes of appearance. Ample applications of ICD have been reported spanning a wide range of fields and it is expected to be ubiquitous in nature. Essentially all the investigations on ICD were for weakly bound systems, like clusters and fluids. We demonstrate that, unexpectedly, ICD can be efficiently active in atomic and molecular gases in spite of the very large distances between the units. We uncover the underlying mechanism, which differs from that prevailing in weakly bound systems. The dynamics of ICD in gases is analyzed. The results considerably broaden the impact of ICD and open the gateway to new kinds of applications.
We define and calculate the triplet-pair population of the $2^1A_g$ dark state in polyenes, as predicted by the $π$-electron Pariser-Parr-Pople (PPP) model, for chains of 8 to 14 C-atoms and Coulomb interaction parameter between 4-14 eV. Our definition of the triplet-pair population is motivated by a two-particle model of the $2^1A_g$ state. We use DMRG to solve the PPP model and we exploit the MPS representation of the DMRG wavefunction to compute the triplet-pair population. Using our results for short chain sizes, we predict a finite-size scaling value of the triplet-pair population of ca. 75% for realistic Coulomb interactions for polyene chains. Our results agree with other theoretical work on the doubly-excited character of polyenes, and represents further evidence that the $2^1A_g$ state is predominantly triplet-pair in character - with implications for singlet fission mechanisms in polyenes.
Machine Learning Interatomic Potentials (MLIPs) enable accurate large-scale atomistic simulations, yet improving their expressive capacity efficiently remains challenging. Here we systematically develop Mixture-of-Experts (MoE) and Mixture-of-Linear-Experts (MoLE) architectures for MLIPs and analyze the effects of routing strategies and expert designs. We show that sparse activation combined with shared experts yields substantial performance gains, and that nonlinear MoE formulations outperform MoLE when shared experts are present, underscoring the importance of nonlinear expert specialization. Furthermore, element-wise routing consistently surpasses configuration-level routing, while global MoE routing often leads to numerical instability. The resulting element-wise MoE model achieves state-of-the-art accuracy across the OMol25, OMat24, and OC20M benchmarks. Analysis of routing patterns reveals chemically interpretable expert specialization aligned with periodic-table trends, indicating that the model effectively captures element-specific chemical characteristics for precise interatomic modeling.
We compute highly accurate first principle based \textit{ab initio} adiabatic potential energy surfaces (PESs) using State-Averaged Multi-Configurational Self-Consistent Field (SA-MCSCF) followed by internally contracted Multi-Reference Configuration Interaction method incorporating fixed-reference Davidson corrections [ic-MRCI(Q)], where a full valence active space of 18 electrons in 12 orbitals and aug-cc-pVQZ basis set are employed for the low-lying four singlet electronic states of ozone ($\tilde{X}^1A'$, $1~^1A''$, $1~^1A'$ and $2^1A''$). It accurately reproduces the dissociation energies of ozone (1.101 eV) as well as the molecular oxygen (5.106 eV) along with vibrational frequencies of O$_3$ in comparison with experimental data. To ensure appropriate accuracy and proper convergence in the interaction as well as asymptotic regions, we (a) extend the number of electronic states in SA-MCSCF calculation (singlet as well as triplet and quintet); (b) systematically expand the active space [(12e,9o) $\rightarrow$ (18e,12o) $\rightarrow$ (24e,15o)] and basis set size (AVDZ $\rightarrow$ AV6Z $\rightarrow$ Complete Basis Set limit); (c) incorporate multi-reference character along with Davidson correction. Conical intersections between the adjacent electronic states (1-2, 2-3 and 3-4) are located at \textit{C}$_{2v}$, \textit{D}$_{3h}$ as well as \textit{C}$_{s}$ geometries through the four-state adiabatic-to-diabatic transformation of non-adiabatic coupling terms (NACTs) computed at Coupled-Perturbed Multi-Configurational Self-Consistent Field (CP-MCSCF) method along the circular contours. Finally, we present: (a) ic-MRCI(Q) calculated minimum energy path of incoming oxygen to the diatom (O$_2$) is devoid of any ``reef'' feature; (b) NACTs and diabatic PES matrix elements as function of hyperangles ($θ$,$φ$) at a fixed hyperradius $ρ= 4$ Bohr for a four state sub-Hilbert space.
Machine learning is rapidly accelerating materials and chemical discovery, but most current models target energies, forces, or selected molecular properties rather than the underlying many-body electronic structure. Learning electronic-structure proxies, such as reduced density matrices, offers a path to surrogates that can predict a broad range of observables from a single ML model. Short of learning the full wavefunction, the two-electron reduced density matrix (2-RDM) is among the most information-rich, minimally lossy targets, providing direct access to expectation values of arbitrary one- and two-electron operators regardless of the strength of the underlying electron correlation. Here we show that learning the 2-RDM is a feasible goal, yielding exceptionally accurate models. We develop surrogates for correlated wavefunction methods (including configuration interaction and coupled cluster) that yield 2-RDMs with sufficient fidelity to provide direct, training-free access to energies and forces for driving energy-conserving molecular dynamics. To tackle realistic molecular condensed phases, we leverage a many-body expansion of the 2-RDM, using our ML models to supply the expansion terms and enabling ML-powered, coupled-cluster-quality electronic structure and energetics for large solvated systems. As a demonstration, we showcase a coupled-cluster-level electronic-structure calculation of glucose solvated by 500 water molecules achieved at Hartree-Fock cost. This work establishes a general framework for learning correlated electronic structure with high fidelity and deploying it to systems beyond the reach of conventional ab initio methods.
Collective light-matter systems host an extensive manifold of dark states whose role in the emergence of thermodynamic behavior remains poorly understood, especially in the presence of disorder and structured environments. Here, we develop a hybrid matrix product state-hierarchical equations of motion (MPS-HEOM) approach that enables numerically exact simulations of polariton dynamics from a few emitters to the thermodynamic limit under both static and dynamic disorder. This allows us, for the first time, to provide a quantitative and operational answer to the long-standing question of what is the minimum system size required to reach the thermodynamic limit in collective polaritonic systems. By introducing a convergence scale, $N_{T}$, i.e., the number of molecules required for the photonic dynamics to reach the thermodynamic limit, we show that dynamic disorder generally poses a greater computational challenge than static disorder. We attribute this behavior to the suppression of collective light-matter dynamics by disorder, which dynamically activates non-collective degrees of freedom. We further find that $N_{T}$ exhibits a turnover behavior as the bath becomes more Markovian, as the bath timescales regulate bright-to-dark energy transfer and the involvement of dark and gray states. Hence, phonon timescales control both the breakdown of collective behavior and the growth of $N_{T}$. Our results establish the suppression of collective behavior as the key mechanism governing thermodynamic convergence in disordered light-matter systems.
Quantitative prediction of thermodynamic properties in solution is essential for translating atomistic simulations into reliable chemical insight. As an exemplar system, the behaviour of CaCO$_3$ in water has been widely studied to understand its mineralization in seawater, with potential implications for carbon-capture strategies. However, making accurate computational predictions has been a long-standing challenge, requiring both highly accurate electronic structure methods and extensive statistical sampling. Here, we combine advances in machine learning and electronic structure theory to fully resolve the ion pairing free energy of CaCO$_3$ with explicit solvation. We show that achieving quantitative agreement with experiment requires going beyond the standard density functional theory up to the "gold-standard" coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We generate a set of systematically improvable models, enabling reliable insights into the initial association mechanism of Ca and CO$_3$ ions prior to nucleation while fully quantifying enthalpic and entropic effects. Our results demonstrate that CCSD(T)-level thermodynamic predictions of complex aqueous systems can now be routinely achieved.
Rare nonadiabatic events play a central role in photochemistry but remain difficult to simulate because excited-state dynamics is computationally demanding and often stochastic. Here we introduce a deterministic and time-reversible implementation of nonadiabatic dynamics that enables the application of transition path sampling (TPS) to excited-state processes. Our approach builds on the Mapping Approach to Surface Hopping (MASH) and establishes the conditions required for path ensemble sampling, in particular time reversibility and detailed balance. Combining this dynamics with the TPS framework yields a new method, termed nonadiabatic transition path sampling (NATPS). Using a model system of electronically coupled potential energy surfaces, we demonstrate that NATPS efficiently generates ensembles of reactive trajectories and provides mechanistic insight into nonadiabatic pathways. Compared with brute-force trajectory simulations and forward-flux sampling approaches, NATPS substantially reduces the computational effort required to obtain reactive trajectories.
We develop a multi-reference perturbation theory for electronic structure calculations based on symmetries of the Hamiltonian. The reference Hamiltonian in the symmetry-based perturbation theory (SBPT) is chosen such that it possesses more symmetries than the original Hamiltonian, leading to a larger reduction of computational resources in terms of both the number of configurations in the configuration interaction expansion and the number of required qubits in quantum computing applications. We provide approximate, scalable solutions for the second-order correction, as well as an application to selected configuration interaction. We show that SBPT is an extension of other existing multi-reference perturbation theories and that it can give better results for some molecular systems in a robust way.
Trions -- Coulomb-bound three-particle excitations composed of two like-charge carriers and one oppositely charged carrier -- are central quasiparticles in two-dimensional semiconductors. Reduced dielectric screening and quantum confinement strongly enhance their binding energies, making them robust and experimentally accessible. This review surveys theoretical and experimental advances in trion physics, emphasizing rigorous few-body approaches and the role of dielectric environment, anisotropy, and external electric and magnetic fields. We analyze computational methods for describing trions in two-dimensional configuration spaces and discuss how reduced dimensionality modifies their structure and stability. Connections to many-body phenomena, including screening, Landau-level mixing, and exciton--polaron crossover, are also highlighted.
Quantum algorithms for electronic-structure simulations are actively being developed, yet many hybrid quantum-classical approaches are bottlenecked by the measurement overhead associated with large molecular Hamiltonians. Here we introduce a diagonalization-driven framework that progressively drives the electronic Hamiltonian toward a (block-)diagonal form in the Slater-determinant basis using sequential Givens rotations. In contrast to Schrödinger-picture methods that variationally optimize a wave function, our approach adopts a Heisenberg-picture viewpoint: the Hamiltonian is iteratively transformed, and rotation angles are determined classically from low-dimensional effective blocks, reducing the quantum workload to a small, fixed set of matrix-element measurements per iteration. Candidate generators are estimated via approximate Baker-Campbell-Hausdorff updates with truncation and cumulant-based approximations that control Hamiltonian growth, complemented by stochastic selection to avoid stagnation. We further introduce an angle-merging procedure that reduces circuit depth by consolidating repeated small-angle rotations. We benchmark the framework on N$_2$ and strongly correlated hydrogen systems, assessing convergence behavior, residual-structure diagnostics, measurement-accuracy trade-offs, circuit costs, and robustness under finite sampling.
Halide segregation in metal halide perovskites limits their bandgap tunability and hinders their adoption in tandem solar cells and light emitting diodes. Here, we reveal the thermodynamic driving force behind halide segregation in mixed halide (Br-I) perovskites. By performing first-principles calculations on slab models with varying bromide and iodide distributions, we demonstrate that bromide ions preferentially occupy surface sites over bulk sites. Our simulations show that the segregation tendency is higher in MAPb(Br$_x$I$_{1-x}$)$_3$ (MA=methylammonium) compared to FA$_{0.8}$Cs$_{0.2}$Pb(Br$_x$I$_{1-x}$)$_3$, highlighting the role of the A-site cation. To quantify this effect, we establish a descriptor for halide segregation: the difference in defect formation energies of Br antisite defects between the bulk and the surface, which confirms the role of the A-site cation at equimolar Br-I concentration. Furthermore, we identify the localization of photo-generated holes near iodide ions, which triggers their oxidation and accelerates the formation of iodide vacancies, thereby promoting segregation. Overall, this work establishes defect thermodynamics as a framework for understanding halide segregation and provides a structural basis for designing stable mixed halide perovskites.
An efficient method of computing the dynamics of a pure quantum state under the time-dependent Schrödinger equation is described: At each timestep, a restricted subspace of the potentially infinite-dimensional total Hilbert space is systematically and naturally constructed via the image of repeated applications of the Hamiltonian operator, and the time evolution is computed exactly within said subspace. The subspace is dynamically recomputed at each timestep such that it co-evolves with the state vector. We benchmark the method using the Holstein model and compare the formal information content of its representation to the matrix-product state formalism. The method is built from the ground up as a parallel algorithm for graphics processing units and is applicable to arbitrary Hamiltonians that are sparse in a given basis. It can be extended to open quantum system dynamics and/or time-dependent generators.
Quantum coherence and other non-classical features are widely discussed in chemical dynamics, yet it remains difficult to quantify when such resources are operationally relevant for a given process and observable. While quantum resource theories provide a comprehensive framework for comparing free and resourceful settings, existing approaches typically rely on resource monotones or on performance bounds under free operations, and do not directly quantify the maximal influence a chosen resource can exert on a fixed chemical dynamics. Here, we introduce task specific, process level quantifiers that upper bound the largest change a quantum resource can induce in a target figure of merit. Central is a resource impact functional $\mathcal{C}_M(Λ)$, defined by comparing a state with its paired resource-free counterpart under the same quantum channel $Λ$, which admits an operational interpretation in binary hypothesis testing. We derive variation and time bounds that constrain how rapidly a resource can modify a target signal, providing resource-aware analogues of quantum speed limits. Moreover, we show that open system dynamics can be decomposed into free and resourceful components such that only the resourceful component contributes to $\mathcal{C}_M(Λ)$, thereby isolating the parts of a generator responsible for resource-induced changes in the observable. We illustrate the framework exemplary for energy transfer in a donor-acceptor dimer in two analytically solvable regimes. Our results provide a general toolbox for diagnosing and benchmarking quantum resource effects in molecular processes.