Showing 1–50 of 186 results for author: Correia, C

Characterisation of the Warm-Jupiter TOI-1130 system with CHEOPS and photo-dynamical approach

Authors: L. Borsato, D. Degen, A. Leleu, M. J. Hooton, J. A. Egger, A. Bekkelien, A. Brandeker, A. Collier Cameron, M. N. Günther, V. Nascimbeni, C. M. Persson, A. Bonfanti, T. G. Wilson, A. C. M. Correia, T. Zingales, T. Guillot, A. H. M. J. Triaud, G. Piotto, D. Gandolfi, L. Abe, Y. Alibert, R. Alonso, T. Bárczy, D. Barrado Navascues, S. C. C. Barros , et al. (71 additional authors not shown)

Abstract: Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as "lonely" and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed… ▽ More Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as "lonely" and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed and well characterised. Within the contest of multi-planet system hosting gas-giant on short orbits, we characterise TOI-1130 system by measuring masses and orbital parameters. This is a 2-transiting planet system with a Jupiter-like planet (c) on a 8.35 days orbit and a Neptune-like planet (b) on an inner (4.07 days) orbit. Both planets show strong anti-correlated transit timing variations (TTVs). Furthermore, radial velocity (RV) analysis showed an additional linear trend, a possible hint of a non-transiting candidate planet on a far outer orbit. Since 2019, extensive transit and radial velocity observations of the TOI-1130 have been acquired using TESS and various ground-based facilities. We present a new photo-dynamical analysis of all available transit and RV data, with the addition of new CHEOPS and ASTEP+ data that achieve the best precision to date on the planetary radii and masses and on the timings of each transit. We were able to model interior structure of planet b constraining the presence of a gaseous envelope of H/He, while it was not possible to assess the possible water content. Furthermore, we analysed the resonant state of the two transiting planets, and we found that they lie just outside the resonant region. This could be the result of the tidal evolution that the system underwent. We obtained both masses of the planets with a precision less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 22 pages, 10 figures, Accepted for publication by Astronomy and Astrophysics

Unveiling the internal structure and formation history of the three planets transiting HIP 29442 (TOI-469) with CHEOPS

Authors: J. A. Egger, H. P. Osborn, D. Kubyshkina, C. Mordasini, Y. Alibert, M. N. Günther, M. Lendl, A. Brandeker, A. Heitzmann, A. Leleu, M. Damasso, A. Bonfanti, T. G. Wilson, S. G. Sousa, J. Haldemann, L. Delrez, M. J. Hooton, T. Zingales, R. Luque, R. Alonso, J. Asquier, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann , et al. (69 additional authors not shown)

Abstract: Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TE… ▽ More Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TESS data to derive planetary radii of $3.410\pm0.046$, $1.551\pm0.045$ and $1.538\pm0.049$ R$_\oplus$ for planets b, c and d, which orbit HIP 29442 with periods of 13.6, 3.5 and 6.4 days. For planet d, this value deviates by more than 3 sigma from the median value reported in the discovery paper, leading us to conclude that caution is required when using TESS photometry to determine the radii of small planets with low per-transit S/N and large gaps between observations. Given the high precision of these new radii, combining them with published RVs from ESPRESSO and HIRES provides us with ideal conditions to investigate the internal structure and formation pathways of the planets in the system. We introduce the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then apply hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identify planetary system analogues in a synthetic population generated with the Bern model for planet formation and evolution. Based on this analysis, we find that the planets likely formed on opposing sides of the water iceline from a protoplanetary disk with an intermediate solid mass. We finally report that the observed parameters of the HIP 29442 system are compatible with both a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope as well as a scenario with water-rich sub-Neptunes. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 30 pages, 17 figures, accepted for publication in A&A

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Evolution of current-carrying string networks

Authors: J. R. C. C. C. Correia, C. J. A. P. Martins, F. C. N. Q. Pimenta

Abstract: Cosmic string networks are expected to form in early Universe phase transitions via the Kibble mechanism and are unavoidable in several Beyond the Standard Model theories. While most predictions of observational signals of string networks assume featureless Abelian-Higgs or Nambu-Goto string networks, in many such extensions the networks can carry additional degrees of freedom, including charges a… ▽ More Cosmic string networks are expected to form in early Universe phase transitions via the Kibble mechanism and are unavoidable in several Beyond the Standard Model theories. While most predictions of observational signals of string networks assume featureless Abelian-Higgs or Nambu-Goto string networks, in many such extensions the networks can carry additional degrees of freedom, including charges and currents; these are often generically known as superconducting strings. All such degrees of freedom can impact the evolution of the networks and therefore their observational signatures. We report on the results of $2048^3$ field theory simulations of the evolution of a current-carrying network of strings, highlighting the different scaling behaviours of the network in the radiation and matter eras. We also report the first numerical measurements of the coherence length scales for the charge and current and of the condensate equation of state, and show that the latter mainly depends on the expansion rate, with chirality occurring for the matter era. Qualitatively, the fact that the matter era is the optimal expansion rate for these networks to reach scaling is in agreement with recent analytic modeling. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 11 pages, 2 figures, 3 tables; to appear in Phys. Lett. B

HIP 41378 observed by CHEOPS: Where is planet d?

Authors: S. Sulis, L. Borsato, S. Grouffal, H. P. Osborn, A. Santerne, A. Brandeker, M. N. Günther, A. Heitzmann, M. Lendl, M. Fridlund, D. Gandolfi, Y. Alibert, R. Alonso, T. Bárczy, D. Barrado Navascues, S. C. Barros, W. Baumjohann, T. Beck, W. Benz, M. Bergomi, N. Billot, A. Bonfanti, C. Broeg, A. Collier Cameron, C. Corral van Damme , et al. (62 additional authors not shown)

Abstract: HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but th… ▽ More HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but the observations show no transit. We find that large ($>22.4$ hours) transit timing variations (TTVs) could explain this non-detection during the CHEOPS observation window. We also investigated the possibility of an incorrect orbital solution, which would have major implications for our knowledge of this system. If $P_\mathrm{d} \neq 278.36$ d, the periods that minimize the eccentricity would be $101.22$ d and $371.14$ d. The shortest orbital period will be tested by TESS, which will observe HIP 41378 in Sector 88 starting in January 2025. Our study shows the importance of a mission like CHEOPS, which today is the only mission able to make long observations (i.e., from space) to track the ephemeris of long-period planets possibly affected by large TTVs. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 686, L18 (2024)

Photo-dynamical characterisation of the TOI-178 resonant chain

Authors: A. Leleu, J. -B. Delisle, L. Delrez, E. M. Bryant, A. Brandeker, H. P. Osborn, N. Hara, T. G. Wilson, N. Billot, M. Lendl, D. Ehrenreich, H. Chakraborty, M. N. Günther, M. J. Hooton, Y. Alibert, R. Alonso, D. R. Alves, D. R. Anderson, I. Apergis, D. Armstrong, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, M. P. Battley, W. Baumjohann , et al. (82 additional authors not shown)

Abstract: The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision ev… ▽ More The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, hence providing important anchors for planet formation models. We aim to improve the characterisation of the architecture of this key system, and in particular the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by each method. We show how stellar activity is preventing us from obtaining a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain. △ Less

Submitted 22 May, 2024; originally announced May 2024.

PASO -- Astronomy and Space Situational Awareness in a Dark Sky Destination

Authors: Domingos Barbosa, Bruno Coelho, Miguel Bergano, Constança Alves, Alexandre C. M. Correia, Luís Cupido, José Freitas, Luís Gonçalves, Bruce Grossan, Anna Guerman, Allan K. de Almeida Jr., Dalmiro Maia, Bruno Morgado, João Pandeirada, Valério Ribeiro, Gonçalo Rosa, George Smoot, Timothée Vaillant, Thyrso Villela, Carlos Alexandre Wuensche

Abstract: The Pampilhosa da Serra Space Observatory (PASO) is located in the center of the continental Portuguese territory, in the heart of a certified Dark Sky destination by the Starlight Foundation (Aldeias do Xisto) and has been an instrumental asset to advance science, education and astrotourism certifications. PASO hosts astronomy and Space Situational Awareness (SSA) activities including a node of t… ▽ More The Pampilhosa da Serra Space Observatory (PASO) is located in the center of the continental Portuguese territory, in the heart of a certified Dark Sky destination by the Starlight Foundation (Aldeias do Xisto) and has been an instrumental asset to advance science, education and astrotourism certifications. PASO hosts astronomy and Space Situational Awareness (SSA) activities including a node of the Portuguese Space Surveillance \& Tracking (SST) infrastructure network, such as a space radar currently in test phase using GEM radiotelescope, a double Wide Field of View Telescope system, a EUSST optical sensor telescope. These instruments allow surveillance of satellite and space debris in LEO, MEO and GEO orbits. The WFOV telescope offers spectroscopy capabilities enabling light curve analysis and cosmic sources monitoring. Instruments for Space Weather are being considered for installation to monitor solar activities and expand the range of SSA services. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 9 pages, 2 figures, Keywords: Dark Sky, radioastronomy, space, space debris, space situational awareness; Draft from Talk presented at Astronomy and Satellite Constellations: Pathways Forward, IAU Symposium 385, October 2023, Ed. C. Walker, D.Turnshek, P.Grimley, D.Galadi-Enriquez & M.Aubé, International Astronomical Union Proceedings Series, Cambridge University Press, 2024

Phasing segmented telescopes via deep learning methods: application to a deployable CubeSat

Authors: Maxime Dumont, Carlos M. Correia, Jean-François Sauvage, Noah Schwartz, Morgan Gray, Jaime Cardoso

