An active galactic nucleus as a massive coaxial pulsar: A massive central object with a Keplerian accretion disk and twisted MHD jet

The similarities in certain empirical relations derived for active galactic nuclei (AGN) and radio pulsars support the idea that the energy release in the jets of AGN is electrodynamical, associated with the presence of a massive central object (MCO) and circumnuclear, magnetized Keplerian disk. For example, in both types of object, the radio luminosities are proportional to characteristic time scales—for precession for AGN or rotation for pulsars—to the 0.8 power. In this approach to studying the nature of the activity, the $\dot P - P$ relation for pulsars can be considered analogous to the P _r-M _MCO relation for AGN, which has the form P _5GHz (W/Hz)=2.73 log(M _MCO/M _⊙)?2.87.     

A comparative analysis of the integrated radio luminosities of normal and millisecond pulsars

We present the results of a comparative statistical analysis of the integrated radio luminosities of millisecond and normal pulsars and their dependences on other parameters of the pulsars. The analysis is based on our own measurements of the flux densities, spectra, and integrated radio luminosities of the millisecond pulsars, as well as data from the literature used to determine the integrated radio luminosities for 545 pulsars, 50 of them millisecond pulsars. Despite large differences in their periods P, period derivatives $\dot P$ , magnetic fields B, and ages τ, the integrated radio luminosities of the millisecond and normal pulsars and their dependences on other parameters are approximately the same. The integrated radio luminosity depends on the parameter B/P ², which is proportional to the potential difference in the polar-cap gap; this may indicate that the radio energy of pulsars is determined by the energy of primary particles accelerated in the polar-cap gap. Secular decreases in the radio luminosities of both normal and millisecond pulsars were also detected.     

Peculiarities of radio pulsars with emission outside the radio range

Parameters of 100 radio pulsars detected outside the radio range (he pulsars) are compared with those of pulsars radiating only in the radio (n pulsars). The periods of he pulsars are, on average, appreciably shorter than those of n pulsars: 〈P〉 = 0.10 and 0.56 s, respectively. The distribution of the magnetic field at the light cylinder is shifted toward higher magnetic fields for the pulsars with high-energy radiation, compared to the distribution for pulsars radiating only in the radio. The magnetic fields at the light cylinder are 〈B _lc〉 = 9×10³ G for he radio pulsars, and 〈B_lc〉 = 56 G formost purely radio pulsars. This suggests the generation of high-energy nonthermal radiation in radio pulsars at the peripheries of their magnetospheres. The distribution of the spin-energy loss rate dE/dt is uniform for he pulsars, and is characterized by a higher average value \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 35.53} \right) \) , compared to n pulsars, \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 32.60} \right) \) . The spatial distribution of he pulsars is nonuniform: they form two well separated clouds.     

The radio pulsar J0205+6449 in the supernova remnant 3C 58

The detection of pulsed radio emission from the recently discovered X-ray pulsar J0205+6449 in the young supernova remnant 3C 58 is reported together with the results of first studies of this emission. The observations were carried out at 111 and 88 MHz on radio telescopes of the Pushchino Observatory. The pulsar period, 65.68 ms, and period derivative, \(\dot P = 1.9 \times 10^{ - 13} \), have been confirmed. The integrated pulse profile at 111 MHz has been obtained and the flux density and spectral index α=2.8 measured. The pulsar dispersion measure DM=141 pc cm^?3 has been confirmed. This dispersion measure yields a distance to the pulsar of d=6.4 kpc, a factor of two or more greater than the previously favored distance to the supernova remnant 3C 58 (2.6 kpc). The problem of the age and distance of the pulsar-SNR system is discussed. If the age of the pulsar J0205+6449 is equal to that of the SNR (820 years), this pulsar is the youngest known radio pulsar. The synchrotron mechanism for the radio and X-ray emission is proposed to explain the lower radio and X-ray luminosity of this new pulsar compared to the Crab pulsar, which is similar to it in many ways. Optical emission with luminosity L_opt=10³¹ erg/s and gamma-ray emission with L_γ=7×10³⁵ erg/s are predicted, and the steep radio spectrum (α≈3) can be explained.     

