Papers

The interplay between strong electron correlation and band topology offers a playground for discovering exotic quantum phases. Here, we predict the emergence of topological electron crystals in a charge-transfer bilayer graphene-Mott insulator heterostructure. In this system, interlayer charge transfer induces a charge-neutral electron-hole bilayer with strong mass asymmetry. While the extreme dilute limit favors a classical triangular dipolar Wigner crystal, we show that increasing the carrier density triggers a critical competition between the Coulomb interaction and the underlying topological band structure of bilayer graphene. This interplay destabilizes the triangular dipolar Wigner crystal and instead stabilizes intrinsic quantum electron crystals with spontaneously formed honeycomb and kagome geometries. Crucially, these new phases can host distinct topological responses including the quantum anomalous and quantum spin Hall effects, which inherit the nonlocal quantum geometry of the bilayer graphene wave functions. Our results establish this artificial heterostructure as a highly tunable platform for mimicking two-dimensional topological solids in a single device.

We systematically investigate the emergence of electron crystal phases in rhombohedral multilayer graphene using comprehensive self-consistent Hartree Fock calculations combined with ab initio tight binding model. As the carrier density increases, we uncover an isospin cascade sequence of phase transitions that gives rise to a rich variety of ordered states, including electron crystal phases with non-zero Chern numbers. We further show the nearly degeneracy of these topological electron crystals hosting extended quantum anomalous Hall effect (EQAH) in the mean field regime and characterize pressure driven phase transitions. Finally, we discuss the thermodynamic signatures, particularly the behavior of the inverse compressibility, in light of recent experimental observations.

Titanium MXenes are two-dimensional inorganic structures composed of titanium and carbon or nitrogen elements, with distinctive electronic, thermal and mechanical properties. Despite the extensive experimental investigation, there is a paucity of computational studies at the level of classical molecular dynamics (MD). As demonstrated in a preceding study, known MD potentials are not capable of fully reproducing the structure and elastic properties of every titanium MXene. In this study, we present a simply trained, but yet efficient, non-density functional theory-based machine learning interatomic potential (MLIP) capable of simulating the structure and elastic properties of titanium MXenes and bulk titanium carbide and nitride with precision comparable to DFT calculations. The training process for the MLIP is delineated herein, in conjunction with a series of dynamical tests. Limitations of the MLIP and steps towards improving its efficacy to simulate titanium MXenes are discussed.

In the present study, we apply shortcut to adiabaticity pulses (time-dependent Rabi frequency and detuning) for the efficient population transfer from the ground to the exciton state in a GaAs/InGaAs quantum dot with phonon-induced dephasing. We use the time-evolving matrix product operator (TEMPO) method to propagate system in time and find that, for temperatures below $ 20 \ \text{K} $ and pulse duration up to $ 10 \ \text{ps} $, a very good transfer efficiency is obtained in general. We explain these results using a Bloch-like equation derived from a generalized Lindblad equation, which adequately describes system dynamics at lower temperatures. For higher temperatures, the transfer efficiency is significantly reduced except for subpicosecond pulses, where the shortcut Rabi frequency reduces to a delta pulse attaining a fast population inversion. The present work is expected to find application in quantum technologies which exploit quantum dots for single-photon generation on demand.

We investigate the superconducting Anderson impurity model for interacting quantum dot Josephson junctions with spin-orbit coupling and a term accounting for tunnelling through higher-energy orbitals. These elements establish the conditions required for spin polarization in the absence of external magnetic field at finite superconducting phase bias. This Hamiltonian has been previously used to model the Andreev spin qubit, where quantum information is encoded in spinful odd-parity subgap states. Here we instead analyse the even-parity sector, i.e., the Andreev pair qubit based on Andreev bound states (ABS). The model is solved using the zero-bandwidth approximation and the numerical renormalization group, with further insight from variational calculations. Electron-electron interaction admixes single-occupancy Yu-Shiba-Rusinov (YSR) components into the ABS states, thereby strongly enhancing spin transitions in the presence of spin-orbit coupling. The ABS states can thus become sensitive to local magnetic field fluctuations, which has implications for decoherence in Andreev pair qubits. For strong interaction $U$, especially in the cross-over region between the ABS and YSR regimes for $U \sim 2\Delta$, charge, spin, and inductive transitions can all become strong, offering avenues for spin control and quantum transduction.

