2026-03-03 | | Total: 158
Heterostructures composed of transition metal dichalcogenides (TMDCs) and organic molecules have been extensively explored for optoelectronic devices. To maximize their application potential, it is essential to investigate the electronic band structures, which govern the charge response of the interfaces to external perturbations. Based on $GW$ calculations, we present a study of organic-inorganic heterostructures with anthracene molecules adsorbed on monolayer MoS2. Building on previous investigations of organic molecule self-assembly at surfaces, we systematically analyze anthracene configurations with various molecular orientations and surface coverages. Partially self-consistent $GW_0$ provides qualitatively different level alignments from those in DFT. Whereas the systems with sparse, horizontally adsorbed anthracenes exhibit type-I alignment, densely packed anthracenes in the head-on position lead to type-II alignment, which indicates the strong dependence of quasiparticle corrections on the interfacial configuration. These findings highlight the importance of level-alignment predictions for both interpreting experiments and guiding the design of organic-inorganic heterostructures.
Investigating the initial conditions that lead many-body quantum systems to an out-of-equilibrium state is fundamental for understanding their thermalization dynamics. In this work we observe the relaxation for two regimes of excitation that can drive the turbulent Bose-Einstein condensate into two distinct final states, and are defined by the amount of energy injected into the system. The subcritical regime is characterized by a lower injection of energy, which can lead to an inverse particle cascade and, consequently, to the BEC mode repopulation during the relaxation process. The supercritical regime is marked by a higher energy injection, that may lead to the BEC dissolution and a final thermal state. In both cases we observe relaxation stages that exhibit the same key features: a direct cascade, a non-thermal fixed point with the same exponents, a prethermalization region and, finally, the thermalization of the system. In the final thermalization stage, universal scaling is observed for both regimes, even though their final states are completely different. By analyzing the coherence length of our turbulent cloud, we clearly visualize the recovery and the loss of the coherence for the subcritical and supercritical regimes after relaxation. These results indicate that the evolution of turbulence occurs independent of its initial conditions and of the final state achieved.
Molecules can form myriad crystalline polymorphs, each with distinct properties affecting their performance across diverse applications, from pharmaceuticals to functional materials and more. Predicting the thermodynamically most stable polymorph from first principles remains a formidable challenge. It requires methods that scale to large, technologically-relevant molecules while achieving very high accuracy (below 1 kJ/mol) on relative lattice energies. Such accuracy, often termed sub-chemical accuracy, is generally beyond the reach of the workhorse density functional theory (DFT). In this work, we introduce a framework, combining advances in correlated wavefunction theory (cWFT) and the many-body expansion, to deliver accurate, cost-effective predictions of complex molecular crystals. For 23 organic molecules and 13 ice polymorphs, we predict crystal lattice energies to within experimental uncertainties at costs comparable to hybrid DFT, while being several orders of magnitude more efficient than previous cWFT approaches. We extend this approach to a set of large, drug-like molecules including axitinib and ROY, previously inaccessible to cWFT and where DFT is insufficient, achieving sub-chemical accuracy on the relative energies between challenging polymorphs. With the reference data generated throughout this work, we have been able to further parametrize a DFT functional with unprecedented accuracy aligning with our predictions. This cWFT framework as well as DFT functional are made openly available, providing new ranking tools to facilitate efficient high-throughput screening of molecular crystal polymorphs.
A profound understanding of terahertz (THz) spin and charge currents in heterostructures involving ferromagnets (FMs) and two-dimensional (2D) materials promises emerging applications in high-speed sensing and data processing. Yet, ultrafast experimental insights remain very limited. Here, we study the efficient photo-generation of THz spin and charge currents in bilayers made from the transition-metal dichalcogenide (TMD) MoS2 and the FM Co. We find that the efficiency of current generation strongly depends on the pump photon energy, as previously reported. Surprisingly, however, we observe that the current dynamics remain identical for pump photon energies above and below the MoS2 band gap. Supported by ab-initio calculations, we conclude that an interfacial hybrid metallic layer forms at the MoS2/Co boundary that has a pronounced photon-energy-dependent absorptance. Thus, the hybrid interfacial layer effectively acts like a pump-energy transducer that increases the spin-current generated in the nearby Co. Our results uncover the vital role of interfacial hybridization as a yet unexplored mechanism for efficient generation of ultrafast photocurrents in 2D-TMD|FM structures.
We introduce a multi-species generalization of the hard-rod gas in which each species has a distinct effective length, and the repulsive scattering shift is set by the smaller of the two colliding rods. We argue that this model shares key quasiparticle and scattering features with the XXZ spin chain. We show that fixing only the functional decay of bare velocities with rod length is sufficient to reproduce the XXZ spin-transport phase diagram: diffusion (with anomalous fluctuations) in the anisotropic regime and superdiffusion at the isotropic point. We then demonstrate that the statistics of charge transfer differs qualitatively from that of particle trajectories. For long rods, trajectories are Gaussian in the diffusive regime and appear to exhibit KPZ statistics at the isotropic point, providing a direct microscopic signature of KPZ physics in integrable quasiparticle motion. In contrast, charge-transfer fluctuations are anomalous in the anisotropic regime, while they cross over to Gaussian statistics at late times at the isotropic point, reconciling non-Gaussian trajectory fluctuations with Gaussian charge-transfer statistics. Our results establish classical hard-rod dynamics as a minimal yet quantitatively faithful framework for anomalous spin and charge transport in integrable systems, and offer new insight into the origin of KPZ fluctuations in isotropic integrable models.
We derive an exact generalization of the well-known Lebowitz--Percus--Verlet (LPV) formula that relates the kinetic energy fluctuations of an isolated system to its specific heat. Our general formula, obtained by the application of expectation identities, is valid for arbitrary steady--state ensembles and system sizes, expressing the relative variance of the kinetic energy in terms of the variance of total energy and the microcanonical specific heat. The usual microcanonical LPV formula can be readily recovered as a particular case where energy fluctuations vanish. We test the validity of the generalized formula by performing Monte Carlo simulations of a superstatistical system of harmonic oscillators, as well as by exact calculation of energy variances in a uniform--energy ensemble, discussing its relevance to systems exhibiting negative heat capacity and ensemble inequivalence, as encountered in finite nuclei and self--gravitating models. Our results may provide useful in the study of non-equilibrium phase transitions in finite systems.
Gyromorphs are a new class of disordered systems that combine an amorphous-like absence of translational order with quasi-long-range rotational order. Gyromorphs can outperform quasicrystals or hyperuniform arrangements in forming isotropic band gaps, suggesting an avenue to realize robust disordered topological phases. However, gyromorphs lack exact rotational symmetry, which is only realized on average, posing an obstacle for existing real-space invariants to correctly diagnose topological gyromorphs. In this work we show that gyromorphs can host higher-order topological insulating (HOTI) phases protected by average rotational symmetry, and we develop and systematically compare tools for diagnosing topological phases protected by such symmetry. We introduce symmetry indicators of the effective Hamiltonian based on average rotational symmetries which, when combined with the spectral localizer and a scattering invariant, draw a consistent topological phase diagram. Our work unlocks gyromorphs as a novel platform to study topological phases beyond crystals, quasicrystals, and amorphous materials.
Upon breaking the integrability, the equations of generalized hydrodynamics (GHD) are supplemented by a Boltzmann collision term. Such terms are typically complicated and stem from a perturbative treatment of integrability-breaking terms in the hamiltonian. In our work, we study a simplified version of the collision operator in a form of relaxation time approximation familiar from kinetic theory. We explicitly compute transport coefficients which characterize the Navier-Stokes (NS) hydrodynamic regime emerging at large space-time scales. We also thoroughly study the crossover between GHD and NS hydrodynamic descriptions, identifying relevant characteristic space-time scales for the transition. In particular, we show how the emergence of NS hydrodynamics is visible in dynamics of conserved and non-conserved charge densities, and in hydrodynamic two-point functions.
We demonstrate how the beyond-mean-field Lee-Huang-Yang (LHY) corrections and its related physics can be naturally incorporated into the representation of an ultracold Bose gas using the truncated Wigner approach without invoking effective energy terms or local density assumptions. By generating a Bogoliubov ground-state representation with appropriately tailored bare interaction strength $g_0$ and condensate density $n_0$, the expected initial energy and densities are obtained while retaining access to quantum effects beyond the reach of the extended Gross-Pitaevskii equation (EGPE) formulation. This approach enables the study of correlations, coherence decay, single realisations, and the onset of quantum fluctuation effects with growing interaction strength. Numerical demonstrations for a weakly interacting single-component Bose gas show that observables deviate significantly from both the plain GPE and the EGPE incorporating LHY corrections. In regimes of strong interaction, many of the interference effects predicted by the GPE and EGPE suppressed, and the EGPE offers no improvement over the plain GPE compared to the full Wigner model. In the weakly interaction limit, the EGPE appears accurate but resolving its deviation from mean-field results requires extensive ensemble averaging.
Anisotropic diffusion governs transport in a wide range of soft and biological materials, where microstructure and molecular interactions jointly shape how matter moves. Here, we quantitatively investigate anisotropic molecular transport in lyotropic chromonic liquid crystals (LCLCs) using fluorescence recovery after photobleaching (FRAP). Disodium cromoglycate (DSCG) serves as a model LCLC system, and diffusion is measured across isotropic, nematic, and columnar phases as concentration and temperature are varied. To disentangle the roles of microstructure and molecular interactions, we employ two fluorescent tracers with distinct affinities for the LCLC aggregates: Acridine Orange (AO), which intercalates into DSCG aggregates, and Bodipy, which interacts weakly and remains largely in the aqueous phase. Fourier-space FRAP analysis independently resolves the parallel and perpendicular diffusion coefficients for both dyes relative to the liquid-crystal alignment. In the nematic phase, diffusion becomes anisotropic, with faster transport along the liquid-crystal director. As the DSCG concentration increases, AO dye molecules that are strongly coupled to the aggregates exhibit a slowdown in all directions, reflecting enhanced packing and steric confinement of the LC microstructure. In contrast, weakly interacting Bodipy dye molecules display enhanced transport along the alignment direction as the DSCG concentration increases in the nematic regime, suggesting the emergence of microscopic channels that guide motion, analogous to transport in oriented porous media. These results reveal how the evolving microstructure of LCLCs controls effective diffusion and provide a quantitative framework for understanding and designing anisotropic transport in aligned soft materials.
In this work we provide an in-depth analysis of the sensing mechanisms of $NO_{2}$ by lead-sulfide nanocrystals (PbS-NCs). A detailed model for the sorption mechanism is proposed, and the correlation is established between experimental sensing characteristics and the surface composition, based on both experimental characterization and ab initio (DFT) simulations. We demonstrated how the sensitivity and the sensing dynamic response can be tuned by a post-deposition multistep dry-thermal process at mild temperature, that alternates vacuum-assisted annealing and heating in open-air. Sensors with different surface compositions were fabricated, and their dynamic response was characterized at low concentration of $NO_{2}$ (0.5 ppm) in air, at ambient temperature. DFT simulations indicate that both surface stoichiometry and oxidation critically govern $NO_{2}$ interaction on PbS, with sulfur-rich terminations favoring weaker binding and faster desorption, while intermediate oxidation enhances interaction and overoxidation leads to surface passivation, in agreement with the measured experimental sensing dynamics. By linking surface composition, adsorption chemistry, and resistance transduction within a single framework, this work provides clear indications to design room-temperature, low-ppm $NO_{2}$ microsensors fabricated through a simple and scalable processes.
We study the topological phase transition in a two-leg Su-Schrieffer-Heeger (SSH) ladder by redefining the unit-cell structure. For both identical hopping dimerization pattern (uniform) and alternate hopping dimerization pattern (staggered), we demonstrate that different unit-cell choices generate different topological phases and phase transitions. In the uniformly dimerized ladder, variation of the inter-leg coupling induces transitions between distinct topological phases mediated by a gapless region. In contrast, the staggered dimerization configuration exhibits a richer phase structure, supporting both topological-topological and trivial-topological transitions occurring through a single gap-closing point, depending on the unit-cell definition. The phases are characterized through bulk-boundary correspondence, edge-state analysis, and bulk invariant. These analyses reveal interesting phenomena, including regimes supporting four localized boundary modes in the weak inter-leg coupling limit and two boundary modes in the strong inter-leg coupling limit. We further perform an experimental implementation using topoelectric circuit simulation by mapping the lattice Hamiltonian to an equivalent LC circuit. Circuit impedance and voltage responses obtained from LTspice simulations confirm the predicted edge-state signatures. These results establish topolectrical platforms as experimentally accessible realizations of tunable lattice topological phases.
This work presents a novel approach to improve phase identification in nanomechanical property maps of multiphase materials, such as those obtained by nanoindentation or atomic force microscopy (AFM). A major challenge in validating clustering strategies for these data is the lack of ground-truth phase labels in experimental measurements, along with the tendency of overly simplistic synthetic datasets to artificially inflate algorithm performance. To address this gap, we construct controlled yet non-trivial synthetic benchmarks with tunable mechanical contrast, graded interfaces, curved boundaries, and diffuse morphologies, enabling rigorous and realistic evaluation of clustering robustness. Conventional clustering methods based only on elastic modulus (E) and hardness (H) often struggle to separate phases when mechanical contrast is low or when diffuse interphase regions are present. We introduce the Kernel-Averaged Mechanical Mismatch (KAMM), a neighborhood-informed feature that quantifies local mechanical heterogeneity by comparing each point to its neighbors in (E, H) space. When incorporated into a three-dimensional clustering space (E, H, KAMM), this framework improves phase separability, enhances interphase detection, and increases robustness to noise. By enabling more reliable segmentation of mechanical domains under realistic contrast conditions, the proposed method facilitates the generation of representative volume elements (RVEs) and supports more accurate extraction of phase-specific properties in heterogeneous microstructures.
The quality, consistency, and information content of training data is often what determines the practical value of machine-learning models for atomistic simulations. Yet, many widely used electronic-structure databases are assembled having materials screening as primary goal rather than robust force-field learning, are limited in their scope to a specific class of chemical compounds, and/or employ inconsistent DFT functionals and settings. Here we introduce MAD-1.5, a highly curated dataset designed explicitly for training broadly applicable atomistic models across the periodic table at high levels of theory. MAD-1.5 extends the MAD dataset with targeted enrichment strategies that improve the coverage of chemical space to 102 elements while keeping the total number of configurations compact. All structures are computed with a single, standardized all-electron DFT workflow using the r$^2$SCAN meta-GGA functional and consistent convergence settings, ensuring uniformity across chemically heterogeneous systems. The dataset encompasses molecules, clusters, bulk crystals, surfaces, and low-dimensional structures, and its quality and consistency are further enhanced by outlier removal using uncertainty quantification. We demonstrate the high accuracy that can be achieved with the proposed dataset by training PET-MAD-1.5, a generally applicable r$^2$SCAN interatomic potential that covers 102 elements in the periodic table and achieves exceptional levels of benchmark accuracy and stability in challenging simulation protocols.
We generalize the scalable tight-binding model for graphene, which allows for efficient quantum transport simulations in the Dirac regime, to account for elastic strain. We show that the original scalable model with scaling factor $s$ is readily applicable to strained graphene, provided that the displacement fields corresponding to the deformed graphene lattice are properly scaled. In particular, we show that the long-wavelength theory remains invariant when the strain tensor is scaled by $s$. This is achieved in practice by scaling the in-plane displacement fields by $s$ while the out-of-plane displacements have to be scaled by $\sqrt{s}$. We confirm these scaling laws by extensive numerical simulations, starting with the pseudomagnetic field and the local density of states for different scaled lattices. The latter allows us to study pseudo-Landau levels as well as hybrid Landau levels in the presence of an external magnetic field. Finally, we consider quantum transport simulations motivated by a recent experiment, where a uniaxial strain barrier is engineered in monolayer graphene by vertically misaligned gates. Our work generalizes the scalable tight-binding model to allow for efficient modeling of quantum transport in large-scale strained graphene devices.
Excitons in anisotropic two-dimensional (2D) materials, defined by direction-dependent effective masses, are of pronounced interest for their roles in excitonic and magneto-optical phenomena. A perpendicular magnetic field complicates the separation of center-of-mass (c.m.) and relative motions, especially when electron and hole masses are comparable. Conventional theories often employ an approximate c.m. separation using factorized wave functions, modifying magnetic Hamiltonian terms and possibly introducing inaccuracies in magnetoexciton energy predictions. This work develops an exact analytical framework for c.m. and relative motion separation in anisotropic 2D magnetoexcitons, without resorting to the stationary-c.m. approximation. Starting from the full electron-hole Hamiltonian in a homogeneous magnetic field, the formalism uses the conserved pseudomomentum to derive a relative-motion Hamiltonian, revealing new anisotropy-dependent couplings and magnetic coefficients absent in approximate models. The resulting Schrödinger equation is treated via the Feranchuk-Komarov operator method and Levi-Civita transformation, allowing non-perturbative, systematically convergent solutions. Application to monolayer black phosphorus and titanium trisulfide, both freestanding and encapsulated in hexagonal boron nitride, yields magnetoexciton energies, diamagnetic coefficients, and probability densities for the ten lowest states across considerable magnetic-field ranges. The results demonstrate the significant influence of anisotropy-dependent coupling on magnetic response in systems with strong mass anisotropy. This formalism is generalizable to other anisotropic 2D semiconductors, establishing a foundation for advanced magneto-optical studies.
The morphology of DNA is strongly influenced by its surrounding environment, including factors such as pH, salt type and valency, and the presence of polymers. Inorganic salts are known to reduce the DNA chain length through mechanisms like electrostatic screening and ion bridging. In contrast, ionic liquids, a new class of organic salts, have previously been found to increase the DNA chain length, indicating a distinct mode of interaction between the ionic liquid and DNA chains. This study utilizes self-assembled DNA-AuNPs as a model system to examine changes in the DNA chain morphology and the nanoscale interaction mechanisms in ionic liquid environment. The DNA chain lengths are measured in solution using X-ray scattering measurements at varying concentrations of two imidazolium ([BMIM] acetate and [EMIM] acetate) based ionic liquids. Additionally, Molecular Dynamics (MD) simulations are performed mimicking the experimental system. Our results suggest an interplay of electrostatic and groove-binding interactions governing the DNA chain morphology, which depends on IL concentration and the composition of the DNA chains. It has been found that for DNA chains with majority ssDNA, electrostatic interaction dominate, however with increasing composition of double strands, the DNA chains exhibits compaction due to non-electrostatic hydrophobic groove-binding mechanism.
Molecular rings of N carbon atoms, that is, cyclo[N]carbons, or $C_N$, can be formed by tip-induced chemistry [1-7]. Because of their monocyclic geometry, cyclocarbons are fundamentally important for testing theoretical models of aromaticity [8-11]. Here, we synthesized large cyclo[N]carbons, with N up to 88, by tip-induced chemistry on a NaCl surface and studied their aromaticity by measuring their transport gaps by scanning tunnelling spectroscopy. We first generated $C_{20}$ and $C_{22}$, and then fused multiple cyclocarbons [5-7] by means of atom manipulation, obtaining $C_{42}$, $C_{44}$, $C_{46}$, $C_{66}$ and $C_{88}$. In line with theory, using a finely tuned density functional approximation [12-15], we observe a substantially smaller transport gap for $C_{20}$ (N = 4n) compared to $C_{22}$ (4n+2), and for $C_{44}$ (4n) compared to $C_{42}$ (4n+2). In larger cyclocarbons, the oscillation of the transport gap between anti-aromatic N = 4n and aromatic N = 4n+2 cyclocarbons becomes smaller, and is expected to eventually vanish with increasing N indicating non-aromaticity. Our experimental results show that aromaticity persists at N = 42, and theory predicts ring currents comparable in magnitude to that of benzene in cyclocarbons of this size. In the future, large cyclocarbons could be used as model systems to study conductance, quantum interference, and the effects of aromaticity in single atomic carbon wires and circuits.
Scalable and symmetry-consistent force-field models are essential for extending quantum-accurate simulations to large spatiotemporal scales. While descriptor-based neural networks can incorporate lattice symmetries through carefully engineered features, we show that graph neural networks (GNNs) provide a conceptually simpler and more unified alternative in which discrete lattice translation and point-group symmetries are enforced directly through local message passing and weight sharing. We develop a GNN-based force-field framework for the adiabatic dynamics of lattice Hamiltonians and demonstrate it for the semiclassical Holstein model. Trained on exact-diagonalization data, the GNN achieves high force accuracy, strict linear scaling with system size, and direct transferability to large lattices. Enabled by this scalability, we perform large-scale Langevin simulations of charge-density-wave ordering following thermal quenches, revealing dynamical scaling and anomalously slow sub--Allen--Cahn coarsening. These results establish GNNs as an elegant and efficient architecture for symmetry-aware, large-scale dynamical simulations of correlated lattice systems.
Computing excitonic spectra in quasicrystal and super-moiré systems constitutes a formidable challenge due to the exceptional size of the excitonic Hilbert space. Here, we demonstrate a tensor-network method for the real-space Bethe-Salpeter Hamiltonian, allowing us to access the spectra of an excitonic $10^{18}$-dimensional Hamiltonian, and enabling the direct computation of bound-exciton spectral functions for systems exceeding one billion lattice sites, several orders of magnitude beyond the capabilities of conventional approaches. Our method combines a tensor-network encoding of the real-space Bethe-Salpeter Hamiltonian with a Chebyshev tensor network algorithm. This strategy bypasses explicit storage of the Hamiltonian while preserving full real-space resolution across widely different length scales. We demonstrate our methodology for one- and two-dimensional super-moiré systems, achieving the simultaneous resolution of atomistic and mesoscopic structures in the excitonic spectra in billion-size systems, showing exciton miniband formation and moiré-induced spatial confinement. Our results establish a real-space methodology enabling the simulation of excitonic physics in large-scale quasicrystal and super-moiré quantum matter.
We investigate the infinite-dimensional limit of nonequilibrium surface growth by numerically integrating stochastic growth equations on a fully connected graph. In particular, we study the Edwards-Wilkinson (EW), Kardar-Parisi-Zhang (KPZ), and tensionless KPZ (TKPZ) equations. Using a network discretization adapted to the all-to-all interaction topology, we analyze the global roughness, height-fluctuation statistics, time power spectra, and two-time correlations. For the EW equation, we obtain an exact expression for the roughness that matches the numerical simulations and shows that the interface becomes flat as $N \to \infty$. We also compute analytically the time power spectrum, show that height fluctuations are Gaussian, and derive an explicit expression for the two-time height autocorrelation function, indicating that the aging behavior is trivial. For the KPZ equation, finite-size and strong-coupling effects can cause deviations from EW behavior at moderate system sizes $N$, often accompanied by numerical instabilities; however, these differences disappear as $N$ increases. In the large-$N$ limit, KPZ dynamics converges to EW behavior, as the four observables analyzed exhibit identical scaling properties. Overall, our results indicate that on a fully connected graph the KPZ nonlinearity is irrelevant as $N\to\infty$, leading to EW-like dynamics with asymptotically flat interfaces. These findings are interpreted in the context of the upper critical dimension of the KPZ universality class.
We fabricate twisted double bilayer graphene devices with zero twist angle and a set of local top and bottom gates aligned perpendicularly to each other. A 1D PN junction can be electrostatically defined when the gate voltages applied to the top gates are the same but different on the bottom gates. Resistance peaks are observed at finite doping instead of at the charge neutrality points, exhibiting an unconventional broken-cross shape that arises from layer polarization of the P and N region, which can be further enhanced with finite magnetic fields. A 0D point junction (PJ) can be electrostatically defined by applying different gate voltages to the top and bottom gates, such that the P and N sides of the device are connected at a single point in the center of the device. As finite magnetic field B increases, the quantum Hall (QH) states are selectively brought into contact or away from each other depending on their layer polarization, leading to unconventional quantum oscillations which characterize the layer-polarized band-crossing. Our work provides new insights into understanding band-structure evolution and layer polarization in twisted bilayers and paves the way for new device functionality based on manipulating layer-polarized electronic states.
We construct a Hubbard model with a nearly flat band whose quantum geometry can be tuned independently of the energy dispersion and the Coulomb interaction. We show that, when the nearly flat band is half-filled, the exact ground state of the model exhibits ferromagnetism and that this ferromagnetism is stabilized by the quantum metric through the spin stiffness. Furthermore, we demonstrate that tuning the quantum geometry alone drives a magnetic phase transition. Our nonperturbative results without resorting to mean-field approximations reveal the quantum-geometric origin of ferromagnetism and the underlying many-body physics in dispersive-band systems.
We perform Discrete Element Method (DEM) simulations of granular particles (polystyrene spheres) vibrated inside a cubic container. The study investigates the evolution of the packing fraction with and without rotational friction at different shaking amplitudes. The mean-squared displacement (MSD) is used to analyze the particles' diffusive, subdiffusive, and superdiffusive behavior. By monitoring both the dynamics and density evolution, one can observe the system's glassification. The comparison with experiments shows that the MSDs from the simulations are significantly higher than the MSDs measured by Diffusing Wave Spectroscopy (DWS) \cite{kunzner2025dynamics}. Following our finding that the rotational MSD is of the same order of magnitude as the MSD measured in DWS experiments, we propose that the experimental signal is not dominated by translational motion but rather by rotational particle dynamics. This provides access to a relevant particle property that has previously been difficult to measure directly. Finally, we conclude that the system reaches a dynamically constrained state well before random close packing, with particle displacements already below the Lindemann length.
The behavior of surface triple junctions (STJ) at emergent grain boundaries on free surfaces are critical contributors to material properties, particularly at the nanoscale. Interest has been stimulated by recent experimental observations of local grain rotation, the restructuring of emergent grain boundaries and the resultant levels of stress generation, and computationally by the introduction of the core-shift method for building wedge disclination structures at STJ and the establishment of relationships between atomic calculation and the continuum elastic theory. Despite this progress, the lowest-energy structure for the surface triple junction is still unknown. Here, we carried out the experimental observation of the local surface deformation at STJ, the local stress analysis and ergodic searches for metastable structures. We identified the zipped Y-shaped notch, rather than the zipped V-shaped notch or wedge disclination, as the most stable structures, which was subsequently verified with machine learning methods for a wide range of boundaries. Local stress analysis explains this energetic preference. Our findings provide fundamental insight into STJ energetics and offer a possible framework for understanding kinetically diffusive deformation and the pathway for engineering related material properties, in particular thin-film properties.