Physical conditions in faint gigahertz-peaked spectrum radio sources

Weak, compact radio sources (～100 mJy peak flux, L～1–10 pc) with their spectral peaks at about a gigahertz are studied, based on the complete sample of 46 radio sources of Snellen, drawn from high-sensitivity surveys, including the low-frequency Westerbork catalog. The physical parameters have been estimated for 14 sources: the magnetic field (H _⊥), the number density of relativistic particles (n _e), the energy of the magnetic field $(E_{H_ \bot } )$ , and the energy of relativistic particles (E _e). Ten sources have $E_{H_ \bot } \ll E_e $ , three have approximate equipartition of the energies $(E_{H_ \bot } \sim E_e )$ , and only one has $E_{H_ \bot } \gg E_e $ . The mean magnetic fields in quasars (10^?3 G) and galaxies (10^?2 G) have been estimated. The magnetic field appears to be related to the sizes of compact features as $H \sim 1/\sqrt L $ .     

Magnetic fields and the $$\dot P - P$$ diagram for radio pulsars

Various mechanisms for the loss of angular momentum of neutron stars are analyzed. Theoretical predictions about the evolution of the period are compared with the observed distribution of pulsars on the log\(\dot P\)log(P) diagram. Pulsars with short periods (P≤0.1 s) cannot be fit well by any of the models considered. Their braking index is n=?1, which requires the development of a new braking mechanism. The evolution of pulsars with P>1.25 s is described by the law \(\dot P \propto P^2\), probably due to processes internal to the neutron stars. The observational data for pulsars with 0.1<P≤1.25 s can be fit with a hybrid model incorporating internal processes and magnetic-dipole losses. The magnetic fields in pulsar catalogs should be recomputed in accordance with the results obtained. For example, the magnetic fields obtained for two magnetars with P=5.16 s and P=7.47 s are B _s =1.7×10¹³ and 2.9×10¹³ G, which are lower than the critical field B_cr=4.4×10¹³ G. For a substantial fraction of pulsars, their characteristic ages \(\tau = P/2\dot P\) cannot serve as measures of their real ages.     

A spindown mechanism for short-period radio pulsars

It is shown that a model with accretion in a “quasi-propeller” mode can explain the observed spindown of pulsars with periods P<0.1 s. The mean accretion rate for 39 selected objects is \(\dot M = 5.6 \times 10^{ - 11} M_ \odot /year\). If \(\dot M\) is constant during the pulsar’s lifetime, the neutron star will stop rotating after 10⁷ years. The mean magnetic field at the neutron-star surface calculated in this model, \(\bar H_0 = 6.8 \times 10^8 G\), is consistent to an order of magnitude with the values of H₀ for millisecond pulsars from known catalogs. However, the actual value of H₀ for particular objects can differ from the catalog values by appreciable factors, and these quantities must be recalculated using more adequate models. The accretion disk around the neutron star should not impede the escape of the pulsar’s radiation, since this radiation is generated near the light cylinder in pulsars with P<0.1 s. Pulsars such as PSR 0531+21 and PSR 0833-45 have probably spun down due to the effect of magnetic-dipole radiation. If the difference in the braking indices for these objects from n=3 is due to the effect of accretion, the accretion rate must be of the order of 10¹⁸ g/s.     

Photometric observations of R Coronae Borealis in the optical and infrared