Density functional theory (DFT) and machine learning potentials (MLPs) are essential for predicting and understanding materials properties, yet preparing, executing, and analyzing these simulations typically requires extensive scripting, multi-step procedures, and significant high-performance computing (HPC) expertise. These challenges hinder reproducibility and slow down discovery. Here, we introduce Masgent, an AI-assisted materials simulation agent that unifies structure manipulation, automated VASP input generation, DFT workflow construction and analysis, fast MLP-based simulations, and lightweight machine learning (ML) utilities within a single platform. Powered by large language models (LLMs), Masgent enables researchers to perform complex simulation tasks through natural-language interaction, eliminating most manual scripting and reducing setup time from hours to seconds. By standardizing protocols and integrating advanced simulation and data-driven tools, Masgent lowers the barrier to performing state-of-the-art computational methodologies, enabling faster hypothesis testing, pre-screening, and exploratory research for both new and experienced practitioners.

We present a theory of the graphene nanoslide, a fundamental device for graphene straintronics that realizes a single pseudogauge barrier. We solve the scattering problem in closed form and demonstrate that the nanoslide gives rise to a hybrid pseudogauge and electrostatic cavity in the bipolar regime, and hosts one-dimensional transverse channels. The latter can be tuned using a bottom gate between valley-chiral or counterpropagating modes, as well as one-dimensional flatbands. Hence, the local density of states near the barrier depends strongly on the gate voltage with a tunable sublattice and electron-hole asymmetry. In the presence of electron-electron interactions, the nanoslide allows for in-situ tuning between a chiral and ordinary Tomonaga-Luttinger liquid.

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.

Inspired by previous works on epidemic-like processes in open quantum systems, we derive an elementary quantum epidemic model that is simple enough to be studied via Quantum Jump Monte Carlo simulations at reasonably large system sizes. We show how some weak symmetries of the Lindblad equation allow us to map the dynamics onto a classical Kinetic Monte Carlo; this simplified, effective dynamics can be described via local stochastic jumps coupled with a local deterministic component. Simulations are then used to reconstruct a phase diagram which displays stationary features completely equivalent to those of completely classical epidemic processes, but richer dynamics with multiple, recurrent waves of infection.

Topological magnetic textures are particle-like spin configurations stabilized by competing interactions. Their formation is commonly attributed to fluctuation-driven, first-order nucleation processes requiring activation over a topological energy barrier. Here, we demonstrate an alternative barrier- and fluctuation-free pathway for nucleating topological magnetic textures, triggered in our experiments by an excitation-induced spin reorientation transition. By combining x-ray imaging, scattering and micromagnetic simulations, we show that the system follows a deterministic cascade of symmetry-breaking phase transitions after excitation. First, the system undergoes a second-order phase transition from a homogeneous state to weak stripe domains, then a first-order transition to topologically trivial bubbles, and finally a topological switching event into skyrmionic textures. Through simulations, we generalize our findings and demonstrate that this pathway is active in a vast range of low-anisotropy materials. This previously unrecognized, spontaneous transition pathway suggests strategies for rapid, low-energy generation of topological spin textures and points to a general role of intrinsic modulational instabilities in phase transitions beyond magnetism.

We study the average bipartite entanglement entropy of Haar-random pure states in quantum many-body systems with global $\mathrm{SU}(2)$ symmetry, constrained to fixed total spin $J$ and magnetization $J_z = 0$. Focusing on spin-$\tfrac12$ lattices and subsystem fractions $f < \frac{1}{2}$, we derive a asymptotic expression for the average entanglement entropy up to constant order in the system volume $V$. In addition to the expected leading volume law term, we prove the existence of a $\frac{1}{2}\log V$ finite-size correction resulting from the scaling of the Clebsch-Gordon coefficients and compute explicitly the $O(1)$ contribution reflecting angular-momentum coupling within magnetization blocks. Our analysis uses features of random matrix ensembles and provides a fully analytical treatment for arbitrary spin densities, thereby extending Page type results to non-Abelian sectors and clarifying how $\mathrm{SU}(2)$ symmetry shapes average entanglement.

