2026-10-19 | | Total: 10
Using the density matrix renormalization group method, we solve the Pariser--Parr--Pople model of $π$-electron systems to determine the character and energies of the low-lying eigenstates of a lycopene dimer comprising two 22-carbon-atom chains. To complement these calculations, we use the molecular-exciton, valence-bond model of polyene systems to derive selection rules governing interchain coupling, which provide a physical framework for interpreting the DMRG results. We find that dimer geometry strongly influences both the character and energetics of the low-lying eigenstates. In vertically stacked dimers, neither the bright-dimer state nor the triplet-pair states contain any interchain charge-transfer (CT) character. In contrast, in slip-stacked dimers, both the bright-dimer state and triplet-pair states do contain interchain CT character. We also find that increasing the interchain coupling preferentially stabilizes states with greater interchain character, leading to a reordering of the low-lying eigenstates. These results have direct implications for energetically possible singlet fission (SF) mechanisms in carotenoid aggregates. For vertically stacked dimers, the only energetically viable pathway identified here proceeds through intermediate intrachain triplet-pair states, whereas slip-stacked geometries allow an additional SF pathway mediated by interchain CT states. Our results establish a direct connection between carotenoid packing geometry, symmetry-controlled electronic coupling, and the energetically accessible pathways for SF.
We derive a vibro-polaritonic Raman activity from third-order cavity Born-Oppenheimer perturbation theory (CBO-PT) accounting for weak interactions between electrons and low-frequency cavity modes in the vibrational strong coupling (VSC) regime. This work extends the recently presented CBO-PT linear response approach [J. Phys. Chem. Lett. 2024, 15, 2262-2269] from infrared (IR) to vibro-polaritonic Raman spectroscopy. In a first step, we discuss cavity-induced corrections to the static dipole polarizability tensor, which manifest as molecular and cavity mode contribution determined by the polarization-projected static hyperpolarizability tensor. We then obtain the corresponding vibro-polaritonic Raman activity in double harmonic approximation, which contains molecular, cavity and mixed contributions similar to previously discussed CBO-PT IR intensities. The cavity Raman activity can be interpreted as effective transmission component, which is exclusively determined by the polaritonic Raman response. The methodology and its relation to vibro-polaritonic IR spectroscopy is illustrated for a formaldehyde molecule under VSC. An implementation of the CBO-PT linear response approach for vibro-polaritonic Raman and IR spectroscopies is provided via the open-source Python package CBOvibSpec.
Companion paper I (Khanal and Arunan, in this issue; arXiv:2610.16629 [physics.chem-ph]) reported reference BEs and F-bond distances for two new databases SEFB-15 (DF$\cdots$CH$_3$) and CSFB-14 (DF$\cdots$NH$_3$). Here, we assess performance of popular ab initio methods (MP2, MP3, MP4(DQ), MP4(SDQ), CCSD, and QCISD) with valence- and full-electron correlation treatments, and sixteen Minnesota DFAs, including M06-2X already assessed in companion paper I. CCSD/AVDZ reproduces reference F-bond distances closely, although BEs predicted by it deviate more. Among the 16 Minnesota DFAs, M08-HX reliably reproduces both reference values for DF$\cdots$CH$_3$ and DF$\cdots$NH$_3$, comparable to M05-2X and M06-2X for DF$\cdots$NH$_3$. Selected meta-GGAs (M06-L, M11-L, and MN15-L) exhibit spurious fractional charge transfer in cationic complexes even at infinite separation, whereas hybrid Minnesota DFAs do not, except M05 and M06 for NH$_3$F$^+\cdots$CH$_3$ and NH$_2$F$_2^+\cdots$CH$_3$. Although M08-HX shows slightly larger monomer delocalization error than M05-2X and M06-2X, it is more reliable for predicting BEs of SEFB complexes, indicating that monomer delocalization error alone does not govern DFA performance. Bader's intermolecular delocalization indices (DIs) further support M08-HX as more reliable than M05-2X and M06-2X for SEFB, while the larger DIs from M05-2X and M06-2X explain the degraded performance of their dispersion-corrected variants for DF$\cdots$CH$_3$, with only a minor impact for DF$\cdots$NH$_3$, highlighting the system-dependent nature of functional delocalization error. Present work reveals that M08-HX is more reliable than M05-2X and M06-2X for open-shell DF$\cdots$CH$_3$ complexes, while all three functionals show comparable reliability for closed-shell DF$\cdots$NH$_3$ complexes.
Electrochemical processes at charged interfaces lie at the heart of modern technologies and their continued improvement depends on insights into the relevant atomistic processes. However, the necessary tools to impose constant-voltage conditions in molecular dynamics (MD) simulations typically directly target the electrode charge, making their extension to methods that go beyond an empirical description of metals more cumbersome. In this work, we propose an MD protocol to realize a thermodynamic ensemble in which voltage control is attained, but which is formulated purely within the framework of finite-field simulations in 3D periodic boundary conditions. We show a derivation and implementation of such a control device, which we coin an electrostat, in the ELECTRODE package within LAMMPS. Our simulations systematically benchmark the derived theory against other established approaches, for example, showing its equivalence to a potentiostat, and quantify the importance of properly accounting for ensemble fluctuations. We find that for systems of several thousands of atoms, one obtains the same equilibrium properties with the electrostat as with a simpler constant-field simulation. The difference between the two methods only manifests itself at smaller length scales, thereby establishing practical limits for the inclusion of voltage fluctuations in atomistic simulations of electrochemical systems. The strength of this framework also lies in its simplicity of extension to first-principles ab initio and reactive machine-learning potentials.
The combinatorial explosion in monomial counts limits the extension of global permutationally invariant polynomial (PIP) potential energy surfaces to larger molecular systems. We use graph-theoretic language to reformulate PIPs as colored multigraphs and develop a series of theoretical and computational tools. First, we introduce an offline PIP generation algorithm whose cost has an upper bound independent of atom count for a fixed degree cutoff and number of element types. Second, we introduce a new online PIP evaluation algorithm, termed tensor-network descent. For the full PIP set through degree $d$ in a system of $M$ atoms, the evaluation complexity for $1\le d\le14$ is $O\!\left(f(d)M^{\left\lfloor(\sqrt{1+8d}+1)/2\right\rfloor}\right)$. At the commonly used cutoff of $d=5$, the online evaluation complexity decreases from the $M^{10}$ scaling of conventional explicit term-by-term computation to $M^3$. Numerical tests using the full PIP set through degree 5 show that our method achieves speedups of several orders of magnitude on most MD17 systems and remains computationally affordable for systems containing more than one thousand atoms. The method has a native tensor form, supports automatic differentiation, and provides a unified theoretical construction of four classes of PIP descriptors defined consistently across systems---global, atom-centered, bond-centered, and cluster-centered---establishing a representation and computational foundation for PIP potentials for large systems and for transfer across systems. Together, these results show that PIPs are neither inherently restricted to small systems nor intrinsically nontransferable.
We have developed a reduced-cost framework for excited-state calculations by combining the state-specific frozen natural orbital (SS-FNO) approach with the third-order algebraic diagrammatic construction [ADC(3)] method and tensor hypercontraction (THC). This strategy enables a systematic and controllable truncation of the virtual orbital space, while significantly lowering the computational scaling through the THC representation of the electron repulsion integrals. The method allows black-box application via adjustable thresholds and includes a perturbative correction that compensates for the truncation error. We have used the CIS(D) method to generate the natural orbitals for the excited states. The performance of the SS-FNO-THC-ADC(3) method is assessed for valence, Rydberg, and charge-transfer excitations. Benchmark calculations demonstrate that the combined approach achieves excellent agreement with the canonical ADC(3) results. The efficiency of the present implementation is demonstrated by calculating the excitation energy of a molecule containing 31 atoms and more than one thousand basis functions.
Transient microscopy has evolved to include a family of well-established techniques for studying the spatiotemporal dynamics of excited states in a variety of materials systems. In these methods, a laser pulse is focused to locally excite a material, and the subsequent evolution of the excited species is visualized to extract transport parameters such as diffusivities. Typically, transport coefficients are extracted from experimental datasets using curve-fitting methods to analyze the change in the mean squared displacement (MSD) of the spatial profiles over time. However, this approach can lead to inaccurate or misleading diffusivity values under many common experimental conditions, and is only formally accurate for the case of a single diffusing species. In this work, we introduce a new approach to the analysis of transient microscopy data based on the evolution of profile shapes in reciprocal space. We use simulated datasets to quantitatively analyze the performance of our proposed methods and compare them to current state-of-the-art analysis techniques. We demonstrate that reciprocal-space-based analysis methods have the highest overall performance in terms of accuracy and precision for normal diffusive transport over a range of realistic initial conditions and signal-to-noise ratios. Furthermore, we highlight the limitations of using MSD alone to characterize anomalous transport in materials, and we show how analysis of spatial frequencies can yield a more comprehensive understanding of complex excited-state dynamics and their underlying physical mechanisms. Our results show that reciprocal-space-based analysis methods are powerful tools for characterizing and quantifying many different types of transport behaviors.
The HCO$^+$/HOC$^+$ abundance ratio has been observed to vary significantly in different astrophysical environments, making it an important diagnostic tool for exploring the interstellar medium. However, uncertainties in the reaction rates and branching ratios of key chemical reactions involving these ions affect the conclusions drawn from observations of either these species or other species that are derived on their basis. Here, we present isomer-resolved ion-molecule reaction rate coefficients of reactions involving HCO$^+$ and HOC$^+$ using a cryogenic ion trap experiment. We measured the rate coefficients of the isomerization reaction of HOC$^+$ to HCO$^+$ with H$_2$ and CO, as well as three different reactions producing both isomers: CO$^+$ + H$_2$ , H$_3^+$ + CO, and ArH$^+$ + CO. Our results show that the reaction of CO$^+$ with H$_2$ predominantly produces HCO$^+$ with a negligible yield of HOC$^+$ . This result is in agreement with recent theoretical predictions, but in major contrast to the generally accepted branching ratio based on previous experimental results. We also find that the isomerization reaction of HOC$^+$ to HCO$^+$ is efficient with both H2 and CO. We were able to assess the impact of our new measured reaction rate coefficients on the predicted HCO$^+$/HOC$^+$ abundance ratio using a state-of-the-art astrochemical model.
We investigate diffusion-controlled reactions on a reactive surface with time-dependent reactivity. Starting from a conventional formulation via diffusion equation with Robin boundary condition for the macroscopic concentration of particles, we aim at quantifying how temporal variations of reactivity affect the survival probability and the associated first-reaction time (FRT) of a single particle on the reactive surface. Using an adiabatic description, we first identify five long-time regimes of the first-reaction time distribution, ranging from exponential and stretched-exponential to power-law decay. In particular, sufficiently rapid decay of the reactivity leads to a nonzero probability of never reacting, despite diffusion in a bounded domain. Second, we provide a probabilistic reformulation of this problem and compare it with a different scenario of encounter-dependent reactivity. Our theoretical predictions are validated via extensive Monte Carlo simulations.
A membrane traversed by a dense array of nanochannels acts as a photonic metamaterial: with channels far smaller than the wavelength of light, the optical response is set by the architecture rather than by the constituent materials. Here each nanochannel is turned into a nonlinear optical meta-atom by filling anodic aluminum oxide (AAO) with benzil (BZL), a noncentrosymmetric molecular crystal, and the crystallization pathway programs their orientation. In free-standing membranes every channel solidifies independently, the orientations remain uncorrelated, and the effective second-harmonic generation (SHG) response is nearly polarization isotropic; its magnitude follows a core - shell model in which an optically inactive interfacial layer of 2.6 nm reduces the active volume of each meta-atom, reproducing on bare alumina the value previously found on alkyl-grafted silica and indicating a host-insensitive limit. Seeding crystallization from a bulk benzil layer instead aligns the meta-atoms over submillimeter domains and produces a macroscopically anisotropic effective susceptibility. Channel diameters from 32 to 270 nm span the transition out of the subwavelength regime, and the anomalous directional SHG emission reported here is confined to the metamaterial limit. Confined crystallization thus offers a route to nonlinear optical metamaterials whose effective susceptibility tensor is written during processing.