The results of long-term photometric observations of R CrB in the UBV JHKLM bands are presented. The temporal and color characteristics of the emission of the star itself and of its extended dust envelope are analyzed in detail. No stable harmonic has been found in the semiregular variations of the optical brightness of R CrB. Two harmonics with periods P≈3.3 and 11.3 yrs have been detected in the brightness variations of the dust envelope; the minima of these variations coincided in 1999, resulting in a record decrease in the LM brightness of the envelope. This by chance coincided in time with a deep minimum of the visual brightness of the star, resulting in a unique decrease in the total brightness of the star and dust envelope. This enabled estimation of the bolometric flux of the hot dust clouds, which made up only a few per cent of the bolometric flux of the dust envelope. The brightness variations of the dust envelope are not accompanied by appreciable color changes and are associated with variations of its optical depth τ(V) in the range 0.2–0.4. The dust envelope forms at a large and fairly constant distance from the star $(r_{in} \approx 110R_* ,T_{gr} \approx 860 K)$ , from material in its stellar wind, whose intensity $(\dot M_{gas} \approx 2.1 \times 10^{ - 7} M_ \odot /year)$ obeys a Reimers law. No variations synchronous with those of the optical depth of the dust envelope, in particular, with the period P≈3.3 yrs, have been found in the optical emission of R CrB, suggesting that the stellar wind is not spherically symmetric. The dust envelope consists of small grains (a _gr≤0.01 µm), while the clouds screening the star from the observer are made up of large grains (a _gr≈0.1 µm). The activity of R CrB, whose nature is unclear, is reflected in variations of the stellar-wind intensity and the appearance of dust clouds in the line of sight: these variations are repeated by corresponding changes in the optical depth of the dust envelope with a delay of ～4 years (the time for a particle moving at V _env≈45 km/s to move from the star to the boundary of the dust envelope).     

The pulsar J1852+0040 in Kes 79 and the nature of central compact objects in supernova remnants

A possible model for the pulsar PSR J1852+0040 associated with the supernova remnant Kes 79 and detected in place of a central compact object in this remnant is discussed. The main observational properties of the pulsar can be understood as consequences of its weak surface magnetic field (B _s < 3 × 10¹¹ G) and short rotational period (P ～ 0.1 s). Its X-ray emission is thermal, and is generated in a small region near the surface of the neutron star due to cooling of the surface as the surface accretes matter from a relict disk surrounding the pulsar. The radio emission is generated in the outer layers of the pulsar magnetosphere by the synchrotron (cyclotron) mechanism. The optical luminosity of J1852+0040 is estimated to be L _opt < 10²⁸ erg/s. If the spectral features in another central compact object, 1E 1207.4+5209, are interpreted as electron cyclotron lines, this provides evidence for a weak surface magnetic field for this neutron star as well (B < 6 × 10¹⁰ G). The hypothesis that all central compact objects have weak surface fields makes it possible to explain the number of detected central compact objects, the absence of pulsar-wind nebulae associated with these objects, and the fact that no pulsar has yet been detected at the position of SN 1987a. We suggest that, after the supernova remnant has dissipated, the central compact object becomes a weak X-ray source (XDINS), whose weak emission is also due to the weakness of its magnetic field.     

Microstructure of pulsar radio pulses measured with a time resolution of 62.5 ns at 1650 MHz

We present an analysis of pulsar observations carried out on two frequency channels at 1634 MHz and 1650 MHz with a time resolution of 62.5 ns on the 70-m radio telescope of the NASA Deep Space Network in Tidbinbilla. The data were recorded using the S2 system, intended primarily for VLBI observations. Microstructure with characteristic timescales of 270, 80, and 150 µs was detected in pulsars B0833-45, B1749-28, and B1933 + 16, respectively. The distribution of microstructure timescales for the Vela pulsar (B0833-45) is characterized by a gradual growth with decreasing timescale to 200 µs; the distribution has a maximum at 20–200 µs and falls off sharply for timescales below 20 µs. The statistical relation between the microstructure modulation index m and the corresponding timescale τ_µ can be approximated by the power law dependence R∝τ _⊙ ^0.5 ; i.e., the intensity is higher for micropulses with longer durations. This contradicts the predictions of nonlinear models for the formation of micropulses by supercompact soliton wave packets. In all the pulsars studied, the time delays of the micropulses between the two frequency channels deviate from the expected dispersion laws for the interstellar plasma. In particular, the micropulses in the low-frequency channel arrive earlier than predicted by the dispersion measures derived previously from the mean pulse profiles. The deviation from the dispersion delay is determined most accurately for B0833-45, and is 4.9±0.2 µs. Such anomalous delays are probably associated with the effects of propagation of the radio emission within the pulsar magnetosphere.     

Differences in the parameters of radio pulsars with short and long periods

