Papers

The asymmetric three-qubit W-class state $|\overline{W_3^L}\rangle$ defines an isosceles entanglement-network geometry, (a) two vertex-base (VB) links form stronger bipartite connections, (b) while the base-base (BB) link is weaker. This suggests that concentrating entanglement into a super-link may be advantageous for quantum-network tasks. Here, we show that this intuition is incomplete. We analytically compare the bipartite concurrence dynamics of the symmetric |W> state and the asymmetric $|\overline{W_3^L}\rangle$ state, which differ both in entanglement-network geometry and excitation sector under standard noise models. In the absence of noise, the concurrence hierarchy is $C_{VB} > C_W > C_{BB}$. Under phase damping, this hierarchy is preserved for all noise strengths and no entanglement sudden death occurs. Under amplitude damping, however, the hierarchy is reordered. The symmetric |W> state becomes the most robust, while the base-base concurrence of $|\overline{W_3^L}\rangle$ vanishes at the finite threshold of parameter $\gamma$. We term this reordering as the Super-Link Fragility Effect. The same structural asymmetry that produces a stronger vertex-base link also makes it more vulnerable to energy dissipation when coupled with multi-excitation amplitudes. Under depolarization, the asymmetry advantage is erased, with $C_W$ and $C_{VB}$ sharing the same sudden-death threshold for some value of the parameter p, while $C_{BB}$ disappears earlier at some other value of the parameter p. The generalized amplitude damping channel continuously connects the damping-dominated regime to the pure-excitation limit, where the initial hierarchy is restored. These results show that entanglement robustness in $W$-class resources is controlled not by initial concurrence alone, but by the joint structure of entanglement-network geometry, excitation sector, and noise symmetry.

The realization of on-chip photonic circuits requires scalable and deterministic single-photon emitters (SPEs) at room temperature, which remain a challenge in van der Waals materials. In this work, we report a novel three-step fabrication process for the generation of spatially controlled SPE arrays in hexagonal boron nitride (hBN). The process comprises site-selective gallium (Ga) focused ion beam milling, nanoscale conformal carbon deposition over the patterned regions, and subsequent thermal annealing. The synergistic combination of these steps resulted in a site-correlated emitter yield of ($\sim 89\%$) across 100 fabrication sites. Second-order autocorrelation measurements revealed pronounced three-level emitter dynamics where the best emitters exhibited high purity ($g^{(2)}(0)=0.15 \pm 0.09$).To the best of our knowledge, this is the first lithography-free, direct-write approach combining Ga-ion milling, selective carbon engineering, and thermal annealing to deterministically generate \hBN{} \SPE{}s. The reproducibility of the method is validated across multiple independently fabricated samples. These results establish a scalable, lithography-free pathway toward on-demand SPE arrays relevant to integrated quantum photonics.

This study investigated how additively manufactured nozzles with body-centered cubic lattice structures reduce the flame flashback propensity in a hydrogen jet flame burner. Five different configurations of a jet flame combustor were investigated, with a focus on mixing duct walls incorporating porous media. The nozzles were manufactured by the powder bed fusion of metals using a laser beam process. The lattice parameters were varied by the volume fraction and the strut diameter. For the experiments, pure hydrogen was used as fuel under atmospheric conditions at various equivalence ratios and Reynolds numbers of 9,000 - 12,000. Flow field measurements, flame imaging, and spectral proper orthogonal decomposition of the flame dynamics were employed to identify possible transition mechanisms from a stable operation to flashback. The flow fields and the flame shapes showed only minor effects from wall modifications, preserving general flow characteristics across configurations. The flow dynamics in the combustion chamber were dominated by large-scale coherent structures in the shear layer, specifically Kelvin-Helmholtz instabilities. The results demonstrated that the nozzle with the coarsest porous wall structure significantly improved the flashback resistance compared to a nozzle with a solid wall. It is concluded that the primary mitigation mechanism was a cooling effect by unburnt mixture flowing through the porous media. The findings confirmed that the integration of lattice structures through additive manufacturing provides a viable strategy for hydrogen flashback mitigation by manipulating the coupled interaction between the flame and the thermal conditions of the wall.

