Engineering topics

Stovall, K.

Publications and source records attributed to Stovall, K..

At least 19 records

NASA NTRS2019

The NANOGrav 11 Yr Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Observations indicate that nearly all galaxies contain supermassive black holes at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays. We have searched the North American Nanohertz Observatory for Gravitational Waves 11 yr data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95% upper limits on the strength of GWs from such sources. At f(gw) = 8 nHz, we placed a sky-averaged upper limit of h(0) < 7.3(3) × 10(exp −15). We also developed a technique to determine the significance of a particular signal in each pulsar using "dropout" parameters as a way of identifying spurious signals. From these upper limits, we ruled out SMBHBs emitting GWs f(gw) = 8 nHz within 120 Mpc for M = 10(exp 9) Solar Mass, and within 5.5 Gpc for M= 10(exp 10) Solar Mass at our most sensitive sky location. We also determined that there are no SMBHBs with M > 1.6 x 10(exp 9) Solar Mass emitting GWs with f(gw) = 2.8–317.8 nHz in the Virgo Cluster. Finally, we compared our strain upper limits to simulated populations of SMBHBs, based on galaxies in the Two Micron All-Sky Survey and merger rates from the Illustris cosmological simulation project, and found that only 34 out of 75,000 realizations of the local universe contained a detectable source.

Aggarwal, K.↗

NASA NTRS2019

The NANOGrav 12.5 yr Data Set: The Frequency Dependence of Pulse Jitter in Precision Millisecond Pulsars

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the moststable millisecond pulsars. Several known sources of noise on short timescales in single radio pulsar observationsare well described by a simple model of three components: template fitting from a finite signal-to-noise ratio, pulsephase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellarscattering variations. Currently template-fitting errors dominate, but as radio telescopes push toward higher signalto-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics ofjitter becomes crucial for properly characterizing arrival time uncertainties. We characterize the radio frequencydependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for GravitationalWaves timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for itsfrequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations ofrms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter istypically 1% of pulse phase. The average frequency dependence for all pulsars using a power-law model hasindex 0.42. We investigate the jitter variations for the interpulse of PSR B1937+21 and find no significantdeviations from the main pulse rms jitter. We also test the time variation of jitter in two pulsars and find thatsystematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly modeljitter as a significant component of the noise within the detector.

Lam, M. T.↗

NASA NTRS2019

The NANOGrav 11 yr Data Set: Solar Wind Sounding through Pulsar Timing

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens ofmillisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM)measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines ofsight cross through the solar wind (SW), which produces correlated DM fluctuations in all pulsars. We develop andapply techniques for extracting the imprint of the SW from the full collection of DM measurements in the recentlyreleased NANOGrav 11 yr data set. We filter out long-timescale DM fluctuations attributable to structure in theinterstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate thecorrelated signature of the wind from signals unique to individual lines of sight. When treating the SW asspherically symmetric and constant in time, we find the electron number density at 1au to be 7.9 ± 0.2 cm(exp -3). Wefind our data to be insensitive to long-term variation in the density of the wind. We argue that our techniques pairedwith a high-cadence, low-radio-frequency observing campaign of near-ecliptic pulsars would be capable ofmapping out large-scale latitudinal structure in the wind.

Madison, D.R.↗

NASA NTRS2018

The NANOGrav 11 Year Data Set: Pulsar-Timing Constraints on the Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints are sensitive to the solar system ephemeris (SSE) model used and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A (sub GWB) < 1.45 × 10 (exp -15) at a frequency of f=1 yr(exp -1) for a fiducial f (exp -2/3) power-law spectrum and with interpulsar correlations modeled. This is a factor of approximately 2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. We use our constraints to characterize the combined influence on the GWB of the stellar mass density in galactic cores, the eccentricity of SMBH binaries, and SMBH-galactic-bulge scaling relationships. We constrain the cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of G (sub mu) < 5.3 × 10(exp -11) - a factor of approximately 2 better than the published NANOGrav nine-year constraints. Our SSE-marginalized 95% upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation universe) is omega (sub GWB)(f) h (exp 2) < 3.4 × 10 (exp -10).

