Search

Alkali-metal optically-pumped magnetometers are prone to inaccuracies arising from the overlap of the average F = I + 1/2 and F = I - 1/2 ground-state Zeeman resonances. We employ density-matrix simulations and experiments to investigate how this hyperfine systematic error varies with spin polarization in a $^{87}$Rb free-induction-decay (FID) magnetometer. At low spin polarizations, ($P \leq 0.5$), this effect causes single-frequency magnetic-field extraction techniques to exhibit inaccuracies up to approximately 3.5 nT. Density-matrix simulations reveal that this bias can be traced to the relative amplitude and phase between the F = I $\pm$ 1/2 hyperfine ground-state manifolds in the FID spin precession signal. We show that this systematic error can be mitigated using either a double-frequency fitting model that accounts for the relative amplitude and phase or synchronous-pulse pumping, that minimizes the F = 1 contribution to the FID signal. Theoretical simulations predict accuracies within 0.5 nT for both techniques across a wide range of spin polarizations, suggesting a sevenfold enhancement over single-frequency extraction methods. Our experiments validate this, showcasing a variation in the extracted field below 1 nT, a 3.5-fold improvement compared to single-frequency extraction methods. Furthermore, the mitigation techniques demonstrate agreement of the extracted magnetic field within 1.5 nT., Comment: 13 pages, 8 figures

The Centaur (60558) Echeclus was discovered on March 03, 2000, orbiting between the orbits of Jupiter and Uranus. After exhibiting frequent outbursts, it also received a comet designation, 174P. If the ejected material can be a source of debris to form additional structures, studying the surroundings of an active body like Echeclus can provide clues about the formation scenarios of rings, jets, or dusty shells around small bodies. Stellar occultation is a handy technique for this kind of investigation, as it can, from Earth-based observations, detect small structures with low opacity around these objects. Stellar occultation by Echeclus was predicted and observed in 2019, 2020, and 2021. We obtain upper detection limits of rings with widths larger than 0.5 km and optical depth of $\tau$ = 0.02. These values are smaller than those of Chariklo's main ring; in other words, a Chariklo-like ring would have been detected. The occultation observed in 2020 provided two positive chords used to derive the triaxial dimensions of Echeclus based on a 3D model and pole orientation available in the literature. We obtained $a = 37.0\pm0.6$ km, $b = 28.4 \pm 0.5$ km, and $c= 24.9 \pm 0.4$ km, resulting in an area-equivalent radius of $30.0 \pm 0.5$ km. Using the projected limb at the occultation epoch and the available absolute magnitude ($\rm{H}_{\rm{v}} = 9.971 \pm 0.031$), we calculate an albedo of $p_{\rm{v}} = 0.050 \pm 0.003$. Constraints on the object's density and internal friction are also proposed., Comment: Accepted for publication on MNRAS

X-ray imaging is a fast, precise and non-invasive method of imaging which, combined with computed tomography, provides detailed 3D rendering of samples. Incorporating propagation-based phase contrast can vastly improve data quality for weakly attenuating samples via material-specific phase retrieval filters, allowing radiation exposure to be reduced. However, applying phase retrieval to multi-material phantoms complicates analysis by requiring a choice of which material boundary to tune the phase retrieval. Filtering for the boundary with strongest phase contrast increases noise suppression, but with the detriment of over-blurring other interfaces, potentially obscuring small or neighbouring features and removing quantitative sample information. Additionally, regions bounded by more than one material type inherently cannot be conventionally filtered to reconstruct the whole boundary. As remedy, we present a computationally-efficient, non-iterative nor AI-mediated method for applying strong phase retrieval, whilst preserving sharp boundaries for all materials within the sample. This technique was tested on phase contrast images of a rabbit kitten brain encased by the surrounding dense skull. Using 24 keV synchrotron radiation with a 5 m propagation distance, our technique provided a 6.9-fold improvement in the signal-to-noise ratio (SNR) of brain tissue compared to the standard phase retrieval procedure, without over-smoothing the images. Simultaneous quantification of edge resolution and SNR gain was performed with an aluminium-water phantom imaged using a microfocus X-ray tube at mean energy 19.58 keV and 0.576 m effective propagation distance. Our method provided a 4.2-fold SNR boost whilst preserving the boundary resolution at 54 $\pm$ 1 $\mu$m, compared to 108 $\pm$ 2 $\mu$m in conventional phase retrieval.