The Ising framework maps the decoding problem in quantum error correction onto ground-state optimization of a classical Hamiltonian, in which $X$-$Z$ error correlations enter as cross terms. Under phenomenological depolarizing noise, the exact joint formulation contains up to 8-body interactions for the toric code and 10-body for the $6.6.6$ color code. These high-order terms degrade solver convergence, inflate runtime, and raise the auxiliary spin overhead when embedding into native 2-body Ising hardware. In this work, we propose the iterative low-order decoding (ILOD) algorithm, which alternates between $X$- and $Z$-type sub-Hamiltonians, approximating cross-type correlations through Bayesian priors that reweight each type's couplings using the other type's inferred error configuration. This halves the maximum body count of interaction terms in the Hamiltonian, accelerating the solver, restoring convergence at larger code distances, and reducing the total spin count for 2-body embedding by a factor of $2.5$. For the toric code, ILOD attains a threshold of $4.73%$ versus $4.83%$ for the joint formulation, with the empirical runtime ratio scaling as $(0.81)^d$. For the $6.6.6$ color code, their thresholds agree within statistical uncertainty for small code distances, and ILOD remains convergent for larger distances where the joint formulation fails to converge despite a larger annealing budget.

We develop a general framework for matching KSVZ-like ultraviolet completions featuring vector-like fermions and a spontaneously broken $\mathrm{U}(1)_{\mathrm{PQ}}$ symmetry onto the Standard Model Effective Field Theory and the low-energy axion-like particle effective theory. The framework applies to arbitrary vector-like fermion representations and PQ-charge assignments, and systematically captures the effective interactions generated in both sectors. We then perform a comprehensive phenomenological analysis of the resulting SMEFT operators, based on global fits to electroweak precision, Higgs, and flavor observables, obtaining robust bounds on the corresponding Wilson coefficients that are largely independent of the details of the ultraviolet realization. These constraints can subsequently be translated into the QCD axion and ALP parameter space, providing indirect probes of ALP couplings. We further investigate several representative examples of the interplay between the SMEFT and ALP sectors, illustrating how direct ALP searches and indirect precision and flavor observables provide complementary information on the same underlying dynamics. We find that, over large regions of parameter space, indirect constraints derived from the SMEFT analysis dominate over direct ALP probes, except in scenarios where the PQ-charge assignment permits mass mixing with Standard Model fermions. Overall, our results establish a unified framework for connecting ultraviolet completions, SMEFT analyses, and ALP searches, enabling both the interpretation of existing constraints and the exploration of future signals within a common theoretical setting.

The origin of cosmic acceleration remains a central problem in cosmology, commonly attributed to a cosmological constant within the $\Lambda$CDM model or to dynamical dark energy. Here, we develop an alternative approach in which acceleration emerges from quantum post-selection, a standard feature of quantum theory that is not usually incorporated into cosmological modelling. While quantum theory admits both pre-selected and post-selected ensembles, quantum cosmological models are almost exclusively formulated in terms of initial conditions. Building on previous work on post-selected quasiclassical dynamics, we construct a minimal predictive cosmological model in which post-selection and coarse-graining generate effective late-time acceleration without introducing a cosmological constant, dark energy, or modifications of general relativity. The resulting expansion history is highly constrained theoretically and depends on at most two parameters beyond standard Friedmann evolution. Confrontation with type Ia supernova and cosmic chronometer data yields statistically competitive fits while naturally avoiding the coincidence problem. The model also reproduces the standard radiation- and matter-dominated behaviour at early times and predicts a present-day jerk parameter significantly different from the $\Lambda$CDM value. These results suggest that cosmic acceleration may arise as a macroscopic quantum cosmological effect rather than from additional cosmological fluids or modified gravitational dynamics.

We investigate a problem, that describes coupling between diffusion and reaction inside a thin circular cylindrical tube. The asymptotic solution of the posed problem is derived by means of the boundary functions method. We perform comparison of this asymptotic solution against corresponding exact solution, which revealed serious methodological drawbacks of known Fick-Jacobs reduction approach. The results obtained may be used to study a wide range of reaction-diffusion problems, when the Fick-Jacobs method cannot be applied.