A comparative analysis of various parameters of pulsars with short (P < 0.1 s) and long (P > 0.1 s) periods is carried out. There is no correlation between the radio and gamma-ray luminosities of the pulsars and their surfacemagnetic fields, but there is a correlation between the X-ray luminosity and the surfacemagnetic field. A dependence of the X-ray and gamma-ray luminosities on the magnetic field at the light cylinder is also found. This result provides evidence for the formation of hard, non-thermal emission at the periphery of the magnetosphere. An appreciable positive correlation between the luminosity and the rate of rotational energy loss by the neutron star is observed, supporting the idea that all radio pulsars have the same basic source of energy. The efficiency of the transformation of rotational energy into radiation is significantly higher in long-period pulsars. The dependence of the pulse width on the pulsar period is steeper for pulsars with short periods than for those with long periods. The results obtained support earlier assertions that there are differences in the processes generating the emission in pulsars with P < 0.1 s and those with P > 0.1 s.     

High-temperature 29Si MAS NMR spectroscopy of anorthite (CaAl2Si2O8) and its P\bar 1-I\bar 1structural phase transition

We present ²⁹Si MAS NMR data for a well-ordered natural anorthite, obtained in situ at temperatures of from 25 to 500° C, which follow the changes in the aluminosilicate framework through the P $\bar 1$ -I $\bar 1$ structural phase transition. Pairs of peaks due to sites offset by approximately 1/2 [111] converge through the P $\bar 1$ phase and only four peaks are present above about 241° C. The variation of the peak positions with temperature and correlations based on structural data for the P $\bar 1$ and I $\bar 1$ phases allow assignment of all the MAS-NMR peaks to crystallographic sites. A Landau-type analysis gives an expression that relates the separation of pairs of con verging peaks to the local order parameter for the P $\bar 1$ -I $\bar 1$ transition, from which we determine its temperature dependence. Data for the best-constrained set of peak positions give for the order parameter critical exponent β = 0.27±0.04, consistent with previous results indicating that the P $\bar 1$ -I $\bar 1$ transition in pure anorthite is tricritical. No significant change in the ²⁹Si spin-lattice relaxation rate occurs across the P $\bar 1$ -I $\bar 1$ transition.     

The P–V–T equation of state of CaPtO3 post-perovskite

Orthorhombic post-perovskite CaPtO₃ is isostructural with post-perovskite MgSiO₃, a deep-Earth phase stable only above 100 GPa. Energy-dispersive X-ray diffraction data (to 9.4 GPa and 1,024 K) for CaPtO₃ have been combined with published isothermal and isobaric measurements to determine its P–V–T equation of state (EoS). A third-order Birch–Murnaghan EoS was used, with the volumetric thermal expansion coefficient (at atmospheric pressure) represented by α(T) = α₀ + α₁(T). The fitted parameters had values: isothermal incompressibility, $ K_{{T_{0} }} $  = 168.4(3) GPa; $ K_{{T_{0} }}^{\prime } $  = 4.48(3) (both at 298 K); $ \partial K_{{T_{0} }} /\partial T $  = ?0.032(3) GPa K^?1; α₀ = 2.32(2) × 10^?5 K^?1; α₁ = 5.7(4) × 10^?9 K^?2. The volumetric isothermal Anderson–Grüneisen parameter, δ _T , is 7.6(7) at 298 K. $ \partial K_{{T_{0} }} /\partial T $ for CaPtO₃ is similar to that recently reported for CaIrO₃, differing significantly from values found at high pressure for MgSiO₃ post-perovskite (?0.0085(11) to ?0.024 GPa K^?1). We also report axial P–V–T EoS of similar form, the first for any post-perovskite. Fitted to the cubes of the axes, these gave $ \partial K_{{aT_{0} }} /\partial T $  = ?0.038(4) GPa K^?1; $ \partial K_{{bT_{0} }} /\partial T $  = ?0.021(2) GPa K^?1; $ \partial K_{{cT_{0} }} /\partial T $  = ?0.026(5) GPa K^?1, with δ _T  = 8.9(9), 7.4(7) and 4.6(9) for a, b and c, respectively. Although $ K_{{T_{0} }} $ is lowest for the b-axis, its incompressibility is the least temperature dependent.     