We present a robust technique for calculating regional volume changes within the lung from X-ray radiograph sequences captured during ventilation, without the use of computed tomography (CT). This technique is based on the change in transmitted X-ray intensity that occurs for each lung region as air displaces the attenuating lung tissue. Lung air volumes calculated from X-ray intensity changes showed a strong correlation ($R^2$=0.98) against the true volumes, measured from high-resolution CT. This correlation enables us to accurately convert projected intensity data into relative changes in lung air volume. We have applied this technique to measure changes in regional lung volumes from X-ray image sequences of mechanically ventilated, recently-deceased newborn rabbits, without the use of CT. This method is suitable for biomedical research studies and shows potential for clinical application., Comment: 12 Pages, 10 Figures

A stellar occultation by Neptune's main satellite, Triton, was observed on 5 October 2017 from Europe, North Africa, and the USA. We derived 90 light curves from this event, 42 of which yielded a central flash detection. We aimed at constraining Triton's atmospheric structure and the seasonal variations of its atmospheric pressure since the Voyager 2 epoch (1989). We also derived the shape of the lower atmosphere from central flash analysis. We used Abel inversions and direct ray-tracing code to provide the density, pressure, and temperature profiles in the altitude range $\sim$8 km to $\sim$190 km, corresponding to pressure levels from 9 {\mu}bar down to a few nanobars. Results. (i) A pressure of 1.18$\pm$0.03 {\mu}bar is found at a reference radius of 1400 km (47 km altitude). (ii) A new analysis of the Voyager 2 radio science occultation shows that this is consistent with an extrapolation of pressure down to the surface pressure obtained in 1989. (iii) A survey of occultations obtained between 1989 and 2017 suggests that an enhancement in surface pressure as reported during the 1990s might be real, but debatable, due to very few high S/N light curves and data accessible for reanalysis. The volatile transport model analysed supports a moderate increase in surface pressure, with a maximum value around 2005-2015 no higher than 23 {\mu}bar. The pressures observed in 1995-1997 and 2017 appear mutually inconsistent with the volatile transport model presented here. (iv) The central flash structure does not show evidence of an atmospheric distortion. We find an upper limit of 0.0011 for the apparent oblateness of the atmosphere near the 8 km altitude., Comment: 52 pages, 26 figures in the main paper, 2 figures in appendix B, 9 figures in appendix C, 1 long table over 5 pages

The rotation state of small asteroids is affected by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, which is a net torque caused by solar radiation directly reflected and thermally reemitted from the surface. Due to this effect, the rotation period slowly changes, which can be most easily measured in light curves because the shift in the rotation phase accumulates over time quadratically. We collected archived light curves and carried out new photometric observations for asteroids (10115) 1992 SK, (1620) Geographos, and (1685) Toro. We applied the method of light curve inversion to fit observations with a convex shape model. The YORP effect was modeled as a linear change of the rotation frequency $\upsilon \equiv \mathrm{d}\omega / \mathrm{d}t$ and optimized together with other spin and shape parameters. We detected the acceleration $\upsilon = (8.3 \pm 0.6) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ of the rotation for asteroid (10115) 1992 SK. This observed value agrees well with the theoretical value of YORP-induced spin-up computed for our shape and spin model. For (1685) Toro, we obtained $\upsilon = (3.3 \pm 0.3) \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$, which confirms an earlier tentative YORP detection. For (1620) Geographos, we confirmed the previously detected YORP acceleration and derived an updated value of $\upsilon$ with a smaller uncertainty. We also included the effect of solar precession into our inversion algorithm, and we show that there are hints of this effect in Geographos' data. The detected change of the spin rate of (10115) 1992 SK has increased the total number of asteroids with YORP detection to ten. In all ten cases, the $\mathrm{d}\omega / \mathrm{d}t$ value is positive, so the rotation of these asteroids is accelerated. It is unlikely to be just a statistical fluke, but it is probably a real feature that needs to be explained.

