How the detailed structure of quantum complexity emerges from quantum dynamics remains a fundamental challenge highlighted by advances in quantum simulators and information processing. The celebrated Small-Incremental-Entangling (SIE) theorem provides a universal constraint on the rate of entanglement generation, yet it leaves open the problem of fully characterizing fine entanglement structures. Here we introduce the concept of Spectral-Entangling strength, which captures the structural entangling power of an operator, and establish a spectral SIE theorem: a universal speed limit for R'enyi entanglement growth at $\alpha \ge 1/2$, revealing a robust $1/s^2$ decay threshold in the entanglement spectrum. Remarkably, our bound at $\alpha=1/2$ is both qualitatively and quantitatively optimal, defining the universal threshold beyond which entanglement growth becomes unbounded. This exposes the detailed structure of Schmidt coefficients and enables rigorous truncation-based error control, linking entanglement structure to computational complexity. Building on this, we derive a generalized entanglement area law under an adiabatic-path condition, extending a central principle of quantum many-body physics to general interactions. As a concrete application, we show that one-dimensional long-range interacting systems admit polynomial bond-dimension approximations for ground, time-evolved, and thermal states, thereby closing the long-standing quasi-polynomial gap and demonstrating that such systems can be simulated efficiently with tensor-network methods. By explicitly controlling R'enyi entanglement, we obtain a rigorous, a priori error guarantee for the time-dependent density-matrix renormalization-group algorithm. Overall, our results extend the SIE theorem to the spectral domain and establish a unified framework that unveils the detailed and universal structure underlying quantum complexity.

The discovery of wide-orbit giant exoplanets has posed a challenge to our conventional understanding of planet formation by coagulation of dust grains and planetesimals, and subsequent accretion of protoplanetary disk gas. As an alternative mechanism, the direct in-situ formation of planets or planetary cores by gravitational instability (GI) in protoplanetary disks has been proposed. However, observational evidence for GI in regions where wide-orbit planets are formed is still lacking. Theoretical studies predict that GI induces spiral arms moving at the local Keplerian speed in a disk. Here, with multiple high angular resolution observations over a seven-year time baseline using the Atacama Large Millimeter/submillimeter Array (ALMA), we report the evidence for spiral arms following the Keplerian rotation in the dust continuum disk around the young star IM Lup. This demonstrates that GI can operate in wide-orbit planet-formation regions, establishing it as a plausible formation mechanism for such planets.

Algol is a well-known eclipsing binary hosting an active and variable star that exhibits frequent stellar flares. Here, we report our pre-planned and coordinated rapid X-ray follow-up observations of an eclipsing flare on Algol. The Monitor of All-sky X-ray Image (MAXI) detected a flare on Algol at 05:52 UT on 2018 July 4. Subsequently, we carried out a prompt X-ray monitoring with the Neutron star Interior Composition Explorer (NICER) starting at 19:45 UT on the same day, and the observation ended at 06:02 UT on 2018 July 6. During the decaying phase of the flare, we successfully detected a 5.8-hour-long eclipse, corresponding to the secondary eclipse in which Algol A blocks the line of sight to Algol B. During the eclipse, the 2--10 keV X-ray flux is decreased to 20\% level from $1.9\times10^{-10}~ \mathrm{erg~cm^{-2}~s^{-1} }$ to $4.5\times10^{-11}~ \mathrm{erg~cm^{-2}~s^{-1} }$. We found a configuration of the flare size and location to explain the X-ray observations; e.g., the flare occurred at the latitude 45°S of the Algol B surface with a flare height of $1.9\times10^{11}~\mathrm{cm}$, corresponding to 0.8 times the stellar radius of Algol B, giving 80% obscuration of the flare loop by Algol A. The apparent absorption increase before the eclipse might originate from coronal mass ejection (CME) in the line of sight ejected during the flare.