Abstract: Capturing high resolution imagery of the Earth's surface often calls for a telescope of considerable size, even from Low Earth Orbits (LEO). A large aperture often requires large and expensive platforms. For instance, achieving a resolution of 1m at visible wavelengths from LEO typically requires an aperture diameter of at least 30cm. Additionally, ensuring high revisit times often prompts the use… ▽ More Capturing high resolution imagery of the Earth's surface often calls for a telescope of considerable size, even from Low Earth Orbits (LEO). A large aperture often requires large and expensive platforms. For instance, achieving a resolution of 1m at visible wavelengths from LEO typically requires an aperture diameter of at least 30cm. Additionally, ensuring high revisit times often prompts the use of multiple satellites. In light of these challenges, a small, segmented, deployable CubeSat telescope was recently proposed creating the additional need of phasing the telescope's mirrors. Phasing methods on compact platforms are constrained by the limited volume and power available, excluding solutions that rely on dedicated hardware or demand substantial computational resources. Neural Network (NN) are known for their computationally efficient inference and reduced on board requirements. Therefore we developed a NN based method to measure co phasing errors inherent to a deployable telescope. The proposed technique demonstrates its ability to detect phasing error at the targeted performance level (typically a wavefront error (WFE) below 15 nm RMS for a visible imager operating at the diffraction limit) using a point source. The robustness of the NN method is verified in presence of high order aberrations or noise and the results are compared against existing state of the art techniques. The developed NN model ensures its feasibility and provides a realistic pathway towards achieving diffraction limited images. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Journal ref: JOSA A 41.3 (2024): 489-499

Dynamical evolution of the Uranian satellite system II. Crossing of the 5/3 Ariel-Umbriel mean motion resonance

Authors: Sérgio R. A. Gomes, Alexandre C. M. Correia

Abstract: At present, the main satellites of Uranus are not involved in any low order mean motion resonance (MMR). However, owing to tides raised in the planet, Ariel and Umbriel most likely crossed the 5/3 MMR in the past. Previous studies on this resonance passage relied on limited time-consuming N-body simulations or simplified models focusing solely on the effects of the eccentricity or the inclination.… ▽ More At present, the main satellites of Uranus are not involved in any low order mean motion resonance (MMR). However, owing to tides raised in the planet, Ariel and Umbriel most likely crossed the 5/3 MMR in the past. Previous studies on this resonance passage relied on limited time-consuming N-body simulations or simplified models focusing solely on the effects of the eccentricity or the inclination. In this paper, we aim to provide a more comprehensive view on how the system evaded capture in the 5/3 MMR. For that purpose, we developed a secular resonant two-satellite model with low eccentricities and low inclinations, including tides using the weak friction model. By performing a large number of numerical simulations, we show that capture in the 5/3 MMR is certain if the initial eccentricities of Ariel, $e_1$, and Umbriel, $e_2$, are related through $(e_1^2 + e_2^2)^{1/2} < 0.007$. Moreover, we observe that the eccentricity of Ariel is the key variable to evade the 5/3 MMR with a high probability. We determine that for $e_1 > 0.015$ and $e_2 < 0.01$, the system avoids capture in at least 60\% of the cases. We also show that, to replicate the currently observed system, the initial inclinations of Ariel and Umbriel must lay within $I_1 \leq 0.05^{\circ}$ and $0.06^{\circ} \leq I_2 \leq 0.11^{\circ}$, respectively. We checked these results using a complete N-body model with the five main satellites and did not observe any significant differences. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: Submitted to Icarus, comments are welcome. arXiv admin note: text overlap with arXiv:2305.08794

Dynamical evolution of the Uranian satellite system I. From the 5/3 Ariel-Umbriel mean motion resonance to the present

Authors: Sérgio R. A. Gomes, Alexandre C. M. Correia

Abstract: Mutual gravitational interactions between the five major Uranian satellites raise small quasi-periodic fluctuations on their orbital elements. At the same time, tidal interactions between the satellites and the planet induce a slow outward drift of the orbits, while damping the eccentricities and the inclinations. In this paper, we revisit the current and near past evolution of this system using a… ▽ More Mutual gravitational interactions between the five major Uranian satellites raise small quasi-periodic fluctuations on their orbital elements. At the same time, tidal interactions between the satellites and the planet induce a slow outward drift of the orbits, while damping the eccentricities and the inclinations. In this paper, we revisit the current and near past evolution of this system using a N-body integrator, including spin evolution and tidal dissipation with the weak friction model. We update the secular eigenmodes of the system and show that it is unlikely that any of the main satellites were recently captured into a high obliquity Cassini state. We rather expect that the Uranian satellites are in a low obliquity Cassini state and compute their values. We also show that the current eccentricities of the satellites are not forced, and estimate the free eccentricities and inclinations. We constrain the quality factor of Uranus to be $Q_U = (8.6 \pm 2.9)\times10^3$, and that of the satellites to be $Q_S \sim 500$. We find that the system most likely encountered the 5/3 mean motion resonance between Ariel and Umbriel in the past, at about ($0.7\pm0.2$) Gyr ago. We additionally determine the eccentricities and inclinations of all satellites just after the resonance passage that comply with the current system. We finally show that, from the crossing of the 5/3 MMR to the present, the evolution of the system is mostly peaceful and dominated by tides raised on Uranus by the satellites. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: Submitted to Icarus

The tidal deformation and atmosphere of WASP-12b from its phase curve

Authors: B. Akinsanmi, S. C. C. Barros, M. Lendl, L. Carone, P. E. Cubillos, A. Bekkelien, A. Fortier, H. -G. Florén, A. Collier Cameron, G. Boué, G. Bruno, B. -O. Demory, A. Brandeker, S. G. Sousa, T. G. Wilson, A. Deline, A. Bonfanti, G. Scandariato, M. J. Hooton, A. C. M. Correia, O. D. S. Demangeon, A. M. S. Smith, V. Singh, Y. Alibert, R. Alonso , et al. (63 additional authors not shown)

Abstract: Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, a… ▽ More Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$σ$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared. △ Less

Submitted 20 February, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: accepted for publication in A&A

Signs of magnetic star-planet interactions in HD 118203. TESS detects stellar variability that matches the orbital period of a close-in eccentric Jupiter-sized companion

Authors: A. Castro-González, J. Lillo-Box, A. C. M. Correia, N. C. Santos, D. Barrado, M. Morales-Calderón, E. L. Shkolnik

Abstract: Planetary systems with close-in giant planets can experience magnetic star-planet interactions that modify the activity levels of their host stars. The induced activity is known to strongly depend on the magnetic moment of the interacting planet. Therefore, such planet-induced activity should be more readily observable in systems with planets in eccentric orbits, since those planets are expected t… ▽ More Planetary systems with close-in giant planets can experience magnetic star-planet interactions that modify the activity levels of their host stars. The induced activity is known to strongly depend on the magnetic moment of the interacting planet. Therefore, such planet-induced activity should be more readily observable in systems with planets in eccentric orbits, since those planets are expected to rotate faster than in circular orbits. However, no evidence of magnetic interactions has been reported in eccentric systems to date. We intend to unveil a possible planet-induced activity in the bright ($V$ = 8.05 $\pm$ 0.03 mag) and slightly evolved star HD 118203, which hosts an eccentric ($e$ = 0.32 $\pm$ 0.02) and close-in ($a$ = 0.0864 $\pm$ 0.0006 au) Jupiter-sized planet. We characterized the system by modelling 56 ELODIE radial velocities and four sectors of TESS photometry. We searched for planet-induced and rotation-related activity signals within the TESS, ELODIE, and ASAS-SN public data. We studied the possible origins of the variability found, analysed its persistence and evolution, and searched for links with the eccentric orbital motion of HD 118203 b. We found evidence of an activity signal within the TESS data that matches the orbital period of HD 118203 b, which suggests the existence of magnetic star-planet interactions. We did not find, however, any additional signal that could be interpreted as the rotation of the star, so we cannot discard stellar rotation as the source of the signal found. Nevertheless, the evolved nature of the star and the orbital eccentricity make the synchronous stellar rotation very unlikely. HD 118203 represents the best evidence that magnetic star-planet interactions can be found in eccentric systems, and it opens the door to future dedicated searches that will allow us to better understand the interplay between close-in planets and their hosts. △ Less

Submitted 31 March, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A; 20 pages, 14 figures, 4 tables

Tangent Velocity constraint for orbital maneuvers with Theory of Functional Connections

Authors: A. K. de Almeida Jr., T. Vaillant, V. M. de Oliveira, D. Barbosa, D. Maia, S. Aljbaae, B. Coelho, M. Bergano, J. Pandeirada, A. F. B. A. Prado, A. Guerman, A. C. M. Correia

Abstract: Maneuvering a spacecraft in the cislunar space is a complex problem, since it is highly perturbed by the gravitational influence of both the Earth and the Moon, and possibly also the Sun. Trajectories minimizing the needed fuel are generally preferred in order to decrease the mass of the payload. A classical method to constrain maneuvers is mathematically modelling them using the Two Point Boundar… ▽ More Maneuvering a spacecraft in the cislunar space is a complex problem, since it is highly perturbed by the gravitational influence of both the Earth and the Moon, and possibly also the Sun. Trajectories minimizing the needed fuel are generally preferred in order to decrease the mass of the payload. A classical method to constrain maneuvers is mathematically modelling them using the Two Point Boundary Value Problem (TPBVP), defining spacecraft positions at the start and end of the trajectory. Solutions to this problem can then be obtained with optimization techniques like the nonlinear least squares conjugated with the Theory of Functional Connections (TFC) to embed the constraints, which recently became an effective method for deducing orbit transfers. In this paper, we propose a tangential velocity (TV) type of constraints to design orbital maneuvers. We show that the technique presented in this paper can be used to transfer a spacecraft (e.g. from the Earth to the Moon) and perform rendezvous maneuvers (e.g. a swing-by with the Moon). In comparison with the TPBVP, solving the TV constraints via TFC offers several advantages, leading to a significant reduction in computational time. Hence, it proves to be an efficient technique to design these maneuvers. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: Submitted to Scientific Reports

Journal ref: Scientific Reports 14, 7479 (2024)

Adaptive Optics Telemetry Standard: Design and specification of a novel data exchange format

Authors: Tiago Gomes, Carlos M. Correia, Lisa Bardou, Sylvain Cetre, Johann Kolb, Caroline Kulcsár, François Leroux, Timothy Morris, Nuno Morujão, Benoît Neichel, Jean-Luc Beuzit, Paulo Garcia