Reference frame mismatch among communication parties introduces errors in quantum cryptographic protocols. As the number of participants increases, aligning reference frames becomes increasingly difficult, complicating multipartite quantum cryptographic implementations. Here, we theoretically and experimentally investigate multipartite reference-frame-independent (RFI) quantum cryptographic communication using Greenberger-Horne-Zeilinger (GHZ) states. We generalize the bipartite RFI security parameter $C$ to an $N$-party parameter $C_N$ and derive the asymptotic secret key rate expressed solely in terms of experimentally accessible quantities. We analyze the key rate under global and local depolarizing noise models and find that increasing the number of parties $N$ enhances robustness against global depolarizing noise while increasing vulnerability to local channel noise. We also present a proof-of-principle experimental demonstration of four-party RFI quantum cryptographic communication using four-photon GHZ states, confirming the reference-frame invariance of both the $C_4$ parameter and the secret key rate under various reference frame rotations.

The quark-gluon plasma (QGP) is a deconfined state of strongly interacting matter formed at extreme temperature and energy density in ultra-relativistic nucleus-nucleus collisions at RHIC and the LHC. High transverse momentum jets, produced in initial hard scatterings, traverse the QGP and lose energy via elastic and radiative processes, an effect known as jet quenching. The nuclear modification factor, $R_{\mathrm{AA}}$, defined as the ratio of the Pb+Pb jet yield to the $pp$ cross section scaled by the nuclear thickness function, is widely used to quantify jet quenching. However, its value depends strongly on both the $pp$ jet spectral shape and the strength of the quenching, complicating comparisons across jet selections. The fractional energy loss, $S_{\text{loss}}$, quantifying the average medium-induced momentum shift of jets, is designed to mitigate this dependence. In central Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02~\mathrm{TeV}$, we compile and compare published ATLAS and ALICE measurements of jet suppression for inclusive single-jet and dijet selections across multiple jet radii, considering (i) the ratio of the nuclear modification factor at a given radius to that at a reference radius of 0.2, and (ii) the fractional energy loss. The radius dependence of this ratio differs between single-jet and dijet measurements, and between ATLAS calorimeter jets and ALICE charged-particle jets, reflecting differences in kinematic event selections and jet constituents. Expressing the results in terms of $S_{\text{loss}}$ allows direct, radius-differential comparisons across experiments with reduced sensitivity to the $pp$ spectral slope. Combining these approaches enables constraints on the radius dependence of jet modification that account for selection biases, and facilitates cross-experiment benchmarking of jet quenching models.

Recently an approximate additivity principle has been successfully applied to molecular chemical potentials in ultrasoft particle systems [R. L. Hendrikse, C. Amador and M. R. Wilson, Phys. Chem. Chem. Phys. 27, 1554 (2025)]. We articulate this additivity principle as one of two ansätze and apply these to dimers in the infinite dilution limit. In this limit, we show that the solvent-mediated potential of mean force (PMF) between solutes, extracted from the hypernetted-chain (HNC) closure of the Ornstein-Zernike equations, can be expressed as a convolution between solute-specific generalised excluded volume functions. In the limit of a structureless solvent of point particles and hard core solutes, this recovers the exact Asakura-Oosawa depletion potential as the overlap between excluded volume spheres. The theory extends to the opposite limit of ultrasoft particle solvents and solutes, such as those encountered in dissipative particle dynamics (DPD), where the solvent-mediated PMF can be recovered with considerable accuracy. These results confirm that in coarse-grained molecular DPD simulations the parametrisation of the solute-solvent interactions is sensitive to the intramolecular bond lengths if they are smaller than the range of the solvent-solvent interaction potential, due to the overlap of the soft excluded volume functions.