The Centaur (10199) Chariklo has the first rings system discovered around a small object. It was first observed using stellar occultation in 2013. Stellar occultations allow the determination of sizes and shapes with kilometre accuracy and obtain characteristics of the occulting object and its vicinity. Using stellar occultations observed between 2017 and 2020, we aim at constraining Chariklo's and its rings physical parameters. We also determine the rings' structure, and obtain precise astrometrical positions of Chariklo. We predicted and organised several observational campaigns of stellar occultations by Chariklo. Occultation light curves were measured from the data sets, from which ingress and egress times, and rings' width and opacity were obtained. These measurements, combined with results from previous works, allow us to obtain significant constraints on Chariklo's shape and rings' structure. We characterise Chariklo's ring system (C1R and C2R), and obtain radii and pole orientations that are consistent with, but more accurate than, results from previous occultations. We confirmed the detection of W-shaped structures within C1R and an evident variation of radial width. The observed width ranges between 4.8 and 9.1 km with a mean value of 6.5 km. One dual observation (visible and red) does not reveal any differences in the C1R opacity profiles, indicating ring particle's size larger than a few microns. The C1R ring eccentricity is found to be smaller than 0.022 (3-sigma), and its width variations may indicate an eccentricity higher than 0.005. We fit a tri-axial shape to Chariklo's detections over eleven occultations and determine that Chariklo is consistent with an ellipsoid with semi-axes of 143.8, 135.2 and 99.1 km. Ultimately, we provided seven astrometric positions at a milliarcseconds accuracy level, based on Gaia EDR3, and use it to improve Chariklo's ephemeris., Comment: 32 pages, 11 Figures in the main text, paper was accepted for publication in Section 10. Planets and planetary systems of Astronomy and Astrophysics on 12/07/2021

The ill-posed problem of phase retrieval in optics, using one or more intensity measurements, has a multitude of applications using electromagnetic or matter waves. Many phase retrieval algorithms are computed on pixel arrays using discrete Fourier transforms due to their high computational efficiency. However, the mathematics underpinning these algorithms is typically formulated using continuous mathematics, which can result in a loss in spatial resolution in the reconstructed images. Herein we investigate how phase retrieval algorithms for propagation-based phase-contrast X-ray imaging can be rederived using discrete mathematics and result in more precise retrieval for single- and multi-material objects and for spectral image decomposition. We validate this theory through experimental measurements of spatial resolution using computed tomography (CT) reconstructions of plastic phantoms and biological tissue, using detectors with a range of imaging system point spread functions (PSFs). We demonstrate that if the PSF substantially suppresses high spatial frequencies, the potential improvement from utilising the discrete derivation is limited. However, with detectors characterised by a single pixel PSF (e.g. direct, photon-counting X-ray detectors), a significant improvement in spatial resolution can be obtained, demonstrated here at up to 17%.

Light-shifts are known to be an important limitation to the mid- and long-term fractional frequency stability of different types of atomic clocks. In this article, we demonstrate the experimental implementation of an advanced anti-light shift interrogation protocol onto a continuous-wave (CW) microcell atomic clock based on coherent population trapping (CPT). The method, inspired by the Auto-Balanced Ramsey (ABR) spectroscopy technique demonstrated in pulsed atomic clocks, consists in the extraction of atomic-based information from two successive light-shifted clock frequencies obtained at two different laser power values. Two error signals, computed from the linear combination of signals acquired along a symmetric sequence, are managed in a dual-loop configuration to generate a clock frequency free from light-shift. Using this method, the sensitivity of the clock frequency to both laser power and microwave power variations can be reduced by more than an order of magnitude compared to normal operation. In the present experiment, the consideration of the non-linear light-shift dependence allowed to enhance light-shift mitigation. The implemented technique allows a clear improvement of the clock Allan deviation for time scales higher than 1000 s. This method could be applied in various kinds of atomic clocks such as CPT-based atomic clocks, double-resonance Rb clocks, or cell-stabilized lasers., Comment: 9 pages with 7 figures