Abstract: The amount of Adaptive Optics (AO) telemetry generated by VIS/NIR ground-based observatories is ever greater, leading to a growing need for a standardised data exchange format to support performance analysis and AO research and development activities that involve large-scale telemetry mining, processing, and curation. This paper introduces the Adaptive Optics Telemetry (AOT) data exchange format a… ▽ More The amount of Adaptive Optics (AO) telemetry generated by VIS/NIR ground-based observatories is ever greater, leading to a growing need for a standardised data exchange format to support performance analysis and AO research and development activities that involve large-scale telemetry mining, processing, and curation. This paper introduces the Adaptive Optics Telemetry (AOT) data exchange format as a standard for sharing AO telemetry from visible/infrared ground-based observatories. AOT is based on the Flexible Image Transport System (FITS) and aims to provide unambiguous and consistent data access across various systems and configurations, including natural and single/multiple laser guide-star AO systems. We designed AOT focused on two key use cases: atmospheric turbulence parameter estimation and point-spread function reconstruction (PSF-R). We prototyped and tested the design using existing AO telemetry datasets from multiple systems: single conjugate with natural and laser guide stars, tomographic systems with multi-channel wavefront sensors, single and multi wavefront correctors in systems featuring either a Shack-Hartmann or Pyramid as main wavefront sensors. The AOT file structure has been thoroughly defined, specifying data fields, descriptions, data types, units, and expected dimensions. To support this format, we have developed a Python package that enables data conversion, reading, writing and exploration of AOT files, which has been made publicly available and compatible with a general-purpose Python package manager. We demonstrate the flexibility of the AOT format by packaging data from five different instruments, installed on different telescopes. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Spin and orbital dynamics of planets undergoing thermal atmospheric tides using a vectorial approach

Authors: Ema F. S. Valente, Alexandre C. M. Correia

Abstract: Earth-mass planets are expected to have atmospheres and experience thermal tides raised by the host star. These tides transfer energy to the planet that can counter the dissipation from bodily tides. Indeed, even a relatively thin atmosphere can drive the rotation of these planets away from the synchronous state. Here we revisit the dynamical evolution of planets undergoing thermal atmospheric tid… ▽ More Earth-mass planets are expected to have atmospheres and experience thermal tides raised by the host star. These tides transfer energy to the planet that can counter the dissipation from bodily tides. Indeed, even a relatively thin atmosphere can drive the rotation of these planets away from the synchronous state. Here we revisit the dynamical evolution of planets undergoing thermal atmospheric tides. We use a novel approach based on a vectorial formalism, which is frame independent and valid for any configuration of the system, including any eccentricity and obliquity values. We provide the secular equations of motion after averaging over the mean anomaly and the argument of the pericenter, which are suitable to model the long-term spin and orbital evolution of the planet. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 15 pages, 1 figure, 2 tables

Journal ref: Astronomy and Astrophysics, vol. 679, A153 (2023)

Simulating METIS SCAO System

Authors: Markus Feldt, Horst Steuer, Carlos Correia, Andreas Obereder, Stefan Raffetseder, Thomas Bertram, Julia Shatokina, Faustine Cantalloube

Abstract: METIS, the Mid-Infrared ELT Imager and Spectrograph, is one of the four first-generation ELT instruments scheduled to see first light in 2028. Its two main science modules are supported by an adaptive optics system featuring a pyramid sensor with 90x90 subapertures working in the H and K bands. During the PDR and FDR phases, extensive simulations were carried out to support the sensing, reconstruc… ▽ More METIS, the Mid-Infrared ELT Imager and Spectrograph, is one of the four first-generation ELT instruments scheduled to see first light in 2028. Its two main science modules are supported by an adaptive optics system featuring a pyramid sensor with 90x90 subapertures working in the H and K bands. During the PDR and FDR phases, extensive simulations were carried out to support the sensing, reconstruction, and control concept of METIS single-conjugate adaptive optics (SCAO) system. We present details on the implementation of the COMPASS-based environment used for the simulations, the metrics used for analyzing our performance expectations, an overview of the main results, and some details on special cases like non-common path aberrations (NCPA) and water vapor seeing, as well as the low-wind effect. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 17 pages, 14 Figures, AO4ELT VII Conference - Avignon - June 2023

A roadmap to the efficient and robust characterization of temperate terrestrial planet atmospheres with JWST

Authors: TRAPPIST-1 JWST Community Initiative, :, Julien de Wit, René Doyon, Benjamin V. Rackham, Olivia Lim, Elsa Ducrot, Laura Kreidberg, Björn Benneke, Ignasi Ribas, David Berardo, Prajwal Niraula, Aishwarya Iyer, Alexander Shapiro, Nadiia Kostogryz, Veronika Witzke, Michaël Gillon, Eric Agol, Victoria Meadows, Adam J. Burgasser, James E. Owen, Jonathan J. Fortney, Franck Selsis, Aaron Bello-Arufe, Emeline Bolmont , et al. (57 additional authors not shown)

Abstract: Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have… ▽ More Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to constrain stellar heterogeneities and gather up to 2$\times$ more transits per JWST hour spent. We conclude that in such systems planets cannot be studied in isolation by small programs, thus large-scale community-supported programs should be supported to enable the efficient and robust exploration of terrestrial exoplanets in the JWST era. △ Less

Submitted 24 October, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

The Adaptive Optics System for the Gemini Infrared Multi-Object Spectrograph: Performance Modeling

Authors: Uriel Conod, Kate Jackson, Paolo Turri, Scott Chapman, Olivier Lardière, Masen Lamb, Carlos Correia, Gaetano Sivo, Suresh Sivanandam, Jean-Pierre Véran

Abstract: The Gemini Infrared Multi-Object Spectrograph (GIRMOS) will be a near-infrared, multi-object, medium spectral resolution, integral field spectrograph (IFS) for Gemini North Telescope, designed to operate behind the future Gemini North Adaptive Optics system (GNAO). In addition to a first ground layer Adaptive Optics (AO) correction in closed loop carried out by GNAO, each of the four GIRMOS IFSs w… ▽ More The Gemini Infrared Multi-Object Spectrograph (GIRMOS) will be a near-infrared, multi-object, medium spectral resolution, integral field spectrograph (IFS) for Gemini North Telescope, designed to operate behind the future Gemini North Adaptive Optics system (GNAO). In addition to a first ground layer Adaptive Optics (AO) correction in closed loop carried out by GNAO, each of the four GIRMOS IFSs will independently perform additional multi-object AO correction in open loop, resulting in an improved image quality that is critical to achieve top level science requirements. We present the baseline parameters and simulated performance of GIRMOS obtained by modeling both the GNAO and GIRMOS AO systems. The image quality requirement for GIRMOS is that 57% of the energy of an unresolved point-spread function ensquared within a 0.1 x 0.1 arcsecond at 2.0 μ m. It was established that GIRMOS will be an order 16 x 16 adaptive optics (AO) system after examining the tradeoffs between performance, risks and costs. The ensquared energy requirement will be met in median atmospheric conditions at Maunakea at 30° from zenith. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 13 pages, 10 figures, Publications of the Astronomical Society of the Pacific

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 135, Number 1052, 2023

Stability maps for the 5/3 mean motion resonance between Ariel and Umbriel with inclination

Authors: Sérgio R. A. Gomes, Alexandre C. M. Correia

Abstract: The evolution of the five largest satellites of Uranus during the crossing of the 5/3 mean motion resonance between Ariel and Umbriel is strongly affected by chaotic motion. Studies with numerical integrations of the equations of motion and analysis of Poincaré surface sections provided helpful insights to the role of chaos on the system. However, they lack of a quantification of this chaos in the… ▽ More The evolution of the five largest satellites of Uranus during the crossing of the 5/3 mean motion resonance between Ariel and Umbriel is strongly affected by chaotic motion. Studies with numerical integrations of the equations of motion and analysis of Poincaré surface sections provided helpful insights to the role of chaos on the system. However, they lack of a quantification of this chaos in the phase-space. Here, we construct stability maps using the frequency analysis method. We determine that for low energies (small eccentricity and/or inclinations), the phase-space is mainly stable. As the energy increases, the chaotic regions replace the stable motion, until only small, localized libration regions remain stable. △ Less

Submitted 9 September, 2023; originally announced September 2023.

Comments: arXiv admin note: text overlap with arXiv:2305.08794

Dynamics of trans-Neptunian objects near the 3/1 mean-motion resonance with Neptune

Authors: Alan J. Alves-Carmo, Timothée Vaillant, Alexandre C. M. Correia

Abstract: The complex classification of trans-Neptunian objects (TNOs) that are captured in mean-motion resonances (MMRs) and the constraint of their multiple origins are two significant open problems concerning the Solar System. The case-by-case study of the different MMRs and their characteristics provide information about their origin and dynamics, which helps us to understand the early stages of the Sol… ▽ More The complex classification of trans-Neptunian objects (TNOs) that are captured in mean-motion resonances (MMRs) and the constraint of their multiple origins are two significant open problems concerning the Solar System. The case-by-case study of the different MMRs and their characteristics provide information about their origin and dynamics, which helps us to understand the early stages of the Solar System evolution. In this paper, we study the dynamics of the detected TNOs close to a 3/1 MMR with Neptune. We initially use a semi-analytic three-body model to investigate the coplanar secular dynamics of these objects and find the stationary points. We then use surface sections and stability maps to analyse the non-averaged dynamics. These methods allow us to isolate the different stability regions and determine the extent of the chaotic regions. We show that stability maps are an extremely powerful tool for studying the resonant dynamics when they are computed in terms of the resonant angle. We then use these maps to study the non-planar three-body problem and the full dynamics in the presence of planetary perturbations. We confirm that TNOs near the 3/1 MMR regions can exist at very high inclinations. In the framework of the three-body problem, many of these objects can also be stable outside the 3/1 MMR owing to a Kozai secular resonance. However, when we take into account the perturbations of the four giant planets, the Kozai regions disappear and only the 3/1 MMR region remains, with eccentricities $e \lesssim 0.5$. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: 17 pages, 17 figures, 1 table

Journal ref: A&A 677, A83 (2023)

TOI-4860 b, a short-period giant planet transiting an M3.5 dwarf

Authors: J. M. Almenara, X. Bonfils, E. M. Bryant, A. Jordán, G. Hébrard, E. Martioli, A. C. M. Correia, N. Astudillo-Defru, C. Cadieux, L. Arnold, É. Artigau, G. Á. Bakos, S. C. C. Barros, D. Bayliss, F. Bouchy, G. Boué, R. Brahm, A. Carmona, D. Charbonneau, D. R. Ciardi, R. Cloutier, M. Cointepas, N. J. Cook, N. B. Cowan, X. Delfosse , et al. (25 additional authors not shown)