Supermassive Black Holes↗

DOE OSTI2017

A Multiwavelength Study of Nearby Millisecond Pulsar PSR J1400−1431: Improved Astrometry and an Optical Detection of Its Cool White Dwarf Companion

In 2012, five high-school students involved in the Pulsar Search Collaboratory discovered the millisecond pulsar (MSP) PSR J1400−1431, and initial timing parameters were published in Rosen et al. a year later. Since then, we have obtained a phase-connected timing solution spanning five years, resolving a significant position discrepancy and measuring P-dot , proper motion, parallax, and a monotonic slope in dispersion measure over time. Due to PSR J1400−1431’s proximity and significant proper motion, we use the Shklovskii effect and other priors to determine a 95% confidence interval for PSR J1400−1431’s distance, d=270{sub −80}{sup +130} pc. With an improved timing position, we present the first detection of the pulsar’s low-mass white dwarf (WD) companion using the Goodman Spectrograph on the 4.1 m SOAR telescope. Deeper imaging suggests that it is a cool DA-type WD with T{sub eff}=3000±100 K and R/R{sub ⊙}=(2.19±0.03)×10{sup −2} (d/270 pc). We show a convincing association between PSR J1400−1431 and a γ-ray point source, 3FGL J1400.5−1437, but only weak (3.3σ) evidence of pulsations after folding γ-ray photons using our radio timing model. We detect an X-ray counterpart with XMM-Newton, but the measured X-ray luminosity (1×10{sup 29} erg s{sup −1}) makes PSR J1400−1431 the least X-ray luminous rotation-powered MSP detected to date. Together, our findings present a consistent picture of a nearby (d≈230 pc) MSP in a 9.5-day orbit around a cool ∼0.3 M {sub ⊙} WD companion, with orbital inclination i≳60{sup ∘}.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗

DOE OSTI2016

Timing of five PALFA-discovered millisecond pulsars

We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the five pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 years. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ∼0.1 M{sub ⊙} companion and possible eclipses that last ∼10% of the orbit. The position of PSR J1957+2516 is also coincident with a near-infrared source. All five MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of γ-ray pulsations in a fold of Fermi Gamma-Ray Space Telescope data. These five MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

OBSERVATIONS OF ROTATING RADIO TRANSIENTS WITH THE FIRST STATION OF THE LONG WAVELENGTH ARRAY

Rotating radio transients (RRATs) are a subclass of pulsars first identified in 2006 that are detected only in searches for single pulses and not through their time averaged emission. Here, we present the results of observations of 19 RRATs using the first station of the Long Wavelength Array (LWA1) at frequencies between 30 and 88 MHz. The RRATs observed here were first detected in higher frequency pulsar surveys. Of the 19 RRATs observed, two sources were detected and their dispersion measures, periods, pulse profiles, and flux densities are reported and compared to previous higher frequency measurements. We find a low detection rate (11%), which could be a combination of the lower sensitivity of LWA1 compared to higher frequency telescopes, and the result of scattering by the interstellar medium or a spectral turnover.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

MEASURING THE MAGNETIC FIELD OF CORONAL MASS EJECTIONS NEAR THE SUN USING PULSARS

The utility of Faraday rotation to measure the magnetic field of the solar corona and large-scale transients within is a small, yet growing field in solar physics. This is largely because it has been recognized as a potentially valuable frontier in space weather studies, because the ability to measure the intrinsic magnetic field within coronal mass ejections (CMEs) when they are close to the Sun is of great interest for understanding a key element of space weather. Such measurements have been attempted over the last few decades using radio signals from artificial sources (i.e., spacecraft on the far side of the Sun), but studies involving natural radio sources are scarce in the literature. We report on a preliminary study involving an attempt to detect the Faraday rotation of a CME that passed in front of a pulsar (PSR B0950+08) in 2015 August. We combine radio measurements with those from a broadband visible light coronagraph, to estimate the upper limit of the magnetic field of the CME when it was in the corona. We find agreement between different approaches for obtaining its density, and values that are consistent with those predicted from prior studies of CME density close to the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