Understanding the mechanisms of hydrogen embrittlement (HE) is essential for advancing next-generation high-strength steels, thereby motivating the development of highly accurate machine-learning interatomic potentials (MLIPs) for the Fe-H binary system. However, the substantial computational expense associated with existing MLIPs has limited their applicability in practical, large-scale simulations. In this study, we construct a new MLIP within the Performant Implementation of the Atomic Cluster Expansion (PACE) framework, trained on a comprehensive HE-related dataset generated through a concurrent-learning strategy. The resulting potential achieves density functional theory-level accuracy in reproducing a wide range of lattice defects in alpha-Fe and their interactions with hydrogen, including both screw and edge dislocations. More importantly, it accurately captures the deformation and fracture behavior of nanopolycrystals containing hydrogen-segregated general grain boundaries-phenomena not explicitly represented in the training data. Despite its high fidelity, the developed potential requires computational resources only several tens of times greater than empirical potentials and is more than an order of magnitude faster than previously reported MLIPs. By delivering both a high-precision and computationally efficient potential, as well as a generalizable methodology for constructing such models, this study significantly advances the atomic-scale understanding of HE across a broad range of metallic materials.

Cellular organization and mechanotransduction pathways are crucial regulators of tissue morphogenesis, whereas their dysregulation contributes to pathologies. Overactive myofibroblasts are key drivers of fibrosis, yet how their presence alters collective cellular ordering remains unclear. Owing to steric interactions, elongated cells exhibit local order. Topological defects, where alignment is disrupted, have been postulated to serve as mechanical centers. In this study, we examine how incorporating slower moving myofibroblast phenotype impacts defect relaxation in monolayers consisting of co-cultured fibroblasts and myofibroblasts. In this system, myofibroblasts act as the less active component. Increasing their fraction increases disorder strength and slows defect recombination. On microgrooved surfaces, higher myofibroblast concentrations lead to worse alignment, suggesting single-cell mechanosensing and cell-cell interactions act jointly. Furthermore, we found that myofibroblasts preferentially localize at negatively charged -1/2 defects, whereas fibroblasts localize at +1/2 defects. Consequently, the slowdown of recombination dynamics can be partially attributed to myofibroblasts' preferential association with the less mobile -1/2 defects, increasing local friction and impeding defect mobility. This localization may also reduce compressive stress on myofibroblasts, as indicated by immunofluorescence of a downstream mechanotransducer. This work provides insights into possible connections between topological defects and cell motility in mixed phenotype monolayers.

Mesoscopic models of the optical response of metals have emerged as fundamental building blocks in quantum plasmonics, in principle overcoming the computational bottlenecks of ab initio techniques by implementing aspects of the atomistic description of the metal in otherwise classical calculations. Nonetheless, even these approaches are eventually hindered by demanding computations due to sophisticated material response. Here, this issue is addressed for the advanced Self-Consistent Hydrodynamic Drude Model (SC-HDM), which captures both nonlocal electron dynamics and electron spill-out, through a Volume Integral Equation (VIE) method. Adopting an IE-based method shifts perspective from the commonly employed Differential Equation (DE)-based ones, demonstrating significant computational efficiency. The VIE approach is a valuable methodological scaffold: It addresses SC-HDM and simpler models, but can also be adapted to more advanced ones. For spherical nanoparticles (NPs), using the inherent symmetries, similar performance for three increasingly complicated material models is achieved, breaking the taboo that increased sophistication in material response requires taxing simulations. Mesoscopic material-response functions can be readily extracted from the VIE implementation, thus circumventing the need for lengthy microscopic calculations. This method opens a new way of modeling quantum hydrodynamic NPs and will serve as essential benchmarking tool for recipes addressing more complicated geometries.