Multipole-Padé approximants provide a compact representation of dynamical response functions in terms of a small number of collective modes. Here, we generalize this framework to incorporate momentum dependence across the full Brillouin zone of 2D metals by constructing a symmetry-conserving, anisotropic representation of the inverse dielectric function. This analytic form enables efficient and accurate evaluation of quantities involving dynamical screening, including spectral features and correlation energies. We construct such compact representations for a set of seven two dimensional metals spanning distinct electronic regimes, and show that a small number of dispersive plasmonic modes suffices to accurately describe the dielectric response across the full Brillouin zone, while also yielding accurate correlation energies. The proposed representation therefore establishes a direct bridge between ab initio calculations and analytical models of screening, opening new avenues for applications in condensed matter systems.

Rydberg excited states are challenging to describe due to their highly diffuse character. Orbital-optimized density functional calculations typically provide more accurate values of the excitation energy of Rydberg states than time-dependent density functional theory approaches, but their performance with respect to excited-state properties such as the dipole moment remains largely unexplored as existing benchmarks are limited to the lowest excited state. Here, orbital-optimized density functional calculations with a plane-wave basis set are used to compute the dipole moment of several Rydberg states of a set of small molecules. Plane waves provide a flexible representation of diffuse Rydberg orbitals, whereas commonly used atomic orbitals basis sets are found to yield significant discrepancies in the magnitude and orientation of the dipole moment, even when diffuse functions are included and the excitation energy agrees well with the plane-wave results. The generalized gradient approximation functional PBE used in combination with the plane-wave representation of the orbitals gives good agreement with higher-level coupled-cluster calculations performed with sufficiently diffuse basis sets, when the latter are available. The hybrid functional PBE0 further improves the results, while PBE with globally scaled explicit Perdew-Zunger self-interaction correction generally leads to larger errors and an overestimation of the dipole moment, despite restoring the correct asymptotic $-1/r$ dependence of the effective Kohn–Sham potential.

Renormalization group (RG) fixed points are commonly associated with scale invariance and a divergent correlation length. We show that this connection can fail in the tricritical $O(N)$ model in three dimensions in the limit $N\to\infty$. Revisiting the line of fixed points identified by Bardeen, Moshe, and Bander, we use the functional renormalization group to clarify the mechanism leading to mass generation at its singular endpoint (the BMB fixed point). We demonstrate that the generated mass is nonuniversal and originates from the nonanalytic structure of the effective potential. We show that the critical exponent $\nu$ which takes the value $\nu = 1/2$ along the regular part of the BMB line, that is, for $0 \leq \lambda < \lambda_{\rm BMB}$, jumps to $\nu = 1/3$ on the singular part of this line with the BMB FP, corresponding to $\lambda = \lambda_{\rm BMB}$, being the pivotal point between these two regimes. We also show how its singular potential emerges dynamically along the renormalization flow.

Quasars can serve as standard candles through an empirical scaling relation between their ultraviolet (UV) and X-ray luminosities. As high-redshift probes, it is critical to test whether this relation evolves with redshift. In this work, we reconstruct the Hubble diagram of the Pantheon+ sample using the deep learning–based LADDER algorithm and use it as a reference to investigate the quasar scaling relation. Our results, which are consistent with those from Gaussian process regression and narrow-bin analyses, show that the potentially contaminated sample at $z<0.7$ differs significantly from the $z>0.7$ sample; thus, it should be further screened or excluded when quasars are used as cosmological probes. We find that the scaling relation exhibits a non-linear redshift dependence that cannot be accounted for by a simple linear correction, and that this behavior is a feature of the current data sample rather than a consequence of cosmological model misspecification. To use quasars as standardizable candles, further modeling of the scaling relation and intrinsic dispersion, or more advanced data processing techniques, is required.

Contemporary theoretical descriptions of nuclear structure rely mainly on microscopic, single-particle frameworks often in competition with collective degrees of freedom, especially when deformation plays a dominant role. Such phenomena are prominent in the rare-earth region, where rotational band structures and enhanced electric quadrupole transitions are systematically examined. The Confined beta-Soft (CBS) rotor model, introduced by N. Pietralla and O.M. Gorbachenko, bridges the gap between the X(5) critical point and the rigid-rotor limit in the region where the R_4/2 = E(4+)/E(2+) ratio lies between 2.904 and 3.333. In the present work, the CBS framework is employed to calculate ground-state band energies, associated B(E2) transition rates, and beta-band excitations of even-even nuclei in the rare-earth region. The theoretical results are systematically compared with available experimental data, and predictions are provided for nuclear observables that have not yet been measured, offering guidance for future experimental investigations.