Abstract: We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03… ▽ More We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this. △ Less

Submitted 12 January, 2024; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: 16 pages, 14 figures, accepted for publication in A&A

Integrated turbulence parameters' estimation from NAOMI adaptive optics telemetry data

Authors: Nuno Morujão, Carlos Correia, Paulo Andrade, Julien Woillez, Paulo Garcia

Abstract: Monitoring turbulence parameters is crucial in high-angular resolution astronomy for various purposes, such as optimising adaptive optics systems or fringe trackers. The former are present at most modern observatories and will remain significant in the future. This makes them a valuable complementary tool for the estimation of turbulence parameters. The feasibility of estimating turbulence param… ▽ More Monitoring turbulence parameters is crucial in high-angular resolution astronomy for various purposes, such as optimising adaptive optics systems or fringe trackers. The former are present at most modern observatories and will remain significant in the future. This makes them a valuable complementary tool for the estimation of turbulence parameters. The feasibility of estimating turbulence parameters from low-resolution sensors remains untested. We perform seeing estimates for both simulated and on-sky telemetry data sourced from the new adaptive optics module installed on the four Auxiliary Telescopes of the Very Large Telescope Interferometer. The seeing estimates are obtained from a modified and optimised algorithm that employs a chi-squared modal fitting approach to the theoretical von Kármán model variances. The algorithm is built to retrieve turbulence parameters while simultaneously estimating and accounting for the remaining and measurement error. A Monte Carlo method is proposed for the estimation of the statistical uncertainty of the algorithm. The algorithm is shown to be able to achieve per cent accuracy in the estimation of the seeing with a temporal horizon of 20s on simulated data. A 0.76" +/- 1.2%$|_\mathrm{stat}$ +/- 1.2%$|_\mathrm{sys}$ median seeing was estimated from on-sky data collected from 2018 and 2020. The spatial distribution of the Auxiliary Telescopes across the Paranal Observatory was found to not play a role in the value of the seeing. △ Less

Submitted 5 October, 2023; v1 submitted 27 July, 2023; originally announced July 2023.

Journal ref: A&A 678, A193 (2023)

Theory of Functional Connections and Nelder-Mead optimization methods applied in satellite characterization

Authors: Allan Kardec de Almeida Junior, Safwan Aljbaae, Timothée Vaillant, Jhonathan M. Piñeros, Bruno Coelho, Domingos Barbosa, Miguel Bergano, João Pandeirada, Francisco C. Carvalho, Leonardo B. T. Santos, Antonio F. B. A. Prado, Anna Guerman, Alexandre C. M. Correia

Abstract: The growing population of man-made objects with the build up of mega-constellations not only increases the potential danger to all space vehicles and in-space infrastructures (including space observatories), but above all poses a serious threat to astronomy and dark skies. Monitoring of this population requires precise satellite characterization, which is is a challenging task that involves analyz… ▽ More The growing population of man-made objects with the build up of mega-constellations not only increases the potential danger to all space vehicles and in-space infrastructures (including space observatories), but above all poses a serious threat to astronomy and dark skies. Monitoring of this population requires precise satellite characterization, which is is a challenging task that involves analyzing observational data such as position, velocity, and light curves using optimization methods. In this study, we propose and analyze the application of two optimization procedures to determine the parameters associated with the dynamics of a satellite: one based on the Theory of Functional Connections (TFC) and another one based on the Nelder-Mead heuristic optimization algorithm. The TFC performs linear functional interpolation to embed the constraints of the problem into a functional. In this paper, we propose to use this functional to analytically embed the observational data of a satellite into its equations of dynamics. After that, any solution will always satisfy the observational data. The second procedure proposed in this research takes advantage of the Nealder-Mead algorithm, that does not require the gradient of the objective function, as alternative solution. The accuracy, efficiency, and dependency on the initial guess of each method is investigated, analyzed, and compared for several dynamical models. These methods can be used to obtain the physical parameters of a satellite from available observational data and for space debris characterization contributing to follow-up monitoring activities in space and astronomical observatories. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: Submitted to Acta Astronautica

Tidal excitation of the obliquity of Earth-like planets in the habitable zone of M-dwarf stars

Authors: Ema F. S. Valente, Alexandre C. M. Correia

Abstract: Close-in planets undergo strong tidal interactions with the parent star that modify their spins and orbits. In the two-body problem, the final stage for tidal evolution is the synchronisation of the rotation and orbital periods, and the alignment of the planet spin axis with the normal to the orbit (zero planet obliquity). The orbital eccentricity is also damped to zero, but over a much longer tim… ▽ More Close-in planets undergo strong tidal interactions with the parent star that modify their spins and orbits. In the two-body problem, the final stage for tidal evolution is the synchronisation of the rotation and orbital periods, and the alignment of the planet spin axis with the normal to the orbit (zero planet obliquity). The orbital eccentricity is also damped to zero, but over a much longer timescale, that may exceed the lifetime of the system. For non-zero eccentricities, the rotation rate can be trapped in spin-orbit resonances that delay the evolution towards the synchronous state. Here we show that capture in some spin-orbit resonances may also excite the obliquity to high values rather than damp it to zero. Depending on the system parameters, obliquities of 60 to 80 degrees can be maintained throughout the entire lifetime of the planet. This unexpected behaviour is particularly important for Earth-like planets in the habitable zone of M-dwarf stars, as it may help to sustain temperate environments and thus more favourable conditions for life. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 12 pages, 9 figures, 2 tables

Journal ref: Astronomy and Astrophysics, vol. 665, A130 (2022)

Tidal evolution for any rheological model using a vectorial approach expressed in Hansen coefficients

Authors: Alexandre C. M. Correia, Ema F. S. Valente

Abstract: We revisit the two body problem, where one body can be deformed under the action of tides raised by the companion. Tidal deformation and consequent dissipation result in spin and orbital evolution of the system. In general, the equations of motion are derived from the tidal potential developed in Fourier series expressed in terms of Keplerian elliptical elements, so that the variation of dissipati… ▽ More We revisit the two body problem, where one body can be deformed under the action of tides raised by the companion. Tidal deformation and consequent dissipation result in spin and orbital evolution of the system. In general, the equations of motion are derived from the tidal potential developed in Fourier series expressed in terms of Keplerian elliptical elements, so that the variation of dissipation with amplitude and frequency can be examined. However, this method introduces multiple index summations and some orbital elements depend on the chosen frame, which is prone to confusion and errors. Here, we develop the quadrupole tidal potential solely in a series of Hansen coefficients, which are widely used in celestial mechanics and depend just on the eccentricity. We derive the secular equations of motion in a vectorial formalism, which is frame independent and valid for any rheological model. We provide expressions for a single average over the mean anomaly and for an additional average over the argument of the pericentre. These equations are suitable to model the long-term evolution of a large variety of systems and configurations, from planet satellite to stellar binaries. We also compute the tidal energy released inside the body for an arbitrary configuration of the system. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 33 pages, 1 figure, 2 tables

Journal ref: Celestial Mechanics and Dynamical Astronomy, vol. 134, article 24 (2022)

Effect of the inclination in the passage through the 5/3 mean motion resonance between Ariel and Umbriel

Authors: Sérgio R. A. Gomes, Alexandre C. M. Correia

Abstract: The orbits of the main satellites of Uranus are expected to slowly drift away owing to tides raised in the planet. As a result, the 5/3 mean motion resonance between Ariel and Umbriel was likely encountered in the past. Previous studies have shown that, in order to prevent entrapment in this resonance, the eccentricities of the satellites must be larger than $\sim 0.01$ at the epoch, which is hard… ▽ More The orbits of the main satellites of Uranus are expected to slowly drift away owing to tides raised in the planet. As a result, the 5/3 mean motion resonance between Ariel and Umbriel was likely encountered in the past. Previous studies have shown that, in order to prevent entrapment in this resonance, the eccentricities of the satellites must be larger than $\sim 0.01$ at the epoch, which is hard to explain. On the other hand, if the satellites experience some temporary capture and then escape, the inclinations rise to high values that are not observed today. We have revisited this problem both analytically and numerically focussing on the inclination, using a secular two-satellite model with circular orbits. We show that if the inclination of Umbriel was around $0.15^{\circ}$ at the time of the 5/3 resonance encounter, capture can be avoided in about $60\%$ of the cases. Moreover, after the resonance crossing, the inclination of Umbriel drops to a mean value around $0.08^{\circ}$, which is close to the presently observed one. The final inclination of Ariel is distributed between $0.01^{\circ}$ and $0.25^{\circ}$ with a nearly equal probability, which includes the present mean value of $0.02^{\circ}$. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: 19 pages, 10 figures

Journal ref: A&A 674, A111 (2023)

The magnetic field and multiple planets of the young dwarf AU~Mic

Authors: J. -F. Donati, P. I. Cristofari, B. Finociety, B. Klein, C. Moutou, E. Gaidos, C. Cadieux, E. Artigau, A. C. M. Correia, G. Boué, N. J. Cook, A. Carmona, L. T. Lehmann, J. Bouvier, E. Martioli, J. Morin, P. Fouqué, X. Delfosse, R. Royon, G. Hébrard, S. H. P. Alencar, J. Laskar, L. Arnold, P. Petit, A. Kospal , et al. (3 additional authors not shown)

Abstract: In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures o… ▽ More In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$σ$ while those of b and c show up at $\simeq$4$σ$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e. △ Less

Submitted 24 April, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS, in press (20 pages and 12 figures + 9 pages of supplementary material)

The ESO's Extremely Large Telescope Working Groups

Authors: Paolo Padovani, Michele Cirasuolo, Remco van der Burg, Faustine Cantalloube, Elizabeth George, Markus Kasper, Kieran Leschinski, Carlos Martins, Julien Milli, Sabine Möhler, Mark Neeser, Benoit Neichel, Angel Otarola, Rubén Sánchez-Janssen, Benoit Serra, Alain Smette, Elena Valenti, Christophe Verinaud, Joël Vernet, Olivier Absil, Guido Agapito, Morten Andersen, Carmelo Arcidiacono, Matej Arko, Pierre Baudoz , et al. (60 additional authors not shown)