We report the influence of Ce$^{3+}$ substitution on the magnetic structures and phonon dynamics in SmCrO$_3$ perovskites. Magnetic landscapes are spanned by long-range canted anti-ferromagnetism, AFM with Neel temperatures $\sim$196 K accompanied by spin-reorientation transitions, T$_{SRPT}$ at 42 K. In Sm$_{0.9}$Ce$_{0.1}$CrO$_3$ (SCCO), Ce$^{3+}$ substitution at Sm$^{3+}$ sites transform the weak ferromagnetic (FM) $\Gamma_4$ state into robust AFM $\Gamma_1$ configuration through a gradual crossover. Such coexistence of magnetic spin configurations $\Gamma_1$(AFM) to $\Gamma_4$ (WFM) results in the enhanced high coercive field and a pronounced exchange bias-field, $H_{\mathrm{EB}}$ $\sim$ 2 kOe. Spin-only driven magneto-crystalline anisotropy of Cr$^{3+}$ and spin-orbit driven magnetic moment in Sm$^{3+}$, and Ce$^{3+}$ exhibits spin-phonon coupling through $A_{1g}(6)$ mode in SCCO are consistent with the temperature dependent spectral features of the isostructural magnetic systems and quite significant in SCCO which is in accordance with the higher structural distortion in SCCO. These results demonstrate that site-specific R$^{3+}$ substitution modulates lattice distortions, spin-phonon coupling, and spin-orbit interactions, offering pathways to optimize perovskites for diverse spintronic applications.

Information flow between subsystems is a central concept in information thermodynamics, which provides the second-law-like inequalities for subsystems. This paper discusses the geometric decomposition of information flow, which was introduced for Markov jump systems [Physical Review Research, 8, 023292 (2026)], and applies it to overdamped Langevin systems. For overdamped Langevin systems, the geometric decomposition of information flow into excess and housekeeping contributions is related to the conventional definition of the 2-Wasserstein distance between marginal distributions in optimal transport theory. This formulation offers an optimal-transport interpretation of subsystem dynamics, and this optimal-transport formulation is simpler for overdamped Langevin systems than for general Markov jump systems. It is also possible to handle features that are specific to overdamped Langevin systems, such as representations based on the Koopman mode decomposition, as well as their relationship with the Fisher information matrix. As with the results for Markov jump systems, we generalize the second law of information thermodynamics using housekeeping and excess information flow, leading to the concepts of excess and housekeeping demons. We also derive a thermodynamic uncertainty relation and an information-thermodynamic speed limit incorporating excess information flow. These results are illustrated for the Gaussian case, and we discuss the conditions under which the excess and housekeeping demons emerge.

Fracton codes have been intensively studied as novel topological states of matter, yet their fault-tolerant properties remain largely unexplored. Here, we investigate the optimal thresholds of self-dual fracton codes, in particular the checkerboard code, against stochastic Pauli noise. By utilizing a statistical-mechanical mapping combined with large-scale parallel tempering Monte Carlo simulations, we calculate the optimal code capacity of the checkerboard code to be $p_{th} \simeq 0.107(3)$. This value is the highest among known three-dimensional codes and nearly saturates the theoretical limit for topological codes. Our results further validate the generalized entropy relation for two mutually dual models, $H(p_{th}) + H(\tilde{p}_{th}) \approx 1$, and extend its applicability beyond standard topological codes. This verification indicates the Haah's code also possesses a code capacity near the theoretical limit $p_{th} \approx 0.11$. These findings highlight fracton codes as highly resilient quantum memory and demonstrate the utility of duality techniques in analyzing intricate quantum error-correcting codes.

We investigate the effect of Ce$^{3+}$ substitution on the magnetic ordering and phonon dynamics of the GdCrO$_3$ orthorhombic perovskite. The Ce doped compound exhibits long-range canted antiferromagnetism with Neel transitions, T$_N$ at $\sim$ 173 K, accompanied by spin-flip, T$_{SF}$ at $\sim$ 10 K. Ce$^{3+}$ incorporation drives a modification of the spin-flip transition from the $\Gamma{(G_x,A_y,F_z)}$ configuration to $\Gamma$ inducing a reorientation of the spin axis between the (001) and (001) crystallographic planes. This spin reorientation is governed by Zeeman energy and produces pronounced field-induced irreversibility between FCC and FCW magnetization processes. The substituted compound Gd$_{0.9}$Ce$_{0.1}$CrO$_3$ (GCCO) exhibits a remarkably large magnetic entropy change, $\Delta$ S $\sim$ 45-40 J/kg-K for $\Delta$ H = 90-70 kOe at 3 K among the highest reported for rare-earth orthochromites. The interplay of spin-only magnetocrystalline anisotropy from Cr$^{3+}$ and spin-orbit-driven magnetic moments of Gd$^{3+}$ and Ce$^{3+}$ results in pronounced spin-phonon coupling, manifested through the A$_{1g}$(6) vibrational mode. The observed temperature-dependent spectral evolution is consistent with behaviour reported in isostructural magnetic perovskites.