We have developed a fabrication technique for Josephson junctions that employs a three-dimensional patterned, low-loss substrate instead of commonly used organic resists. The technique enables the fabrication of high-quality trilayer junctions from a wide range of geometries and materials, including high-melting-point superconductors such as tantalum or niobium. The junction electrodes are free from intentionally introduced oxides and organic materials, which are known sources of decoherence. We fabricate and characterize underdamped Nb/AlOx/Nb junctions of different sizes in several geometries. Such junctions enable manufacturing of quantum circuits operating at higher speeds and elevated temperatures.

Given a finite collection of Lagrangian submanifolds $\mathscr L$ in a compact symplectic manifold $X$, we construct a cyclic, filtered, strictly unital curved $A_{\infty}$ category $\mathcal L$ and develop Floer theory of closed-open maps and open-closed maps. Using them, we prove that, whenever the map from the quantum cohomology of $X$ to the Hochschild cohomology of the Fukaya category $\mathcal L$ with objects $\mathscr L$ is injective, the following consequences follow: (1) any other Lagrangian submanifold equipped with a weak bounding cochain lies in the category split-generated by $\mathscr L$, and (2) the Hochschild homology and cohomology of the Fukaya category are isomorphic to quantum cohomology. In the exact case a similar result was obtained in [Ab]. We also provide some applications.

Diffuse-interface models are a widely used framework for interfacial dynamics in complex fluids, in which interfaces are represented through smooth transition layers and capillary effects are encoded by a free-energy functional. For incompressible mixtures with more than two phases, however, robust computation is substantially more difficult because the numerical method should preserve the balance structure of the continuum model, maintain the saturation constraint, dissipate energy, and treat all phases symmetrically even when density ratios are arbitrary. Existing structure-preserving methods are largely developed for binary flows or for formulations that distinguish a reference phase, so a genuinely symmetric N-phase discretization remains lacking. The practical problem is therefore to construct a fully-discrete method for N-phase incompressible Navier–Stokes–Cahn–Hilliard mixture models that retains the key thermodynamic and conservation properties of the continuum equations for arbitrary density ratios. Here we propose a symmetric fully-discrete method for the N-phase incompressible Navier–Stokes–Cahn–Hilliard mixture model with arbitrary density ratios. The method yields a fully-discrete problem in which every solution satisfies exact phase volume conservation, phase mass conservation, total volume conservation, total mass conservation, and a discrete energy-dissipation law. In addition, if the volume-saturation constraint holds for the initial data, then it is preserved at every time step. We numerically verify these structure-preserving properties and demonstrate the robustness of the method in representative multiphase flow problems. The resulting scheme provides a computational framework for incompressible N-phase mixture flows with complex interfacial dynamics and arbitrary density contrasts.

Nonlinear structure growth generates higher-order correlations and morphological features in the cosmic density field that cannot be fully characterized by two-point statistics. Upcoming surveys will measure these features with greater precision, making it essential to develop methods capable of extracting as much cosmological information as possible from them. Field-level inference (FLI) is one such approach, in which cosmological parameters are constrained by comparing observed maps to forward-modeled maps, either directly or through learned summaries that retain map-level information. In this work, we compare FLI with power-spectrum-based inference using the same forward-modeling pipeline for generating weak lensing maps, with the goal of quantifying the gain from map-level analysis relative to two-point statistics. We perform this comparison with both implicit and explicit inference methods, using 8-million-parameter forward models based on Lagrangian perturbation theory and particle-mesh (PM) N-body evolution. The two FLI approaches yield closely consistent posteriors; this agreement, together with coverage tests confirming the calibration of the implicit analyses, gives us confidence in the recovered field-level constraints. Relative to the power-spectrum-based analyses, these results show significant gains in cosmological information, especially when small scales are included in the PM-based forward model. We then discuss the remaining challenges that must be addressed before PM-based explicit FLI can be applied to observational datasets.