Water is essential to our understanding of the planet-formation process and habitability on Earth. Although trace amounts of water are seen across all phases of star and planet formation, the bulk of the water reservoir often goes undetected, hiding crucial parts of its journey from giant molecular clouds to planets. This raises the question of whether water molecules in comets and (exo-)planets is largely inherited from the interstellar medium or if the water molecules are destroyed and then reformed in the disk. Water isotopologue ratios involving doubly deuterated water (D$_2$O) are a sensitive tracer to answer this question. We present strong evidence of inheritance through an enhancement of D$_2$O in the outbursting V883 Ori disk. The high D$_2$O/H$_2$O ratio of $(3.2 \pm 1.2) \times 10^{-5}$ is consistent with values seen in protostellar envelopes and a comet and is two orders of magnitude higher than expected if water is reprocessed. The high deuteration of the heaviest isotopologues D$_2$O/HDO = $(2.3 \pm 1.0) \times$HDO/H$_2$O further establishes the inheritance of water. We conclude that water ice in disks originates from the earliest phases of star formation, providing the missing link between cold dark clouds and (exo-)comets.

We investigate whether a novel method of quantum machine learning (QML) can identify anomalous events in X-ray light curves as transient events and apply it to detect such events from the XMM-Newton 4XMM-DR14 catalog. The architecture we adopt is a quantum version of the long-short term memory (LSTM) where some fully connected layers are replaced with quantum circuits. The LSTM, making predictions based on preceding data, allows identification of anomalies by comparing predicted and actual time-series data. The necessary training data are generated by simulating active galactic nucleus-like light curves as the species would be a significant population in the XMM-Newton catalog. Additional anomaly data used to assess trained quantum LSTM (QLSTM) models are produced by adding flares like quasi-periodic eruptions to the training data. Comparing various aspects of the performances of the quantum and classical LSTM models, we find that QLSTM models incorporating quantum superposition and entanglement slightly outperform the classical LSTM (CLSTM) model in expressive power, accuracy, and true-positive rate. The highest-performance QLSTM model is then used to identify transient events in 4XMM-DR14. Out of 40154 light curves in the 0.2--12 keV band, we detect 113 light curves with anomalies, or transient event candidates. This number is $\approx$ 1.3 times that of anomalies detectable with the CLSTM model. By utilizing SIMBAD and four wide-field survey catalogs made by ROSAT, SkyMapper, Pan-STARRS, and WISE, no possible counterparts are found for 12 detected anomalies.

The interaction of carbon atoms with solid carbon monoxide (CO) is a fundamental process in astrochemistry, influencing the formation of complex organic molecules in interstellar environments. This study investigates the adsorption and reaction mechanisms of carbon atoms on solid CO under cryogenic conditions, employing a combination of experimental techniques, including photostimulated desorption and resonance-enhanced multiphoton ionization (PSD-REMPI) and infrared spectroscopy, alongside quantum chemical calculations. The results reveal the formation of oxygenated carbon chains, such as CCO, C$_3$O$_2$, and C$_5$O$_2$, as well as CO$_2$. The findings highlight the role of chemisorption and subsequent reactions in driving molecular complexity on solid CO, with implications for the chemical evolution of interstellar ices and the potential formation of prebiotic molecules.

The timing (cross-)calibration of astronomical instruments is often done by comparing pulsar times-of-arrival (TOAs) to a reference timing model. In high-energy astronomy, the choice of solar system ephemerides and source positions used to barycenter the photon arrival times has a significant impact on the procedure, requiring a full reprocessing the data each time a new convention is used. Our method, developed as part of the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC), adapts an existing pulsar solution to arbitrary JPL ephemerides and source positions by simulating geocentric TOAs and refitting timing models (implemented with PINT). We validate the procedure and apply it to thousands of observations of the Crab pulsar from 14 missions spanning 2002--2025, demonstrating inter-ephemeris TOA consistency at the $\lesssim5\,\mu$s level, using the DE200/FK5-based Jodrell Bank Monthly Ephemeris as a reference. We release open-source tools (TOAextractor) and a TOA database to support future calibration and scientific studies. Instrument timing performance is broadly consistent with mission specifications; the X-ray-to-radio phase offset varies with energy and time at a level that is marginally compatible with the uncertainties of the radio ephemeris, motivating coordinated multiwavelength follow-up.