A deep campaign to characterize the synchronous radio/X-ray mode switching of PSR B0943+10

Observations of PSR B0943+10 with XMM-Newton and the LOFAR, LWA, and Arecibo radio telescopes in 2014 November confirm the synchronous X-ray/radio switching between a radio-bright (B) mode and a radio-quiet (Q) mode, in which the X-ray flux is a factor ∼2.4 higher than in the B-mode. We discovered X-ray pulsations during the B-mode (0.5–2 keV pulsed fraction of (38 ± 5)%) and confirm their presence in the Q-mode, where the pulsed fraction increases with energy from ∼20% to ∼65% at 2 keV. We found marginal evidence for an increase in the X-ray pulsed fraction during the B-mode on a timescale of hours. The X-ray spectrum during the Q-mode requires a fit with either a power law plus blackbody or the sum of two blackbodies, while in the B-mode it is well fit by a single blackbody (a single power law is rejected). In the Q-mode, the pulsed emission has a blackbody spectrum with temperature ∼3.4×10{sup 6} K and the unpulsed emission is a power law with photon index ∼2.5, while during the B-mode both the pulsed and unpulsed emission can be fit by either a blackbody or a power law with similar values of temperature and index. A Chandra image does not show diffuse X-ray emission. These results support a scenario in which unpulsed non-thermal emission, likely magnetospheric, and pulsed thermal emission from a small polar cap (∼1500 m{sup 2}) with a non-dipolar field (∼10{sup 14} G) are present during both modes and vary in a correlated way. This is broadly consistent with the partially screened gap model and does not necessarily imply global magnetospheric rearrangements to explain the mode switching.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

A LOW FREQUENCY SURVEY OF GIANT PULSES FROM THE CRAB PULSAR

We present a large survey of giant pulses from the Crab Pulsar as observed with the first station of the Long Wavelength Array. Automated methods for detecting giant pulses at low frequencies where scattering becomes prevalent are also explored. More than 1400 pulses were detected across four frequency bands between 20 and 84 MHz over a seven-month period beginning in 2013, with additional followup observations in late 2014 and early 2015. A handful of these pulses were detected simultaneously across all four frequency bands. We examine pulse characteristics, including pulse broadening and power law indices for amplitude distributions. We find that the flux density increases toward shorter wavelengths, consistent with a spectral turnover at 100 MHz. Our observations uniquely span multiple scattering epochs, manifesting as a notable trend in the number of detections per observation. These results are characteristic of the variable interface between the synchrotron nebula and the surrounding interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

NEW DISCOVERIES FROM THE ARECIBO 327 MHz DRIFT PULSAR SURVEY RADIO TRANSIENT SEARCH

We present Clusterrank, a new algorithm for identifying dispersed astrophysical pulses. Such pulses are commonly detected from Galactic pulsars and rotating radio transients (RRATs), which are neutron stars with sporadic radio emission. More recently, isolated, highly dispersed pulses dubbed fast radio bursts (FRBs) have been identified as the potential signature of an extragalactic cataclysmic radio source distinct from pulsars and RRATs. Clusterrank helped us discover 14 pulsars and 8 RRATs in data from the Arecibo 327 MHz Drift Pulsar Survey (AO327). The new RRATs have DMs in the range 23.5–86.6 pc cm{sup −3} and periods in the range 0.172–3.901 s. The new pulsars have DMs in the range 23.6–133.3 pc cm{sup −3} and periods in the range 1.249–5.012 s, and include two nullers and a mode-switching object. We estimate an upper limit on the all-sky FRB rate of 10{sup 5} day{sup −1} for bursts with a width of 10 ms and flux density ≳83 mJy. The DMs of all new discoveries are consistent with a Galactic origin. In comparing statistics of the new RRATs with sources from the RRATalog, we find that both sets are drawn from the same period distribution. In contrast, we find that the period distribution of the new pulsars is different from the period distributions of canonical pulsars in the ATNF catalog or pulsars found in AO327 data by a periodicity search. This indicates that Clusterrank is a powerful complement to periodicity searches and uncovers a subset of the pulsar population that has so far been underrepresented in survey results and therefore in Galactic pulsar population models.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