Studying young stellar populations within galaxies can help refine our understanding of recent star formation in galaxies and their evolution. With this motivation, we present a catalog of ~25,000 recently formed (within 400 Myr) star-forming clumps (SFCs) in 17 morphologically diverse nearby galaxies, including 8 massive, classic spirals, 6 intermediate-mass, flocculent spirals, and 3 dwarf irregulars. We used far- and near-UV observations from the UltraViolet Imaging Telescope (UVIT), whose ~1.5" angular resolution and 28' field-of-view allow us to probe SFCs at a mean physical scale of ~54 parsec, within the full extent of our galaxies. We adopted a homogeneous SFC detection criterion, corrected for spatially varying dust attenuation (using 6" resolution A_V maps, made by combining FUV with archival infrared observations), and estimated the SFC ages by comparing the observed UV color-magnitude diagrams with Starburst99 simple stellar population models. Using our SFC catalog, we studied the age demographic of the recently formed stellar populations across different galaxy morphologies and observed age trends consistent with several well-known phenomena, such as the inside-out formation of disc galaxies, local gravitational instabilities leading to flocculent spiral arms, and the stochastic nature of star formation in dwarf galaxies. Leveraging full galaxy coverage and FUV data, our catalog complements existing optically-identified star cluster catalogs in the literature towards improving our understanding of star formation across a wide range of galaxy morphologies, masses, and environments. We make the SFC catalog and A_V maps of our 17 galaxies publicly available with this paper.

In this review, we discuss the most recent developments and applications of the Subtracted Second RPA (SSRPA), an extension of the Second RPA (SRPA), which overcomes its pathological issues encountered within the Energy Density Functional theory. After recalling the formal properties of the SRPA and SSRPA, the anomalous behavior of SRPA is shown and discussed by presenting several applications with different kinds of nuclear interactions. The most recent pathology-free SSRPA studies are then presented both for charge-conserving and charge-exchange nuclear excitations. The comparison with experimental data is presented to assess and quantify the improvement introduced by the SSRPA with respect to the RPA and SRPA. The impact of beyond-mean-field correlations induced in SSRPA is also qualitatively estimated in connection with the modeling of the nuclear equation of state. We conclude by discussing the future perspectives of the SSRPA, focusing on its potential connections with some current experimental challenges and outlining necessary theoretical extensions and numerical developments.

An accurate determination of the dust attenuation within galaxies is essential to derive key physical properties such as the star formation rate (SFR). We present an analysis using the JWST/NIRISS data from the GLASS-JWST ERS programme to investigate and characterise the stellar and nebular reddening of galaxies at $1.0<z<2.4$, down to the sub-kpc scale. We use a multiregion fitting method to extract high-quality H$\alpha$ and H$\beta$ emission line maps for 99 individual galaxies across a stellar mass range $7.0<\log_{10}(M_*/\mathrm{M}_{\odot})<10.5$. We find no evidence for ratios of the Balmer decrement (H$\alpha$/H$\beta$) below the intrinsic limit for Case B recombination, beyond the expected variation from observational uncertainties. We reproduce the local correlation between the Balmer decrement and total stellar mass, and find no measurable difference when splitting the sample by redshift, with negligible attenuation below $\log_{10}(M_*/\mathrm{M}_{\odot})\lesssim8.5$. Similarly, the best-fit relation between the nebular and continuum reddening follows the same relation as in local starburst galaxies, $E(B-V)_{\mathrm{SED}} = (0.46\pm0.02)E(B-V)_{\mathrm{neb}}$, together indicating no significant evolution in the dust geometry within galaxies out to $z\lesssim2.4$. We derive best-fit linear relations between the differential nebular-stellar reddening and the SED-derived star formation rate (SFR) and stellar mass, finding statistically significant relations for both quantities. We use our spatially-resolved measurements to derive an empirical calibration between the resolved differential reddening, and the SFR surface density. These will enable crucial dust attenuation corrections for spatially-resolved science at higher redshifts where the Balmer lines are inaccessible, such as with future Roman grism observations.