We show that the light shift in atomic clocks can be suppressed using time variation of the interrogation field intensity. By measuring the clock output at two intensity levels, error signals can be generated that simultaneously stabilize a local oscillator to an atomic transition and correct for the shift of this transition caused by the interrogating optical field. These methods are suitable for optical clocks using one- and two-photon transitions, as well as for microwave clocks based on coherent population trapping or direct interrogation. The proposed methods can be widely used both for high-precision scientific instruments and for a wide range of commercial clocks, including chip-scale atomic clocks., Comment: 7 pages, 6 figures

Coherent dark states in atoms, created by simultaneous interaction of two coherent light fields with a 3-level system, are of prime importance in quantum state manipulation. They are used extensively in quantum sensing and quantum information applications to build atomic clocks, magnetometers, atomic interferometers and more. Here we study the formation and decay of coherent dark-states in an ensemble of laser-cooled free-falling atoms. We measure the optical-pumping rate into the dark-state in the $\sigma ^+-\sigma ^-$ polarization configuration. We find that the pumping rate is linear with the optical field intensity, but about an order-of-magnitude slower than the rate predicted by the commonly used, but simplistic, three-level analytic formula. Using a numerical model we demonstrate that this discrepancy is due to the multi-level Zeeman manifold. Taking into account the slower pumping rate we explain quantitatively the relation between the light-shift and the duration of pumping into dark-state in Ramsey spectroscopy. We also measure the decay of the dark-state coherence and find that in our apparatus it is dominated by the mechanical motion of the atoms out of the probing region, while the atomic decoherence is negligible.

We present results from three world-wide campaigns that resulted in the detections of two single-chord and one multi-chord stellar occultations by the Plutino object (84922) 2003~VS$_2$. From the single-chord occultations in 2013 and 2014 we obtained accurate astrometric positions for the object, while from the multi-chord occultation on November 7th, 2014, we obtained the parameters of the best-fitting ellipse to the limb of the body at the time of occultation. We also obtained short-term photometry data for the body in order to derive its rotational phase during the occultation. The rotational light curve present a peak-to-peak amplitude of 0.141 $\pm$ 0.009 mag. This allows us to reconstruct the three-dimensional shape of the body, with principal semi-axes $a = 313.8 \pm 7.1$ km, $b = 265.5^{+8.8}_{-9.8}$ km, and $c = 247.3^{+26.6}_{-43.6}$ km, which is not consistent with a Jacobi triaxial equilibrium figure. The derived spherical volume equivalent diameter of $548.3 ^{+29.5}_{-44.6}$ km is about 5\% larger than the radiometric diameter of 2003~VS$_2$ derived from Herschel data of $523 \pm 35$ km, but still compatible with it within error bars. From those results we can also derive the geometric albedo ($0.123 ^{+0.015}_{-0.014}$) and, under the assumption that the object is a Maclaurin spheroid, the density $\rho = 1400^{+1000}_{-300}$ for the plutino. The disappearances and reappearances of the star during the occultations do not show any compelling evidence for a global atmosphere considering a pressure upper limit of about 1 microbar for a pure nitrogen atmosphere, nor secondary features (e.g. rings or satellite) around the main body., Comment: 22 pages, 14 figures, 1 appendix

We present a new non-convex model of the binary asteroid (809) Lundia. A SAGE (Shaping Asteroids with Genetic Evolution) method using disc-integrated photometry only was used for deriving physical parameters of this binary system. The model of (809) Lundia improves former system's pole solution and gives the ecliptic coordinates of the orbit pole - $\lambda=122^{\circ}$, $\beta=22^{\circ}$, $\sigma=\pm5^{\circ}$ - and the orbital period of $15.41574 \pm 0.00001$ h. For scaling our results we used effective diameter of $D_{eff} = 9.6 \pm 1.1$ km obtained from Spitzer observations. The non-convex shape description of the components permitted a refined calculation of the components' volumes, leading to a density estimation of $2.5\pm0.2$ g/cm$^3$ and macroporosity of 13-23\%. The intermediate-scale features of the model may also offer new clues on the components' origin and evolution.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere., Comment: 21 pages, 11 figures