Based on a cocycle structure, we identify a new derivation of the Boltzmann distribution for finite energy-level systems from the maximal entropy principle (MEP). Our approach does not rely on the method of the Lagrange multiplier, and it provides a more transparent way to understand the dependence on the energy levels of the temperature $T = 1/\beta$ for the equilibrium distribution. Finally, we make two curious observations associated with our derivations.

We continue the variational study of the discrete-to-continuum evolution of lattice systems of Blume-Emery-Griffith type which model two immiscible phases in the presence of a surfactant. In our previous work \cite{CFS}, we analyzed the case of a completely wetted crystal and described how the interplay between surfactant evaporation and mass conservation leads to a transition between crystalline mean curvature flow and pinned evolutions. In the present paper, we extend the analysis to the regime of partial wetting, where the surfactant occupies only a portion of the interface. Within the minimizing-movements scheme, we rigorously derive the continuum evolution and show how partial wetting introduces a complex coupling between interfacial motion and redistribution of surfactant. The resulting evolution exhibits new features absent in the fully wetted case, including the coexistence of moving and pinned facets or the emergence and long-lived metastable states. This provides, to our knowledge, the first discrete-to-continuum variational description of partially wetted crystalline interfaces, bridging the gap between microscopic lattice models and experimentally observed surfactant-induced pinning phenomena in immiscible systems.

$SrVO_3$ (SVO), a model correlated metal and a promising transparent conducting oxide, develops a several-nanometer-thick near-surface region (NSR), rich in $V^{5+}$ species under ambient conditions. This oxidized layer obscures the intrinsic correlated-metallic $V^{4+}$ character and limits both fundamental studies of the physics and the material's integration into electronic devices. Here, we demonstrate a direct and controllable approach for recovering the metallic SVO surface by thermally reducing the NSR under ultra-high vacuum. Real-time in-situ X-ray photoelectron spectroscopy (XPS) reveals a sharp transformation from a $V^{5+}$-dominated surface to mixed valence states, dominated by $V^{4+}$, and a recovery of its metallic character. Ex-situ X-ray diffraction (XRD), atomic force microscopy (AFM), and high-resolution scanning electron microscopy (HR-SEM) suggest that this transformation is accompanied by mass redistribution and partial oxygen loss, leading to nanoscale surface reorganization and modest lattice expansion. While thermodynamic considerations motivate evaluation of a $V_2O_5$ volatilization pathway, the combined experimental evidence instead points toward a predominantly structural surface reorganization. These findings establish a practical method for obtaining predominantly $V^{4+}$ SVO surfaces without protective capping layers, a capability that expands the utility of SVO for advanced spectroscopies, interface engineering, and oxide-electronics device integration.

In this work, we investigate a simple nonequilibrium system with many interconnected, open subsystems, each exchanging a globally conserved resource with an external reserve. The system is represented by a random graph, where nodes represent the subsystems connected through edges. At each time step, a randomly selected node gains a token (i.e, a resource) from the reserve with probability (1-p) or loses a token to the reserve with probability p. When a node loses a token, its neighbors also lose a token each. This asymmetric token exchange breaks the detailed balance. We investigate the steady state behavior of our model for different types of random graphs: graphs without edges, regular graphs, Erdős-Rényi, and Barabási-Albert graphs. In all cases, the system exhibits a sharp, switch-like transition between a token-saturated state and an empty state. When the control parameter p is below a critical threshold, almost all tokens accumulate on the graph. Furthermore, in a non-regular graph, most tokens accumulate or condense on nodes of minimum degree. A slight increase in p beyond the threshold drains almost all the tokens from the graph. This switching transition results from the interplay between drift and the conservation of tokens. However, the position of the critical threshold and the behavior at the transition zone depend on graph topology.