Core-collapse supernovae, occurring at the end of massive star evolution, produce heavy elements, including those in the iron peak. Although the explosion mechanism is not yet fully understood, theoretical models can reproduce optical observations and observed elemental abundances. However, many nuclear reaction rates involved in explosive nucleosynthesis have large uncertainties, impacting the reliability of abundance predictions. To address this, we have previously developed a Monte Carlo-based nucleosynthesis code that accounts for reaction rate uncertainties and has been applied to nucleosynthesis processes beyond iron. Our framework is also well suited for studying explosive nucleosynthesis in supernovae. In this paper, we investigate 1D explosion models using the "PUSH method", focusing on progenitors with varying metallicities and initial masses around $M_\mathrm{ZAMS} = 16 M_{\odot}$. Detailed post-process nucleosynthesis calculations and Monte Carlo analyses are used to explore the effects of reaction rate uncertainties and to identify key reaction rates in explosive nucleosynthesis. We find that many reactions have little impact on the production of iron-group nuclei, as these elements are primarily synthesized in the nuclear statistical equilibrium. However, we identify a few "key reactions" that significantly influence the production of radioactive nuclei, which may affect astrophysical observables. In particular, for the production of ${}^{44}$Ti, we confirm that several traditionally studied nuclear reactions have a strong impact. However, determining a single reaction rate is insufficient to draw a definitive conclusion.

Isotopic measurements of Solar System bodies provide a primary paradigm within which to understand the origins and histories of planetary materials. The D/H ratio in particular, helps reveal the relationship between (and heritage of) different H$_2$O reservoirs within the Solar System. Here we present interferometric maps of water (H$_2$O) and semiheavy water (HDO) in the gas-phase coma of a comet (Halley-type comet 12P/Pons-Brooks), obtained using the Atacama Large Millimeter/submillimeter Array (ALMA). The maps are consistent with outgassing of both H$_2$O and HDO directly from the nucleus, and imply a coma D/H ratio (for water) of $(1.71 \pm 0.44)\times10^{-4}$. This is at the lower end of the range of previously-observed values in comets, and is consistent with D/H in Earth's ocean water. Our results suggest a possible common heritage between a component of the Oort cloud's water ice reservoir, and the water that was delivered to the young Earth during the early history of the Solar System.

Trillium lattices, where magnetic ions form a three-dimensional chiral network of corner-sharing equilateral triangular motifs, offer a prominent platform to explore exotic quantum states. In this work, we report ground-state properties of the $S$ = 5/2 trillium lattice compound K$_{2}$FeSn(PO$_{4}$)$_{3}$ through thermodynamic, electron spin resonance (ESR), and muon spin relaxation (${\mu}$SR) experiments. Thermodynamic and ESR measurements reveal the two-step evolution of magnetic correlations across $T^{*}$ = 11 K, which results from an interplay between dominant antiferromagnetic Heisenberg interactions and subleading interactions. Below $T^{*}$, \textit{dc} and \textit{ac} magnetic susceptibilities indicate weak \textcolor{black}{magnetic ordering} at $T_{\rm N} \approx 2$ K under low fields, which is suppressed for $\mu_{0}H \geq 2$ T, consistent with a power-law dependence of magnetic specific heat at low temperatures. $\mu$SR experiments confirm the dominance of persistent spin dynamics and the absence of conventional spin freezing, supporting the subtle nature of weak magnetic ordering coexisting with spin-liquid-like fluctuations. These findings underscore the potential for realizing a classical spin-liquid ground state with exotic excitations in high-spin trillium lattice systems.

