This is the list of the papers for the past 5 days that include local authors affiliated with Princeton University's Astrophysical Sciences department.
We present JWST MIRI and NIRCam imaging of the inner ~5 kpc of the M82 superwind at 0.05-0.375'' (~0.9-6.5 pc) resolution. Targeted filters probe emission from polycyclic aromatic hydrocarbons (PAHs; F335M, F360M, F770W, F1130W) and continuum (F250M, F360M). The images reveal a network of cool wind filaments traced by PAHs. PAH surface brightness declines with the inverse square of distance to the midplane, suggesting that the incident radiation field from the starburst drives the observed PAH intensity out to 2.5 kpc. The 3.3/11.3 and 3.3/7.7 band ratios show uniformity with distance from the starburst, though comparisons with mid-IR dust emission models indicate a modest shift toward larger PAHs. Outside the disk, 11.3/7.7 increases moderately, reflecting that PAHs become more neutral with distance from the starburst as they are exposed to a declining radiation field and ionization parameter. Overall, PAHs in the wind are consistent with standard-to-large sizes and standard-to-high ionization states. Including Spitzer and Herschel data, PAH abundance (qPAH) is set at ~1% in the starburst and remains unchanging out to 5 kpc off the disk. This flat qPAH profile suggests that PAHs are shielded from the hot wind, perhaps residing in the surface layers of cool clouds, with possible replenishment from cloud interiors and enrichment of the halo from previous bursts. In this picture, clouds are not dense enough to promote PAH growth, and they likely undergo radiative cooling and mixing with the hot phase to survive the gauntlet for at least ~20 Myr.
Protostellar winds can theoretically lift solids from the planet-forming disks, but direct evidence for launched dust has been scarce so far. Numerous atomic lines that are unique to mid-infrared (IR) wavelengths reveal refractories eroded from dust grains and provide information on wind properties in the earliest stages of the star formation process. We present JWST/MIRI-MRS spectral imaging of the inner 2000 au of the BHR71-IRS1 blueshifted side of the outflow. Atomic line intensities are compared to shock models to constrain the physical conditions and elemental abundances of the outflowing gas. Dust continuum maps are constructed from PSF-subtracted cubes, and the dust spectral energy distribution is analyzed. The ionized central jet of BHR71-IRS1 is spatially resolved and imaged for the first time, revealing a unique inventory of refractory, volatile, and noble-gas fine-structure lines (Fe, Ni, Co, Cl, S, Ne, Ar). The emission is concentrated along four bright knots that wiggle along the jet axis. PSF-subtracted continuum maps reveal extended mid-IR continuum emission co-spatial with the jet bullets and within the H$_2$-traced outflow cone. Spectral energy distributions along the jet are fit together with the extinction, revealing a warm (200-400 K) and a cold (70-90 K) dust component. Shock modeling constrained by the mid-IR lines indicates a decline in shock velocity from 70 to 35 km s$^{-1}$ and pre-shock density from $>$10$^5$ to $ 4\times 10^4$ cm$^{-3}$ with distance from the protostar. Gas-phase Fe and Ni are measurably depleted relative to Solar abundances, consistent with a substantial fraction of refractories remaining locked in grains in spite of the shocks. These JWST observations provide direct evidence that dust is launched in a Class 0 jet and at least partly survives shock processing.
The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization universe as a means of measuring large-scale cosmic structure. Using 9 nights of observations targeting the North Celestial Pole (NCP) field, totaling approximately 107 hours of integration time, we analyze data in the frequency range 700-800 MHz (corresponding to redshift $z \sim 0.9$). We do the data format conversion, radio frequency interference (RFI) flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition (SVD). The spherically averaged power spectrum $\Delta^2(k)$ is obtained. This work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.
While photoionization rates of Ganymede's O$_2$ dominated atmosphere are well constrained, the contribution of electron-impact ionization is rather uncertain. Previous quantitative estimates have relied on assumptions about densities and energy distributions of precipitating electrons, or on rare spacecraft measurements that cannot be unambiguously mapped to the regions of ionization. In this study, we present a novel approach to quantify electron-impact ionization rates directly through OI 1356 Å emission brightness observations. The analysis of measured cross sections reveals that the ionization-to-excitation ratio is limited to 10-60 over all electron energies, reducing the uncertainty of estimating ionization rates to a factor less than 6. We apply this method to Juno UVS observations of Ganymede's aurora. We find that the OI 1356 Å brightness of the auroral ovals is well described by 3-5° latitude wide Gaussian distributions centered on the open-closed field line boundary, with an average peak of 120 R. The average brightness outside the ovals in the polar and equatorial background regions is ~8 R. From these observations, we derive a global map of electron-impact ionization rates, which are at least an order of magnitude higher than photoionization rates. The estimated total global ionization rate is 1.3-7.6$\times$10$^{26}$ s$^{-1}$, with average column rates of ~5$\times$10$^{9}$ cm$^{-2}$s$^{-1}$ in the ovals and ~3$\times$10$^{8}$ cm$^{-2}$s$^{-1}$ in the background regions. Comparison of radio occultation measurements with predicted electron densities indicates that transport processes are the dominant loss mechanism in Ganymede's ionosphere. The rate of ionospheric outflow of O$_2^+$ is 0.1-2$\times$10$^{26}$ s$^{-1}$ or 0.5-11 kg s$^{-1}$, indicating 0.03-0.5 cm Myr$^{-1}$ erosion of Ganymede's surface ice.
The temporal and spatial coincidence between the gravitational wave (GW) event GW190425 and the fast radio burst (FRB) event FRB 20190425A raises the intriguing possibility of a physical connection between the two. The widely discussed possibility invoking the collapse of a supermassive neutron star as the merger product suffers the inconsistency between the model prediction and the measured inclination angle of the system. Here, we propose a novel physical mechanism to account for the association. We envisage a magnetar located at about 2.5 light hours away from the binary neutron star merger site. The kiloherz GWs generated by the merger are converted into kiloherz electromagnetic (EM) radiation via the Gertsenshtein-Zeldovich (GZ) effect near the magnetar. Subsequent inverse Compton scattering off the kilohertz EM waves by relativistic particles generates the observed gigahertz FRB emission. Our calculation reveals that, with appropriate parameter choices, the properties of FRB 20190425A can be reproduced.
Little Red Dots (LRDs) are a population of compact, red sources that have emerged as one of the most puzzling findings of JWST. Variability provides a direct probe of their central engines. Here we present the first joint spectroscopic and photometric time-domain study of LRDs undertaken with the JWST TWINKLE slitless spectroscopy program. Surveying the FRESCO GOODS-North legacy field, TWINKLE monitors a complete, H$\alpha$-flux-limited sample of 18 LRDs at z = 3.9-6.8, achieving a rest-frame baseline of $\sim$140-220 days. We detect no variability in photometry, H$\alpha$ line flux, or line shape across the sample. If LRDs resembled AGN in reverberation mapping samples -- the foundation for black hole mass calibrations and luminosity scaling relations -- we would expect >10 sources to show measurable fluctuations. Observing none implies a 5.9$\sigma$ deficit. The non-detections hold across all broad H$\alpha$ emitters within TWINKLE's field of view -- the 18 V-shaped LRDs as well as 9 non-LRDs. Comparison with simulated light curves disfavors sub-Eddington accretion and is instead consistent with super-Eddington accretion, other mechanisms that suppress variability, or perhaps no AGN whatsoever. If LRDs do harbor black holes, calibrations derived from sub-Eddington systems may not apply, thereby explaining JWST's apparently "overmassive" black holes. These observations provide unique constraints on the physics of one of the most enigmatic populations discovered by JWST.
Reddened Type-1 quasars offer a unique window into the structure and evolution of active galactic nuclei (AGN), yet their physical origin and the source of their reddening remain uncertain. Optical surveys often miss these dust-obscured objects, resulting in an incomplete view of the quasar population. In this work, we construct a sample of 6,600 Type-1 quasars at redshifts $0.5 \leq z \leq 1.2$ by combining deep optical imaging from HSC with mid-infrared photometry from WISE, enabling a more complete selection that is not biased against reddened objects. We perform detailed SED modeling using the CIGALE code, enhanced by synthetic photometry derived from SDSS spectra to better constrain the optical continuum. We classify quasars into blue and reddened Type-1 populations based on their continuum slopes and compare their SEDs and emission line properties. As expected from this definition, reddened Type-1 AGN show higher dust extinction, with a median $A_V = 0.60^{+0.32}_{-0.19}$ mag, compared to $A_V = 0.06^{+0.10}_{-0.03}$ mag for blue objects. But they also exhibit smaller torus half-opening angles, with a median of $25.7^{+10.1}_{-8.7}$ deg, compared to $33.3^{+11.1}_{-5.9}$ deg for blue objects. While such extinction could arise on either galaxy or nuclear scales, the systematically stronger narrow-line equivalent widths and weaker Balmer broad lines in reddened Type-1s indicate that the obscuration acts on nuclear scales, likely from dust concentrated near the polar axis. We discuss the possibility that these structural differences may be linked to a sub-pc outflow, that carries dusty gas into the polar region and evacuates the torus region.
The YORP effect plays an important role in the spin evolution of asteroids. Although craters are ubiquitous surface features, their influence on YORP torque has received limited attention. In this paper, we investigate the YORP torque of a circular crater on a spherical asteroid, focusing specifically on how lateral thermal conduction breaks symmetry to produce a net torque. Using three-dimensional finite element simulations, we calculate the resulting spin and obliquity accelerations and examine their dependence on the crater's location, depth, and thermal parameters. Our results show that the crater-induced spin torque is consistently positive, and craters at different latitudes drive the spin axis toward obliquity equilibria at 0, 90 or 180 degree. We demonstrate that the spin torque arises primarily from the lateral heat conduction inside the asteroid that occurs only in 3D model, while the contributions from self-heating and shadowing effects are negligible. While the YORP effect induced by internal heat conduction may be overtaken by torque components arising from shadowing and crater orientation, particularly on large asteroids, our numerical results show that for small craters, this spin torque amounts to approximately 10 to 100 percent of the normal YORP torque. Its persistent positivity may help explain the observed prevalence of positive spin accelerations in asteroids.
The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB~170817A, heralded the birth of multi-messenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV--TeV GW follow-up strategies for CTAO using a multi-step simulation pipeline and to estimate the expected rate of joint GW-GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV--TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star-black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.
Supernova (SN) 1987A provides a unique laboratory for investigating many aspects of SN physics and evolution. An observation at Day 12927 (35.4 yr) since the explosion with the Mid-Infrared Instrument (MIRI) Medium Resolution Spectrometer (MRS) on the James Webb Space Telescope (JWST) provided the first spatially resolved spectroscopic study of SN 1987A in the mid-IR, yielding insights into the evolution of dust, the ejecta, the equatorial ring (ER), and shocks in the system. Here we present a second epoch with MIRI/MRS at Day 13311 (36.4 yr) allowing the mid-IR spatially resolved spectroscopic temporal evolution of SN 1987A to be probed for the first time. Analysis of the ER-dominated dust continuum showed little evolution between Days 12927 and 13311. However, a spatial analysis reveals the inner ER to be fading while the outermost regions are brightening. Broad ejecta emission lines detected at Day 12927 are evolving rapidly, driven by the recent onset of the ejecta/equatorial ring interaction in the northeast and southwest of the ER. Most lines from the ER show no change during the 384 days between the epochs, though some such as [Ne II] and [Ar II] have faded. We identify mid-IR H2 emission associated with the ejecta for the first time. Using the near- and mid-IR [Fe II] lines as density and temperature diagnostics of the ejecta in the interaction region we find it likely that the dense inner Fe-rich ejecta has now reached the reverse shock. Continued monitoring of SN 1987A is essential to observe the evolving ejecta/ER interaction and dust components.
We provide a comprehensive analysis of GW190711_030756 and GW200114_020818, two of the most significant binary black hole merger candidates in the IAS catalog, with probabilities of astrophysical origin $p_{\rm astro}=0.99$ and $0.71$, respectively, and signal-to-noise ratios of approximately $10.0$ and $13.4$. We employ numerical relativity surrogate models to infer both the source properties and the remnant properties of these two candidates. We find that both GW190711_030756 and GW200114_020818 are asymmetric-mass binaries, with inferred mass ratios of $0.35^{+0.32}_{-0.15}$ and $\leq 0.20$. In addition, GW200114_020818 is inferred to have a source-frame total mass of approximately $220M_{\odot}$ and highly spinning black holes, with primary (secondary) dimensionless spin magnitudes of $0.96^{+0.03}_{-0.07}$ ($0.84^{+0.13}_{-0.34}$), closely resembling GW231123_135430. We further find that GW200114\_020818 has a confidently negative effective inspiral spin of $\chi_{\rm eff}=-0.60^{+0.22}_{-0.13}$ and exhibits strong spin precession, characterized by an effective precession parameter of $\chi_{\rm p}=0.60^{+0.21}_{-0.19}$. GW200114_020818 (when considered alongside GW231123_135430) points towards an emerging population of massive, rapidly spinning BBH mergers. While GW231123_135430 is consistent with mergers in globular clusters, producing systems like GW200114_020818 in such environments remains difficult even under hierarchical merger scenarios. The probability that the remnant black hole of GW190711_030756 (GW200114_020818) is retained in its host environment is $0.079$ ($0.0002$), $0.62$ ($0.965$), and $0.997$ ($1$) if the merger occurred in a globular cluster, a nuclear star cluster, or an elliptical galaxy, respectively.
The Cosmological Collider (CC) provides a unique opportunity to probe the particle spectrum and fundamental interactions at extremely high energies. Massive particles, via their decay into inflaton quanta, can induce a non-analytic, oscillatory, primordial non-Gaussianity (NG), including the bispectrum. At tree level, three classes of such processes contribute to the bispectrum: 'single exchange', 'double exchange', and 'triple exchange', depending on the number of massive particle propagators. We provide a unified evaluation of all three diagrams and derive the explicit shape functions for the bispectrum, valid across the entire kinematic space. We perform a search for these three processes with the Planck data, finding no evidence for NG. We also consider simple extensions of the minimal scenario that can counter the exponential suppression of the non-analytic signature, and produce on-shell particles with masses $M\gg H$, the Hubble scale during inflation. In particular, we focus on the 'scalar chemical potential' mechanism and extend our previous search to a wider range of chemical potential ($\omega$) and $M$, finding global 1.5$\sigma$ evidence for non-zero NG for the parameter space $\omega - M \simeq 3H$.
The pulsar timing array (PTA) is a powerful technique for detecting nanohertz gravitational wave backgrounds (GWBs). However, conventional PTAs lack sensitivity to parity violation in the GWB. In this work, we propose a dipole pulsar timing array system (dPTA). By deriving the overlap reduction functions (ORFs) from the cross-correlation of timing signals, we find that this system exhibits sensitivity to chiral GWBs in the nanohertz regime. Furthermore, through numerical calculations of its sensitivity curves, we demonstrate that the dPTA extends the detectable frequency range of PTAs for GWBs from the nanohertz to the microhertz regime.
Supernova remnants (SNRs) can strongly affect the chemical composition of the interstellar dust. In this paper we investigate to what degree the dust and ices are modified by observing four stars expected to be absorbed by a giant molecular cloud interacting with SNR W44, using medium-resolution spectroscopy in 2-5 $\mu$m. Absorption from H2O ice around 3.0 $\mu$m and aliphatic hydrocarbon dust around 3.4 $\mu$m were detected towards two stars, while probable CO ice at 4.67 $\mu$m towards one of them. Millimeter gas-phase CO J = 1-0 lines and three-dimensional dust extinction maps show that the dense molecular gas associated with W44 dominates (> 60%) the total interstellar extinction (A_K ~ 2.6) along these two sightlines. The H2O ice column densities are a factor of 1.5-3 lower than nearby MCs at similar extinctions, possibly because of the destruction of ice by shocks and cosmic rays (CRs) from W44, consistent with the low CO ice abundance relative to H2O (< 12%). One of the sightlines shows an unusually strong 3.4 $\mu$m aliphatic hydrocarbon absorption. If the carriers are located in diffuse dust along the sightline, unrelated to W44, its strength is ~ 4 times larger than those typically observed for diffuse dust clouds. Alternatively, the carriers may be enhanced in the W44 environment. We discuss several possible explanations, including shock formation of aliphatic hydrocarbons in diffuse clouds associated with W44, contribution from aliphatic hydrocarbons in shocked and CR-bombarded molecular clouds, and changes in the extinction law due to the SNR interaction.
The next nearby core-collapse supernova (SN) promises to yield a treasure of scientific information through multi-messenger astronomy. Early observations of the shock breakout (SBO) emissions are especially critical to understand the SN explosive mechanism as well as the properties of the progenitor star. Neutrino observatories are able to provide an early alert of a SN before the arrival of the SBO radiation. Super-Kamiokande (SK) has the unique capability to independently reconstruct an accurate SN pointing direction as part of its real-time monitoring system, ``SNWATCH.'' Recent upgrades to SK by adding gadolinium (Gd) to the detection volume have been accompanied by efforts to improve the speed and accuracy of SN direction reconstruction. A new, novel HEALPix-based approach (``HP-Fitter'') can calculate the SN direction from the reconstructed burst event directions in less than one second. As well, the previous maximum-likelihood direction fitter (``ML-Fitter'') was upgraded by incorporating event information from Gd neutron-capture as well as using the HP-Fitter for the initial fit parameters and from code refactoring and optimization. The improved ML-Fitter has better angular resolution but direction reconstruction time is $\mathcal{O}$(sec). Together with improvements in burst detection and event reconstruction times, SNWATCH is now able to generate an SN alert with pointing information in about 90 seconds. These upgrades have been implemented at SK and integrated into a new automated system to provide GCN notices.
We report the discovery of six faint millisecond pulsars (MSPs) in the globular clusters NGC 6517 and NGC 7078 (M15) using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These discoveries were enabled by stacking power spectra from multiple observations, a method that effectively boosts the signal-to-noise ratio of faint sources. In NGC 6517, we identified four new MSPs (NGC 6517S-V) with spin periods ranging from 3.68 to 6.02 ms and dispersion measures (DMs) between 182.45 and 182.85 pc cm^-3. In M15, two additional MSPs (M15M and M15N) were discovered, with spin periods of 4.83 and 9.28 ms, and DMs of 67.89 and 66.65 pc cm^-3, respectively. A phase-coherent timing solution has been obtained for M15M; however, sparse detection rates currently preclude phase-connected solutions for the remaining five pulsars. Current timing parameters suggest all six MSPs are isolated, which is consistent with the expected pulsar populations in core-collapsed globular clusters. Notably, pulsars M15N, NGC 6517U, and NGC 6517V eluded detection by standard frequency-domain searches (e.g., PRESTO-based) and the Fast Folding Algorithm, demonstrating that the stack search technique significantly enhances detection sensitivity to inherently faint pulsar signals.
Semi-analytic models (SAMs) have been treating galaxy populations as dynamical systems for $\gtrsim50$ years, but their evolution equations remain poorly constrained. We introduce sapphire, a modular, automatically differentiable, GPU-accelerated SAM written from scratch in JAX. For the first time, we compute exact Jacobian matrices of our nonlinear differential equations and show that they have interpretable, non-random structures, using the Pandya et al. (2023) physical model as an initial example. Both local and global sensitivity analyses reveal that supernova energy loading is a key astrophysical parameter for galaxy evolution. We use gradient descent and Hamiltonian Monte Carlo (HMC) to perform comprehensive mock parameter recovery tests. These indicate that the z=0 stellar-to-halo-mass relation alone does not contain enough information to infer many astrophysical parameters. Using observations of star-forming galaxies from the MaNGA survey and the Behroozi et al. (2019) empirical model as one baseline, we derive multiple posteriors assuming different combinations of data, including z=0 interstellar medium gas fractions and metallicities. The inferred physical parameters suggest that galaxies self-regulate their star formation primarily through preventative rather than ejective feedback. Both Fisher and HMC forecasts demonstrate the potential of sapphire to enable precision inference for galaxy formation, but more work is needed to expand its library of models. We discuss how our unique blend of differentiability, massive GPU parallelization, numerical robustness and principled Bayesian methods sets the stage for hybrid physics-informed, data-driven discovery of galaxy formation astrophysics and cosmology. We make sapphire publicly available at this https URL.
The structures and dynamics of the interstellar medium are governed by a combination of self-gravity, external gravity, and various sources of ordered and random motions on different spatial scales. This paper uses ALMA CO (3-2) observations at 0.1" $\approx$ 5 pc resolution to examine the scale dependence of molecular gas structure and dynamics in the central molecular zone (CMZ) of a nearby galaxy, NGC 3351. We use the dendrogram technique to characterize hierarchical molecular gas structures spanning two decades in spatial scales and measure their size, gas mass, and velocity dispersion. Their size-linewidth relation shows a power-law slope of 0.58, comparable to measurements for CMZs in other galaxies and suggestive of significant contribution from ordered motion on large scales. We further decompose the observed velocity dispersion in each gas structure into ordered versus random motions. The former appears stronger in gas structures at $\gtrsim$ 30 pc while the latter becomes more dominant at $\lesssim$ 30 pc. Modulo uncertainties with the CO-to-H$_2$ conversion factor, the estimated gravitational free-fall time is comparable to the crossing time of ordered motions for structures on all spatial scales, and both becomes longer than the crossing time of random motions at small, $\lesssim$ 10 pc scales. Our results highlight the varying sources and drivers of gas motions on different spatial scales in the CMZ of a Milky Way-like galaxy.
Improving predictions of the geomagnetic impact of coronal mass ejections (CMEs) requires understanding how solar source properties relate to in-situ measurements at Earth. However, major geomagnetic storms frequently arise from interacting CMEs, complicating the link back to their solar origins. We analyze a CME interaction event that caused a major geomagnetic storm in 2024 October 10-11 (D$_{st}$ $\sim$-333 nT). Multiviewpoint observations reveal that the storm was related to a sympathetic eruption involving a quiescent filament and an active-region CME. The coronagraph on board the Advanced Space-based Solar Observatory clearly shows that this sympathetic eruption resulted in two distinct CMEs. Due to the overlap of the CMEs in the coronagraph field of view (FOV), a spheroid shock model was used to fit the observed shock. Kinematic analysis indicates that the interacting CMEs had completed their impulsive acceleration phase before entering the coronagraph FOV, with a slow deceleration continuing beyond 100 R$_\odot$. In-situ measurements indicate that the enhanced southward magnetic fields, arising from compression during CME interactions, were the primary driver of the storm. Compared to photospheric fields, the in-situ magnetic fields suggest that the trailing CME maintained flux-rope-like signatures consistent with the source region. In contrast, the compressed leading CME displayed varying magnetic configurations between Wind and STEREO-A, featuring distorted flux-rope signatures and inconsistent inferred axis orientations. Our study bridges solar source dynamics to in-situ multipoint measurements, providing key insights for space weather prediction. Nevertheless, the direct linkage between source-region magnetic field configurations and these measurements remains tentative and requires further investigation.
The chemical abundance of host stars plays a pivotal role in shaping the formation history of planetary systems, yet the influence of elements beyond iron remains poorly understood. Here, we investigate the relationship between the carbon-to-oxygen (C/O) ratio of host stars and the orbital periods of giant planets. By analyzing high-resolution spectroscopic data from 598 planet-hosting stars (hosting 929 planets) across SDSS, Keck, and HARPS surveys, we identify a correlation: stars with higher C/O ratios are more likely to host longer-period giant planets. Theoretical models of pebble-driven planet formation and migration further support this observation, demonstrating that elevated C/O ratios enhance solid material availability at outer disk regions, promoting giant planet formation at larger distances and subsequent moderate inward migration. Our findings establish stellar C/O as a critical factor in shaping the orbital architecture of giant planets, bridging disk chemistry to planetary system evolution.
We incorporate an arbitrarily high-order method for the Laplacian operator into the Spectral Difference method (SD). The resulting method is capable of capturing shocks thanks to its a-posteriori limiting methodology, and therefore it is able to survive scenarios in which the dissipative scales (viscous and diffusive) are not properly described. Moreover, it is capable of capturing these scales at lower resolution compared to lower-order methods and therefore attains convergence at lower resolution. We show that the method at hand has exponential convergence when describing smooth solutions and is able to recover a high-order solution when solving the dissipative scales.