We describe $\texttt{stistools.barycentric_correction}$, a new Python utility for calculating barycentric timing corrections for HST/STIS observations. This tool replaces the deprecated $\texttt{stsdas.hst_calib.stis.odelaytime}$ IRAF function that was previously used for HST barycentric corrections. Our new utility uses $\texttt{astropy}$ for conversion between time formats and standards and introduces a new way to calculate HST's position through JPL Horizons, replacing the need to download separate HST orbital ephemeris files. Here, we describe the methods used in the new utility, the tests that were carried out to verify its accuracy, and explain some of the complexities involved in determining light travel times to accuracies down to a millisecond for HST. We also summarize the current understanding of the absolute accuracy of STIS time stamps.

The location of a star along the horizontal branch (HB) during core-helium burning is primarily determined by the amount of mass lost by its progenitor. We investigate the formation and properties of stripped core-helium-burning stars, focusing on how the residual hydrogen-envelope mass ($M_{\mathrm{env}}$) and the timing of envelope removal shape their properties. We used the MESA stellar evolution code to model stars that lose their hydrogen envelopes on the first giant branch. We explored two limiting cases for the timing of stripping, corresponding to the minimum and maximum core masses for helium ignition, for progenitors with initial masses below $\sim$6 $M_{\odot}$ at two metallicities ($Z=0.02$ and $Z=0.004$), while systematically varying $M_{\mathrm{env}}$. As expected, the effective temperature along the HB decreases as $M_{\mathrm{env}}$ increases. We determined the maximum $M_{\mathrm{env}}$ required to avoid subsequent evolution through the thermally pulsing asymptotic giant branch, which ranges from $\sim0.05$ $M_{\odot}$, for low-mass progenitors to $\sim0.30$ $M_{\odot}$ for intermediate-mass progenitors. In low-mass progenitors, early envelope removal triggers a late hot flash, naturally explaining the hottest blue hook stars. In intermediate-mass systems, partial envelope stripping can produce extended pre-HB configurations consistent with puffed-up stripped stars observed in binaries with Be companions. Our post-stripping evolutionary tracks are publicly available for use in binary evolution and population synthesis studies.

The interplay among superconductivity, magnetism, and nontrivial band topology represents one of the most compelling frontiers in condensed matter physics. The exploration of novel superconductivity in 3d transition-metal compounds, particularly the rare Cr-based systems containing strongly magnetic Cr ions, has long attracted attention owing to their unconventional pairing mechanisms that challenge conventional wisdom. Yet, Cr-based superconductors remain scarce, especially those possessing nontrivial topological character, underscoring the urgent need to uncover new members. Here we report the observation of superconductivity in pressure-amorphized Cr-based topological insulator CrP$_4$. Upon compression, CrP$_4$ undergoes an anomalous quantum phase transition from a metallic to a semiconducting-like state at around 15 GPa, driven by significant changes in the electronic structure. At approximately 70 GPa, re-metallization with superconductivity occurs alongside an irreversible amorphization. The superconducting transition temperature Tc increases monotonically with pressure, reaching 4.8 K at 141.3 GPa. Furthermore, theoretical calculations predict multiple topological phase transitions from a strong topological insulator to a trivial state and finally back to a strong topological state under pressure. Our study not only establishes CrP$_4$ as the first Cr-based amorphous superconductor but also opens a new paradigm for exploring superconducting and topological properties in amorphous materials.