Methanol is a seed species of complex organic molecules that is of fundamental importance in astrochemistry. Although various isotopologues of CH$_3$OH have been detected in the interstellar medium (ISM), CH$_{3}$$^{17}$OH is only tentatively detected in Sgr~B2. To confirm the presence of CH$_{3}$$^{17}$OH in the ISM and to investigate its abundance, we search for its emission lines in the Orion~KL region. We have obtained image cubes covering the frequency ranges 236.40~GHz-236.65~GHz and 231.68~GHz-231.88~GHz using ALMA archival data observed toward the Orion~KL region. The column densities of CH$_3$$^{17}$OH and CH$_3$$^{18}$OH are estimated under the assumption of local thermodynamic equilibrium condition with fixed excitation temperatures at the two CH$_3$$^{18}$OH peaks, MeOH1 and MeOH2,. We have identified six emission lines of CH$_{3}$$^{17}$OH in MeOH1 and MeOH2 and confirmed that the line profiles and spatial distributions are consistent with those of CH$_3$$^{18}$OH. The abundance ratios of CH$_3$$^{18}$OH/CH$_3$$^{17}$OH are evaluated to be $\sim 3.4-3.5$ and are similar to the canonical value of $^{18}$O/$^{17}$O $\sim 3-4$ derived from CO observations in the Orion~KL region. We have compared the results with the previous study of CH$_3$OH and evaluated CH$_3$$^{16}$OH/CH$_3$$^{17}$OH ratios to be $\sim 2300-2500$ at a resolution of $\sim 4$~arcsec. The ratios are close to the $^{16}$O/$^{17}$O ratio in the local ISM. This result indicates that the CH$_3$OH isotopologues can serve as new tracers of oxygen isotope ratios in star-forming regions because the opacity of CH$_3$OH can be evaluated using transition lines spanning a wide range of line intensities. Moreover, this method enables us to study the star-formation history of our Galaxy with the aid of the Galactic chemical evolution models.

Protostellar multiplicity is a common outcome of the star formation process. To fully understand the formation and evolution of these systems, the physical parameters of the molecular gas together with the dust must be systematically characterized. Using observations of molecular gas tracers, we characterize the physical properties of cloud cores in the Perseus molecular cloud (average distance of 295 pc) at envelope scales (5000-8000 AU). We used Atacama Pathfinder EXperiment (APEX) and Nobeyama 45m Radio Observatory (NRO) observations of DCO$^+$, H$_2$CO and c-C$_3$H$_2$ in several transitions to derive the physical parameters of the gas toward 31 protostellar systems in Perseus. Gas kinetic temperature was obtained from DCO$^+$, H$_2$CO and c-C$_3$H$_2$ line ratios. Column densities and gas masses were then calculated for each species and transition. Gas kinetic temperature and gas masses were compared with bolometric luminosity, envelope dust mass, and multiplicity to search for statistically significant correlations. Gas kinetic temperature derived from DCO$^+$, H$_2$CO and c-C$_3$H$_2$ line ratios have average values of 14 K, 26 and 16 K, respectively, with a range of 10-26 K for DCO$^+$ and c-C$_3$H$_2$. The gas kinetic temperature obtained from H$_2$CO line ratios have a range of 13-82 K. Column densities of all three molecular species are on the order of 10$^{11}$ to 10$^{14}$ cm$^{-2}$, resulting in gas masses of 10$^{-11}$ to 10$^{-9}$ M$_{\odot}$. Statistical analysis of the physical parameters finds: i) similar envelope gas and dust masses for single and binary protostellar systems; ii) multiple (>2 components) protostellar systems tend to have slightly higher gas and dust masses than binaries and single protostars; iii) a continuous distribution of gas and dust masses is observed regardless of separation between components in protostellar systems.

We report the first statistically significant detection of X-ray polarization from the high-mass X-ray binary (HMXB) 4U 1700-377, observed using the Imaging X-ray Polarimetry Explorer (IXPE). A polarization degree exceeding 10% was detected above 5 keV, placing it among the highest polarization observed in HMXBs to date. The observation was conducted over a full orbital period of the binary system, during which several sporadic and instantaneous flares were detected. We identify a clear correlation between the polarization degree and orbital phase, with the highest polarization occurring just before and after the eclipse, reaching over 20% for a few tens of ks. These results suggest that the scattering medium responsible for the observed polarization is spatially localized between the compact object and the O-type companion star, likely created by large-scale inhomogeneities in the stellar wind and its interaction with the compact star's emission. We also explore the roles of disk winds and orbital reflection in the observed polarization variability. While both mechanisms contribute to the polarization, the substantial increase in polarization before and after the eclipse cannot be fully explained by these models alone, suggesting that the involvement of additional factors. The properties of the X-ray polarization observed by IXPE provide new insights into the accretion processes, X-ray emission, and wind structure in 4U 1700-377, advancing our understanding of their complex environments and the nature of the compact objects within.