Abstract: Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations an… ▽ More Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations and calibrations of the telescope, thirteen specific working groups were created to coordinate the effort between ESO, the instrument consortia, and the wider community. We describe here the goals of these working groups as well as their achievements so far. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 8 pages, 4 figures

Journal ref: The Messenger, vol. 189, (December 2022) p. 23-30

Improving circumbinary planet detections by fitting their binary's apsidal precession

Authors: Thomas A. Baycroft, Amaury H. M. J. Triaud, João Faria, Alexandre C. M. Correia, Matthew R. Standing

Abstract: Apsidal precession in stellar binaries is the main non-Keplerian dynamical effect impacting the radial-velocities of a binary star system. Its presence can notably hide the presence of orbiting circumbinary planets because many fitting algorithms assume perfectly Keplerian motion. To first order, apsidal precession ($\dotω$) can be accounted for by adding a linear term to the usual Keplerian model… ▽ More Apsidal precession in stellar binaries is the main non-Keplerian dynamical effect impacting the radial-velocities of a binary star system. Its presence can notably hide the presence of orbiting circumbinary planets because many fitting algorithms assume perfectly Keplerian motion. To first order, apsidal precession ($\dotω$) can be accounted for by adding a linear term to the usual Keplerian model. We include apsidal precession in the kima package, an orbital fitter designed to detect and characterise planets from radial velocity data. In this paper, we detail this and other additions to kima that improve fitting for stellar binaries and circumbinary planets including corrections from general relativity. We then demonstrate that fitting for $\dotω$ can improve the detection sensitivity to circumbinary exoplanets by up to an order of magnitude in some circumstances, particularly in the case of multi-planetary systems. In addition, we apply the algorithm to several real systems, producing a new measurement of aspidal precession in KOI-126 (a tight triple system), and a detection of $\dotω$ in the Kepler-16 circumbinary system. Although apsidal precession is detected for Kepler-16, it does not have a large effect on the detection limit or the planetary parameters. We also derive an expression for the precession an outer planet would induce on the inner binary and compare the value this predicts with the one we detect. △ Less

Submitted 6 April, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: 10 pages, 5 figures Re-submitted to MNRAS after reviewer comments

Radial-velocity discovery of a second planet in the TOI-1338/BEBOP-1 circumbinary system

Authors: Matthew R. Standing, Lalitha Sairam, David V. Martin, Amaury H. M. J. Triaud, Alexandre C. M. Correia, Gavin A. L. Coleman, Thomas A. Baycroft, Vedad Kunovac, Isabelle Boisse, Andrew Collier Cameron, Georgina Dransfield, João P. Faria, Michaël Gillon, Nathan C. Hara, Coel Hellier, Jonathan Howard, Ellie Lane, Rosemary Mardling, Pierre F. L. Maxted, Nicola J. Miller, Richard P. Nelson, Jerome A. Orosz, Franscesco Pepe, Alexandre Santerne, Daniel Sebastian , et al. (2 additional authors not shown)

Abstract: We report the detection of a gas-giant planet in orbit around both stars of an eclipsing binary star system that also contains the smaller, inner transiting planet TOI-1338b. The new planet, called TOI-1338/BEBOP-1c, was discovered using radial-velocity data collected with the HARPS and ESPRESSO spectrographs. Our analysis reveals it is a $65.2~\rm{M_{\oplus}}$ circumbinary planet with a period of… ▽ More We report the detection of a gas-giant planet in orbit around both stars of an eclipsing binary star system that also contains the smaller, inner transiting planet TOI-1338b. The new planet, called TOI-1338/BEBOP-1c, was discovered using radial-velocity data collected with the HARPS and ESPRESSO spectrographs. Our analysis reveals it is a $65.2~\rm{M_{\oplus}}$ circumbinary planet with a period of $215.5~$days. This is the first detection of a circumbinary planet using radial-velocity observations alone, and makes TOI-1338/BEBOP-1 only the second confirmed multiplanet circumbinary system to date. We do not detect the smaller inner transiting planet with radial-velocity data, and can place an upper limit on the inner planet's mass at $21.8~\mathrm{M}_\oplus$ with $99\%$ confidence. The inner planet is the first circumbinary planet amenable for atmospheric characterisation, using the James Webb Space Telescope. △ Less

Submitted 12 June, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-023-01948-4 Main: 14-pages, 4-Figures, 1-table. Supporting material: 48-pages, 12-figures, 4-tables

Reversible time-step adaptation for the integration of few-body systems

Authors: Tjarda C. N. Boekholt, Timothee Vaillant, Alexandre C. M. Correia

Abstract: The time step criterion plays a crucial role in direct N-body codes. If not chosen carefully, it will cause a secular drift in the energy error. Shared, adaptive time step criteria commonly adopt the minimum pairwise time step, which suffers from discontinuities in the time evolution of the time step. This has a large impact on the functioning of time step symmetrisation algorithms. We provide new… ▽ More The time step criterion plays a crucial role in direct N-body codes. If not chosen carefully, it will cause a secular drift in the energy error. Shared, adaptive time step criteria commonly adopt the minimum pairwise time step, which suffers from discontinuities in the time evolution of the time step. This has a large impact on the functioning of time step symmetrisation algorithms. We provide new demonstrations of previous findings that a smooth and weighted average over all pairwise time steps in the N-body system, improves the level of energy conservation. Furthermore, we compare the performance of 27 different time step criteria, by considering 3 methods for weighting time steps and 9 symmetrisation methods. We present performance tests for strongly chaotic few-body systems, including unstable triples, giant planets in a resonant chain, and the current Solar System. We find that the harmonic symmetrisation methods (methods A3 and B3 in our notation) are the most robust, in the sense that the symmetrised time step remains close to the time step function. Furthermore, based on our Solar System experiment, we find that our new weighting method based on direct pairwise averaging (method W2 in our notation), is slightly preferred over the other methods. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted by MNRAS. 13 pages, 6 figures

Developing a data fusion concept for radar and optical ground based SST station

Authors: Bruno Coelho, Domingos Barbosa, Miguel Bergano, João Pandeirada, Paulo Marques, Alexandre C. M. Correia, José Matias de Freitas

Abstract: As part of the Portuguese Space Surveillance and Tracking (SST) program, a tracking radar and a double Wide Field of View Telescope system (4.3° x 2.3°) are being installed at the Pampilhosa da Serra Space Observatory (PASO) in the centre of continental Portugal, complementing an already installed deployable optical sensor for MEO and GEO surveillance. The tracking radar will track space debris in… ▽ More As part of the Portuguese Space Surveillance and Tracking (SST) program, a tracking radar and a double Wide Field of View Telescope system (4.3° x 2.3°) are being installed at the Pampilhosa da Serra Space Observatory (PASO) in the centre of continental Portugal, complementing an already installed deployable optical sensor for MEO and GEO surveillance. The tracking radar will track space debris in Low Earth Orbit (LEO) up to 1000 km and at the same time the telescope will also have LEO tracking capabilities. This article intends to discuss possible ways to take advantage of having these two sensors at the same location. Using both types of sensors takes advantage of the radar measurements which give precise radial velocity and distance to the objects, while the telescope gives better sky coordinates measurements. With the installation of radar and optical sensors, PASO can extend observation time of space debris and correlate information from optical and radar provenances in real time. During twilight periods both sensors can be used simultaneously to rapidly compute new TLEs for LEO objects, eliminating the time delays involved in data exchange between sites in a large SST network. This concept will not replace the need for a SST network with sensors in multiple locations around the globe, but will provide a more complete set of measurements from a given object passage, and therefore increase the added value for initial orbit determination, or monitoring of reentry campaigns of a given location. PASO will contribute to the development of new solutions to better characterize the objects improving the overall SST capabilities and constitute a perfect site for the development and testing of new radar and optical data fusion algorithms and techniques for space debris monitoring. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 3 pages, 3 figures, Draft accepted for publication at 73rd International Astronautical Congress, Paris 2022

Chasing extreme planetary architectures: I- HD196885Ab, a super-Jupiter dancing with two stars?

Authors: G. Chauvin, M. Videla, H. Beust, R. Mendez, A. C. M. Correia, S. Lacour, A. Tokovinin, J. Hagelberg, F. Bouchy, I. Boisse, C. Villegas, M. Bonavita, S. Desidera, V. Faramaz, T. Forveille, A. Gallenne, X. Haubois, J. S. Jenkins, P. Kervella, A. -M. Lagrange, C. Melo, P. Thebault, S. Udry, D. Segransan

Abstract: Planet(s) in binaries are unique architectures for testing predictions of planetary formation and evolution theories in very hostile environments. We used the IRDIS dual-band imager of SPHERE at VLT, and the speckle interferometric camera HRCAM of SOAR, to acquire high-angular resolution images of HD 196885 AB between 2015 and 2020. Radial velocity observations have been extended over almost 40 yr… ▽ More Planet(s) in binaries are unique architectures for testing predictions of planetary formation and evolution theories in very hostile environments. We used the IRDIS dual-band imager of SPHERE at VLT, and the speckle interferometric camera HRCAM of SOAR, to acquire high-angular resolution images of HD 196885 AB between 2015 and 2020. Radial velocity observations have been extended over almost 40 yr extending the radial velocity measurements HD 196885 A and resolving both the binary companion and the inner giant planet HD 196885 Ab. Finally, we took advantage of the exquisite astrometric precision of the dual-field mode of VLTI/GRAVITY (down to 30 μas) to monitor the relative position of HD 196885 A and B to search for the 3.6 yr astrometric wobble of the circumprimary planet Ab imprinted on the binary separation. Our observations enable to accurately constrain the orbital properties of the binary HD 196885 AB, seen on an inclined and retrograde orbit (iAB = 120.43 deg) with a semi-major axis of 19.78 au, and an eccentricity of 0.417. The GRAVITY measurements confirm for the first time the nature of the inner planet HD 196885 Ab by rejecting all families of pole-on solutions in the stellar or brown dwarf masses. The most favored island of solutions is associated with a Jupiter-like planet (MAb = 3.39 MJup), with moderate eccentricity (eAaAb = 0.44), and inclination close to 143.04 deg. This results points toward a significant mutual inclination (Phi = 24.36 deg) between the orbital planes (relative to the star) of the binary companion B and the planet Ab. Our dynamical simulations indicate that the system is dynamically stable over time. Eccentricity and mutual inclination variations could be expected for moderate von Zipele Kozai Lidov cycles that may affect the inner planet. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 12 pages, 6 figures, accepted in A&A

