2026-07-17 | | Total: 8
We report on the ground and numerous excited electronic states. In the ground state the C4n rings are closed-shell systems possessing polyynic structures and can be classified as double anti-aromatic molecules. In their energetically lowest lying triplet state the rings exhibit aromatic cumulenic structures. The overall change in the electronic structures is rather dramatic upon the found moderate geometric changes from polyynic to cumulenic structure. Among others, Hund's rule is violated in both C8 and C12 in their cumulenic structures. We mention that until now, graphene is the only carbon allotrope reported to violate Hund's rule. The reasons for the violation are analyzed. Much effort has been invested to understand the relaxation pathways of the low-lying states leading the C8 from polyynic to cumulenic geometry and vice versa. On its minimum energy path, the first singlet excited state changes from open-shell character in the polyynic structure to a closed-shell state in the cumulenic structure. The cumulenic state lowest in energy is an open-shell singlet which relaxes to the closed-shell polyynic ground state.
The kinetics of protein-ligand binding systems are increasingly recognized as a key determinant of drug efficacy, yet remain far harder to compute than binding affinities. Existing kinetics methods either bias the dynamics along a collective variable (CV), demanding careful system-specific CV design, or use path sampling, which keeps the dynamics unbiased but can struggle to converge rates out of deep free-energy wells and often relies on hand-engineered descriptors. By combining the `best of both worlds', we propose a method to compute accurate kinetics for general ligand-unbinding problems at modest computational expense and minimal fine tuning, building on the AI for Molecular Mechanism Discovery (AIMMD) path sampling framework. To avoid the need for feature engineering, we opt for modelling the committor with a single descriptor-free, equivariant graph neural network shared across all systems. We also partially flatten deep bound-state wells with a static, basin-restricted bias potential. This improves convergence by lifting the path sampling state boundary out of regions, where the committor is hard to learn, while leaving the reactive region strictly unbiased. Across host-guest and protein-ligand systems spanning roughly 17 orders of magnitude in residence time, the method robustly recovers rates in line with reference and experimental values. Simultaneously, and without further sampling, it also reconstructs the underlying unbinding mechanisms. We additionally find that accurate rates do not require globally accurate committor models, allowing for efficient kinetics estimation even in a low-data training regime. Requiring little system-specific setup, our approach offers an efficient and broadly generalizable route to binding kinetics, and its shared committor architecture lays crucial groundwork for probing structure-kinetics relationships across ligand series in drug discovery.
Dissociative chemisorption (DC) of H$_2$ on the Cu(111) surface is a prototypical problem for understanding elements of heterogeneous catalysis [Science 326, 832 (2009)]. The challenge lies in modeling the reaction dynamics that in turn reflects a classical potential for atomic deformations, friction, and inelastic scattering. Here, I test the use of a set of range-separated hybrid (RSH) van der Waals density functionals (vdW-DFs) [JPCM 37, 211501 (2025)] on their ability to describe the classical barrier for dynamics in this H$_2$+Cu(111) DC problem. I furthermore document use of a variant for fast accurate predictions of the molecular quasi-particles (QPs), finding excellent performance across a set of small molecules that are often studied in catalysis. Finally, I suggest and implement a way to use that QP focus to identify what I consider a best-possible non-empirical (yet adsorbate specific) RSH vdW-DF version, denoted AHBR($γ^*$) for H$_2$ DC modeling, navigating what are partly conflicting requirements on the molecule and metal sides. I find that the AHBR($γ^*$) can determine the classical H$_2$+Cu(111) DC barrier height close to chemical accuracy. I suggest that DC modeling can test broader relevance of the physics underpinning these RSH vdW-DFs.
Machine-learning force fields (MLFFs) are reliable only near their training distribution, making efficient construction of diverse training sets a major bottleneck for both train-from-scratch and foundation fine-tuning workflows. Active learning can reduce this cost, but standard model-committee uncertainty is impractical for foundation MLFFs because each committee member requires a separate fine-tuning run. We present an active-learning workflow based on last-layer-projection regression (LLPR), a forward-pass-cheap per-configuration uncertainty estimator. Across molecular, condensed-phase, and electrolyte systems, LLPR identifies compact, high-value training sets that recover full-data accuracy using only a small fraction of electronic-structure labels. In foundation-model fine-tuning, LLPR-selected configurations reach the full-pool fine-tuning ceiling with substantially fewer labels than random selection. In iterative electrolyte fine-tuning, LLPR detects unphysical local coordination before DFT labelling, provides an absolute force-error threshold, and enables automatic termination of the learning loop. The resulting models reproduce reference density and ion-coordination structure, providing a scalable uncertainty-quantification strategy across MLFF training regimes.
Hydrated electrons of opposite spins pair to form dielectrons at sufficiently high concentrations that can be achieved by dissolution of alkali metals in water. While experimental investigations of these systems are challenging due to their vigorous, even explosive, reactivity, simulations open the possibility to characterize the structure and reactivity of hydrated dielectrons and the effects of alkali cations thereon. Here, we present ab initio molecular dynamics simulations of a hydrated dielectron without or with explicit Li$^+$ or Cs$^+$ counterions. While the overall solvation structure is preserved in all these systems, the presence of cations has a distinct effect of increasing the dielectron gyration radius by about 10% and forming cation-specific cation-dielectron arrangements. Moreover, analysis of water bond lengths reveals in all studied systems a substantial elongation of first-shell O-H bonds oriented toward the dielectron, providing a structural explanation for vibrational red-shifts observed in resonance Raman measurements. Finally, at the 10 ps simulation timescale, rare reactive events were observed, albeit only for the system without metal cations, where hydride intermediates stable on picosecond timescales were identified. These observations also suggest that on the investigated timescales, metal cations may suppress hydrated dielectron reactivity.
We present a theory of a flow of salt solutions in long axisymmetric channels induced by concentration and pressure drops between their ends. The consideration is restricted to thin, compared to the local radius, electrostatic diffuse layers, but remains valid even when the concentration drop is quite large. We show that the magnitude of the diffusio-osmotic fluid flow rate $Q_{DO}$ in the cylinder is the same as in the slit of equal to its diameter thickness, but channels of variable cross-sections could either retard or enhance it, depending on their geometry. The application of the pressure drop $Δp \neq 0$ results in an extra contribution $Q_{P}$ to the total flow rate of fluid $Q$, but does not affect $Q_{DO}$. We calculate the curves $Δp (Q)$ for several axisymmetric channels and conclude that they are nearly linear, with the sensitive to the shape slopes. This leads to the possibility of introducing a simple, but rather accurate, cylinder approximation, where the radius of the imaginary cylinder is related to a hydrodynamic resistivity of the real channel and can be easily determined, if its geometry is known. We also derive an equation relating the ionic flux with the total flow rate of fluid and demonstrate that both the sign and magnitude of ionic flux could be tuned by using the appropriate channel shape. Our analysis provides a framework for interpreting experimental and numerical data, as well as may guide the design of micro- and nanofluidic devices.
An accurate estimation of the molecular abundances of isomers in the interstellar medium (ISM) is necessary to unravel the underlying chemistry and physics. After the recent detections of both isomers of formic acid ($cis-$ and $trans-$HCOOH) in dense dark cold clouds, their accurate molecular line modeling became of interest. The conditions of these environments do not necessarily follow the local thermodynamic equilibrium, thus taking into account the competition between the radiative and collisional processes is required. This involves the knowledge of the rotational excitation data for collisions with the most abundant interstellar species \textemdash He and H$_2$. In this paper, the first potential energy surfaces (PES) for the interaction of the two rotamers of formic acid with He atoms are computed using the explicitly correlated coupled-cluster theory [CCSD(T)-F12]. The obtained PESs demonstrate qualitative similarities and high anisotropy. The global minima are found with $V=-53.0$ cm$^{-1}$ and $V=-46.0$ cm$^{-1}$ for $cis-$HCOOH and $trans-$HCOOH respectively. Collisional excitation cross sections calculated for total energies up to 100 cm$^{-1}$ demonstrate similar propensity rules for both isomers. Quantitative differences of the cross sections associated with the two rotamers are also discussed.
Reentrant coil-globule-coil transitions, in which a polymer collapses and then reexpands as a single parameter is varied, have been observed across diverse soft matter systems, yet the minimal ingredients required to produce them remain unclear. Using molecular dynamics simulations of coarse-grained polymers interacting with a single species of attractive crowder, we show that crowder volume fraction $φ_c$ alone is sufficient to drive a complete reentrant transition. At low $φ_c$, crowders bridge distant monomers and drive cooperative collapse; at high $φ_c$, saturation of monomer binding sites suppresses bridging connectivity and produces reentrant expansion. This density-driven transition is absent with purely repulsive crowders, which produce only monotonic compaction while preserving self-avoiding walk (SAW) chain statistics. In contrast, bridging breaks SAW universality: the rescaled size distributions no longer collapse onto a universal curve, and the conformational distributions trace the full coil-globule-coil trajectory as $φ_c$ is varied. For charged polymers with explicit counterions, electrostatics amplifies rather than suppresses reentrance: bridging crowders displace counterions from the chain, and upon saturation the unscreened backbone charges drive expansion well beyond the original chain size. Saturable geometric bridging thus emerges as a minimal mechanism linking reentrant phenomena across neutral and charged polymers in crowded environments.