2024-11-06 | | Total: 8

In materials and pharmaceutical development, rapidly and accurately determining the similarity between X-ray powder diffraction (XRPD) measurements is crucial for efficient solid form screening and analysis. We present SMolNet, a classifier based on a Siamese network architecture, designed to automate the comparison of XRPD patterns. Our results show that training SMolNet on loss functions from the self-supervised learning domain yields a substantial boost in performance with respect to class separability and precision, specifically when classifying phases of previously unseen compounds. The application of SMolNet demonstrates significant improvements in screening efficiency across multiple active pharmaceutical ingredients, providing a powerful tool for scientists to discover and categorize measurements with reliable accuracy.

Gas-phase molecules capable of repeatable, narrow-band spontaneous photon scattering are prized for direct laser cooling and quantum state detection. Recently, large molecules incorporating phenyl rings have been shown to exhibit similar vibrational closure to the small molecules demonstrated so far, and it is not yet known if the high vibrational-mode density of even larger species will eventually compromise optical cycling. Here, we systematically increase the size of hydrocarbon ligands attached to single alkaline-earth-phenoxides from (-H) to -C$_{14}$H$_{19}$ while measuring the vibrational branching fractions of the optical transition. We find that varying the ligand size from 1 to more than 30 atoms does not systematically reduce the cycle closure, which remains around 90%. Theoretical extensions to larger diamondoids and bulk diamond surface suggest that alkaline earth phenoxides may maintain the desirable scattering behavior as the system size grows further, with no indication of an upper limit.

We explore the applicability of the transcorrelated method to the elements in the second row of the periodic table. We use transcorrelated Hamiltonians in conjunction with full configuration interaction quantum Monte Carlo and coupled cluster techniques to obtain total energies and ionisation potentials, investigating their dependence on the nature and size of the basis sets used. Transcorrelation accelerates convergence to the complete basis set limit relative to conventional approaches, and chemically accurate results can generally be obtained with the cc-pVTZ basis, even with a frozen Ne core in the post-Hartree--Fock treatment.

Porphyrins offer a malleable and cost-efficient platform to sculpt bioinspired technologies with tunable charge transfer, energy conversion, and photocatalytic properties. Yet, despite decades of research, the physical mechanisms that determine their electronic spectra remain elusive. Even for metal-free porphyrins, no consensus exists on the origin of the splitting of their $Q$-bands, an energetic region critical in photosynthesis. We leverage our recent statistical treatment of molecular motions in the condensed phase to predict their linear absorption spectra. By bridging exact quantum-dynamical expressions with atomistic simulations, our theory is the first to capture the spectra of representative porphyrins in solution: porphine, tetraphenylporphyrin, and tetraphenylporpholactone. Our work reveals that $Q$-band splitting arises from extreme timescale separation of nuclear motions that turn symmetry-forbidden transitions bright. We exploit these insights to propose and confirm how chemical modifications tune the optical properties of porphyrins, demonstrating the potential for developing design principles to tailor their optoelectronic behavior.

Recent experiments have realized exciton condensation in bilayer materials such as graphene double layers and the van der Waals heterostructure MoSe$_2$-WSe$_2$ with the potential for nearly frictionless energy transport. Here we computationally observe the microscopic beginnings of exciton condensation in a molecular-scale fragment of MoSe$_2$-WSe$_2$, using advanced electronic structure methods based on reduced density matrices. We establish a connection between the signature of exciton condensation -- the presence of a large eigenvalue in the particle-hole reduced density matrix -- and experimental evidence of exciton condensation in the material. The presence of a "critical seed" of exciton condensation in a molecular-scale fragment of a heterostructure bilayer provides insight into how local short-range strongly correlated effects may give rise to macroscopic exciton condensation. We find that molecular-scale properties such as layer alignment and interlayer distance can impact the formation of nonclassical long-range order in heterostructure bilayers, demonstrating the importance of geometric considerations for the rational design of exciton condensate material. Mechanistic insights into the microscopic origins of exciton condensation have potential implications for the design and development of new materials with enhanced energy transport properties.

The structure of the low-lying rotational levels of the X$^+$ $ ^2 \Sigma_g ^+$ ($v^+=0$) vibronic ground state of para-H$_2^+$ has been determined by combining frequency-comb calibrated continuous-wave high-resolution laser spectroscopy of $n$f Rydberg series in the range of principal quantum number $n$ between 28 and 115 and Rydberg-series extrapolation using multichannel-quantum-defect theory (MQDT). The use of accurate quantum-defect parameters obtained from new ab initio calculations enabled the experimental determination of the pure rotational term values of the $N^+= 2$, 4 and 6 rotational levels of H$_2^+$ with sub-MHz accuracy (174.236\,744\,6(77), 575.455\,632\,5(86) and 1191.385\,571(240) cm$^{-1}$, respectively), and of the corresponding spin-rotational coupling constants with an accuracy of better than 100 kHz (42.21(4), 41.26(8) and 40.04(8) MHz, respectively). These values are in agreement with the results of first-principles calculations that include high-order relativistic and quantum-electrodynamics corrections to the level energies. To reach the reported accuracy in the Rydberg series extrapolation, it was necessary to correct for artificial level shifts arising in the MQDT calculations in the vicinity of local perturbations of high-$n$ Rydberg states with a $v^+=0$ H$_2^+$ ion core caused by low-$n$ core-excited Rydberg states, and resulting from approximations in the treatment of the Rydberg-electron energy in the interacting channels.

The mechanisms underlying chemifriction, i.e. the contribution of interfacial bonding to friction in defected twisted graphene interfaces are revealed using fully atomistic machine-learning molecular dynamics simulations. This involves stochastic events of consecutive bond formation and rupture, that are spatially separated but not necessarily independent. A unique shear-induced interlayer atomic transfer healing mechanism is discovered that can be harnessed to design a run-in procedure to restore superlubric sliding. This mechanism should be manifested as negative differential friction coefficients that are expected to emerge under moderate normal loads. A physically motivated phenomenological model is developed to predict the effects of chemifriction in experimentally relevant sliding velocity regimes. This allows us to identify a distinct transition between logarithmic increase and logarithmic decrease of frictional stress with increasing sliding velocity. While demonstrated for homogeneous graphene interfaces, a similar mechanism is expected to occur in other homogeneous or heterogeneous defected two-dimensional material interfaces.

We present KinetiX, a software toolkit to generate computationally efficient fuel-specific routines for the chemical source term, thermodynamic and mixture-averaged transport properties for use in combustion simulation codes. The C++ routines are designed for high-performance execution on both CPU and GPU architectures. On CPUs, chemical kinetics computations are optimized by eliminating redundant operations and using data alignment and loops with trivial access patterns that enable auto-vectorization, reducing the latency of complex mathematical operations. On GPUs, performance is improved by loop unrolling, reducing the number of costly exponential evaluations and limiting the number of live variables for better register usage. The accuracy of the generated routines is checked against reference values computed using Cantera and the maximum relative errors are below $10^{-5}\%$. We evaluate the performance of the kernels on some of the latest CPU and GPU architectures from AMD and NVIDIA, i.e., AMD EPYC 9653, AMD MI250X, and NVIDIA H100. The routines generated by KinetiX outperform the general-purpose Cantera library, achieving speedups of up to 2.4x for species production rates and 3.2x for mixture-averaged transport properties on CPUs. Compared to the routines generated by PelePhysics (CEPTR), KinetiX achieves speedups of up to 2.6x on CPUs and 1.7x on GPUs for the species production rates kernel on a single-threaded basis.