The NANOGrav Nine-Year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr(exp -1) for a GWB from supermassive black hole binaries of A(sub gw) less than 1.5 x 10(exp -15). We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, OMEGA(sub gw) (f) h squared less than 4.2 x 10(exp -10). Our limit on the cosmic string GWB, OMEGA(sub gw) (f) h squared less than 2.2 x 10(exp -10), translates to a conservative limit on the cosmic string tension with G mu less than 3.3 x 10(exp -8), a factor of four better than the joint Planck and high-l‚ cosmic microwave background data from other experiments.

Arzoumanian, Z.↗

DOE OSTI2015

Pulsar J0453+1559: A double neutron star system with a large mass asymmetry

To understand the nature of supernovae and neutron star (NS) formation, as well as binary stellar evolution and their interactions, it is important to probe the distribution of NS masses. Until now, all double NS (DNS) systems have been measured as having a mass ratio close to unity (q ≥ 0.91). Here, we report the measurement of the individual masses of the 4.07-day binary pulsar J0453+1559 from measurements of the rate of advance of periastron and Shapiro delay: the mass of the pulsar is M{sub p} = 1.559 ± 0.005 M{sub ⊙} and that of its companion is M{sub c}=1.174±0.004 M{sub ⊙}; q = 0.75. If this companion is also an NS, as indicated by the orbital eccentricity of the system (e = 0.11), then its mass is the smallest precisely measured for any such object. The pulsar has a spin period of 45.7 ms and a spin period derivative of P-dot = (1.8616±0.0007)×10{sup −19} s s{sup −1}; from these, we derive a characteristic age of ∼ 4.1×10{sup 9} years and a magnetic field of ∼ 2.9×10{sup 9} G, i.e., this pulsar was mildly recycled by the accretion of matter from the progenitor of the companion star. This suggests that it was formed with (very approximately) its current mass. Thus, NSs form with a wide range of masses, which is important for understanding their formation in supernovae. It is also important for the search for gravitational waves released during an NS–NS merger: it is now evident that we should not assume that all DNS systems are symmetric.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Discovery and follow-up of rotating radio transients with the Green Bank and LOFAR telescopes

We have discovered 21 Rotating Radio Transients (RRATs) in data from the Green Bank Telescope (GBT) 350 MHz Drift-scan and the Green Bank North Celestial Cap pulsar surveys using a new candidate sifting algorithm. RRATs are pulsars with sporadic emission that are detected through their bright single pulses rather than Fourier domain searches. We have developed RRATtrap, a single-pulse sifting algorithm that can be integrated into pulsar survey data analysis pipelines in order to find RRATs and Fast Radio Bursts. We have conducted follow-up observations of our newly discovered sources at several radio frequencies using the GBT and Low Frequency Array, yielding improved positions and measurements of their periods, dispersion measures (DMs), and burst rates, as well as phase-coherent timing solutions for four of them. The new RRATs have DMs ranging from 15 to 97 pc cm{sup −3}, periods of 240 ms to 3.4 s, and estimated burst rates of 20 to 400 pulses hr{sup −1} at 350 MHz. We use this new sample of RRATs to perform statistical comparisons between RRATs and canonical pulsars in order to shed light on the relationship between the two populations. We find that the DM and spatial distributions of the RRATs agree with those of the pulsars found in the same survey. We find evidence that slower pulsars (i.e., P>200 ms) are preferentially more likely to emit bright single pulses than are faster pulsars (P<200 ms), although this conclusion is tentative. Our results are consistent with the proposed link between RRATs, transient pulsars, and canonical pulsars as sources in various parts of the pulse activity spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