Fitting matrix-valued variogram models by simultaneous diagonalization (Part I: Theory)

Suppose that ¯(x₁),...,¯Z(x_n). are observations of vector-valued random function ¯(x). In the isotropic situation, the sample variogram γ^*(h) for a given lag h is $$\bar \gamma ^ * (h) = \frac{1}{{2N(h)}}\mathop \sum \limits_{s(h)} (\overline Z (x_1 ) - \overline Z (x_1 )) \overline {(Z} (x_1 ) - \overline Z (x_1 ))^T $$ where s(h) is a set of paired points with distance h and N(h) is the number of pairs in s(h).. For a selection of lags h₁, h₂, .... h_k such that N (h₁) > O. we obtain a ktuple of (semi) positive definite matrices $\bar \gamma ^ * (h_{ 1} ),. . . ., \bar \gamma ^ * (h_{ k} )$ . We want to determine an orthonormal matrix B which simultaneously diagonalizes the $\bar \gamma ^ * (h_{ 1} ),. . . ., \bar \gamma ^ * (h_{ k} )$ or nearly diagonalizes them in the sense that the sum of squares of offdiagonal elements is small compared to the sum of squares of diagonal elements. If such a B exists, we linearly transform $\overline Z (x)$ by $\overline Y (x) = B\overline Z (x)$ . Then, the resulting vector function $\overline Y (x)$ has less spatial correlation among its components than $\overline Z (x)$ does. The components of $\overline Y (x)$ with little contribution to the variogram structure may be dropped, and small crossvariograms fitted by straightlines. Variogram models obtained by this scheme preserve the negative definiteness property of variograms (in the matrix-valued function sense). A simplified analysis and computation in cokriging can be carried out. The principles of this scheme arc presented in this paper.     

The dilatant behaviour of sand–pile interface subjected to loading and stress relief

Property and behaviour of sand–pile interface are crucial to shaft resistance of piles. Dilation or contraction of the interface soil induces change in normal stress, which in turn influences the shear stress mobilised at the interface. Although previous studies have demonstrated this mechanism by laboratory tests and numerical simulations, the interface responses are not analysed systematically in terms of soil state (i.e. density and stress level). The objective of this study is to understand and quantify any increase in normal stress of different pile–soil interfaces when they are subjected to loading and stress relief. Distinct element modelling was carried out. Input parameters and modelling procedure were verified by experimental data from laboratory element tests. Parametric simulations of shearbox tests were conducted under the constant normal stiffness, constant normal load and constant volume boundary conditions. Key parameters including initial normal stress ( $ \sigma_{{{\text{n}}0}}^{\prime } $ ), initial void ratio (e ₀), normal stiffness constraining the interface and loading–unloading stress history were investigated. It is shown that mobilised stress ratio ( $ \tau /\sigma_{\text{n}}^{\prime } $ ) and normal stress increment ( $ \Updelta \sigma_{\text{n}}^{\prime } $ ) on a given interface are governed by $ \sigma_{{{\text{n}}0}}^{\prime } $ and e ₀. An increase in $ \sigma_{{{\text{n}}0}}^{\prime } $ from 100 to 400 kPa leads to a 30 % reduction in $ \Updelta \sigma_{\text{n}}^{\prime } $ . An increase in e ₀ from 0.18 to 0.30 reduces $ \Updelta \sigma_{\text{n}}^{\prime } $ by more than 90 %, and therefore, shaft resistance is much lower for piles in loose sands. A unique relationship between $ \Updelta \sigma_{\text{n}}^{\prime } $ and normal stiffness is established for different soil states. It can be applied to assess the shaft resistance of piles in soils with different densities and subjected to loading and stress relief. Fairly good agreement is obtained between the calculated shaft resistance based on the proposed relationship and the measured results in centrifuge model tests.     

Biradical states of oxygen-vacancy defects in ��-quartz: centers E_{2}^{\prime \prime } and E_{4}^{\prime \prime }