We explore a pulse-based variational quantum eigensolver (VQE) algorithm for Rydberg atoms in optical tweezer arrays and evaluate its performance on prototypical quantum spin models. We numerically demonstrate that the ground states of the one-dimensional antiferromagnetic Heisenberg model and the mixed-field Ising model can be accurately prepared using an adaptive update algorithm that randomly segments pulse sequences, for systems of up to ten qubits. Furthermore, we propose and validate a hybrid scheme that integrates this pulse-level analog quantum algorithm with a variational quantum gate approach, where digital quantum gates are approximated by optimized analog pulses. This enables efficient measurement of the cost function for target many-body Hamiltonians.

Cosmic rays are often modeled as charged particles. This allows their non-ballistic propagation in magnetized structures to be captured. In certain situations, a neutral cosmic ray component can arise. For example, cosmic ray neutrons are produced in considerable numbers through hadronic pp and p$\gamma$ interactions. At ultrahigh energies, the decay timescales of these neutrons is dilated, allowing them to traverse distances on the scale of galactic and cosmological structures. Unlike charged cosmic rays, neutrons are not deflected by magnetic fields. They propagate ballistically at the speed of light in straight lines. The presence of a neutral baryonic cosmic ray component formed in galaxies, clusters and cosmological filaments can facilitate the escape and leakage of cosmic rays from magnetic structures that would otherwise confine them. We show that, by allowing confinement breaking, the formation of cosmic-ray neutrons by high-energy hadronic interactions in large scale astrophysical structures can modify the exchange of ultra high-energy particles across magnetic interfaces between galaxies, clusters, cosmological filaments and voids.

We present new two-dimensional radiation hydrodynamic simulations of supernova shock breakout from red supergiants using the $\texttt{CASTRO}$ code. Our progenitors are 20 and 25 M$_{\odot}$ solar-metallicity stars evolved from the zero-age main sequence with $\texttt{MESA}$ and exploded in one dimension using $\texttt{FLASH}$. We consider a range of circumstellar media (CSM) produced by stellar winds to investigate how pre-explosion mass-loss affects shock breakout. The multigroup flux-limited diffusion scheme in $\texttt{CASTRO}$ captures the interaction between the explosion shock, its radiation precursor, and the surrounding CSM. We find that strong radiation precursors, generated by radiation leakage behind the shock, can drive fluid instabilities and move the effective photosphere outward before the shock reaches the stellar surface. The resulting breakout emissions reach peak luminosities of ${\sim}10^{44}$ erg s$^{-1}$ with full-width half-maximum durations of 1-3 hr, which are much dimmer and longer than those from blue supergiants. The light-curve colors gradually evolve from blue to red after the peak. The 25 M$_{\odot}$ model with explosion energy $E \sim 1.69\times10^{51}$ erg produces ${\sim}$10-30\% higher maximum luminosity than the 20 M$_{\odot}$ model with $E \sim 1.09\times10^{51}$ erg. The dense CSM further extends the breakout rise time by increasing the photon diffusion. These results provide new constraints on red supergiant atmospheres and mass-loss histories prior to core collapse.

Supernova (SN) 1987A provides a unique window into the aftermath of a massive stellar explosion, offering key insights into the ejecta's morphology, composition, explosion mechanism, progenitor system, and circumstellar medium (CSM) interaction. We investigate large-scale ejecta asymmetries in SN 1987A. By comparing the simulations with JWST observations and making predictions for XRISM, we aim to refine our understanding of the explosion mechanism and the remnant's evolution. We performed 3D MHD simulations that trace the evolution of SN 1987A from the SN to the SNR, extending our predictions up to 5000 years into the future and considering the Ni-bubble effects. The simulation results are compared with JWST observations and used to predict XRISM spectra, to evaluate the accuracy of the modeled ejecta structure. Our simulations reproduce the large-scale Fe-rich ejecta morphology seen by JWST, revealing two clumps suggestive of a bipolar explosion. Ni-bubble effects in the first year boost Fe-rich ejecta expansion and their interaction with the reverse shock. However, discrepancies with JWST observations in clump velocities and spatial distribution suggest stronger explosion asymmetries than modeled. Since 2021, our models predict that shocked ejecta have contributed increasingly to X-ray emission, now rivaling shocked CSM and soon dominating as the latter fades. Future XRISM observations will trace the evolution of these ejecta structures, refining constraints on explosion geometry. Early remnant asymmetries from CSM interaction may persist for at least 100 years. Our results underscore the role of asymmetric core-collapse mechanisms in shaping SN 1987A's ejecta and constraining its explosion geometry. Future studies should explore more extreme asymmetries, in neutrino-driven core collapse or magneto-rotational SN models, to identify the origin of its bipolar Fe-rich structure.