We investigate the electronic structure, spin textures, and charge to spin/orbital transport in graphene encapsulated by 1T-TaS$_{2}$ monolayers in the charge density wave phase. Using first-principles calculations, tight-binding modeling, and the Kubo formalism, we show that the encapsulation stacking dictates fundamentally distinct transport regimes. In the asymmetrical (AA) stacking, proximity fields from both interfaces constructively interfere, yielding a cumulative Rashba phase of nearly $\pi/2$. This pure radial Rashba spin pattern leads to the unconventional Rashba-Edelstein effect, which robustly dominates over the conventional response by a factor of 35 across a wide energy range. Conversely, the symmetrical (AA') stacking preserves a horizontal mirror symmetry, establishing a stable, purely out-of-plane persistent spin texture. Furthermore, the computed orbital Hall effect is exceptionally efficient, surpassing the spin Hall effect by three orders of magnitude. Within the proximity-induced spectral gaps, the orbital Hall conductivity exhibits a finite plateau, whereas the spin Hall conductivity vanishes. Our findings establish graphene encapsulated heterostructures as a promising system for realizing distinct charge to spin and charge to orbital interconversion regimes determined by the choice of stacking order.

We propose a phase field model able to simulate hydrogen-assisted cracking in polycrystalline materials. Within a variational framework, the model simultaneously describes crack propagation and hydrogen segregation on crack surfaces and grain boundaries together with the associated reduction in interfacial energies. In the context of hydrogen-enhanced decohesion (HEDE) mechanisms, we demonstrate the ability of the model to capture the transition from transgranular cracking to hydrogen-assisted intergranular cracking.

In this paper, one constructs a nonlocal extension of the Rarita-Schwinger theory for spin-$3/2$ fermions. Two classes of analytic form factors are considered: scalar form factors $f(\Box)$ and Dirac-operator form factors $f(\slashed{\partial})$. The massless theory is treated together with a covariant nonlocal gauge fixing, which allows the propagator to be written directly in terms of the spin-$3/2$ projector. In the massive theory, we show that the free Rarita-Schwinger constraints remain intact for analytic form factors, so that the unphysical spin-$1/2$ sector does not become dynamical. For $f(\Box)$ the tensor-spinor structure of the propagator is the same as in the local theory, while the pole equation is deformed by the scalar form factor. For $f(\slashed{\partial})$ the physical modes obey a nonlocal Dirac-type equation, leading to modified dispersion relations that can be written explicitly for exponential form factors. We discuss the conditions under which the construction is ghost-free at the free level and identify the natural limitations that must be addressed before interactions are introduced.

The Stokes parameters are three real parameters that completely characterize partially coherent optical fields spanned by two modes – whether a pair of polarization or spatial modes – and their use is thus ubiquitous in optics. Because the Stokes parameters are defined through an expansion of the $2\times2$ coherence matrix in terms of the Pauli matrices, they cannot be applied to optical fields comprising three modes, which are described by a $3\times3$ coherence matrix. Examples of such fields include the polarization of non-paraxial fields (spanned by three orthogonal polarization modes), and fields comprising three spatial or temporal modes. It has long been theorized that the $3\times3$ Gell-Mann matrices – developed in high-energy particle physics – can serve as a basis for $3\times3$ optical coherence matrices, with 8~expansion coefficients known as the Stokes-Gell-Mann (SGM) parameters, but the measurement procedure is daunting, and the SGM parameters have not been measured directly to date in optics. Here we present the first measurements of the SGM parameters for partially coherent three-mode light in a photonic integrated platform comprising a hexagonal mesh of Mach-Zehnder interferometers. Measuring the SGM parameters on chip, from which we reconstruct the $3\times3$ coherence matrix facilitates exploring the full space of iso-entropy fields that can be inter-converted into each other unitarily, and those that share the same value of entropy and yet cannot be inter-converted unitarily. These results pave the way to utilizing multimode partially coherent light in applications involving optical communications, sensing, and information processing.

Origami-inspired metamaterials exploit the interplay between geometry and elasticity to achieve programmable mechanical responses. Yet the origin and tunability of snap-through instabilities in non-rigidly foldable patterns remain poorly understood. Here we show that the Mars tessellation, a degree-4 vertex origami pattern composed of alternating square and rhombic faces, is not rigidly foldable because the folding-speed ratios required for vertex compatibility cannot be propagated consistently across neighboring units. This geometric incompatibility forces the facets to bend during folding, giving rise to a reproducible snap-through discontinuity in the force-displacement curve with a mean force drop of about 92.6 +/- 5.5 %, marking a transition between metastable states. Laser scoring of additional diagonal creases, guided by strain-field simulations, enables continuous tuning of the snap magnitude. These results reveal a general mechanism by which geometric frustration can be harnessed to program multistability in thin-sheet metamaterials.