Journal ref: A&A 675, A114 (2023)

Three-sided pyramid wavefront sensor. II. Preliminary demonstration on the new CACTI testbed

Authors: Lauren Schatz, Johanan Codona, Joseph D. Long, Jared R. Males, Weslin Pullen, Jennifer Lumbres, Kyle Van Gorkom, Vincent Chambouleyron, Laird M. Close, Carlos Correia, Olivier Fauvarque, Thierry Fusco, Olivier Guyon, Michael Hart, Pierre Janin-Potiron, Robert Johnson, Nemanja Jovanovic, Mala Mateen, Jean-Francois Sauvage, Benoit Neichel

Abstract: The next generation of giant ground and space telescopes will have the light-collecting power to detect and characterize potentially habitable terrestrial exoplanets using high-contrast imaging for the first time. This will only be achievable if the performance of Giant Segmented Mirror Telescopes (GSMTs) extreme adaptive optics (ExAO) systems are optimized to their full potential. A key component… ▽ More The next generation of giant ground and space telescopes will have the light-collecting power to detect and characterize potentially habitable terrestrial exoplanets using high-contrast imaging for the first time. This will only be achievable if the performance of Giant Segmented Mirror Telescopes (GSMTs) extreme adaptive optics (ExAO) systems are optimized to their full potential. A key component of an ExAO system is the wavefront sensor (WFS), which measures aberrations from atmospheric turbulence. A common choice in current and next-generation instruments is the pyramid wavefront sensor (PWFS). ExAO systems require high spatial and temporal sampling of wavefronts to optimize performance, and as a result, require large detectors for the WFS. We present a closed-loop testbed demonstration of a three-sided pyramid wavefront sensor (3PWFS) as an alternative to the conventional four-sided pyramid wavefront (4PWFS) sensor for GSMT-ExAO applications on the new Comprehensive Adaptive Optics and Coronagraph Test Instrument (CACTI). The 3PWFS is less sensitive to read noise than the 4PWFS because it uses fewer detector pixels. The 3PWFS has further benefits: a high-quality three-sided pyramid optic is easier to manufacture than a four-sided pyramid. We detail the design of the two components of the CACTI system, the adaptive optics simulator and the PWFS testbed that includes both a 3PWFS and 4PWFS. A preliminary experiment was performed on CACTI to study the performance of the 3PWFS to the 4PWFS in varying strengths of turbulence using both the Raw Intensity and Slopes Map signal processing methods. This experiment was repeated for a modulation radius of 1.6 lambda/D and 3.25 lambda/D. We found that the performance of the two wavefront sensors is comparable if modal loop gains are tuned. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 28 Pages, 15 Figures, and 4 Tables

A direct N-body integrator for modelling the chaotic, tidal dynamics of multi-body extrasolar systems: TIDYMESS

Authors: Tjarda C. N. Boekholt, Alexandre C. M. Correia

Abstract: Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop ne… ▽ More Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop new N-body codes, which allow for experimentation with various tidal models and numerical implementations. We present the open-source N-body code TIDYMESS, which stands for ``TIdal DYnamics of Multi-body ExtraSolar Systems''. This code implements a creep deformation law for the bodies, parametrized by their fluid Love numbers and fluid relaxation times. Due to tidal and centrifugal deformations, we approximate the general shape of a body to be an ellipsoid. We calculate the associated gravitational field to quadruple order, from which we derive the gravitational accelerations and torques. The equations of motion for the orbits, spins and deformations are integrated directly using a fourth-order integration method based on a symplectic composition. We implement a novel integration method for the deformations, which allows for a time step solely dependent on the orbits, and not on the spin periods or fluid relaxation times. This feature greatly speeds up the calculations, while also improving the consistency when comparing different tidal regimes. We demonstrate the capabilities and performance of TIDYMESS, particularly in the niche regime of parameter space where orbits are chaotic and tides become non-linear. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: 17 pages, 6 figures

Super-resolution wavefront reconstruction

Authors: Sylvain Oberti, Carlos Correia, Thierry Fusco, Benoit Neichel, Pierre Guiraud

Abstract: Super-Resolution (SR) is a technique that seeks to upscale the resolution of a set of measured signals. SR retrieves higher-frequency signal content by combining multiple lower resolution sampled data sets. SR is well known both in the temporal and spatial domains. It is widely used in imaging to reduce aliasing and enhance the resolution of coarsely sampled images.This paper applies the SR techni… ▽ More Super-Resolution (SR) is a technique that seeks to upscale the resolution of a set of measured signals. SR retrieves higher-frequency signal content by combining multiple lower resolution sampled data sets. SR is well known both in the temporal and spatial domains. It is widely used in imaging to reduce aliasing and enhance the resolution of coarsely sampled images.This paper applies the SR technique to the bi-dimensional wavefront reconstruction. In particular, we show how SR is intrinsically suited for tomographic multi WaveFront Sensor (WFS) AO systems revealing many of its advantages with minimal design effort. This paper provides a direct space and Fourier-optics description of the wavefront sensing operation and demonstrate how SR can be exploited through signal reconstruction, especially in the framework of Periodic Nonuniform Sampling. Both meta uniform and nonuniform sampling schemes are investigated. Then, the SR bi-dimensional model for a Shack Hartmann (SH) WFS is provided and the characteristics of the sensitivity function are analyzed. The SR concept is finally validated with numerical simulations of representative multi WFS SH AO systems. Our results show that combining several WFS samples in a SR framework grants access to a larger number of modes than the native one offered by a single WFS and that despite the fixed sub-aperture size across samples. Furthermore, we show that the associated noise propagation is not degraded under SR. Finally, the concept is extended to the signal produced by single Pyramid WFS. In conclusion, SR applied to wavefront reconstruction offers a new parameter space to explore as it decouples the size of the subaperture from the desired wavefront sampling resolution. By cutting short with old assumptions, new, more flexible and better performing AO designs become now possible. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: Astronomy & Astrophysics Section: 13. Astronomical instrumentation AA/2022/43954

Journal ref: A&A 667, A48 (2022)

Multitension strings in high-resolution U(1)$\times$U(1) simulations

Authors: J. R. C. C. C. Correia, C. J. A. P. Martins

Abstract: Topological defects are a fossil relic of early Universe phase transitions, with cosmic strings being the best motivated example. While in most cases one studies Nambu-Goto or Abelian-Higgs strings, one also expects that cosmologically realistic strings should have additional degrees of freedom in their worldsheets, one specific example being superstrings from Type IIB superstring theory. Here we… ▽ More Topological defects are a fossil relic of early Universe phase transitions, with cosmic strings being the best motivated example. While in most cases one studies Nambu-Goto or Abelian-Higgs strings, one also expects that cosmologically realistic strings should have additional degrees of freedom in their worldsheets, one specific example being superstrings from Type IIB superstring theory. Here we continue the scientific exploitation of our recently developed multi-GPU field theory cosmic strings code to study the evolution of U(1)$\times$U(1) multitension networks, which are a numerically convenient proxy: these contain two lowest-tension strings networks able to interact and form bound states, providing a convenient first approximation to the behaviour expected from cosmic superstrings. (...) We rely on the largest field theory simulations of this model so far, specifically $4096^3$, $Δx = 0.5$ boxes. We present robust evidence of scaling for the lightest strings, measured through a complete and self-consistent set of correlation length and velocity diagnostics. We also find a linearly growing average length of the bound state segments, consistent with a scaling behaviour. (In previously reported lower-resolution simulations, such behaviour had only been identified with carefully engineered initial conditions, rich in those segments.) Finally, while we see no evidence of a large population of bound states forming at early stages of the network evolution, we do present tentative evidence for an asymptotic constant value of the fraction of bound states, with this value being different in the radiation and the matter eras. Our work demonstrates that our GPU-accelerated field theory code can by successfully extended beyond the simple Abelian-Higgs approximation, and enables future detailed studies of realistic string networks and of their observational signatures. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 23 pages, 11 figures, 5 tables; Phys. Rev. D (in press). Abstract abridged due to arXiv's 1920 character limit

Daytime calibration and testing of the Keck All sky Precision Adaptive Optics Tomography System

Authors: Avinash Surendran, Jacques R. Delorme, Carlos M. Correia, Steve Doyle, Sam Ragland, Paul Richards, Peter Wizinowich, Philip M. Hinz, Daren Dillon, Cesar Laguna, Sylvain Cetre, Scott Lilley, Ed Wetherell, Jason C. Y. Chin, Eduardo Marin

Abstract: The development of the Keck All sky Precision Adaptive optics (KAPA) project was initiated in September 2018 to upgrade the Keck I adaptive optics (AO) system to enable laser tomography adaptive optics (LTAO) with a four laser guide star (LGS) asterism. The project includes the replacement of the existing LMCT laser with a Toptica laser, the implementation of a new real-time controller (RTC) and w… ▽ More The development of the Keck All sky Precision Adaptive optics (KAPA) project was initiated in September 2018 to upgrade the Keck I adaptive optics (AO) system to enable laser tomography adaptive optics (LTAO) with a four laser guide star (LGS) asterism. The project includes the replacement of the existing LMCT laser with a Toptica laser, the implementation of a new real-time controller (RTC) and wavefront sensor optics and camera, and a new daytime calibration and test platform to provide the required infrastructure for laser tomography. The work presented here describes the new daytime calibration infrastructure to test the performance for the KAPA tomographic algorithms. This paper outlines the hardware infrastructure for daytime calibration and performance assessment of tomographic algorithms. This includes the implementation of an asterism simulator having fiber-coupled light sources simulating four Laser Guide Stars (LGS) and two Natural Guide Stars (NGS) at the AO bench focus, as well as the upgrade of the existing TelSim on the AO bench to simulate focal anisoplanatism and wind driven atmospheric turbulence. A phase screen, that can be adjusted in effective altitude, is used to simulate wind speeds up to 10 m/s for a duration of upto 3 s. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 15 pages, 13 figures. Publishd as part of SPIE Astronomical telescopes + Instrumentation 2022 proceedings

Towards virtual access of adaptive optics telemetry data