Light-induced frequency shifts can be a key limiting contribution to the mid and long-term frequency instability in atomic clocks. In this letter, we demonstrate the experimental implementation of the combined error signal interrogation protocol to a cold-atom clock based on coherent population trapping (CPT) and Ramsey spectroscopy. The method uses a single error signal that results from the normalized combination of two error signals extracted from two Ramsey sequences of different dark periods. The single combined error signal is used to stabilize the atomic clock frequency. Compared to the standard Ramsey-CPT interrogation, this method reduces the clock frequency sensitivity to light-shift variations by more than one order of magnitude. This method can be applied in various kinds of Ramsey-based atomic clocks, sensors and instruments.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process., Comment: Submitted to Icarus on 2019 January 7. The Electronic Supplementary Information to this paper is at http://www.asu.cas.cz/~asteroid/astpairscomplex_si.pdf

Sophisticated Ramsey-based interrogation protocols using composite laser pulse sequences have been recently proposed to provide next-generation high-precision atomic clocks with a near perfect elimination of frequency shifts induced during the atom-probing field interaction. We propose here a simple alternative approach to the auto-balanced Ramsey (ABR) interrogation protocol and demonstrate its application to a cold-atom microwave clock based on coherent population trapping. The main originality of the method, based on two consecutive Ramsey sequences with different dark periods, is to sample the central Ramsey fringes with frequency-jumps finely-adjusted by an additional frequency-displacement concomitant parameter, scaling as the inverse of the dark-period. The advantage of this displaced frequency-jumps Ramsey (DFJR) method is that the local oscillator frequency is used as a single physical variable to control both servo loops of the sequence, simplifying its implementation and avoiding noise associated with controlling the LO phase. Compared to the usual Ramsey-CPT technique, the DFJR scheme reduces the sensitivity of the clock frequency to variations of the CPT sideband ratio and to the one-photon laser detuning by more than an order of magnitude. This simple method could be applied in a wide variety of Ramsey-spectroscopy based applications including frequency metrology with CPT-based and optical atomic clocks, mass spectrometry, and precision spectroscopy.

We have developed a universal method to form the reference signal for the stabilization of arbitrary atomic clocks based on Ramsey spectroscopy. Our approach uses an interrogation scheme of the atomic system with two different Ramsey periods and a specially constructed combined error signal (CES) computed by subtracting two error signals with the appropriate calibration factor. CES spectroscopy allows for perfect elimination of probe-induced light shifts and does not suffer from the effects of relaxation, time-dependent pulse fluctuations and phase-jump modulation errors and other imperfections of the interrogation procedure. The method is simpler than recently developed auto-balanced Ramsey spectroscopy techniques [Ch. Sanner, et al., Phys. Rev. Lett. 120, 053602 (2018); V. I. Yudin, et al., Phys. Rev. Appl. 9, 054034 (2018)], because it uses a single error signal that feeds back on the clock frequency. CES universal technique can be applied to many applications of precision spectroscopy., Comment: 12 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:1712.03365

The off-resonant AC Stark shift for coherent population trapping (CPT) resonances probed with Ramsey spectroscopy is investigated experimentally and theoretically. Measurements with laser-cooled $^{87}$Rb atoms show excellent quantitative agreement with a simple theory. The shift depends on the relative intensity of the two CPT light fields, but depends only weakly on the total intensity. Since the origin of the shift is through couplings of the interrogation light to off-resonant excited state hyperfine levels, the size and sign of the shift depend on the specific interrogation scheme. However, the simple theory shows that for several commonly used interrogation schemes the off-resonant shift goes to zero at specific values of this relative intensity., Comment: 7 pages, 7 figures