The Parallel Minority Game (PMG) is a synchronous adaptive multi-agent model that exhibits active-absorbing transitions characteristic of non-equilibrium statistical systems. We perform a comprehensive numerical study of the PMG under two families of microscopic decision rules: (i) agents update their choices based on instantaneous population in their alternative choices, and (ii) threshold-based activation that activates agents movement only after overcrowding density crossing a threshold. We measure time-dependent and steady state limits of activity $A(t)$, overcrowding fraction $F(t)$ as functions of the control parameter $g=N/D$, where $N$ is the number of agents and $D$ is the total number of sites. Instantaneous rules display mean-field directed-percolation (MF-DP) scaling with $\beta\approx1.00$, $\delta\approx0.5$, and $\nu_{\parallel}\approx2.0$. Threshold rules, however, produce a distinct non-mean-field universality class with $\beta\approx0.75$ and a systematic failure of dynamical scaling. We show that thresholding acts as a relevant perturbation to the critical behavior of the model. The results highlight how minimal cognitive features at the agent level fundamentally alter large-scale critical behavior in socio-economic and active systems.

Chromatin is repeatedly deformed in vivo during transcription, nuclear remodeling, and confined migration - yet how mechanical response varies from locus to locus, and how it relates to epigenetic state, remains unclear. We develop a theory to infer locus-specific viscoelasticity from three-dimensional genome organization. Using chromatin structures derived from contact maps, we calculate frequency-dependent storage and loss moduli for individual loci and establish that the mechanical properties are determined both by chromatin epigenetic marks and organization. On large length scales, chromatin exhibits Rouse-like viscoelastic scaling, but this coarse behavior masks extensive heterogeneity at the single-locus level. Loci segregate into two mechanical subpopulations with distinct longest relaxation times: one characterized by single-timescale and another by multi-timescale relaxation. The multi-timescale loci are strongly enriched in active marks, and the longest relaxation time for individual loci correlates inversely with effective local stiffness. Pull-release simulations further predict a time-dependent susceptibility: H3K27ac-rich loci deform more under sustained forcing yet can resist brief, large impulses. At finer genomic scales, promoters, enhancers, and gene bodies emerge as "viscoelastic islands" aligned with their focal interactions. Together, these results suggest that chromatin viscoelasticity is an organized, epigenetically coupled property of the 3D genome, providing a mechanistic layer that may influence enhancer-promoter communication, condensate-mediated organization, and response to cellular mechanical stress. The prediction that locus-specific mechanics in chromatin are controlled by 3D structures as well as the epigenetic states is amenable to experimental test.

Deformable and flexible electronics have garnered significant attention due to their distinctive properties; however, their current applications are primarily limited to the thermoelectric domain. Expanding the range of these electronics and their application scope represents a pivotal trend in their development. In this work, a plastic inorganic semiconductor material, Sn2BiS2I3, with a band gap of 1.2 eV was synthesized and fabricated into a three-electrode flexible and portable electronic tongue capable of detecting heavy metal elements. The electronic tongue device exhibits exceptional linearity and demonstrates resistance against interference from impurity ions. The linear regression equation is expressed as Y=0.24+19.06X, yielding a linear coefficient of approximately 0.96, and the detectable limit stands at around 1.1 ppb, surpassing the 2.0 ppb limit of the ICP-AES instrument. Furthermore, mechanical testing reveals the favorable plasticity of Sn2BiS2I3, as evidenced by the absence of cracks during nanoindentation. The indentation hardness of Sn2BiS2I3 is approximately 1.67 GPa, slightly exceeding the hardness of Cu, which is 1.25 GPa. This study expands the repertoire of deformable and flexible electronics, offering a new and exceptional choice for biomimetic tongue sensor materials.

We study energy dissipation and propagation of information encoded by Josephson vortices in two types of circular shift register: a) a uniform register composed of sections of discrete Josephson transmission lines (JTL) forming a closed loop with a flux pump allowing to change the number of moving fluxon; b) a nonuniform register composed of sections of the regular JTL and sections of JTLs utilizing nSQUIDs - dc-SQUIDs with negative inductance between their arms - instead of single Josephson junctions. nSQUIDs are parametric devices with a flexible double-well potential that were proposed as components for reversible computing. For the uniform register, we demonstrate the energy dissipation per bit-shift operation below the Landauer's thermodynamic limit $E_T=k_BTln2$ up to propagation delays of ~0.7 ns, corresponding to the circular information motion with frequencies up to ~1.4 GHz. This does not contradict Landauer's minimum energy requirement for computations since information is not destroyed. For the nonuniform register, we find the minimum energy dissipation per bit-shift of about 16$E_T$ and attribute this to a nonuniform movement of vortices and energy barriers between the regular JTL and nSQUID sections. Differences of Josephson vortex propagation in both types of circular registers are discussed based on the measured current-voltage characteristics, extracted effective resistance and the terminal speed of Josephson vortices, and their dependences on the number of moving vorticies. nSQUID inductance connecting JJs to the ground leads to an unusual type of lossless discrete transmission line with frequency-dependent impedance and propagation speed, both different from the regular JTLs.