Authors: Tiago Gomes, Carlos Correia, Lisa Bardou, Olivier Beltramo-Martin, Thierry Fusco, Caroline Kulcsár, Timothy Morris, Nuno Morujão, Benoît Neichel, James Osborn, Paulo Garcia

Abstract: Large amounts of Adaptive-Optics (AO) control loop data and telemetry are currently inaccessible to end-users. Broadening access to those data has the potential to change the AO landscape on many fronts, addressing several use-cases such as derivation of the system's PSF, turbulence characterisation and optimisation of system control. We address one of the biggest obstacles to sharing these data:… ▽ More Large amounts of Adaptive-Optics (AO) control loop data and telemetry are currently inaccessible to end-users. Broadening access to those data has the potential to change the AO landscape on many fronts, addressing several use-cases such as derivation of the system's PSF, turbulence characterisation and optimisation of system control. We address one of the biggest obstacles to sharing these data: the lack of standardisation, which hinders access. We propose an object-oriented Python package for AO telemetry, whose data model abstracts the user from an underlining archive-ready data exchange standard based on the Flexible Image Transport System (FITS). Its design supports data from a wide range of existing and future AO systems, either in raw format or abstracted from actual instrument details. We exemplify its usage with data from active AO systems on 10m-class observatories, of which two are currently supported (AOF and Keck), with plans for more. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Key wavefront sensors features for laser-assisted tomographic adaptive optics systems on the Extremely Large Telescope

Authors: Thierry Fusco, Guido Agapito, Benoit Neichel, Sylvain Oberti, Carlos Correia, Pierre Haguenauer, Cédric Plantet, Felipe Pedreros, Zibo Ke, Anne Costille, Pierre Jouve, Lorenzo Busoni, Simone Esposito

Abstract: Laser guide star (LGS) wave-front sensing (LGSWFS) is a key element of tomographic adaptive optics system. However, when considering Extremely Large Telescope (ELT) scales, the LGS spot elongation becomes so large that it challenges the standard recipes to design LGSWFS. For classical Shack-Hartmann wave-front sensor (SHWFS), which is the current baseline for all ELT LGS-assisted instruments, a tr… ▽ More Laser guide star (LGS) wave-front sensing (LGSWFS) is a key element of tomographic adaptive optics system. However, when considering Extremely Large Telescope (ELT) scales, the LGS spot elongation becomes so large that it challenges the standard recipes to design LGSWFS. For classical Shack-Hartmann wave-front sensor (SHWFS), which is the current baseline for all ELT LGS-assisted instruments, a trade-off between the pupil spatial sampling [number of sub-apertures (SAs)], the SA field-of-view (FoV) and the pixel sampling within each SA is required. For ELT scales, this trade-off is also driven by strong technical constraints, especially concerning the available detectors and in particular their number of pixels. For SHWFS, a larger field of view per SA allows mitigating the LGS spot truncation, which represents a severe loss of performance due to measurement biases. For a given number of available detectors pixels, the SA FoV is competing with the proper sampling of the LGS spots, and/or the total number of SAs. We proposed a sensitivity analysis, and we explore how these parameters impacts the final performance. In particular, we introduce the concept of super resolution, which allows one to reduce the pupil sampling per WFS and opens an opportunity to propose potential LGSWFS designs providing the best performance for ELT scales. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 21 pages, 11 figures, 1 table, part of the 2022 JATIS Special Section on Extremely Large Telescopes

Journal ref: J. of Astronomical Telescopes, Instruments, and Systems, 8(2), 021514 (2022)

Dynamics of co-orbital exoplanets in a first order resonance chain with tidal dissipation

Authors: Jérémy Couturier, Philippe Robutel, Alexandre C. M. Correia

Abstract: Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time… ▽ More Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time-lag model, which extends the Hamiltonian formalism of the point-mass case. We show the existence of several families of equilibria, and how these equilibria extend to the complete system. In one family, which we call the main branch, a secular resonance between the libration frequency of the co-orbitals and the precession frequency of the pericentres has unexpected dynamical consequences when tidal dissipation is added. We report the existence of two distinct mechanisms that make co-orbital planets much more stable within the $p:p:p+1$ resonance chain rather than outside it. The first one is due to negative real parts of the eigenvalues of the linearised system with tides, in the region of the secular resonance mentioned above. The second one comes from non-linear contributions of the vector field and it is due to eccentricity damping. These two stabilising mechanisms increase the chances of a still-to-come detection of exoplanets in the co-orbital configuration. △ Less

Submitted 17 April, 2022; originally announced April 2022.

Journal ref: A&A 664, A1 (2022)

Photodynamical analysis of the nearly resonant planetary system WASP-148: Accurate transit-timing variations and mutual orbital inclination

Authors: J. M. Almenara, G. Hébrard, R. F. Díaz, J. Laskar, A. C. M. Correia, D. R. Anderson, I. Boisse, X. Bonfils, D. J. A. Brown, V. Casanova, A. Collier Cameron, M. Fernández, J. M. Jenkins, F. Kiefer, A. Lecavelier des Étangs, J. J Lissauer, G. Maciejewski, J. McCormac, H. Osborn, D. Pollacco, G. Ricker, J. Sánchez, S. Seager, S. Udry, D. Verilhac , et al. (1 additional authors not shown)

Abstract: WASP-148 is a recently announced extra-solar system harbouring at least two giant planets. The inner planet transits its host star. The planets travel on eccentric orbits and are near the 4:1 mean-motion resonance, which implies significant mutual gravitational interactions. In particular, this causes transit-timing variations of a few minutes, which were detected based on ground-based photometry.… ▽ More WASP-148 is a recently announced extra-solar system harbouring at least two giant planets. The inner planet transits its host star. The planets travel on eccentric orbits and are near the 4:1 mean-motion resonance, which implies significant mutual gravitational interactions. In particular, this causes transit-timing variations of a few minutes, which were detected based on ground-based photometry. This made WASP-148 one of the few cases where such a phenomenon was detected without space-based photometry. Here, we present a self-consistent model of WASP-148 that takes into account the gravitational interactions between all known bodies in the system. Our analysis simultaneously fits the available radial velocities and transit light curves. In particular, we used the photometry secured by the TESS space telescope and made public after the WASP-148 discovery announcement. The TESS data confirm the transit-timing variations, but only in combination with previously measured transit times. The system parameters we derived agree with those previously reported and have a significantly improved precision, including the mass of the non-transiting planet. We found a significant mutual inclination between the orbital planes of the two planets: I=41.0 +6.2 -7.6 deg based on the modelling of the observations, although we found I=20.8 +/- 4.6 deg when we imposed a constraint on the model enforcing long-term dynamical stability. When a third planet was added to the model - based on a candidate signal in the radial velocity - the mutual inclination between planets b and c changed significantly allowing solutions closer to coplanar. We conclude that more data are needed to establish the true architecture of the system. If the significant mutual inclination is confirmed, WASP-148 would become one of the only few candidate non-coplanar planetary systems. We discuss possible origins for this misalignment. △ Less

Submitted 16 September, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 26 pages, 20 figures, accepted for publication in A&A

Journal ref: A&A 663, A134 (2022)

Eviction-like resonances for satellite orbits. Application to Phobos, the main satellite of Mars

Authors: Timothée Vaillant, Alexandre C. M. Correia

Abstract: The motion of a satellite can experience secular resonances between the precession frequencies of its orbit and the mean motion of the host planet around the star. Some of these resonances can significantly modify the eccentricity (evection resonance) and the inclination (eviction resonance) of the satellite. In this paper, we study in detail the secular resonances that can disturb the orbit of a… ▽ More The motion of a satellite can experience secular resonances between the precession frequencies of its orbit and the mean motion of the host planet around the star. Some of these resonances can significantly modify the eccentricity (evection resonance) and the inclination (eviction resonance) of the satellite. In this paper, we study in detail the secular resonances that can disturb the orbit of a satellite, in particular the eviction-like ones. Although the inclination is always disturbed while crossing one eviction-like resonance, capture can only occur when the semi-major axis is decreasing. This is, for instance, the case of Phobos, the largest satellite of Mars, that will cross some of these resonances in the future because its orbit is shrinking owing to tidal effects. We estimate the impact of resonance crossing in the orbit of the satellite, including the capture probabilities, as a function of several parameters, such as the eccentricity and the inclination of the satellite, and the obliquity of the planet. Finally, we use the method of the frequency map analysis to study the resonant dynamics based on stability maps, and we show that some of the secular resonances may overlap, which leads to chaotic motion for the inclination of the satellite. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 657, A103 (2022)

Spin-orbit coupling for close-in planets

Authors: Alexandre C. M. Correia, Jean-Baptiste Delisle

Abstract: We study the spin evolution of close-in planets in multi-body systems and present a very general formulation of the spin-orbit problem. This includes a simple way to probe the spin dynamics from the orbital perturbations, a new method for computing forced librations and tidal deformation, and general expressions for the tidal torque and capture probabilities in resonance. We show that planet-plane… ▽ More We study the spin evolution of close-in planets in multi-body systems and present a very general formulation of the spin-orbit problem. This includes a simple way to probe the spin dynamics from the orbital perturbations, a new method for computing forced librations and tidal deformation, and general expressions for the tidal torque and capture probabilities in resonance. We show that planet-planet perturbations can drive the spin of Earth-size planets into asynchronous or chaotic states, even for nearly circular orbits. We apply our results to Mercury and to the KOI-1599 system of two super-Earths in a 3/2 mean motion resonance. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: 12 pages, 9 figures

Journal ref: Astronomy & Astrophysics, vol. 630, A102 (2019)

BEBOP III. Observations and an independent mass measurement of Kepler-16 (AB) b -- the first circumbinary planet detected with radial velocities

Authors: Amaury H. M. J. Triaud, Matthew R. Standing, Neda Heidari, David V. Martin, Isabelle Boisse, Alexandre Santerne, Alexandre C. M. Correia, Lorana Acuña, Matthew Battley, Xavier Bonfils, Andrés Carmona, Andrew Collier Cameron, Pía Cortés-Zuleta, Georgina Dransfield, Shweta Dalal, Magali Deleuil, Xavier Delfosse, João Faria, Thierry Forveille, Nathan C. Hara, Guillaume Hébrard, Sergio Hoyer, Flavien Kiefer, Vedad Kunovac, Pierre F. L. Maxted , et al. (8 additional authors not shown)

