More Like this

1. The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

   https://doi.org/10.3847/1538-4357/ad5291  

   Larsen, Bjorn; Mingarelli, Chiara_M_F; Hazboun, Jeffrey_S; Chalumeau, Aurélien; Good, Deborah_C; Simon, Joseph; Agazie, Gabriella; Anumarlapudi, Akash; Archibald, Anne_M; Arzoumanian, Zaven; et al (August 2024, The Astrophysical Journal)

   Abstract Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from ${\log }_{10}{A}_{\mathrm{RN}}=-{14.1}_{-0.1}^{+0.1}$to $-{14.7}_{-0.5}^{+0.3}$, and the spectral index broadens from ${\gamma }_{\mathrm{RN}}={2.6}_{-0.4}^{+0.5}$to ${\gamma }_{\mathrm{RN}}={3.5}_{-0.9}^{+1.2}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. 

   more » « less
2. HSTPROMO Internal Proper-motion Kinematics of Dwarf Spheroidal Galaxies. I. Velocity Anisotropy and Dark Matter Cusp Slope of Draco

   https://doi.org/10.3847/1538-4357/ad571c  

   Vitral, Eduardo; van_der_Marel, Roeland_P; Sohn, Sangmo_Tony; Libralato, Mattia; del_Pino, Andrés; Watkins, Laura_L; Bellini, Andrea; Walker, Matthew_G; Besla, Gurtina; Pawlowski, Marcel_S; et al (July 2024, The Astrophysical Journal)

   Abstract We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value $\overline{{\beta }_{\mathrm{B}}}=-{0.20}_{-0.53}^{+0.28}$, (where 1 – ${\beta }_{\mathrm{B}}\equiv 〈v_{\tan }^2〉/〈v_{\mathrm{rad}}^2〉$in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γ_dark, is $-{0.83}_{-0.37}^{+0.32}$, consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco’s low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1σ, 2σ, and 3σconfidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is ${24.19}_{-2.97}^{+6.31}\mathrm{km}{\mathrm{s}}^{-1}$. We infer a dynamical distance of ${75.37}_{-4.00}^{+4.73}$kpc and show that Draco has a modest LOS rotation, with $\left(v/\sigma \right)=0.22\pm 0.09$. Our results provide a new stringent test of the so-called “cusp–core” problem that can be readily extended to other dwarfs. 

   more » « less
3. Divergence beneath the Brillouin sphere and the phenomenology of prediction error in spherical harmonic series approximations of the gravitational field

   https://doi.org/10.1088/1361-6633/ad44d5  

   Bevis, M; Ogle, C; Costin, O; Jekeli, C; Costin, R D; Guo, J; Fowler, J; Dunne, G V; Shum, C K; Snow, K (June 2024, Reports on Progress in Physics)

   The Brillouin sphere is defined as the smallest sphere, centered at the origin of the geocentric coordinate system, that incorporates all the condensed matter composing the planet. The Brillouin sphere touches the Earth at a single point, and the radial line that begins at the origin and passes through that point is called the singular radial line. For about 60 years there has been a persistent anxiety about whether or not a spherical harmonic (SH) expansion of the external gravitational potential,V, will converge beneath the Brillouin sphere. Recently, it was proven that the probability of such convergence is zero. One of these proofs provided an asymptotic relation, called Costin’s formula, for the upper bound,E_N, on the absolute value of the prediction error,e_N, of a SH series model, $V_N(\theta ,\lambda ,r)$, truncated at some maximum degree, $N=n_{max}$. When the SH series is restricted to (or projected onto) a particular radial line, it reduces to a Taylor series (TS) in $1/r$. Costin’s formula is $E_N\simeq B{N}^{-b}(R/r{)}^N$, whereRis the radius of the Brillouin sphere. This formula depends on two positive parameters:b, which controls the decay of error amplitude as a function ofNwhenris fixed, and a scale factorB. We show here that Costin’s formula derives from a similar asymptotic relation for the upper bound,A_non the absolute value of the TS coefficients,a_n, for the same radial line. This formula, $A_n\simeq K{n}^{-k}$, depends on degree,n, and two positive parameters,kandK, that are analogous tobandB. We use synthetic planets, for which we can compute the potential,V, and also the radial component of gravitational acceleration, $g_r=\partial V/\partial r$, to hundreds of significant digits, to validate both of these asymptotic formulas. Let superscriptVrefer to asymptotic parameters associated with the coefficients and prediction errors for gravitational potential, and superscriptgto the coefficients and predictions errors associated withg_r. For polyhedral planets of uniform density we show that $b^V=k^V=7/2$and $b^g=k^g=5/2$almost everywhere. We show that the frequency of oscillation (around zero) of the TS coefficients and the series prediction errors, for a given radial line, is controlled by the geocentric angle,α, between that radial line and the singular radial line. We also derive useful identities connecting $K^V,B^V,K^g$, andB^g. These identities are expressed in terms of quotients of the various scale factors. The only other quantities involved in these identities areαandR. The phenomenology of ‘series divergence’ and prediction error (whenr < R) can be described as a function of the truncation degree,N, or the depth,d, beneath the Brillouin sphere. For a fixed $r\text{⩽}R$, asNincreases from very low values, the upper error boundE_Nshrinks until it reaches its minimum (best) value whenNreaches some particular or optimum value, $N_{\mathrm{opt}}$. When $N>N_{\mathrm{opt}}$, prediction error grows asNcontinues to increase. Eventually, when $N\gg N_{\mathrm{opt}}$, prediction errors increase exponentially with risingN. If we fix the value ofNand allow $R/r$to vary, then we find that prediction error in free space beneath the Brillouin sphere increases exponentially with depth,d, beneath the Brillouin sphere. Because $b^g=b^V-1$everywhere, divergence driven prediction error intensifies more rapidly forg_rthan forV, both in terms of its dependence onNandd. If we fix bothNandd, and focus on the ‘lateral’ variations in prediction error, we observe that divergence and prediction error tend to increase (as doesB) as we approach high-amplitude topography. 

   more » « less
4. TOI-5344 b: A Saturn-like Planet Orbiting a Super-solar Metallicity M0 Dwarf

   https://doi.org/10.3847/1538-3881/ad09c2  

   Han, Te; Robertson, Paul; Kanodia, Shubham; Cañas, Caleb; Lin, Andrea S. J.; Stefánsson, Gumundur; Libby-Roberts, Jessica E.; Larsen, Alexander; Kobulnicky, Henry A.; Mahadevan, Suvrath; et al (December 2023, The Astronomical Journal)

   Abstract We confirm the planetary nature of TOI-5344 b as a transiting giant exoplanet around an M0-dwarf star. TOI-5344 b was discovered with the Transiting Exoplanet Survey Satellite photometry and confirmed with ground-based photometry (the Red Buttes Observatory 0.6 m telescope), radial velocity (the Habitable-zone Planet Finder), and speckle imaging (the NN-Explore Exoplanet Stellar Speckle Imager). TOI-5344 b is a Saturn-like giant planet (ρ= 0.80 ${}_{-0.15}^{+0.17}$g cm⁻³) with a planetary radius of 9.7 ± 0.5R_⊕(0.87 ± 0.04R_Jup) and a planetary mass of ${135}_{-18}^{+17}{M}_{\oplus }$(0.42 ${}_{-0.06}^{+0.05}{M}_{\mathrm{Jup}}$). It has an orbital period of ${3.792622}_{-0.000010}^{+0.000010}$days and an orbital eccentricity of ${0.06}_{-0.04}^{+0.07}$. We measure a high metallicity for TOI-5344 of [Fe/H] = 0.48 ± 0.12, where the high metallicity is consistent with expectations from formation through core accretion. We compare the metallicity of the M-dwarf hosts of giant exoplanets to that of M-dwarf hosts of nongiants (≲8R_⊕). While the two populations appear to show different metallicity distributions, quantitative tests are prohibited by various sample caveats. 

   more » « less
5. Magnetogram-matching Biot–Savart Law and Decomposition of Vector Magnetograms

   https://doi.org/10.3847/1538-4357/add895  

   Titov, Viacheslav S; Downs, Cooper; Török, Tibor; Linker, Jon A; Prazak, Michael; Qiu, Jiong A (July 2025, The Astrophysical Journal)

   Abstract We generalize a magnetogram-matching Biot–Savart law (BSl) from planar to spherical geometry. For a given coronal current densityJ, this law determines the magnetic field $\tilde{\mathit{B}}$whose radial component vanishes at the surface. The superposition of $\tilde{\mathit{B}}$with a potential field defined by a given surface radial field,B_r, provides the entire configuration whereB_rremains unchanged by the currents. Using this approach, we (1) upgrade our regularized BSls for constructing coronal magnetic flux ropes (MFRs) and (2) propose a new method for decomposing a measured photospheric magnetic field as $\mathit{B}={\mathit{B}}_{\mathrm{pot}}+{\mathit{B}}_T+{\mathit{B}}_{\tilde{S}}$, where the potential,B_pot, toroidal,B_T, and poloidal, ${\mathit{B}}_{\tilde{S}}$, fields are determined byB_r,J_r, and the surface divergence ofB–B_pot, respectively, all derived from magnetic data. OurB_Tis identical to the one in the alternative Gaussian decomposition by P. W. Schuck et al., whileB_potand ${\mathit{B}}_{\tilde{S}}$are different from their poloidal fields ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$, which arepotentialin the infinitesimal proximity to the upper and lower side of the surface, respectively. In contrast, our ${\mathit{B}}_{\tilde{S}}$has no such constraints and, asB_potandB_T, refers to thesameupper side of the surface. In spite of these differences, for a continuousJdistribution across the surface,B_potand ${\mathit{B}}_{\tilde{S}}$are linear combinations of ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$. We demonstrate that, similar to the Gaussian method, our decomposition allows one to identify the footprints and projected surface-location of MFRs in the solar corona, as well as the direction and connectivity of their currents. 

   more » « less