PULSAR OBSERVATIONS USING THE FIRST STATION OF THE LONG WAVELENGTH ARRAY AND THE LWA PULSAR DATA ARCHIVE

We present initial pulsar results from the first station of the Long Wavelength Array (LWA1) obtained during the commissioning period of LWA1 and in early science results. We present detections of periodic emission from 44 previously known pulsars, including 3 millisecond pulsars. The effects of the interstellar medium (ISM) on pulsar emission are significantly enhanced at the low frequencies of the LWA1 band (10–88 MHz), making LWA1 a very sensitive instrument for characterizing changes in the dispersion measure (DM) and other effects from the ISM. Pulsars also often have significant evolution in their pulse profile at low frequency and a break in their spectral index. We report DM measurements for 44 pulsars, mean flux density measurements for 36 pulsars, and multi-frequency component spacing and widths for 15 pulsars with more than one profile component. For 27 pulsars, we report spectral index measurements within our frequency range. We also introduce the LWA1 Pulsar Data Archive, which stores reduced data products from LWA1 pulsar observations. Reduced data products for the observations presented here can be found in the archive. Reduced data products from future LWA1 pulsar observations will also be made available through the archive.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The green bank northern celestial cap pulsar survey. I. Survey description, data analysis, and initial results

We describe an ongoing search for pulsars and dispersed pulses of radio emission, such as those from rotating radio transients (RRATs) and fast radio bursts, at 350 MHz using the Green Bank Telescope. With the Green Bank Ultimate Pulsar Processing Instrument, we record 100 MHz of bandwidth divided into 4096 channels every 81.92 μs. This survey will cover the entire sky visible to the Green Bank Telescope (δ > –40°, or 82% of the sky) and outside of the Galactic Plane will be sensitive enough to detect slow pulsars and low dispersion measure (<30 pc cm{sup –3}) millisecond pulsars (MSPs) with a 0.08 duty cycle down to 1.1 mJy. For pulsars with a spectral index of –1.6, we will be 2.5 times more sensitive than previous and ongoing surveys over much of our survey region. Here we describe the survey, the data analysis pipeline, initial discovery parameters for 62 pulsars, and timing solutions for 5 new pulsars. PSR J0214+5222 is an MSP in a long-period (512 days) orbit and has an optical counterpart identified in archival data. PSR J0636+5129 is an MSP in a very short-period (96 minutes) orbit with a very low mass companion (8 M{sub J}). PSR J0645+5158 is an isolated MSP with a timing residual RMS of 500 ns and has been added to pulsar timing array experiments. PSR J1434+7257 is an isolated, intermediate-period pulsar that has been partially recycled. PSR J1816+4510 is an eclipsing MSP in a short-period orbit (8.7 hr) and may have recently completed its spin-up phase.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

DETECTION OF RADIO EMISSION FROM FIREBALLS

We present the findings from the Prototype All-Sky Imager, a back end correlator of the first station of the Long Wavelength Array, which has recorded over 11,000 hr of all-sky images at frequencies between 25 and 75 MHz. In a search of this data for radio transients, we have found 49 long-duration (10 s of seconds) transients. Ten of these transients correlate both spatially and temporally with large meteors (fireballs), and their signatures suggest that fireballs emit a previously undiscovered low frequency, non-thermal pulse. This emission provides a new probe into the physics of meteors and identifies a new form of naturally occurring radio transient foreground.

79 ASTRONOMY AND ASTROPHYSICS↗