Quantum light sources using single-walled carbon nanotubes show promise for quantum technologies but face challenges in achieving precise control over color center formation. Here we present a novel technique for deterministic creation of single organic color centers in carbon nanotubes using \textit{in-situ} photochemical reaction. By monitoring discrete intensity changes in photoluminescence spectra, we achieve precise control over the formation of individual color centers. Furthermore, our method allows for position-controlled formation of color centers as validated through photoluminescence imaging. We also demonstrate photon antibunching from a color center, confirming the quantum nature of the defects formed. This technique represents a significant step forward in the precise engineering of atomically defined quantum emitters in carbon nanotubes, facilitating their integration into advanced quantum photonic devices and systems.

Phase separation, the spontaneous segregation of density, is a ubiquitous phenomenon observed across diverse physical and biological systems. Within a crowd of self-propelled elements, active phase separation emerges from the interplay of activity and density interactions. A striking feature of active phase separation is the persistent formation of dilute-phase bubbles within the dense phase, which has been explored in theoretical models. However, the fundamental parameters that systematically control bubble formation remain unclear in conventional active particle models. Here, we introduce an active binary mixture model, where active solutes and solvents dynamically exchange their positions, and find that spontaneous bubble formation can be tuned by the asymmetry in their activities. Numerical simulations reveal that moderate solvent activity enhances bubble formation, while larger solvent activity, comparable to solute activity, suppresses it. By employing mean-field theory, which captures essential phase behaviors, we consider the mechanism for the enhancement of bubble formation induced by solvent activity. Beyond these findings, when solute and solvent activities are equal, numerical simulations indicate that critical phenomena of active phase separation under the suppression of bubbles belong to the Ising universality class. Our findings establish activity asymmetry as a key control parameter for active matter phase transitions, offering new insights into universality in nonequilibrium systems.

Context. Deriving the mass and large-scale asymmetries of radioactive isotopes offers valuable insights into the complex phases of a supernova explosion. Important examples are $^{56}$Ni, with its decay products $^{56}$Co and $^{56}$Fe, and $^{44}$Ti, which are studied through their X-rays emission lines and provide a powerful diagnostic tool to probe the explosive nucleosynthesis processes in the inner layers of the exploding star. Aims. In this framework, SN 1987A provides a privileged laboratory being the youngest supernova remnant from which the mass of Ti has been estimated. However, some tension exists in determining the initial mass of $^{44}$Ti. Previous analysis, relying on \textit{NuSTAR} and \textit{INTEGRAL} data, report $M_{44} = (1.5 \pm 0.3) \times 10^{-4}$ $M_\odot$ and $M_{44}=(3.1 \pm 0.8) \times 10^{-4} M_\odot$, respectively. In this paper we estimate the initial mass of $^{44}$Ti via its decay product, the $^{44}$Sc emission line at 4.09 keV, using \textit{Chandra} observations. Methods. We perform multi-epoch spectral analysis focusing on the inner part of the remnant, to minimize the contamination from the X-ray emission stemming from the shocked plasma. As a result, we provide the detection of $^{44}$Sc emission line in the central part of SN 1987A with a $\sim$99.7\% (3 $\sigma$) significance. Results. The simultaneous fit of the spectra extracted from observations between 2016 and 2021 provides a line flux of $6.8^{+2.2}_{-2.3}\times 10^{-7}$ photons s$^{-1}$ cm$^{-2}$ corresponding to a $^{44}$Ti mass $M_{44}=(1.6\pm0.5) \times 10^{-4} M_\odot$ (errors at the $90\%$ confidence level). The results obtained with our spectral analysis seem to align with those derived with NuSTAR.