When performing precision measurements, the quantity being measured is often perturbed by the measurement process itself. This includes precision frequency measurements for atomic clock applications carried out with Ramsey spectroscopy. With the aim of eliminating probe-induced perturbations, a method of generalized auto-balanced Ramsey spectroscopy (GABRS) is presented and rigorously substantiated. Here, the usual local oscillator frequency control loop is augmented with a second control loop derived from secondary Ramsey sequences interspersed with the primary sequences and with a different Ramsey period. This second loop feeds back to a secondary clock variable and ultimately compensates for the perturbation of the clock frequency caused by the measurements in the first loop. We show that such a two-loop scheme can lead to perfect compensation of measurement-induced light shifts and does not suffer from the effects of relaxation, time-dependent pulse fluctuations and phase-jump modulation errors that are typical of other hyper-Ramsey schemes. Several variants of GABRS are explored based on different secondary variables including added relative phase shifts between Ramsey pulses, external frequency-step compensation, and variable second-pulse duration. We demonstrate that a universal anti-symmetric error signal, and hence perfect compensation at finite modulation amplitude, is generated only if an additional frequency-step applied during both Ramsey pulses is used as the concomitant variable parameter. This universal technique can be applied to the fields of atomic clocks, high-resolution molecular spectroscopy, magnetically induced and two-photon probing schemes, Ramsey-type mass spectrometry, and to the field of precision measurements. Some variants of GABRS can also be applied for rf atomic clocks using CPT-based Ramsey spectroscopy of the two-photon dark resonance., Comment: 11 pages, 7 figures

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

Two narrow and dense rings (called C1R and C2R) were discovered around the Centaur object (10199) Chariklo during a stellar occultation observed on 2013 June 3. Following this discovery, we planned observations of several occultations by Chariklo's system in order to better characterize the physical properties of the ring and main body. Here, we use 12 successful occulations by Chariklo observed between 2014 and 2016. They provide ring profiles (physical width, opacity, edge structure) and constraints on the radii and pole position. Our new observations are currently consistent with the circular ring solution and pole position, to within the $\pm 3.3$ km formal uncertainty for the ring radii derived by Braga-Ribas et al. The six resolved C1R profiles reveal significant width variations from $\sim 5$ to 7.5 km. The width of the fainter ring C2R is less constrained, and may vary between 0.1 and 1 km. The inner and outer edges of C1R are consistent with infinitely sharp boundaries, with typical upper limits of one kilometer for the transition zone between the ring and empty space. No constraint on the sharpness of C2R's edges is available. A 1$\sigma$ upper limit of $\sim 20$ m is derived for the equivalent width of narrow (physical width <4 km) rings up to distances of 12,000 km, counted in the ring plane.

Photometric lightcurves were also obtained for 14 of them during the Spitzer observations to provide the context of the observations and reliable estimates of their absolute magnitudes. The extracted mid-IR spectra were analyzed using a modified standard thermal model (STM) and a thermophysical model (TPM) that takes into account the shape and geometry of the large primary at the time of the Spitzer observation. We derived a reliable estimate of the size, albedo, and beaming factor for each of these asteroids, representing three main taxonomic groups: C, S, and X. For large (volume-equivalent system diameter Deq $\lt$ 130 km) binary asteroids, the TPM analysis indicates a low thermal inertia ($\Gamma$ < $\sim$100 J s-1/2K-1m-2) and their emissivity spectra display strong mineral features, implying that they are covered with a thick layer of thermally insulating regolith. The smaller (surface-equivalent system diameter Deff $\lt$17 km) asteroids also show some emission lines of minerals, but they are significantly weaker, consistent with regoliths with coarser grains, than those of the large binary asteroids. The average bulk densities of these multiple asteroids vary from 0.7-1.7 g/cm3 (P-, C- type) to $\sim$2 g/cm3 (S-type). The highest density is estimated for the M-type (22) Kalliope (3.2 $\pm$ 0.9 g/cm3). The spectral energy distributions (SED) and emissivity spectra, made available as a supplement document, could help to constrain the surface compositions of these asteroids., Comment: 69 pages, 5 Figures, 8 Tables

