Ultra-hot Jupiters exhibit day-to-night temperature contrasts upwards of 1000 K due to competing effects of strong winds, short radiative timescales, magnetic drag, and H2 dissociation/recombination. Spectroscopic phase curves provide critical insights into these processes by mapping temperature distributions and constraining the planet's energy budget across different pressure levels. Here, we present the first NIRISS/SOSS phase curve of an ultra-hot Jupiter, WASP-121 b. The instrument's bandpass [0.6 - 2.85 micron] captures an estimated 50-83% of the planet's bolometric flux, depending on orbital phase, allowing for unprecedented constraints on the planet's global energy budget; previous measurements with HST/WFC3 and JWST/NIRSpec/G395H captured roughly 20% of the planetary flux. Accounting for the unobserved regions of the spectrum, we estimate effective day and nightside temperatures of T_day = 2717 +/- 17 K and T_night = 1562 +/- 19 K corresponding to a Bond albedo of A_B = 0.277 +/- 0.016 and a heat recirculation efficiency of epsilon = 0.246 +/- 0.014. Matching the phase-dependent effective temperature with energy balance models yields a similar Bond albedo of 0.3 and a mixed layer pressure of 1 bar consistent with photospheric pressures, but unexpectedly slow winds of 0.2 km/s, indicative of inefficient heat redistribution. The shorter optical wavelengths of the NIRISS/SOSS Order 2 yield a geometric albedo of A_g = 0.093 +/- 0.029 (3 sigma upper limit of 0.175), reinforcing the unexplained trend of hot Jupiters exhibiting larger Bond albedos than geometric albedos. We also detect near-zero phase curve offsets for wavelengths above 1.5 micron, consistent with inefficient heat transport, while shorter wavelengths potentially sensitive to reflected light show eastward offsets.

As Radial velocity (RV) spectrographs reach unprecedented precision and stability below 1 m/s, the challenge of granulation in the context of exoplanet detection has intensified. Despite promising advancements in post-processing tools, granulation remains a significant concern for the EPRV community. We present a pilot study to detect and characterise granulation using the High-Accuracy Radial-velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. We observed HD166620, a K2 star in the Maunder Minimum phase, intensely for two successive nights, expecting granulation to be the dominant nightly noise source in the absence of strong magnetic activity. Following the correction for a newly identified instrumental signature arising from illumination variations across the CCD, we detected the granulation signal using structure functions and a one-component Gaussian Process (GP) model. The granulation signal exhibits a characteristic timescale of 43.65$\pm$15.8 minutes, within one $\sigma$, and a standard deviation of 22.9$\pm$0.77 cm/s, with in three $\sigma$ of the predicted value. By examining spectra and RVs as a function of line formation temperature , we investigated the sensitivity of granulation-induced RV variations across different photospheric layers. We extracted RVs from various photospheric depths using both the line-by-line (LBL) and cross-correlation function (CCF) methods to mitigate any extraction method biases. Our findings indicate that granulation variability is detectable in both temperature bins, with the cooler bins, corresponding to the shallower layers of the photosphere, aligning more closely with predicted values.

GJ 536 is a low-mass star, located 10 pc away from the Sun, that hosts a low-mass planet orbiting with a period of 8.71 days. Based on an analysis of the radial velocity (RV) time series obtained from the available data of the spectrographs HARPS, HARPS-N, CARMENES and HIRES, we announce the discovery of a second low-mass planet orbiting the star. We performed a RV global analysis on RV, spectroscopic activity indicators, and ASAS photometry, within the multidimensional Gaussian process framework, updated the parameters of GJ 536 b, and found significant evidence of the presence of a second planet. GJ 536 c is a low-mass planet ($m_{p} \sin i$ = 5.89 $\pm$ 0.70 M$_{\oplus}$), orbiting with a period of 32.761 $\pm$ 0.015 days, at a distance of 0.1617 $\pm$ 0.0028 au from its parent star. It induces an RV semi-amplitude of 1.80 $\pm$ 0.20 m$\cdot$s$^{-1}$. Given its distance to the star, it receives a flux of 1.692 $\pm$ 0.069 F$_{\oplus}$, for an equilibrium temperature of 290.5 $\pm$ 9.5 K. We update the mass of the planet GJ 536 b to $m_{p} \sin i$ = 6.37 $\pm$ 0.38 M$_{\oplus}$. The orbits of both planets are consistent with circular. We explored the use of statistical Doppler imaging on the photometric and RV data, and find a tentative projected obliquity of the stellar rotation axis of 58$^{+16}_{-19}$ deg. Current evidence does not support the presence of additional planets with masses > 5 M$_{\oplus}$for orbital periods up to 100 days, or > 10 M$_{\oplus}$ for periods up to 1000 days.