The recent discovery of high-temperature superconductivity in pressurized nickelates has renewed interest in the broken-symmetry states of their ambient-pressure parent phases, where a density-wave (DW) order emerges and competes with superconductivity, but its microscopic origin remains unresolved. Using ultrafast optical spectroscopy, we track quasiparticle relaxation dynamics across the DW transition at $T_{\rm DW} \approx$ 136 K in trilayer nickelate {\LNO} single crystals, revealing the opening of an energy gap of $\sim$52 meV. Multiple coherent phonons, including $A_g$ modes near 3.88, 5.28, and 2.09 THz, display pronounced mode-selective anomalies across the transition, indicating that the DW is strongly coupled to lattice degrees of freedom and suggesting an important role of electron-phonon coupling. At higher excitation densities, the DW is nonthermally suppressed, producing a temperature-fluence phase diagram that parallels pressure-tuned behavior. These results establish the DW in {\LNO} as a lattice-entangled instability involving multiple phonon modes, and highlight ultrafast optical excitation as a nonequilibrium tuning parameter for suppressing density-wave order in nickelates.

The impact of solid intruders into granular media is commonly described by a combination of quasi-static resistance and an inertial drag force proportional to the square of the impact speed. While intruder geometry is known to influence force magnitudes, its role in controlling the onset of inertial drag has remained largely unexplored. Here we present systematic impact experiments using conical intruders spanning a wide range of apex angles. By measuring the peak acceleration during impact, we show that the emergence of a well-defined inertial response depends sensitively on cone geometry. Blunt cones exhibit quadratic scaling with impact speed over the full range of velocities studied, whereas sharper cones display a delayed transition to inertial behavior at higher speeds. We define a geometry-dependent crossover speed marking the onset of the inertial regime and find that it scales approximately linearly with the cone angle through $\tan\phi$. Once the inertial regime is established, the peak force collapses when rescaled by $\tan\phi$, indicating that cone geometry controls the effective momentum transfer to the grains. These results demonstrate that intruder geometry governs not only the magnitude of inertial drag, but also the impact speed at which it becomes dominant.

Emerging spin-orbit-torque devices based on spin valves require a thin magnetic free layer to maximize the torque per moment. However, reducing the free-layer thickness to $\lesssim 2$ nm deteriorates the giant magnetoresistance (GMR) signal for electrical readout. Here, we demonstrate that the addition of a 1-nm Cu seed layer promotes sharp interfaces in simple polycrystalline Co-based spin valves, enabling high GMR ratios of 5-7% at sub-2-nm free-layer thicknesses. Our work offers a pathway for engineering high-signal GMR readout in spin-orbit-torque digital memories and neuromorphic computers.

We analyze the spin susceptibility of spin-triplet superconductors from the zero-field to finite-field regimes, with emphasis on its implications for Knight-shift measurements. In the zero-field limit, we review the general expression for the static spin susceptibility and highlight the universal zero-temperature sum rule, $\sum_i \chi_{ii}(T=0)=2\chi^N$, which constrains the residual susceptibility components for any triplet state. Using representative isotropic, helical, and chiral $\vec{d}$-vectors, we illustrate how the Knight shift encodes the spin configuration of the order parameter and show that the sum rule remains robust even for anisotropic Fermi surfaces. We then incorporate magnetic field effects through a semiclassical Doppler shift of quasiparticle energies in the vortex state. The resulting field dependence of the susceptibility - including both longitudinal (Knight-shift) and transverse magnetic susceptibility components - provides a sensitive probe of nodal directions and the momentum dependence of the $\vec{d}$-vector. Applying this framework to UTe$_2$, we demonstrate how the distinct irreducible representations allowed by orthorhombic symmetry can be differentiated by their field-dependent susceptibility.