We analyze two multi-chord stellar occultations by Pluto observed on July 18th, 2012 and May 4th, 2013, and monitored respectively from five and six sites. They provide a total of fifteen light-curves, twelve of them being used for a simultaneous fit that uses a unique temperature profile, assuming a clear (no-haze) and pure N_2 atmosphere, but allowing for a possible pressure variation between the two dates. We find a solution that fits satisfactorily (i.e. within the noise level) all the twelve light-curves, providing atmospheric constraints between ~1,190 km (pressure ~ 11 \mubar) and ~ 1,450 km (pressure ~0.1 \mubar) from Pluto's center. Our main results are: (1) the best-fitting temperature profile shows a stratosphere with strong positive gradient between 1,190 km (at 36 K, 11 \mubar) and r = 1,215 km (6.0 \mubar), where a temperature maximum of 110 K is reached; above it is a mesosphere with negative thermal gradient of -0.2 K/km up to ~ 1,390 km (0.25 \mubar), where, the mesosphere connects itself to a more isothermal upper branch around 81 K; (2) the pressure shows a small (6 %) but significant increase (6-\sigma level) between the two dates; (3) without troposphere, Pluto's radius is found to be R_P = 1,190 +/- 5km. Allowing for a troposphere, R_P is constrained to lie between 1,168 and 1,195 km; (4) the currently measured CO abundance is too small to explain the mesospheric negative thermal gradient. Cooling by HCN is possible, but only if this species is largely saturated; Alternative explanations like zonal winds or vertical compositional variations of the atmosphere are unable to explain the observed mesospheric trend., Comment: 41 pages, 13 figures, 4 tables

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

Using our photometric observations taken between 1996 and 2013 and other published data, we derived properties of the binary near-Earth asteroid (175706) 1996 FG3 including new measurements constraining evolution of the mutual orbit with potential consequences for the entire binary asteroid population. We also refined previously determined values of parameters of both components, making 1996 FG3 one of the most well understood binary asteroid systems. We determined the orbital vector with a substantially greater accuracy than before and we also placed constraints on a stability of the orbit. Specifically, the ecliptic longitude and latitude of the orbital pole are 266{\deg} and -83{\deg}, respectively, with the mean radius of the uncertainty area of 4{\deg}, and the orbital period is 16.1508 +/- 0.0002 h (all quoted uncertainties correspond to 3sigma). We looked for a quadratic drift of the mean anomaly of the satellite and obtained a value of 0.04 +/- 0.20 deg/yr^2, i.e., consistent with zero. The drift is substantially lower than predicted by the pure binary YORP (BYORP) theory of McMahon and Scheeres (McMahon, J., Scheeres, D. [2010]. Icarus 209, 494-509) and it is consistent with the theory of an equilibrium between BYORP and tidal torques for synchronous binary asteroids as proposed by Jacobson and Scheeres (Jacobson, S.A., Scheeres, D. [2011]. ApJ Letters, 736, L19). Based on the assumption of equilibrium, we derived a ratio of the quality factor and tidal Love number of Q/k = 2.4 x 10^5 uncertain by a factor of five. We also derived a product of the rigidity and quality factor of mu Q = 1.3 x 10^7 Pa using the theory that assumes an elastic response of the asteroid material to the tidal forces. This very low value indicates that the primary of 1996 FG3 is a 'rubble pile', and it also calls for a re-thinking of the tidal energy dissipation in close asteroid binary systems., Comment: Many changes based on referees comments

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

To take full advantage of the September 2008 opposition passage of the M-type asteroid (216) Kleopatra, we have used near-infrared adaptive optics (AO) imaging with the W.M. Keck II telescope to capture unprecedented high resolution images of this unusual asteroid. Our AO observations with the W.M. Keck II telescope, combined with Spitzer/IRS spectroscopic observations and past stellar occultations, confirm the value of its IRAS radiometric radius of 67.5 km as well as its dog-bone shape suggested by earlier radar observations. Our Keck AO observations revealed the presence of two small satellites in orbit about Kleopatra (see Marchis et al., 2008). Accurate measurements of the satellite orbits over a full month enabled us to determine the total mass of the system to be 4.64+/-0.02 10^18 Kg. This translates into a bulk density of 3.6 +/-0.4 g/cm3, which implies a macroscopic porosity for Kleopatra of ~ 30-50%, typical of a rubble-pile asteroid. From these physical characteristics we measured its specific angular momentum, very close to that of a spinning equilibrium dumbbell., Comment: 35 pages, 3 Tables, 9 Figures. In press to Icarus

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation., Comment: 12 pages, 2 figures, 1 table + Supplementary Information