Thousands of tight (<1 AU) main sequence binaries have been discovered, but it is uncertain how they formed. There is likely too much angular momentum in a collapsing, fragmenting protostellar cloud to form such binaries in situ, suggesting some post processing. One probe of a binary's dynamical history is the angle between the stellar spin and orbital axes -- its obliquity. The classical method for determining stellar obliquity is the Rossiter-McLaughlin effect. It has been applied to over 100 hot Jupiters, but to less than a dozen stellar binaries. In this paper, we present the Rossiter-McLaughlin measurement of EBLM J0021-16, a $0.19M_\odot$ M-dwarf eclipsing a $1.05M_\odot$ G-dwarf on a 5.97 day, almost-circular orbit. We combine CORALIE spectroscopy with TESS photometry of primary and secondary eclipses and star spot modulation. We show that the orbital axis is well-aligned with the primary star's spin axis, with a true 3D obliquity of $\psi=4.01\pm0.54^{\circ}$. EBLM J0021-16 becomes one of only a handful of eclipsing binaries where a true obliquity has been measured. Finally, we derive the M-dwarf's mass and radius to a fractional precision better than 1\%. The radius of the M dwarf is inflated by 6\% ($7.4\sigma$) with respect to stellar models, consistent with many other M-dwarfs in the literature.

This paper provides a comprehensive overview of the subsystems of the NIGHT instrument. NIGHT (the Near Infrared Gatherer of Helium Transits) is a narrowband, high-resolution spectrograph, marking the first dedicated survey instrument for exoplanetary atmosphere observations. Developed through a collaboration between the Observatory of Geneva and the Universite de Montreal, NIGHT aims to conduct an extensive statistical survey of helium atmospheres around 100+ exoplanets over several years. The instrument will report new detections of helium in exoplanet atmospheres and perform temporal monitoring of a subset of these. NIGHT measures absorption from the metastable helium state during exoplanet transits, observable in a triplet of lines around 1083nm. The instrument comprises a vacuum enclosure housing the spectrograph, a front end unit for fiber injection at the telescope's focal plane, and a calibration and control rack containing calibration light sources and control hardware. The spectrograph is optimized for efficiency, achieving a uniform throughput of approximately 71%. The primary disperser employs a VPH grating in a unique double-pass configuration, enabling a spectral resolution of 75,000 while maintaining high throughput. The detector is a HAWAII-1 infrared array, cooled to 85K, with the spectrograph operating at room temperature. Thanks to its relatively high throughput, NIGHT on a 2m class telescope is predicted to be as sensitive as existing instruments on 4m class telescopes. The front end unit injects starlight and sky background into two separate fibers leading to the spectrograph. It also performs near-infrared guiding and includes a mechanism for injecting calibration light. The assembly and optical alignment of NIGHT's spectrograph and front end unit are scheduled for July to September 2024, with the first light anticipated before early 2025.

High quality wavelength calibration is crucial for science cases like radial-velocity studies of exoplanets, the search for a possible variation of fundamental constants, and the redshift drift experiment. However, for state-of-the-art spectrographs it has become difficult to verify the wavelength calibration on sky, because no astrophysical source provides spectra with sufficiently stable or accurate wavelength information. We therefore propose to use iodine absorption cells to validate the wavelength calibration. Observing a bright and featureless star through the iodine cell emulates an astrophysical target with exactly known spectral features that can be analyzed like any other science target, allowing to verify the wavelength calibration derived from the internal calibration sources and to identify systematics in the data processing. As demonstration, we temporarily installed an I$_2$ absorption cell at ESPRESSO. Employing a full forward modeling approach of the I$_2$ spectrum, including the instrumental line-spread function, we demonstrate wavelength calibration accuracy at the level of a few m/s. We also show that wavelength measurements do depend on the geometry of the light-injection into the spectrograph fibers. This highlights the importance of probing exactly the same light path as science targets, something not possible with internal calibration sources alone. We also demonstrate excellent radial-velocity stability at the <20 cm/s level in a full end-to-end fashion, from sky to data product. Our study therefore showcases the great potential of absorption cells for the verification and long-term monitoring of the wavelength calibration as well as the unique insights they can provide.