Abstract: The radial velocity method is amongst the most robust and most established means of detecting exoplanets. Yet, it has so far failed to detect circumbinary planets despite their relatively high occurrence rates. Here, we report velocimetric measurements of Kepler-16A, obtained with the SOPHIE spectrograph, at the Observatoire de Haute-Provence's 193cm telescope, collected during the BEBOP survey fo… ▽ More The radial velocity method is amongst the most robust and most established means of detecting exoplanets. Yet, it has so far failed to detect circumbinary planets despite their relatively high occurrence rates. Here, we report velocimetric measurements of Kepler-16A, obtained with the SOPHIE spectrograph, at the Observatoire de Haute-Provence's 193cm telescope, collected during the BEBOP survey for circumbinary planets. Our measurements mark the first radial velocity detection of a circumbinary planet, independently determining the mass of Kepler-16~(AB)~b to be $0.313 \pm 0.039\,{\rm M}_{\rm Jup}$, a value in agreement with eclipse timing variations. Our observations demonstrate the capability to achieve photon-noise precision and accuracy on single-lined binaries, with our final precision reaching $\rm 1.5~m\,s^{-1}$ on the binary and planetary signals. Our analysis paves the way for more circumbinary planet detections using radial velocities which will increase the relatively small sample of currently known systems to statistically relevant numbers, using a method that also provides weaker detection biases. Our data also contain a long-term radial velocity signal, which we associate with the magnetic cycle of the primary star. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: under review at MNRAS, 7 pages, 4 figures, all RV data in appendix

BEBOP II: Sensitivity to sub-Saturn circumbinary planets using radial-velocities

Authors: Matthew R. Standing, Amaury H. M. J. Triaud, João P. Faria, David V. Martin, Isabelle Boisse, Alexandre C. M. Correia, Magali Deleuil, Georgina Dransfield, Michaël Gillon, Guillaume Hébrard, Coel Hellier, Vedad Kunovac, Pierre F. L. Maxted, Rosemary Mardling, Alexandre Santerne, Lalitha Sairam, Stéphane Udry

Abstract: BEBOP is a radial-velocity survey that monitors a sample of single-lined eclipsing binaries, in search of circumbinary planets by using high-resolution spectrographs. Here, we describe and test the methods we use to identify planetary signals within the BEBOP data, and establish how we quantify our sensitivity to circumbinary planets by producing detection limits. This process is made easier and m… ▽ More BEBOP is a radial-velocity survey that monitors a sample of single-lined eclipsing binaries, in search of circumbinary planets by using high-resolution spectrographs. Here, we describe and test the methods we use to identify planetary signals within the BEBOP data, and establish how we quantify our sensitivity to circumbinary planets by producing detection limits. This process is made easier and more robust by using a diffusive nested sampler. In the process of testing our methods, we notice that contrary to popular wisdom, assuming circular orbits in calculating detection limits for a radial velocity survey provides over-optimistic detection limits by up to $40\%$ in semi-amplitude with implications for all radial-velocity surveys. We perform example analyses using three BEBOP targets from our Southern HARPS survey. We demonstrate for the first time a repeated ability to reach a residual root mean squared scatter of $3~\rm m.s^{-1}$ (after removing the binary signal), and find we are sensitive to circumbinary planets with masses down to that of Neptune and Saturn, for orbital periods up to $1000~\rm days$. △ Less

Submitted 10 December, 2021; originally announced December 2021.

Comments: 17 pages, 11 figures

Cyber-Cosmos: A New Citizen Science Concept in A Dark Sky Destination

Authors: Domingos Barbosa, Bruno Coelho, Miguel Bergano, Catarina Magalhães, David Mendonça, Daniela Silva, Alexandre C. M. Correia, João Pandeirada, Valério Ribeiro, Thomas Esposito, Franck Marchis

Abstract: Astrotourism and related citizen science activities are becoming a major trend of a sustainable, high-quality tourism segment, core elements to the protection of Dark skies in many countries. In the Summer of 2020, in the middle of COVID pandemics, we started an initiative to train young students - Cyber-Cosmos - using an Unistellar eVscope, a smart, compact and user-friendly digital telescope tha… ▽ More Astrotourism and related citizen science activities are becoming a major trend of a sustainable, high-quality tourism segment, core elements to the protection of Dark skies in many countries. In the Summer of 2020, in the middle of COVID pandemics, we started an initiative to train young students - Cyber-Cosmos - using an Unistellar eVscope, a smart, compact and user-friendly digital telescope that offers unprecedented opportunities for deep-sky observation and citizen science campaigns. Sponsored by the Ciência Viva Summer program, this was probably the first continuous application of this equipment in a pedagogical and citizen-science context, and in a pandemic context. Pampilhosa da Serra, home to a certified Dark Sky destination (Aldeias do Xisto) in central Portugal, was the chosen location for this project, where we expect astrotourism and citizen science to flourish and contribute to space sciences education. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: Manuscript presented at the International Astronautical Congress, IAC 2021, Dubai, United Arab Emirates, 25 - 29 October 2021. Copyright by IAF

High resolution calibration of string network evolution

Authors: J. R. C. C. C. Correia, C. J. A. P. Martins

Abstract: The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution contains a number of free parameters which cannot be obtained ab initio. Therefore it must be calibrated using high resolution numerical simulations. We exploit our state of the art graphically accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a statistically robust calibr… ▽ More The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution contains a number of free parameters which cannot be obtained ab initio. Therefore it must be calibrated using high resolution numerical simulations. We exploit our state of the art graphically accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a statistically robust calibration of this model. In order to do so, we will make use of the largest set of high resolution simulations carried out to date, for a variety of cosmological expansion rates, and explore the impact of key numerical choices on model calibration, including the dynamic range, lattice spacing, and the choice of numerical estimators for the mean string velocity. This sensitivity exploration shows that certain numerical choices will indeed have consequences for observationally crucial parameters, such as the loop chopping parameter. To conclude, we will also briefly illustrate how our results impact observational constraints on cosmic strings. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: Summary of a talk given at the From Cosmic Strings to Superstrings parallel session of the Sixteenth Marcel Grossmann Meeting, partially summarizing work previously reported in arXiv:2108.07513. To appear in the proceedings

The Three-Sided PyramidWavefront Sensor. I. Simulations and Analysis for Astronomical Adaptive Optics

Authors: Lauren Schatz, Jared R. Males, Carlos Correia, Benoit Neichel, Vincent Chambouleyron, Johanan Codona, Olivier Fauvarque, Jean-François Sauvage, Thierry Fusco, Michael Hart, Pierre Janin-Potiron, Robert Johnson, Joseph Long, Mala Mateen

Abstract: For ExAO instruments for the Giant Segmented Mirror Telescopes (GSMTs), alternative architectures of WFS are under consideration because there is a tradeoff between detector size, speed, and noise that reduces the performance of GSMT-ExAO wavefront control. One option under consideration for a GSMT-ExAO wavefront sensor is a three-sided PWFS (3PWFS). The 3PWFS creates three copies of the telescope… ▽ More For ExAO instruments for the Giant Segmented Mirror Telescopes (GSMTs), alternative architectures of WFS are under consideration because there is a tradeoff between detector size, speed, and noise that reduces the performance of GSMT-ExAO wavefront control. One option under consideration for a GSMT-ExAO wavefront sensor is a three-sided PWFS (3PWFS). The 3PWFS creates three copies of the telescope pupil for wavefront sensing, compared to the conventional four-sided PWFS (4PWFS) which uses four pupils. The 3PWFS uses fewer detector pixels than the 4PWFS and should therefore be less sensitive to read noise. Here we develop a mathematical formalism based on the diffraction theory description of the Foucault knife edge test that predicts the intensity pattern after the PWFS. Our formalism allows us to calculate the intensity in the pupil images formed by the PWFS in the presence of phase errors corresponding to arbitrary Fourier modes. We then use the Object Oriented MATLAB Adaptive Optics toolbox (OOMAO) to simulate an end-to-end model of an adaptive optics system using a PWFS with modulation and compare the performance of the 3PWFS to the 4PWFS. In the case of a low read noise detector, the Strehl ratios of the 3PWFS and 4PWFS are within 0.01. When we included higher read noise in the simulation, we found a Strehl ratio gain of 0.036 for the 3PWFS using Raw Intensity over the 4PWFS using Slopes Maps at a stellar magnitude of 10. At the same magnitude, the 4PWFS RI also outperformed the 4PWFS SM, but the gain was only 0.012 Strehl. This is significant because 4PWFS using Slopes Maps is how the PWFS is conventionally used for AO wavefront sensing. We have found that the 3PWFS is a viable wavefront sensor that can fully reconstruct a wavefront and produce a stable closed-loop with correction comparable to that of a 4PWFS, with modestly better performance for high read-noise detectors. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 39 pages, 15 figures

High resolution calibration of the cosmic strings velocity dependent one-scale model

Authors: J. R. C. C. C. Correia, C. J. A. P. Martins

Abstract: The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution must be calibrated using high resolution numerical simulations, We exploit our state of the art graphics processing unit accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a detailed and statistically robust calibration of the VOS model. We rely on the largest set of high r… ▽ More The canonical velocity-dependent one-scale (VOS) model for cosmic string evolution must be calibrated using high resolution numerical simulations, We exploit our state of the art graphics processing unit accelerated implementation of the evolution of local Abelian-Higgs string networks to provide a detailed and statistically robust calibration of the VOS model. We rely on the largest set of high resolution simulations carried out to date, with a wide range of cosmological expansion rates, and explore the impact of key numerical parameters, including the dynamic range (comparing box sizes from $1024^3$ to $4096^3$), the lattice spacing, and the choice of numerical estimators for the string velocity. We explore the sensitivity of the VOS model parameters to these numerical parameters, with a particular emphasis on the observationally crucial loop chopping efficiency, and also identify key differences between the equation of state and conjugate momentum estimators for the string velocities, showing that the latter one is more reliable for fast expansion rates (while in Minkowski space the opposite has been previously shown). Finally, we briefly illustrate how our results impact observational constraints on cosmic strings. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 16 pages, 9 figures, 2 tables; Phys. Rev. D (in press)