Atomic structure, cold atoms, and molecular spectroscopy
The intermolecular interaction between a water molecule and the electrons in aromatic π systems--the water-π bond--lies at the heart of many chemical processes, yet its properties remain challenging to measure experimentally and model computationally. Infrared spectroscopy of pyrene anions hydrated by a single water molecule reveals vibrational and electronic motions that are often hidden in condensed phase measurements. Results from new machine-learning approaches to potentials and dipole moments show that the electron dynamics of the aromatic π cloud quench signals from some of water's vibrations and amplify others. The observed interplay between electronic and vibrational motions has general implications for modeling intermolecular interactions between water and aromatic systems in clusters, solutions, and at interfaces.
The ultrafast nonlinear optical response of molecular ensembles is fundamentally altered under strong light-matter coupling. To rigorously isolate the genuine many-body contributions, an exact time-domain field-subtraction protocol is developed within a fully non-perturbative Maxwell-Liouville framework explicitly incorporating the two-exciton manifold in real space and time. This approach reveals that while collective cavity delocalization drives the macroscopic nonlinear signal toward a severe harmonic cancellation (an effect termed "spectral starvation"), intrinsic many-body molecular interactions robustly resurrect genuine polaritonic double-quantum coherences (DQCs). This many-body resurrection is governed by a universal two-photon matching rule, $Δ_B + 4J = Ω_R$, linking molecular anharmonicity ($Δ_B$) to the macroscopic Rabi splitting ($Ω_R$) and excitonic coupling ($J$). Crucially, this dictates that J-aggregates ($J < 0$) uniquely isolate the resonant many-body state below the dense two-exciton scattering continuum, protecting the macroscopic coherence from spatial fragmentation. This predictive framework establishes a direct phase diagram to engineer and protect optical nonlinearities across diverse strongly coupled platforms.
We present a molecular extension of our recently proposed Green's function embedding method, interacting-bath dynamical embedding theory (ibDET), for computing charged excitation energies at the $GW$ and EOM-CCSD levels. Starting from atom-centered impurities, we construct bath representations that capture the frequency-dependent entanglement between the impurity and its environment and can be systematically improved via the construction of cluster-specific natural orbitals. Utilizing a $GW$ or coupled-cluster Green's function solver, the self-energy of the full system is assembled from all embedding problems to obtain the interacting Green's function. We show that ibDET provides accurate spectral properties with much reduced cost for a broad range of systems, including conjugated molecules and nanoclusters. Compared with full-system results, the errors in the predicted ionization potentials and electron affinities are around 0.1 eV or smaller, while each embedding problem includes only a small fraction of the total orbital space. This work provides an efficient and scalable framework for computing spectral properties of molecular systems.
Machine learning force field (MLFF) has emerged as a powerful data-driven tool for atomistic simulations, enabling large-scale and complex atomic systems to be simulated with accuracy comparable to \textit{ab initio} methods. However, MLFFs often suffer from low training efficiency in the phase transition regime, where structural fluctuations are significantly elevated. To address this challenge, we propose a Central-Peripheral Distillation (CPD) algorithm for training dataset distillation. By strategically integrating representative samples with critical corner cases, the CPD algorithm ensures that the distilled dataset retains maximum structural diversity. We validated the efficacy of the CPD method on the liquid-liquid phase transition of dense hydrogen. Results show that, with the CPD approach, only 200 configurations are sufficient to train a MLFF that can fully reproduce the structural and dynamical properties of liquid hydrogen in the vicinity of its phase transition regime. This work paves the way for high-fidelity labeling of the MLFF training datasets, for instance by adopting high-level \textit{ab initio} calculations beyond the standard density functional theory, thereby enhancing the predictive accuracy of MLFFs.
Aza-BODIPY dimers represent promising molecular systems for efficient triplet-state generation through either intramolecular-singlet fission (iSF) or spin-orbit charge transfer intersystem crossing (SOCT-ISC). In this work, we investigate the role of molecular geometry in governing these mechanisms across four regioisomeric aza-BODIPY dimers (D[1,1], D[1,3], D[3,3], and D[2,2]) using multireference quantum-chemical calculations. Ground- and excited-state properties were analyzed at the MP2 and SA-XMCQDPT levels of theory, while diabatic couplings and spin-orbit matrix elements were evaluated to estimate iSF and SOCT-ISC rate constants, respectively. Our results reveal that triplet formation is strongly governed by the torsional angle (Φ) between monomer units, with regio-connectivity exerting a secondary influence. Dimers D[1,1] and D[1,3] exhibit favorable iSF energetics and coupling magnitudes, whereas D[2,2] displays low iSF rate constant (kSF ) but enhanced SOCT-ISC activity. The D[3,3] dimer shows exothermic multiexciton formation but reduced iSF efficiency due to destructive coupling interactions. The dominant ISC channel proceeds through the S1-T3 transition with large spin-orbit coupling and a small energy gap. These findings provide critical mechanistic insights into geometry-dependent triplet generation in aza-BODIPY dimers.
Optically pumped magnetometers (OPMs) have demonstrated significant potential in weak magnetic field detection due to their high sensitivity. In this study, we developed an Mz-type optically pumped rubidium magnetometer using a paraffin-coated anti-relaxation vapor cell. The system optimization and performance characterization were conducted inside a magnetic shield. Specifically, the pump light intensity and radio-frequency (RF) magnetic field were jointly optimized by using the linewidth-amplitude ratio as the core metric. Based on the frequency-domain noise spectrum, the sensitivity in open-loop mode was measured to be approximately 30.8 pT/Hz^{1/2}. Furthermore, a closed-loop feedback locking technique was applied, reducing the measured noise floor under the tested conditions and improving the sensitivity to 22.9 pT/Hz^{1/2}, with a measured -3 dB bandwidth of 123 Hz. The dynamic characteristics were evaluated via magnetic-field step response, showing that the system could track magnetic-field changes stably under closed-loop operation. Finally, by using tri-axial modulation and frequency-domain demodulation, we overcame the scalar measurement limitation of traditional Mz magnetometers. This work realizes vector magnetic field detection and provides a technical basis for applications such as geomagnetic navigation and magnetic anomaly detection.
The Bell-Bloom-type optically pumped atomic magnetometers are well suited for weak geomagnetic field detection. However, conventional single-beam pumping introduces an atomic spin polarization gradient, which limits the measurement accuracy and sensitivity. To address this issue, this paper proposes and experimentally demonstrates a Bell-Bloom-type rubidium FID magnetometer scheme integrating orthogonally polarized counter-propagating pumping and multi-pass probe detection. This design homogenizes the atomic spin polarization distribution and suppresses light shifts and power broadening effects induced by the pump beam. Meanwhile, the five-pass probe configuration significantly enhances the signal amplitude. Experimental results reveal that, compared with the traditional single-beam pumping and single-pass detection scheme, the proposed magnetometer achieves a remarkable improvement in magnetic field measurement accuracy, and the magnetic field sensitivity is improved from 18.9 pT/\sqrt{Hz} to 3.1 pT/\sqrt{Hz}. This work provides an effective technical approach and reference for optimizing the performance of atomic magnetometers and extending their applications in integrated arrays.
We report a study of the diazabicyclo[2.2.2]octane (DABCO) molecule photoionized using VUV synchrotron radiation in combination with an ion--electron coincidence spectrometer. We determine accurately the adiabatic ionization energy to $7.199\pm0.006$~eV. Two vibrational progressions of DABCO cation ground state are resolved at $847~\text{cm}^{-1}\pm27~\text{cm}^{-1}$ and $1257~\text{cm}^{-1}\pm67~\text{cm}^{-1}$, which we assign to modes of $e'$ symmetry. Analysis of the photoelectron angular distribution shows that the anisotropy parameter depends on the vibrational excitation. This dependence of the $β$ parameter with the vibrational excitation is attributed to the scattering of the outgoing wavefunction mediated by high-lying Rydberg states.
The correlation discrete variable representation (CDVR) enables efficient quantum dynamics calculation with the multi-layer multi-configurational time-dependent Hartree (MCTDH) approach on general potential energy surfaces. It employs a time-dependent quadrature to compute potential energy matrix elements, thereby eliminating the need to refit the potential to a sum of products form. The non-hierarchical CDVR conserves the inherent symmetry properties of tree-shaped wavefunction representations and drastically reduces the number of grid points compared to the original hierarchical CDVR. However, it requires projection on the space spanned by the single-hole functions (SHFs) at each node of the tree, which can introduce unphysical couplings for unconverged basis sets. In this work, the non-hierarchical CDVR is revisited and a revised approach that avoids explicit projection on the single-hole space is introduced. The computational costs of the revised approach scale favorably with the number of single-particle functions (SPFs): for a tree with three edges at each node and $n$ SPFs at each edge, a n^4 scaling is achieved. Furthermore, a revised scheme that uses artificial SPFs to systematically increase the accuracy of the CDVR quadrature is presented. Computations studying the photodissociation of NOCl, the vibrational states of methyl, and the non-adiabatic quantum dynamics of photoexcited pyrazine demonstrate the accuracy and efficiency of the revised non-hierarchical CDVR. Notably, for the 24-dimensional pyrazine system the use of the CDVR does not increase the required CPU time compared to calculations utilizing the sum of products form of the vibronic coupling model.
2604.02144Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this work, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.
We demonstrate that strong-field ionization of atoms in circularly polarized laser fields generates a photoelectron spin texture with toroidal topology in momentum space. Using time-dependent Schrödinger equation simulations, spin-resolved classical-trajectory Monte Carlo calculations, and an extended spin-resolved strong-field approximation including intermediate excitation pathways, we show that the rotation angle of this spin torus provides access to attosecond relative time delays associated with photoelectron wave packets released by tunneling from the counter-rotating and co-rotating \(p\)-orbital channels. When intermediate-state dynamics become significant, the torus develops a clear splitting. These results establish photoelectron spin textures as a complementary source of dynamical information beyond conventional momentum spectroscopy, and identify spin polarization as a robust internal degree of freedom for self-referenced attosecond metrology.
Auxiliary Field Quantum Monte Carlo (AFQMC) has emerged as a powerful framework for treating strongly correlated electronic systems, offering a favorable balance between computational cost and accuracy. In this paper, we present a novel AFQMC method that uses the isometric tensor hypercontraction (ITHC) technique to diagonalize the two-body Coulomb interaction of molecular electronic Hamiltonians by introducing additional fictitious fermionic modes. Our method shows reduced theoretical complexity and better practical performance for both propagation and local energy evaluation compared to the standard AFQMC method. We demonstrate the efficacy of this approach by computing the ground-state energies of a linear $\ce{H10}$-chain and the benzene molecule. Our results show that the extended-basis AFQMC recovers many-body correlations with a precision comparable to that of high-level wavefunction methods such as Coupled Clusters (CC) or Density Matrix Renormalization Group (DMRG), while offering significantly improved scaling.
In this work, we develop a new framework for computing atom-resolved contributions to optical properties based on atoms-in-molecules (AIM) schemes. The formalism is independent of the specific AIM method and is made rigorous by partitioning momentum matrix elements into atomic contributions while exactly satisfying the relevant sum rules. We apply it to second-harmonic generation (SHG) in six representative UV and deep-UV nonlinear-optical crystals, namely $β$-\ce{BaB2O4} (BBO), \ce{LiB3O5} (LBO), \ce{CsB3O5} (CBO), \ce{CsLiB6O10} (CLBO), \ce{KBe2BO3F2} (KBBF), and \ce{LiCs2PO4} (LCPO). The atom-triplet decomposition reveals a clear hierarchy for the largest SHG component of each crystal. In general, two-center terms provide the leading contribution, one-center terms remain comparatively small, and fully three-center terms supply an important secondary contribution. A motif-triplet decomposition further indicates behavior dominated by the anionic framework in KBBF and LBO. In BBO, CBO, and CLBO, contributions from the anionic framework and the cation sublattice act cooperatively, although the cation contribution is crystal dependent. Moreover, cooperative contributions from the phosphate framework and the Cs sublattice are also observed in LCPO, where the O-Cs contribution is particularly significant. These results may provide a new perspective for understanding the microscopic origin of SHG in nonlinear-optical materials.
We present a method for performing low frequency electric field sensing via ionization detection of Rydberg atoms in a collimated atomic beam. A collimated beam avoids much of the electric field screening effects that are common in warm vapor cells due to the accumulation of alkali-metal atoms on glass surfaces. Further, a beam facilitates a spatially separated region for high signal-to-noise readout via ionization detection. Using this approach, we measure DC Stark shifts from external fields with frequencies as low as 1 Hz. The sensor demonstrates a sensitivity of better than 1 mV/m$\sqrt{\rm {Hz}}$ for frequencies above 20 Hz and $0.14(4)$ mV/m$\sqrt{\rm {Hz}}$ above 500 Hz with a linear dynamic range of over 50 dB.
We present TUNA, an open-source quantum chemistry program specifically designed for atoms and diatomic molecules. Within this narrow molecular domain, a broad and consistent set of electronic structure methods and calculation types is available. Energies, optimisations, vibrational frequencies, response properties, coordinate scans and ab initio molecular dynamics trajectories can be accessed through an intuitive command-line interface. A single principle underlies TUNA: once a method can be used to evaluate the energy, all properties follow from numerical differentiation. This makes the program both a transparent teaching platform and a compact environment for benchmarking methods on diatomics $\unicode{x2014}$ among the most simple yet instructive systems in quantum chemistry. Reference implementations including density functional theory, many-body perturbation theory and coupled cluster theory, supported by detailed theoretical documentation, make TUNA an accessible foundation for developing improved methods and algorithms in electronic structure.
Predicting the perceived intensity of odorants remains a fundamental challenge in sensory science due to the complex, non-linear behavior of their response, as well as the difficulty in correlating molecular structure with human perception. While traditional deep learning models, such as Graph Convolutional Networks (GCNs), excel at capturing molecular topology, they often fail to account for the biological and perceptual context of olfaction. This study introduces VIANA, a novel "tri-pillar" framework that integrates structural graph theory, character value embeddings, and phenomenological behavior. This methodology systematically evaluates knowledge transfer across three distinct domains: molecular structure via GCNs, semantic odor character values via Principal Odor Map (POM) embeddings, and biological dose-response logic via Hill's law. We demonstrate that knowledge transfer is not inherently positive; rather, a balance must be maintained in the volume of information provided to the model. While raw semantic data led to "information overload" in domain-informed models, applying Principal Component Analysis (PCA) to distill the 95% most impactful semantic variance yielded a superior "signal distillation" effect. Results indicate that the synthesis of these three knowledge transfer pillars significantly outperforms baseline structural models, with VIANA achieving a peak R^2 of 0.996 and a test Mean Squared Error (MSE) of 0.19. In this context, VIANA successfully captures the physical ceiling of saturation, the sensitivity of detection thresholds, and the nuance of odor character value expression, providing a domain grounded simulation of the human olfactory experience. This research provides a robust framework for digital olfaction, effectively bridging the gap between molecular informatics and sensory perception.
Computational chemistry has become an indispensable tool for generating data and insights, pervading all branches of experimental chemistry. Its most central concept is the potential energy hypersurface, key to all chemistry and materials science, as it assigns an energy to a molecular structure, the necessary ingredient for reaction mechanism elucidation and reaction rate calculation. Density functional theory (DFT) has been the most important method in practice for obtaining such energies, which is mirrored in the use of high-performance computing hardware. In the last two decades, a new class of surrogate potential energy functions has been evolving with remarkable properties: quantum accuracy combined with force-field speed. Until very recently, their application was hampered by the fact that they needed to be trained on truly large system-specific data sets, generated before a computational chemistry study could be started (in sharp contrast to DFT, which, as a first-principles method, works out of the box, but at a far higher price of computational cost). Very recently, this roadblock has been overcome by so-called foundation machine learning interatomic potentials, which are poised to completely change the way we do computational chemistry, likely prompting us to abandon DFT as the prime method of choice for this purpose in less than a decade.
Electric field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the use of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter their orientations. In this work, we propose two molecular reference frames -- the principal axis frame and the local reference frame -- to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis- and trans-formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric field-controlled chemistry.
Optical atomic clocks have been rapidly developing in recent decades, resulting in major improvements in both precision and accuracy. As a result, they have become instrumental in multiple areas of applied and fundamental research. Despite all atomic frequency references having more than two energy-levels, the commonly used model for evaluating their ultimate limits assumes a two-level atom. This leads to frequency interrogation protocols and theoretical stability bounds that are suboptimal for a true multi-level atom. The most fundamental stability bound assumes two noise sources - quantum projection noise and spontaneous decay from the excited state. In this work, we analyze a model that includes these noise types and is generalized beyond the two-level assumption, where spontaneous decay can branch to more than a single ground state. This model allows for detection and exclusion of atomic frequency interrogations in which the atom decayed, leading to a frequency stability improvement of up to $\approx 4.5 \text{ dB}$ compared with the two-level model. Furthermore, we identify an even greater stability enhancement of $\approx 5.4 \text{ dB}$ for frequency comparisons between atoms in an odd parity Bell state. These enhancements are particularly relevant for the numerous trapped-ion optical clock species that operate close to lifetime-limited stability. We calculate new stability limits for those cases and provide a detailed experimental protocol for frequency interrogation with an $^{27}\text{Al}^{+}$ optical ion clock.
Accurate prediction of acid dissociation constants (p$K_{\rm a}$) and the determination of dominant protonation states is critical in drug discovery, influencing molecular properties such as solubility, permeability, and protein-ligand binding. We present Acep$K_{\rm a}$, an advanced application integrated into the PlayMolecule AI platform. Acep$K_{\rm a}$ is built upon the theoretically rigorous Uni-p$K_{\rm a}$ framework, which unifies statistical mechanics with representation learning. By modeling the complete protonation ensemble rather than treating p$K_a$ as a scalar regression target, Acep$K_{\rm a}$ ensures thermodynamic consistency across coupled ionization sites. We describe the application's enhanced architecture, which features a retrained Uni-Mol backbone achieving state-of-the-art performance on standard benchmarks. Furthermore, we detail critical engineering advancements. These include AceConfgen, a proprietary GPU-accelerated conformer generator that achieves a ~40x speed-up compared to NVIDIA's nvmolkit, a streamlined inference engine to directly protonate molecules, and a 3D-aware modality for applying protonation states to bound ligand poses. Finally, we discuss the integration of Acep$K_{\rm a}$ into the PlayMolecule AI ecosystem, a modern AI-assisted environment for molecular modelling and drug discovery.