Several new radiation defects with total electron spin S?=?1 occurring in electron-irradiated, synthetic ??-quartz have been observed by using electron paramagnetic resonance spectroscopy. These defects are considered to be biradicals, i.e., pairs of S?=?1/2 species. The concentration of these centers depends on the condition of the fast-electron irradiation. They have different decay behaviors that allow measurements of any individual species especially when it predominates over the others. The primary spin Hamiltonian parameter matrices g ₁, g ₂, D have now been determined for two similar defects, which herein are labeled $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ . Inter-electron distances estimated by using the magnetic dipole model, suggest that the structures of centers $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ both involve the unpaired electrons each located in orbitals of two silicon atoms next to a common oxygen vacancy but which have slightly different Si?CSi distances at 0.90 and 0.79?nm, respectively. This model is consistent with previous DFT calculations of the triplet configurations with local energetic minima. Observed decay behaviors suggest a transformation of centers $ E_{2,4}^{\prime \prime } $ to the analogous $ E_{1}^{\prime \prime } $ center. These triplet centers in quartz provide new insights into the structures of analogous defects in amorphous silica.     

Vigrishinite, Zn2Ti4 − x Si4O14(OH,H2O,□)8, a new mineral from the Lovozero alkaline complex, Kola Peninsula, Russia

A new mineral vigrishinite, epistolite-group member and first layer titanosilicate with species-defining Zn, was found at Mt. Malyi Punkaruaiv, in the Lovozero alkaline complex, Kola Peninsula, Russia. It occurs in a hydrothermally altered peralkaline pegmatite and is associated with microcline, ussingite, aegirine, analcime, gmelinite-Na, and chabazite-Ca. Vigrishinite forms rectangular or irregularly shaped lamellae up to 0.05 × 2 × 3 cm flattened on [001]. They are typically slightly split and show blocky character. The mineral is translucent to transparent and pale pink, yellowish-pinkish or colorless. The luster is vitreous. The Mohs’ hardness is 2.5–3. Vigrishinite is brittle. Cleavage is {001} perfect. D _meas = 3.03(2), D _calc = 2.97 g/cm³. The mineral is optically biaxial (?), α = 1.755(5), β = 1.82(1), γ = 1.835(8), 2V _meas = 45(10)°, 2V _calc = 50°. IR spectrum is given. The chemical composition (wt %; average of 9 point analyses, H₂O is determined by modified Penfield method) is as follows: 0.98 Na₂O, 0.30 K₂O, 0.56 CaO, 0.05 SrO, 0.44 BaO, 0.36 MgO, 2.09 MnO, 14.39 ZnO, 2.00 Fe₂O₃, 0.36 Al₂O₃, 32.29 SiO₂, 29.14 TiO₂, 2.08 ZrO₂, 7.34 Nb₂O₅, 0.46 F, 9.1 H₂O, ?0.19 O=F₂, total is 101.75. The empirical formula calculated on the basis of Si + Al = 4 is: H_7.42(Zn_1.30Na_0.23Mn_0.22Ca_0.07Mg_0.07K_0.05Ba_0.02)_Σ1.96(Ti_2.68Nb_0.41Fe _0.18 ³⁺ Zr_0.12)_Σ3.39(Si_3.95Al_0.05)_{Σ4 20.31}F_0.18. The simplified formula is: Zn₂Ti_4?x Si₄O₁₄(OH,H₂O,□)₈ (x < 1). Vigrishinite is triclinic, space group P $\bar 1$ , a = 8.743(9), b = 8.698(9), c = 11.581(11)Å, α = 91.54(8)°, β = 98.29(8)°, γ = 105.65(8)°, V = 837.2(1.5) ų, Z = 2. The strongest reflections in the X-ray powder pattern (d, Å, ?I[hkl]) are: 11.7-67[001], 8.27-50[100], 6.94-43[0 $\bar 1$ 1, $\bar 1$ 10], 5.73–54[1 $\bar 1$ 1, 002], 4.17-65[020, $\bar 1$ $\bar 1$ 2, 200], and 2.861-100[3 $\bar 1$ 0, 2 $\bar 2$ 2, 004, 1 $\bar 3$ 1]. The crystal structure model was obtained on a single crystal, R = 0.171. Vigrishinite and murmanite are close in the structure of the TiSiO motif, but strongly differ from each other in part of large cations and H-bearing groups. Vigrishinite is named in honor of Viktor G. Grishin (b. 1953), a Russian amateur mineralogist and mineral collector, to pay tribute to his contribution to the mineralogy of the Lovozero Complex. The type specimen is deposited in the Fersman Mineralogical Museum of Russian Academy of Sciences, Moscow.     

