2025-12-30 | | Total: 11
H2O is a unique substance with exceptional thermal properties arising from the subtle interplay between its electronic, phononic, and structural degrees of freedom. Of particular interest in H2O are the negative thermal expansion (NTE) phenomena, observed in its solid phase (ice) at low temperature, and in its liquid phase (water) near the freezing temperature. Furthermore, ice and water exhibit the abnormal volume isotope effect (VIE), where volume expansions occur when replacing H with its heavier isotope, deuterium (D). In order to capture more conceptual and intuitive understanding of intriguing NTE and VIE phenomena in ice and water, we have explored isotope effects in their NTE and melting properties by employing a type of Born-Oppenheimer-approximation approach and the Lindemann criterion. Our findings demonstrate that unusual isotope effects in these phenomena stem from competition between zero-point-energy phonons, thermal phonons, and the hydrogen bonding in H2O. All these components originate from nuclear quantum mechanical (QM) processes, revealing that QM physics plays a crucial role in the seemingly classical ice/water systems.
Polymers-macromolecular systems composed of repeating chemical units-constitute the molecular foundation of living organisms, while their synthetic counterparts drive transformative advances across medicine, consumer products, and energy technologies. While machine learning (ML) models have been trained on millions of quantum chemical atomistic simulations for materials and/or small molecular structures to enable efficient, accurate, and transferable predictions of chemical properties, polymers have largely not been included in prior datasets due to the computational expense of high quality electronic structure calculations on representative polymeric structures. Here, we address this shortcoming with the creation of the Open Polymers 2026 (OPoly26) dataset, which contains more than 6.57 million density functional theory (DFT) calculations on up to 360 atom clusters derived from polymeric systems, comprising over 1.2 billion total atoms. OPoly26 captures the chemical diversity that makes polymers intrinsically tunable and versatile materials, encompassing variations in monomer composition, degree of polymerization, chain architectures, and solvation environments. We show that augmenting ML model training with the OPoly26 dataset improves model performance for polymer prediction tasks. We also publicly release the OPoly26 dataset to help further the development of ML models for polymers, and more broadly, strive towards universal atomistic models.
The discovery of novel odorant molecules is key for the fragrance and flavor industries, yet efficiently navigating the vast chemical space to identify structures with desirable olfactory properties remains a significant challenge. Generative artificial intelligence offers a promising approach for \textit{de novo} molecular design but typically requires large sets of molecules to learn from. To address this problem, we present a framework combining a variational autoencoder (VAE) with a quantitative structure-activity relationship (QSAR) model to generate novel odorants from limited training sets of odor molecules. The self-supervised learning capabilities of the VAE allow it to learn SMILES grammar from ChemBL database, while its training objective is augmented with a loss term derived from an external QSAR model to structure the latent representation according to odor probability. While the VAE demonstrated high internal consistency in learning the QSAR supervision signal, validation against an external, unseen ground truth dataset (Unique Good Scents) confirms the model generates syntactically valid structures (100\% validity achieved via rejection sampling) and 94.8\% unique structures. The latent space is effectively structured by odor likelihood, evidenced by a Fréchet ChemNet Distance (FCD) of $\approx$ 6.96 between generated molecules and known odorants, compared to $\approx$ 21.6 for the ChemBL baseline. Structural analysis via Bemis-Murcko scaffolds reveals that 74.4\% of candidates possess novel core frameworks distinct from the training data, indicating the model performs extensive chemical space exploration beyond simple derivatization of known odorants. Generated candidates display physicochemical properties ....
Studies of quantum thermal effects on molecular excitation dynamics have often relied on oversimplified models, such as energy eigenstates or low-dimensional potentials, which fail to capture the complexity of real chemical systems. In reality, molecules are spatially extended and embedded in anisotropic environments, where molecular orbitals (MOs) play a central role in determining quantum behavior. To advance beyond these limitations, we propose a three-dimensional rotationally invariant system-bath (3D-RISB) model within the MO framework, with explicit inclusion of intramolecular vibrational motion. From this MO foundation, we derive numerically ``exact'' hierarchical equations of motion (MO-HEOM). As a demonstration, we analyze hydrogen molecules and hydrogen molecular ions with vibrational degrees of freedom, revealing their linear absorption spectra.
How quantum electron and nuclei motions affect biomolecular chemical reactions remains a central challengeable question at the interface of quantum chemistry and biology. Ultrafast charge migration in deoxyribonucleic acid (DNA) has long been hypothesized to play a critical role in photochemistry, genome stability, and long-range biomolecular signaling, however, direct real-time observation of these electronic processes has remained elusive. Here, we present a theoretical investigation and propose the concept of future experimental measurements of laser-driven charge dynamics in the canonical DNA nucleobase pairs thymine_adenine and cytosine_guanine. Attosecond-resolved simulations employing high-level ab initio methods reveal base-dependent ionization mechanisms, directional charge migration pathways, and electronic coherences that govern sub-femtosecond redistribution of electron density across hydrogen-bonded nucleobase interfaces. Accordingly, we propose the concept of a quantum attosecond scanning electron microscope, termed the quantum attomicroscope (Q-attomicroscope), a capable of imaging photoinduced quantum chemistry reactions in attosecond temporal resolution and sub-nanometer spatial precision. As a proof of principle, we propose to image the charge migrations dynamics in DNA which we studied theoretically. Together, our preceptive bridges theory, instrumentation, and control, outlining a pathway toward laser mediated manipulation of DNA structure with implications for repair processes, chemical reactivity, and future personalized medicine.
Understanding the intricate interplay between soot dynamics and chemical reactions within catalytic diesel particulate filters (CDPF) is crucial for enhancing both filtration efficiency and regeneration performance. In this paper, we establish a unified pore-scale multiphysics model based on the Eulerian-Lagrangian framework to comprehensively resolve the transport, deposition, and oxidation of soot. Distinguishing itself from conventional empirical correlations and stochastic-based approximations, the system models soot deposition through fundamental physical principles, integrating elastic deformation and surface adhesion mechanics at the particle-wall interface. Simultaneously, it incorporates a robust oxidation model that accounts for the competitive kinetics of both $\textrm{O}_2$ and $\textrm{NO}_2$ pathways, enabling comprehensive coverage of all CDPF operating regimes. Validated against three classical benchmark cases, the model demonstrates superior accuracy in capturing interfacial mass transfer and particle-wall interactions. Simulation under a typical CDPF low-temperature operating condition emphasizes the pivotal role of $\textrm{NO}_2$ and catalyst in promoting regeneration and reveals complex synergistic and competitive effects between distinct reaction pathways. Notably, the reaction rate of direct $\textrm{O}_2$ pathway is accelerated by a factor of 87 in the presence of the catalyst. For ultra-fine soot particles ($50~\mathrm{nm}$), the Brownian motion and thermophoretic forces directly dictate the deposition efficiency. Their strong thermal sensitivity also underscores the necessity for an integrated soot transport-deposition-oxidation framework. To support further research, the model implementation can be accessed at https://github.com/zhangyujing2001.
Water adopts many different crystal structures in its solid form. These provide insight into potential structures of water even in its liquid phase, and they can be used to calibrate pair potentials used for simulation of water. In crowded biological environments, water may behave more like ice than bulk water. The different ice structures have different dielectric properties. This brief primer is intended to facilitate further research.
We propose a simple and accurate approach to estimate the random close packing (RCP) fraction of binary hard-disk mixtures. By introducing a parameter based on the reduced third virial coefficient of the mixture, we show that the RCP fraction depends nearly linearly on this parameter, leading to a universal collapse of simulation data across a wide range of size ratios and compositions. Comparisons with previous models by Brouwers and Zaccone demonstrate that our approach provides the most consistent and accurate predictions. The method can be naturally extended to polydisperse mixtures with continuous size distributions, offering a robust framework for understanding the universality of RCP in hard-disk systems.
Generating low-energy conformer ensembles and identifying ground-state conformations from molecular graphs remain computationally demanding with physics-based pipelines. Current learning-based approaches often suffer from a fragmented paradigm: generative models capture diversity but lack reliable energy calibration, whereas deterministic predictors target a single structure and fail to represent ensemble variability. Here we present EnFlow, a unified framework that couples flow matching (FM) with an explicitly learned energy model through an energy-guided sampling scheme defined along a non-Gaussian FM path. By incorporating energy-gradient guidance during sampling, our method steers trajectories toward lower-energy regions, substantially improving conformational fidelity, particularly in the few-step regime. The learned energy function further enables efficient energy-based ranking of generated ensembles for accurate ground-state identification. Extensive experiments on GEOM-QM9 and GEOM-Drugs demonstrate that EnFlow simultaneously improves generation metrics with 1--2 ODE-steps and reduces ground-state prediction errors compared with state-of-the-art methods.
Superionic conduction in solid-state materials is governed not only by static factors, such as structure and composition, but also by dynamic interactions between the mobile ion and the crystal lattice. Specifically, the dynamics of lattice vibrations, or phonons, have attracted interest because of their hypothesized ability to facilitate superionic conduction. However, direct experimental measurement of the role of phonons in ionic conduction is challenging due to the fast intrinsic timescales of ion hopping and the difficulty of driving relevant phonon modes, which often lie in the low-energy THz regime. To overcome these limitations, we use laser-driven ultrafast impedance spectroscopy (LUIS). LUIS resonantly excites phonons using a THz field and probes ion hopping with picosecond time resolution. We apply LUIS to understand the dynamical role of phonons in $Li_7La_3Zr_2O_{12}$ (LLZO). When in its cubic phase (c-LLZO), this garnet-type solid electrolyte has an ionic conductivity two orders of magnitude greater than its tetragonal phase (t-LLZO). T-LLZO is characterized by an ordered and filled $Li^+$ sublattice necessitating synchronous ion hopping. In contrast, c-LLZO is characterized by a disordered and vacancy-rich $Li^+$ sublattice, and has a conduction mechanism dominated by single hops. We find that, upon excitation of phonons in the 0.5-7.5 THz range, the impedance of t-LLZO experiences a longer ion hopping decay signal in comparison to c-LLZO. The results suggest that phonon-mediated ionic conduction by THz modes may lead to larger ion displacements in ordered and fully occupied mobile ion sublattices as opposed to those that are disordered and vacancy-rich. Overall, this work highlights the interplay between static and dynamic factors that enables improved ionic conductivity in otherwise poorly conducting inorganic solids.
Reducing the formation temperature of single-phase multioxides is one of the classic challenges in ceramic processing, including wet-chemical synthesis routes. Toward pursuing this aim for diopside (MgCaSi2O6), the merit of different sol-gel and coprecipitation processes using the related chloride precursors followed by calcination was compared from the viewpoints of crystallinity and homogeneity. In accordance to the results, the use of the sol-gel techniques, directed with/without an alkaline catalyst, gave rise to the unfavorable creation of multiphase and low-crystallinity structures. Regarding the coprecipitation methods, the one-step addition of a precipitant agent is accompanied by an indirect low-temperature formation of nano-diopside, while a direct crystallization into this phase was explored in the dropwise condition, albeit with a lower crystallinity. Thus, by employing a suitable synthesis processing, it is feasible to take control of a wide range of nanoparticulate diopside-based structures achieved after a low-temperature calcination.