Mutual event observations between the two components of 90 Antiope were carried out in 2007-2008. The pole position was refined to lambda0 = 199.5+/-0.5 eg and beta0 = 39.8+/-5 deg in J2000 ecliptic coordinates, leaving intact the physical solution for the components, assimilated to two perfect Roche ellipsoids, and derived after the 2005 mutual event season (Descamps et al., 2007). Furthermore, a large-scale geological depression, located on one of the components, was introduced to better match the observed lightcurves. This vast geological feature of about 68 km in diameter, which could be postulated as a bowl-shaped impact crater, is indeed responsible of the photometric asymmetries seen on the "shoulders" of the lightcurves. The bulk density was then recomputed to 1.28+/-0.04 gcm-3 to take into account this large-scale non-convexity. This giant crater could be the aftermath of a tremendous collision of a 100-km sized proto-Antiope with another Themis family member. This statement is supported by the fact that Antiope is sufficiently porous (~50%) to survive such an impact without being wholly destroyed. This violent shock would have then imparted enough angular momentum for fissioning of proto-Antiope into two equisized bodies. We calculated that the impactor must have a diameter greater than ~17 km, for an impact velocity ranging between 1 and 4 km/s. With such a projectile, this event has a substantial 50% probability to have occurred over the age of the Themis family., Comment: 30 pages, 3 Tables, 8 Figures. Accepted for publication in Icarus

We report on the results of a six-month photometric study of the main-belt binary C-type asteroid 121 Hermione, performed during its 2007 opposition. We took advantage of the rare observational opportunity afforded by one of the annual equinoxes of Hermione occurring close to its opposition in June 2007. The equinox provides an edge-on aspect for an Earth-based observer, which is well suited to a thorough study of Hermione's physical characteristics. The catalog of observations carried out with small telescopes is presented in this work, together with new adaptive optics (AO) imaging obtained between 2005 and 2008 with the Yepun 8-m VLT telescope and the 10-m Keck telescope. The most striking result is confirmation that Hermione is a bifurcated and elongated body, as suggested by Marchis et al., (2005). A new effective diameter of 187 +/- 6 km was calculated from the combination of AO, photometric and thermal observations. The new diameter is some 10% smaller than the hitherto accepted radiometric diameter based on IRAS data. The reason for the discrepancy is that IRAS viewed the system almost pole-on. New thermal observations with the Spitzer Space Telescope agree with the diameter derived from AO and lightcurve observations. On the basis of the new AO astrometric observations of the small 32-km diameter satellite we have refined the orbit solution and derived a new value of the bulk density of Hermione of 1.4 +0.5/-0.2 g cm-3. We infer a macroscopic porosity of ~33 +5/-20%., Comment: 39 pages, 4 Tables, 8 figures

In 2007, the M-type binary asteroid 22 Kalliope reached one of its annual equinoxes. As a consequence, the orbit plane of its small moon, Linus, was aligned closely to the Sun's line of sight, giving rise to a mutual eclipse season. A dedicated international campaign of photometric observations, based on amateur-professional collaboration, was organized and coordinated by the IMCCE in order to catch several of these events. The set of the compiled observations is released in this work. We developed a relevant model of these events, including a topographic shape model of Kalliope refined in the present work, the orbit solution of Linus as well as the photometric effect of the shadow of one component falling on the other. By fitting this model to the only two full recorded events, we derived a new estimation of the equivalent diameter of Kalliope of 166.2+/-2.8km, 8% smaller than its IRAS diameter. As to the diameter of Linus, considered as purely spherical, it is estimated to 28+/-2 km. This substantial "shortening" of Kalliope gives a bulk density of 3.35+/-0.33g/cm3, significantly higher than past determinations but more consistent with its taxonomic type. Some constraints can be inferred on the composition., Comment: 59 pages, 8 tables, 15 figures submitted to Icarus