The origin of the rotational and spatial velocities of radio pulsars

We analyze possible origins of the observed high rotational and spatial velocities of radio pulsars. In particular, these can be understood if all radio pulsars originate in close binary systems with orbital periods of 0.1–100 days, with the neutron star being formed by a type Ib,c supernova. The high spatial velocities of pulsars (v _p up to 1000 km/s) reflect the high Keplerian velocities of the components of these binaries, while their short periods of rotation (P _p < 4 s) are due to the rapid rotation of the presupernova helium-star components with masses of 2.5–10 M_⊙, which is synchronous with their orbital rotation. Single massive stars or components in wide binaries are likely to produce only slowly rotating (P _p > 4 s) neutron stars or black holes, which cannot be radio pulsars. As a result, the rate of formation of radio pulsars should be a factor of a few lower than the rate of type II and type Ib,c supernovae estimated from observations. This scenario for the formation of radio pulsars is supported by (i) the bimodal spatial velocity distribution of radio pulsars; (ii) the coincidence of the observed spatial velocities of radio pulsars with the orbital velocities of the components of close binaries with nondegenerate helium presupernovae; (iii) the correlation between the orbital and rotational periods for 22 observed radio pulsars in binaries with elliptical orbits; and (iv) the similarity of the observed rate of formation of radio pulsars and the rate of type Ib,c supernovae.     

Thermal equation of state of synthetic orthoferrosilite at lunar pressures and temperatures

Iron-rich orthopyroxene plays an important role in models of the thermal and magmatic evolution of the Moon, but its density at high pressure and high temperature is not well-constrained. We present in situ measurements of the unit-cell volume of a synthetic polycrystalline end-member orthoferrosilite (FeSiO₃, fs) at simultaneous high pressures (3.4–4.8 GPa) and high temperatures (1,148–1,448 K), to improve constraints on the density of orthopyroxene in the lunar interior. Unit-cell volumes were determined through in situ energy-dispersive synchrotron X-ray diffraction in a multi-anvil press, using MgO as a pressure marker. Our volume data were fitted to a high-temperature Birch–Murnaghan equation of state (EoS). Experimental data are reproduced accurately, with a  $\varDelta P$ Δ P  standard deviation of 0.20 GPa. The resulting thermoelastic parameters of fs are: V ₀ = 875.8 ± 1.4 Å³, K ₀ = 74.4 ± 5.3 GPa, and $\frac{{\text d}K}{{\text d}T} = -0.032 \pm 0.005\,\hbox{GPa K}^{-1}$ d K d T = - 0.032 ± 0.005 GPa K - 1 , assuming ${K}^{\prime}_{0} = 10 $ K 0 ′ = 10 . We also determined the thermal equation of state of a natural Fe-rich orthopyroxene from Hidra (Norway) to assess the effect of magnesium on the EoS of iron-rich orthopyroxene. Comparison between our two data sets and literature studies shows good agreement for room-temperature, room-pressure unit-cell volumes. Preliminary thermodynamic analyses of orthoferrosilite, FeSiO₃, and orthopyroxene solid solutions, (Mg_1?x Fe _x ) SiO₃, using vibrational models show that our volume measurements in pressure–temperature space are consistent with previous heat capacity and one-bar volume–temperature measurements. The isothermal bulk modulus at ambient conditions derived from our measurements is smaller than values presented in the literature. This new simultaneous high-pressure, high-temperature data are specifically useful for calculations of the orthopyroxene density in the Moon.