The measurement of the temporal evolution in the redshift of distant objects, the redshift drift, is a probe of universal expansion and cosmology. We perform the first steps towards a measurement of such effect using the Lyman-$\alpha$ forest in the spectra of bright quasars as a tracer of cosmological expansion. Our goal is to determine to which precision a velocity shift measurement can be carried out with the signal-to-noise (S/N) level currently available and whether this precision aligns with previous theoretical expectations. A precise assessment of the achievable measurement precision is fundamental for estimating the time required to carry out the whole project. We acquire 12 hours of ESPRESSO observations distributed over 0.875 years of the brightest quasar known, J052915.80-435152.0 (z=3.962), to obtain high-resolution spectra of the Lyman-$\alpha$ forest, with median S/N of ~86 per 1 km/s pixel at the continuum. We divide the observations into two epochs and analyse them using both a pixel-by-pixel method and a model-based approach. This comparison allows us to estimate the velocity shift between the epochs, as well as the velocity precision that can be achieved at this S/N. The model-based method is calibrated using high-resolution simulations of the intergalactic medium, and it provides greater accuracy compared to the pixel-by-pixel approach. We measure a velocity drift of the Lyman-$\alpha$ forest consistent with zero: $\Delta v = -1.25\pm 4.45 {\rm ms^{-1}}$, equivalent to a cosmological drift of $\dot{v}=-1.43\pm 5.09 {\rm ms^{-1}yr^{-1}}$ or $\dot{z}= (-2.19\pm7.77) \times 10^{-8}{\rm yr^{-1}}$. The measurement uncertainties are on par with the expected precision. We estimate that reaching a 99% detection of the cosmic drift requires a monitoring campaign of 5400 hours of integration time over 54 years with an ELT and an ANDES-like high-resolution spectrograph.

(Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$\delta$. The derived values vary from 26.29 days (H$\gamma$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$\alpha$, H$\beta$, H$\gamma$, H$\epsilon$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$\alpha$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$\alpha$ activity indicators and between the S index and radial velocity variations.

We analyse spectroscopic and photometric transits of the hot Jupiters WASP-52b and HAT-P30b obtained with ESPRESSO, Eulercam and NGTS for both targets, and additional TESS data for HAT-P-30. Our goal is to update the system parameters and refine our knowledge of the host star surfaces. For WASP-52, the companion planet has occulted starspots in the past, and as such our aim was to use the reloaded Rossiter-McLaughlin technique to directly probe its starspot properties. Unfortunately, we find no evidence for starspot occultations in the datasets herein. Additionally, we searched for stellar surface differential rotation (DR) and any centre-to-limb variation (CLV) due to convection, but return a null detection of both. This is unsurprising for WASP-52, given its relatively cool temperature, high magnetic activity (which leads to lower CLV), and projected obliquity near 0 degrees (meaning the transit chord is less likely to cross several stellar latitudes). For HAT-P-30, this result was more surprising given its hotter effective temperature, lower magnetic field, and high projected obliquity (near 70 degrees). To explore the reasons behind the null DR and CLV detection for HAT-P-30, we simulated a variety of scenarios. We find that either the CLV present on HAT-P-30 is below the solar level or the presence of DR prevents a CLV detection given the precision of the data herein. A careful treatment of both DR and CLV is required, especially for systems with high impact factors, due to potential degeneracies between the two. Future observations and/or a sophisticated treatment of the red noise present in the data (likely due to granulation) is required to refine the DR and CLV for these particular systems; such observations would also present another opportunity to try to examine starspots on WASP-52.

Doppler Imaging (DI) is a well-established technique to map a physical field at a stellar surface from a time series of high-resolution spectra. In this proof-of-concept study, we aim to show that traditional DI algorithms, originally designed for rapidly-rotating stars, have also the ability to model the activity of Sun-like stars, when observed with new-generation highly-stable spectrographs, and search for low-mass planets around them. We used DI to retrieve the relative brightness distribution at the surface of the Sun from radial velocity (RV) observations collected by HARPS-N between 2022 and 2024. The brightness maps obtained with DI have a typical angular resolution of about 36 degrees and are a good match to low-resolution disc-resolved Dopplergrams of the Sun at epochs when the absolute, disc-integrated RV exceeds ~2 m/s. The RV residuals after DI correction exhibit a dispersion of about 0.6 m/s, comparable with existing state-of-the-art activity correction techniques. Using planet injection-recovery tests, we also show that DI can be a powerful tool for blind planet searches, so long as the orbital period is larger than ~100days (i.e. 3 to 4 stellar rotation periods), and that it yields planetary mass estimates with an accuracy comparable to, for example, multi-dimensional Gaussian process regression. Finally, we highlight some limitations of traditional DI algorithms, which should be addressed to make DI a reliable alternative to state-of-the-art RV-based planet search techniques.

We re-determine planetary occurrences around M dwarfs using 20 years of observations from HARPS on 197 targets. The first aim of this study is to propose more precise occurrence rates using the large volume of the sample but also variations to previous calculations, particularly by considering multiplicity, which is now an integral part of planetary occurrence calculations. The second aim is to exploit the extreme longevity of HARPS to determine occurrence rates in the unexplored domain of very long periods. This work relies entirely on the 197 radial velocity time series obtained and analysed in our previous study. By considering they are cleaned of any detectable signal, we convert them into detection limits. We use these 197 limits to produce a detectability map and combine it with confirmed planet detections to establish our occurrence rates. Finally, we also convert the detection limits from orbital period to insolation in order to construct an occurrence statistics for the temperate zone. We find a strong prevalence of low-mass planets around M dwarfs, with an occurrence rate of 120% for planets with a mass between 0.75 and 3 Me. In addition, we compute an occurrence rate of 45.3% +20-16% for temperate zone planets around M dwarfs. We obtain an occurrence rate of a few percent for giant planets with wide separations. In our sample these giant planets with wide separations are only detected around the most massive M dwarfs.

We search for excess in-transit absorption of neutral helium at 1.083 $\mu$m in the atmospheres of the young (<800 Myr) sub-Jovian (0.2-0.5 $\rm R_{J}$) planets HD 63433b, K2-100b, and V1298 Tau c using high-resolution (R~25,000) transit observations taken with Keck II/NIRSPEC. Our observations do not show evidence of helium absorption for any of the planets in our sample. We calculate 3$\sigma$ upper limits on the planets' excess helium absorption of <0.47% for HD 63433b, <0.56% for K2-100b, and <1.13% for V1298 Tau c. In terms of equivalent width, we constrain these to <2.52, <4.44, and <8.49 mA for HD 63433b, K2-100b, and V1298 Tau c, respectively. We fit our transmission spectra with one-dimensional Parker wind models to determine upper limits on the planets' mass-loss rates of <7.9$\times10^{10}$, <1.25$\times10^{11}$, and <$7.9\times10^{11}$g s$^{-1}$. Our non-detections align with expectations from one-dimensional hydrodynamic escape models, magnetic fields, and stellar wind confinement. The upper limits we measure for these planets are consistent with predicted trends in system age and He equivalent width from 1D hydrodynamic models.

Stellar activity is the main limitation to the detection of Earth-twins using the RV technique. Despite many efforts in trying to mitigate the effect of stellar activity using empirical and statistical techniques, it seems that we are facing an obstacle that will be extremely difficult to overcome using current techniques. In this paper, we investigate a novel approach to derive precise RVs considering the wealth of information present in high-resolution spectra. This new method consists in building a master spectrum from all observations and measure the RVs of each spectral line in a spectrum relative to it. When analysing several spectra, the final product is the RVs of each line as a function of time. We demonstrate on three stars intensively observed with HARPS that our new method gives RVs that are extremely similar to the ones derived from the HARPS data reduction software. Our new approach to derive RVs demonstrates that the non-stability of daily HARPS wavelength solution induces night-to-night RV offsets with an standard deviation of 0.4 m/s, and we propose a solution to correct for this systematic. Finally, and this is probably the most astrophysically relevant result of this paper, we demonstrate that some spectral lines are strongly affected by stellar activity while others are not. By measuring the RVs on two carefully selected subsample of spectral lines, we demonstrate that we can boost by a factor of 2 or mitigate by a factor of 1.6 the red noise induced by stellar activity in the 2010 RVs of Alpha Cen B. By measuring the RVs of each spectral line, we are able to reach the same RV precision as other approved techniques. In addition, this new approach allows to demonstrate that each line is differently affected by stellar activity. Preliminary results show that studying in details the behaviour of each spectral line is probably the key to overcome stellar activity.

Close-in planets evolve under extreme conditions, raising questions about their origins and current nature. Two predominant mechanisms are orbital migration, which brings them close to their star, and atmospheric escape under the resulting increased irradiation. Yet, their relative roles remain unclear because we lack models that couple the two mechanisms with high precision on secular timescales. To address this need, we developed the JADE code, which simulates the secular atmospheric and dynamical evolution of a planet around its star, and can include the perturbation induced by a distant third body. On the dynamical side, the 3D evolution of the orbit is modeled under stellar and planetary tidal forces, a relativistic correction, and the action of the distant perturber. On the atmospheric side, the vertical structure of the atmosphere is integrated over time based on its thermodynamical properties, inner heating, and the evolving stellar irradiation, which results, in particular, in photo-evaporation. The JADE code is benchmarked on GJ436 b, prototype of evaporating giants on eccentric, misaligned orbits at the edge of the hot Neptunes desert. We confirm that its orbital architecture is well explained by Kozai migration and unveil a strong interplay between its atmospheric and orbital evolution. During the resonance phase, the atmosphere pulsates in tune with the Kozai cycles, which leads to stronger tides and an earlier migration. This triggers a strong evaporation several Gyr after the planet formed, refining the paradigm that mass loss is dominant in the early age of close-in planets. This suggests that the edge of the desert could be formed of warm Neptunes whose evaporation was delayed by migration. It strengthens the importance of coupling atmospheric and dynamical evolution over secular timescales, which the JADE code will allow simulating for a wide range of systems.

WASP-52b is a low density hot Jupiter orbiting a moderately active K2V star. Previous low-resolution studies have revealed a cloudy atmosphere and found atomic Na above the cloud deck. Here we report the detection of excess absorption at the Na doublet, the H$\alpha$ line, and the K D$_1$ line. We derived a high-resolution transmission spectrum based on three transits of WASP-52b, observed with the ultra stable, high-resolution spectrograph ESPRESSO at the VLT. We measure a line contrast of $1.09\pm 0.16$% for Na D$_1$, $1.31\pm 0.13$% for Na D$_2$, $0.86\pm 0.13$% for H$\alpha$, and $0.46\pm 0.13$$_1$, with a line FWHM range of 11-22 km s$^{-1}$. We also find that the velocity shift of these detected lines during the transit is consistent with the planet orbital motion, thus confirming their planetary origin. We do not observe any significant net blueshift or redshift that can be attributed to planetary winds. We use activity indicator lines as control but find no excess absorption. However, we do notice signatures arising from the Center-to-Limb variation (CLV) and the Rossiter-McLaughlin (RM) effect at these control lines. This highlights the importance of the CLV+RM correction in correctly deriving the transmission spectrum, which, if not corrected, could resemble or cancel out planetary absorption in certain cases. WASP-52b is the second non-ultra-hot Jupiter to show excess H$\alpha$ absorption, after HD 189733b. Future observations targeting non-ultra-hot Jupiters that show H$\alpha$ could help reveal the relation between stellar activity and the heating processes in the planetary upper atmosphere.

K2-18 b is a transiting mini-Neptune that orbits a nearby (38 pc) cool M3 dwarf and is located inside its region of temperate irradiation. We report on the search for hydrogen escape from the atmosphere K2-18 b using Lyman-$\alpha$ transit spectroscopy with the Space Telescope Imaging Spectrograph (STIS) instrument installed on the Hubble Space Telescope (HST). We analyzed the time-series of the fluxes of the stellar Lyman-$\alpha$ emission of K2-18 in both its blue- and redshifted wings. We found that the average blueshifted emission of K2-18 decreases by $67\% \pm 18\%$ during the transit of the planet compared to the pre-transit emission, tentatively indicating the presence of H atoms escaping vigorously and being blown away by radiation pressure. This interpretation is not definitive because it relies on one partial transit. Based on the reconstructed Lyman-$\alpha$ emission of K2-18, we estimate an EUV irradiation between $10^1-10^2$ erg s$^{-1}$ cm$^{-2}$ and a total escape rate in the order of $10^8$ g s$^{-1}$. The inferred escape rate suggests that the planet will lose only a small fraction (< 1%) of its mass and retain its volatile-rich atmosphere during its lifetime. More observations are needed to rule out stellar variability effects, confirm the in-transit absorption and better assess the atmospheric escape and high-energy environment of K2-18 b.

(Abridged) The quiet M2.5 star GJ 436 hosts a warm Neptune that displays an extended atmosphere that dwarfs its own host star. Predictions of atmospheric escape in such planets state that H atoms escape from the upper atmosphere in a collisional regime and that the flow can drag heavier atoms to the upper atmosphere. It is unclear, however, what astrophysical mechanisms drive the process. Our objective is to leverage the extensive coverage of HST/COS observations of the far-ultraviolet (FUV) spectrum of GJ 436 to search for signals of metallic ions in the upper atmosphere of GJ 436 b. We analyzed flux time-series of species present in the FUV spectrum of GJ 436, as well as the Lyman-$\alpha$ line. GJ 436 displays FUV flaring events with a rate of $\sim$10 d$^{-1}$. There is evidence for a possibly long-lived active region or longitude that modulates the FUV metallic lines of the star with amplitudes up to 20%. Despite the strong geocoronal contamination in the COS spectra, we detected in-transit excess absorption signals of $\sim$50% and $\sim$30% in the blue and red wings, respectively, of the Lyman-$\alpha$ line. We rule out a wide range of excess absorption levels in the metallic lines of the star during the transit. The large atmospheric loss of GJ 436 b observed in Lyman-$\alpha$ transmission spectra is stable over the timescale of a few years, and the red wing signal supports the presence of a variable hydrogen absorption source besides the stable exosphere. The previously claimed in-transit absorption in the Si III line is likely an artifact resulting from the stellar magnetic cycle. The non-detection of metallic ions in absorption could indicate that the escape is not hydrodynamic or that the atmospheric mixing is not efficient in dragging metals high enough for sublimation to produce a detectable escape rate of ions to the exosphere.

Radial velocity perturbations induced by stellar surface inhomogeneities including spots, plages and granules currently limit the detection of Earth-twins using Doppler spectroscopy. Such stellar noise is poorly understood for stars other than the Sun because their surface is unresolved. In particular, the effects of stellar surface inhomogeneities on observed stellar radial velocities are extremely difficult to characterize, and thus developing optimal correction techniques to extract true stellar radial velocities is extremely challenging. In this paper, we present preliminary results of a solar telescope built to feed full-disk sunlight into the HARPS-N spectrograph, which is in turn calibrated with an astro-comb. This setup enables long-term observation of the Sun as a star with state-of-the-art sensitivity to radial velocity changes. Over seven days of observing in 2014, we show an average 50\cms radial velocity rms over a few hours of observation. After correcting observed radial velocities for spot and plage perturbations using full-disk photometry of the Sun, we lower by a factor of two the weekly radial velocity rms to 60\cms. The solar telescope is now entering routine operation, and will observe the Sun every clear day for several hours. We will use these radial velocities combined with data from solar satellites to improve our understanding of stellar noise and develop optimal correction methods. If successful, these new methods should enable the detection of Venus over the next two to three years, thus demonstrating the possibility of detecting Earth-twins around other solar-like stars using the radial velocity technique.