Researcher profile

C. J. A. P. Martins

C. J. A. P. Martins contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
43works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid&#39;s Flagship simulation to investigate its impact on Euclid&#39;s weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock&#39;) galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid&#39;s Data Release 1 (DR1).

preprint2025arXiv

Unveiling the nature of HE 0107-5240

The vast majority of the most iron-poor stars in the Galaxy exhibit a strong carbon enhancement, with C/H ratios only about two orders of magnitude below solar. This unusual chemical composition likely reflects the properties of the gas cloud from which these stars formed, having been enriched by one, or at most a few, supernovae. A remarkable member of this stellar class, HE 0107-5240 with [Fe/H]=-5.56, has been identified as part of a binary system. To constrain its orbital parameters, radial velocity monitoring has been carried out using the ESPRESSO spectrograph. Radial velocities were derived using cross-correlation with a template, taking advantage of the strong G-band feature. Combining all observations yielded a high signal-to-noise spectrum, which has been used to refine our understanding of the stellar chemical composition. Additionally, a co-added UVES spectrum in the blue was used to complement the wavelength coverage of ESPRESSO. Observations of HE 0107-5240 over a span of more than four years have yielded a revised orbital period of about 29 years. Updated elemental abundances have been determined for Sc, Cr, Co, and, tentatively, Al, along with a new upper limit for Be. The iron abundance has been derived from ionised Fe lines. Significant upper limits have been established for Li, Si, and Sr. The star is confirmed to be a long-period binary. Iron abundances derived from neutral and ionised lines are consistent with local thermodynamical equilibrium (LTE) assumption, casting doubt on published deviation from LTE corrections for Fe for this star. The heavy elements Sr and Ba remain undetected, confirming the classification of HE 0107-5240 as a carbon enhanced metal-poor and non enhanced in heavy elements (CEMP-no) star and supporting the absence of an n-capture element plateau at the lowest metallicities.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2023arXiv

The pristine nature of SMSS 1605$-$1443 revealed by ESPRESSO

SMSS J160540.18$-$144323.1 is the carbon-enhanced metal-poor (CEMP) star with the lowest iron abundance ever measured, [Fe/H]=-6.2, which was first reported with the SkyMapper telescope. The carbon abundance is A(C)~6.1 in the low-C band, as the majority of the stars in this metallicity range. Yet, constraining the isotopic ratio of key species, such as carbon, sheds light on the properties and origin of these elusive stars. We performed high-resolution observations of SMSS1605$-$1443 with the ESPRESSO spectrograph to look for variations in the radial velocity ($v_{rad}$) with time. These data have been combined with older MIKE and UVES archival observations to enlarge the temporal baseline. The $^{12}$C/$^{13}$C isotopic ratio is also studied to explore the possibility of mass transfer from a binary companion. A cross-correlation function against a natural template was applied to detect $v_{rad}$ variability and a spectral synthesis technique was used to derive $^{12}$C/$^{13}$C in the stellar atmosphere. We confirm previous indications of binarity in SMSS1605$-$1443 and measured a lower limit $^{12}$C/$^{13}$C$>60$ at more than a 3$σ$ confidence level, proving that this system is chemically unmixed and that no mass transfer from the unseen companion has happened so far. Thus, we confirm the CEMP-no nature of SMSS1605$-$1443 and show that the pristine chemical composition of the cloud from which it formed is currently imprinted in its stellar atmosphere free of contamination.

preprint2022arXiv

A candidate short-period sub-Earth orbiting Proxima Centauri

Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough evaluation of the presence of a third low-mass planetary companion, which started emerging during a previous campaign. Radial velocities (RVs) were calculated using both a cross-correlation function (CCF) and a template matching approach. The RV analysis includes a component to model Proxima&#39;s activity using a Gaussian process (GP). We use the CCF&#39;s full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit.

preprint2022arXiv

A novel framework for semi-Bayesian radial velocities through template matching

The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alternative for M-dwarf stars. In this paper, we describe a new implementation of template matching within a semi-Bayesian framework, providing a more statistically principled characterization of the RV measurements. In this context, a common RV shift is used to describe the difference between each spectral order of a given stellar spectrum and a template built from the available observations. Posterior probability distributions are obtained for the relative RV associated with each spectrum, after marginalizing with respect to the continuum. This methodology was named S-BART: Semi-Bayesian Approach for RVs with Template-matching, and it can be applied to HARPS and ESPRESSO. The application of our method to HARPS archival observations of Barnard&#39;s star allowed us to validate our implementation against HARPS-TERRA and SERVAL. Then, we applied it to 33 ESPRESSO targets, evaluating its performance and comparing it with the CCF method. We found a decrease in the median RV scatter of \sim 10\% and \sim 4\% for M- and K-type stars, respectively. S-BART yields more precise RV estimates than the CCF method, particularly in the case of M-type stars where a median uncertainty of \sim 15 cm/s is achieved over 309 observations. Further, we estimated the nightly zero point (NZP) of ESPRESSO, finding a weighted NZP scatter below \sim 0.7 m/s. As this includes stellar variability, photon noise, and potential planetary signals, it should be taken as an upper limit of the RV precision attainable with ESPRESSO data.

preprint2022arXiv

CaRM: Exploring the chromatic Rossiter-McLaughlin effect. The cases of HD 189733b and WASP-127b

In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO.

preprint2022arXiv

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

The standard cosmological model is based on the fundamental assumptions of a spatially homogeneous and isotropic universe on large scales. An observational detection of a violation of these assumptions at any redshift would immediately indicate the presence of new physics. We quantify the ability of the Euclid mission, together with contemporary surveys, to improve the current sensitivity of null tests of the canonical cosmological constant $Λ$ and the cold dark matter (LCDM) model in the redshift range $0<z<1.8$. We considered both currently available data and simulated Euclid and external data products based on a LCDM fiducial model, an evolving dark energy model assuming the Chevallier-Polarski-Linder (CPL) parameterization or an inhomogeneous Lema\^ıtre-Tolman-Bondi model with a cosmological constant $Λ$ (LLTB), and carried out two separate but complementary analyses: a machine learning reconstruction of the null tests based on genetic algorithms, and a theory-agnostic parametric approach based on Taylor expansion and binning of the data, in order to avoid assumptions about any particular model. We find that in combination with external probes, Euclid can improve current constraints on null tests of the LCDM by approximately a factor of three when using the machine learning approach and by a further factor of two in the case of the parametric approach. However, we also find that in certain cases, the parametric approach may be biased against or missing some features of models far from LCDM. Our analysis highlights the importance of synergies between Euclid and other surveys. These synergies are crucial for providing tighter constraints over an extended redshift range for a plethora of different consistency tests of some of the main assumptions of the current cosmological paradigm.

preprint2022arXiv

Fundamental physics with ESPRESSO: Constraints on Bekenstein and dark energy models from astrophysical and local probes

Dynamical scalar fields in an effective four-dimensional field theory are naturally expected to couple to the rest of the theory&#39;s degrees of freedom, unless some new symmetry is postulated to suppress these couplings. In particular, a coupling to the electromagnetic sector will lead to spacetime variations of the fine-structure constant, $α$. Astrophysical tests of the space-time stability of $α$ are therefore a powerful probe of new physics. Here we use ESPRESSO and other contemporary measurements of $α$, together with background cosmology data, local laboratory atomic clock and Weak Equivalence Principle measurements, to place stringent constraints on the simplest examples of the two broad classes of varying $α$ models: Bekenstein models and quintessence-type dark energy models, both of which are parametric extensions of the canonical $Λ$CDM model. In both cases, previously reported constraints are improved by more than a factor of ten. This improvement is largely due to the very strong local constraints, but astrophysical measurements can help to break degeneracies between cosmology and fundamental physics parameters.

preprint2022arXiv

Low redshift constraints on scale-covariant models

The search for a physical model which explains the observed recent acceleration of the universe is a compelling task of modern fundamental cosmology. Recently Fernandes \textit{et al.} presented low redshift observational constraints on a scale invariant model by Maeder. Phenomenologically this can be interpreted as a bimetric theory with a time-dependent cosmological constant. It was shown that a matter density $Ω_m\sim0.3$ is a poor fit to the data, and the best-fit model would require a fluid with a much smaller density and a significantly positive equation of state parameter. This model is a particular case of an earlier and broader class of models by Canuto \textit{et al.}, which we study here. Specifically, we consider it in two distinct scenarios: as a genuine alternative to $Λ$CDM (i.e., without any cosmological constant) and as a parametric extension thereof (where both a cosmological constant and the new mechanism can coexist, and the relative contributions of both are determined by the data). We find that the first scenario can in principle fit the low-redshift data (but a good fit would require values of model parameters, such as the matter equation of state, in conflict with other data), while in the second one the deviation from $Λ$CDM is constrained to be small.

preprint2022arXiv

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

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.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Observational constraints on nonlinear matter extensions of general relativity

We present a phenomenological analysis of current observational constraints on classes of FLRW cosmological models in which the matter side of Einstein&#39;s equations includes, in addition to the canonical term, a term proportional to some function of the energy-momentum tensor ($T^2=T_{αβ}T^{αβ}=ρ^2+3p^2$), or of its trace ($T=ρ-3p$). Qualitatively, one may think of these models as extensions of general relativity with a nonlinear matter Lagrangian. As such they are somewhat different from the usual dynamical dark energy or modified gravity models: in the former class of models one adds further dynamical degrees of freedom to the Lagrangian (often in the form of scalar fields), while in the latter the gravitational part of the Lagrangian is changed. We study both of these models under two different scenarios: (1) as phenomenological two-parameter or three-parameter extensions of the standard $Λ$CDM, in which case the model still has a cosmological constant but the nonlinear matter Lagrangian leads to additional terms in Einstein&#39;s equations, which cosmological observations tightly constrain, and (2) as alternatives to $Λ$CDM, where there is no cosmological constant, and the nonlinear matter term would have to provide the acceleration (which would be somewhat closer in spirit to the usual modified gravity models). A comparative analysis of the observational constraints obtained in the various cases provides some insight on the level of robustness of the $Λ$ model and on the parameter space still available for phenomenological alternatives.

preprint2022arXiv

Observational constraints on nonlinear matter extensions of general relativity: Separable trace power models

The search for the physical mechanism underlying the observational evidence for the acceleration of the recent universe is a compelling goal of modern fundamental cosmology. Here we quantitatively study a class of homogeneous and isotropic cosmological models in which the matter side of Einstein&#39;s equations includes, in addition to the canonical term, a term proportional to the trace of the energy-momentum tensor, $T=ρ-3p$, and constrain these models using low redshift background cosmology data. One may think of these models as extensions of general relativity with a nonlinear matter Lagrangian, and they can be studied either as phenomenological extensions of the standard $Λ$CDM model, containing both matter and a cosmological constant, or as direct alternatives to it, where there is no cosmological constant but the additional terms would have to be responsible for accelerating the universe. Overall, our main finding is that parametric extensions of $Λ$CDM are tightly constrained, with additional model parameters being constrained to their canonical behaviours to within one standard deviation, while alternative models in this class (which do not have a $Λ$CDM limit) are ruled out. This provides some insight on the level of robustness of the $Λ$CDM model and on the parameter space still available for phenomenological alternatives and extensions.

preprint2022arXiv

The polar orbit of the warm Neptune GJ436b seen with VLT/ESPRESSO

GJ436b might be the prototype of warm Neptunes that have undergone late migration induced by an outer companion. Precise determination of the orbital architecture of such systems is critical to constraining their dynamical history and evaluating the role of delayed migration in the exoplanet population. To this purpose we analyzed the Rossiter-McLaughlin (RM) signal of GJ436 b in two transits - recently observed with ESPRESSO - using three different techniques. The high level of precision achieved in radial velocity (RV) measurements allows us to detect the deviation from the Keplerian orbit, despite the slow rotation of the M dwarf host (vsini = 272.0+40.0-34.0 m/s), and to measure the sky-projected obliquity ($λ$ = 102.5+17.2-18.5$^{\circ}$). The Reloaded RM technique, which allows the stellar RV field along the transit chord to be analyzed, yields $λ$ = 107.5+23.6-19.3$^{\circ}$ and vsini = 292.9+41.9-49.9 m/s. The RM Revolutions technique, which allows us to fit the spectral profiles from all planet-occulted regions together, yields $λ$ = 114.1+22.8-17.8$^{\circ}$ and vsini = 300.5+45.9-57.0 m/s. The consistent results between these three techniques, and with published results from HARPS/HARPS-N data, confirm the polar orbit of GJ436b and support the hypothesis that its origin lies in Kozai migration. Results from a joint RM Revolutions analysis of the ESPRESSO, HARPS, and HARPS-N datasets ($λ$ = 113.5+23.3-17.3$^{\circ}$; vsini = 293.5+43.7-52.2 m/s) combined with a revised stellar inclination ($i_\star$ = 35.7+5.9-7.6$^{\circ}$ or 144.2+7.6-5.9$^{\circ}$) lead us to constrain the 3D obliquity $Ψ$ to 103.2+12.8-11.5$^{\circ}$.

preprint2022arXiv

Transmission spectroscopy of MASCARA-1b with ESPRESSO: Challenges of overlapping orbital and Doppler tracks

Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$α$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b.

preprint2022arXiv

Varying fine-structure constant cosmography

Cosmography is a phenomenological and relatively model-independent approach to cosmology, where physical quantities are expanded as a Taylor series in the cosmological redshift, or in related variables. Here we apply this methodology to constrain possible cosmological variations of the fine-structure constant, $α$. Two peculiarities of this case are the existence of high-redshift data, and the fact that one term in the series is directly and tightly constraint by local laboratory tests with atomic clocks. We use this atomic clock data, together with direct model-independent high-resolution astrophysical spectroscopy measurements of $α$ up to redshift $z\sim7$ and additional model-dependent constraints on $α$ from the cosmic microwave background and big bang nucleosynthesis, to place stringent (parts per million level) constraints on the first two terms in the $α$ cosmographic series.

preprint2021arXiv

Euclid: Forecast constraints on the cosmic distance duality relation with complementary external probes

In metric theories of gravity with photon number conservation, the luminosity and angular diameter distances are related via the Etherington relation, also known as the distance-duality relation (DDR). A violation of this relation would rule out the standard cosmological paradigm and point at the presence of new physics. We quantify the ability of Euclid, in combination with contemporary surveys, to improve the current constraints on deviations from the DDR in the redshift range $0<z<1.6$. We start by an analysis of the latest available data, improving previously reported constraints by a factor of 2.5. We then present a detailed analysis of simulated Euclid and external data products, using both standard parametric methods (relying on phenomenological descriptions of possible DDR violations) and a machine learning reconstruction using Genetic Algorithms. We find that for parametric methods Euclid can (in combination with external probes) improve current constraints by approximately a factor of six, while for non-parametric methods Euclid can improve current constraints by a factor of three. Our results highlight the importance of surveys like Euclid in accurately testing the pillars of the current cosmological paradigm and constraining physics beyond the standard cosmological model.

preprint2021arXiv

Retrieving the transmission spectrum of HD 209458b using CHOCOLATE: A new chromatic Doppler tomography technique

Multiband photometric transit observations or low-resolution spectroscopy (spectro-photometry) are normally used to retrieve the broadband transmission spectra of transiting exoplanets in order to assess the chemical composition of their atmospheres. In this paper, we present an alternative approach for recovering the broadband transmission spectra using chromatic Doppler tomography based on physical modeling through the SOAP tool: CHOCOLATE (CHrOmatiC line prOfiLe tomogrAphy TEchnique). To validate the method and examine its performance, we use observational data recently obtained with the ESPRESSO instrument to retrieve the transmission spectra of the archetypal hot Jupiter HD 209458b. Our findings indicate that the recovered transmission spectrum is in good agreement with the results presented in previous studies, which used different methodologies to extract the spectrum, achieving similar precision. We explored several atmospheric models and inferred from spectral retrieval that a model containing H2O and NH3 is the preferred scenario. The CHOCOLATE methodology is particularly interesting for future studies of exoplanets around young and active stars or moderate to fast rotating stars, considering SOAP&#39;s ability to model stellar active regions and the fact that the rotational broadening of spectral lines favors its application. Furthermore, CHOCOLATE will allow the broad transmission spectrum of a planet to be retrieved using high signal-to-noise ratio, high-resolution spectroscopy with the next generation of Extremely Large Telescopes (ELTs), where low-resolution spectroscopy will not always be accessible.

preprint2021arXiv

Six transiting planets and a chain of Laplace resonances in TOI-178

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth&#39;s density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.

preprint2021arXiv

The atmosphere of HD 209458b seen with ESPRESSO. No detectable planetary absorptions at high resolution

We observed two transits of the iconic gas giant HD 209458b between 380 and 780 nm, using the high-resolution ESPRESSO spectrograph. The derived planetary transmission spectrum exhibits features at all wavelengths where the parent star shows strong absorption lines, for example, NaI, MgI, FeI, FeII, CaI, VI, H$α$, and KI. We interpreted these features as the signature of the deformation of the stellar line profiles due to the Rossiter-McLaughlin effect, combined with the centre-to-limb effects on the stellar surface, which is in agreement with similar reports recently presented in the literature. We also searched for species that might be present in the planetary atmosphere but not in the stellar spectra, such as TiO and VO, and obtained a negative result. Thus, we find no evidence of any planetary absorption, including previously reported NaI, in the atmosphere of HD 209458b. The high signal-to-noise ratio in the transmission spectrum allows us to compare the modelled deformation of the stellar lines in assuming different one-dimensional stellar atmospheric models. We conclude that the differences among various models and observations remain within the precision limits of the data. However, the transmission light curves are better explained when the centre-to-limb variation is not included in the computation and only the Rossiter-McLaughlin deformation is considered. This demonstrates that ESPRESSO is currently the best facility for spatially resolving the stellar surface spectrum in the optical range using transit observations and carrying out empirical validations of stellar models.

preprint2021arXiv

Varying Alpha Generalized Dirac-Born-Infeld Models

We study the cosmological consequences of a class of Dirac-Born-Infeld models, and assess their viability as a candidate for the recent acceleration of the Universe. The model includes both the rolling tachyon field and the generalized Chaplygin gas models as particular limits, and phenomenologically each of these provides a possible mechanism for a deviation of the value of the dark energy equation of state from its canonical (cosmological constant) value. The field-dependent potential that is characteristic of the rolling tachyon also leads to variations of the fine-structure constant $α$, implying that the model can be constrained both by standard cosmological probes and by astrophysical measurements of $α$. Our analysis, using the latest available low-redshfit data and local constraints from atomic clock and weak equivalence principle experiments, shows that the two possible deviations of the dark energy equation of state are constrained to be $\log_{10}{(1+w_0)_V}<-7.85$ and $\log_{10}{(1+w_0)_C}<-0.85$, respectively for the rolling tachyon and Chaplygin components, both being at the $95.4\%$ confidence level (although the latter depends on the choice of priors, in a way that we quantify). Alternatively, the $95.4\%$ confidence level bound on the dimensionless slope of the potential is $\log_{10}λ<-5.36$. This confirms previous analyses indicating that in these models the potential needs to be extremely flat.

preprint2020arXiv

Abelian-Higgs Cosmic String Evolution with CUDA

Topological defects form at cosmological phase transitions by the Kibble mechanism, with cosmic strings and superstrings having the most interesting phenomenology. A rigorous analysis of their astrophysical consequences is limited by the availability of accurate numerical simulations, and therefore by hardware resources and computation time. Improving the speed and efficiency of existing codes is therefore essential. All current cosmic string simulations were performed on Central Processing Units. In previous work we presented a General Purpose Graphics Processing Unit implementation of the evolution of cosmological domain wall networks. Here we continue this paradigm shift and discuss an analogous implementation for local Abelian-Higgs strings networks. We discuss the implementation algorithm (including the discretization used and how to calculate network averaged quantities) and then showcase its performance and current bottlenecks. We validate the code by directly comparing our results for the canonical scaling properties of the networks in the radiation and matter eras with those in the literature, finding very good agreement. We finally highlight possible directions for improving the scalability of the code.

preprint2020arXiv

Atmospheric Rossiter-McLaughlin effect and transmission spectroscopy of WASP-121b with ESPRESSO

WASP-121b is one of the most studied Ultra-hot Jupiters: many recent analyses of its atmosphere report interesting features at different wavelength ranges. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in 1-telescope mode, and one partial transit taken during the commissioning of the instrument in 4-telescope mode. We investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. By analysing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter-McLaughlin effect. We measured the height of the planetary atmospheric layer that correlates with the stellar mask (mainly Fe) to be 1.052$\pm$0.015 Rp and we also confirmed the blueshift of the planetary atmosphere. By examining the planetary absorption signal on the stellar cross-correlation functions we confirmed the presence of a temporal variation of its blueshift during transit, which could be investigated spectrum-by-spectrum. We detected significant absorption in the transmission spectrum for Na, H, K, Li, CaII, Mg, and we certified their planetary nature by using the 2D tomographic technique. Particularly remarkable is the detection of Li, with a line contrast of $\sim$0.2% detected at the 6$σ$ level. With the cross-correlation technique we confirmed the presence of FeI, FeII, CrI and VI. H$α$ and CaII are present up to very high altitudes in the atmosphere ($\sim$1.44 Rp and $\sim$2 Rp, respectively), and also extend beyond the transit-equivalent Roche lobe radius of the planet. These layers of the atmosphere have a large line broadening that is not compatible with being caused by the tidally-locked rotation of the planet alone, and could arise from vertical winds or high-altitude jets in the evaporating atmosphere.

preprint2020arXiv

Broadband transmission spectroscopy of HD209458b with ESPRESSO: Evidence for Na, TiO, or both

The detection and characterization of exoplanet atmospheres is currently one of the main drivers pushing the development of new observing facilities. In this context, high-resolution spectrographs are {proving} their potential and showing that high-resolution spectroscopy will be paramount in this field. We aim to make use of ESPRESSO high-resolution spectra, which cover two transits of HD209458b, to probe the broadband transmission optical spectrum of the planet. We applied the chromatic Rossiter-McLaughin method to derive the transmission spectrum of HD209458b. We compared the results with previous HST observations and with synthetic spectra. We recover a transmission spectrum of HD209458b similar to the one obtained with HST data. The models suggest that the observed signal can be explained by only Na, only TiO, or both Na and TiO, even though none is fully capable of explaining our observed transmission spectrum. Extra absorbers may be needed to explain the full dataset, though modeling approximations and observational errors can also be responsible for the observed mismatch. Using the chromatic Rossiter-McLaughlin technique, ESPRESSO is able to provide broadband transmission spectra of exoplanets from the ground, in conjunction with space-based facilities, opening good perspectives for similar studies of other planets.

preprint2020arXiv

Distinguishing freezing and thawing dark energy models through measurements of the fine-structure constant

Mapping the behaviour of dark energy is a pressing task for observational cosmology. Phenomenological classification divides dynamical dark energy models into freezing and thawing, depending on whether the dark energy equation of state is approaching or moving away from $w=p/ρ=-1$. Moreover, in realistic dynamical dark energy models the dynamical degree of freedom is expected to couple to the electromagnetic sector, leading to variations of the fine-structure constant $α$. We discuss the feasibility of distinguishing between the freezing and thawing classes of models with current and forthcoming observational facilities and using a parametrisation of the dark energy equation of state, which can have either behaviour, introduced by Mukhanov as fiducial paradigm. We illustrate how freezing and thawing models lead to different redshift dependencies of $α$, and use a combination of current astrophysical observations and local experiments to constrain this class of models, improving the constraints on the key coupling parameter by more than a factor of two, despite considering a more extended parameter space than the one used in previous studies. We also briefly discuss the improvements expected from future facilities and comment on the practical limitations of this class of parametrisations. In particular, we show that sufficiently sensitive data can distinguish between freezing and thawing models, at least if one assumes that the relevant parameter space does not include phantom dark energy models.

preprint2020arXiv

ESPRESSO high resolution transmission spectroscopy of WASP-76b

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b.

preprint2020arXiv

Generalised velocity-dependent one-scale model for current-carrying strings

We develop an analytic model to quantitatively describe the evolution of superconducting cosmic string networks. Specifically, we extend the velocity-dependent one-scale (VOS) model to incorporate arbitrary currents and charges on cosmic string worldsheets under two main assumptions, the validity of which we also discuss. We derive equations that describe the string network evolution in terms of four macroscopic parameters: the mean string separation (or alternatively the string correlation length) and the root mean square (RMS) velocity which are the cornerstones of the VOS model, together with parameters describing the averaged timelike and spacelike current contributions. We show that our extended description reproduces the particular cases of wiggly and chiral cosmic strings, previously studied in the literature. This VOS model enables investigation of the evolution and possible observational signatures of superconducting cosmic string networks for more general equations of state, and these opportunities will be exploited in a companion paper.

preprint2020arXiv

Nightside condensation of iron in an ultra-hot giant exoplanet

Ultra-hot giant exoplanets receive thousands of times Earth&#39;s insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night (&#34;evening&#34;) and night-to-day (&#34;morning&#34;) terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

preprint2020arXiv

Primordial nucleosynthesis with varying fundamental constants: Degeneracies with cosmological parameters

The success of primordial nucleosynthesis as a cornerstone of the hot Big Bang model has been limited by the long-standing lithium problem. Recent work presented a self-consistent perturbative analysis of the effects of variations in nature&#39;s fundamental constants on primordial nucleosynthesis for a broad class of grand unified theory models, showing that such models provide a possible solution to the lithium problem, provided the value of the fine-structure constant $α$ at the nucleosynthesis epoch is larger than the current laboratory one by a few parts per million of relative variation. Here we extend the earlier analysis, focusing on how this preferred value of $α$ is affected if relevant cosmological parameters are also allowed to vary--specifically focusing on the baryon-to-photon ratio, the number of neutrinos, and the neutron lifetime. We rephrase the lithium problem in terms of the values of these parameters that would be needed to solve it within this class of grand unified theories, thus obtaining values that would disagree with the results of other experiments by several standard deviations. Using these experimental results as priors in the analysis, we find that a larger value of $α$ is still preferred, confirming our previous results. By excluding lithium from the analysis, we also obtain upper limits on possible variations of $α$ at the primordial nucleosynthesis epoch. At the two-sigma level, these are $|Δα/α|<50$ ppm without nuclear physics, cosmology, or atomic clocks priors, or alternatively $|Δα/α|<5$ ppm if these priors are used. While the simplest solution to the lithium problem is likely to be found within observational astrophysics, our work shows that varying fundamental constants remain a viable alternative.

preprint2020arXiv

Primordial nucleosynthesis with varying fundamental constants: Improved constraints and a possible solution to the Lithium problem

Primordial nucleosynthesis is an observational cornerstone of the Hot Big Bang model and a sensitive probe of physics beyond the standard model. Its success has been limited by the so-called Lithium problem, for which many solutions have been proposed. We report on a self-consistent perturbative analysis of the effects of variations in nature&#39;s fundamental constants, which are unavoidable in most extensions of the standard model, on primordial nucleosynthesis, focusing on a broad class of Grand Unified Theory models. A statistical comparison between theoretical predictions and observational measurements of ${}^4$He, D, ${}^3$He and, ${}^7$Li consistently yields a preferred value of the fine-structure constant $α$ at the nucleosynthesis epoch that is larger than the current laboratory one. The level of statistical significance and the preferred extent of variation depend on model assumptions but the former can be more than four standard deviations, while the latter is always compatible with constraints at lower redshifts. If Lithium is not included in the analysis, the preference for a variation of $α$ is not statistically significant. The abundance of ${}^3$He is relatively insensitive to such variations. Our analysis highlights a viable and physically motivated solution to the Lithium problem, which warrants further study.

preprint2020arXiv

Prospects for Fundamental Physics with LISA

In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a &#34;science-first&#34; approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.

preprint2020arXiv

Quantifying the effect of cooled initial conditions on cosmic string network evolution

Quantitative studies of the evolution and cosmological consequences of networks of cosmic strings (or other topological defects) require a combination of numerical simulations and analytic modeling with the velocity-dependent one-scale (VOS) model. In previous work, we demonstrated that a GPU-accelerated code for local Abelian-Higgs string networks enables a statistical separation of key dynamical processes affecting the evolution of the string networks and thus a precise calibration of the VOS model. Here we further exploit this code in a detailed study of two important aspects connecting the simulations with the VOS model. First, we study the sensitivity of the model calibration to the presence (or absence) of thermal oscillations due to high gradients in the initial conditions. This is relevant since in some Abelian-Higgs simulations described in the literature a period of artificial (unphysical) dissipation---usually known as cooling---is introduced with the goal of suppressing these oscillations and accelerating the convergence to scaling. We show that a small amount of cooling has no statistically significant impact on the VOS model calibration, while a longer dissipation period does have a noticeable effect. Second, in doing this analysis we also introduce an improved Markov Chain Monte Carlo based pipeline for calibrating the VOS model, Comparison to our previous bootstrap based pipeline shows that the latter accurately determined the best-fit values of the VOS model parameter, but underestimated the uncertainties in some of the parameters. Overall, our analysis shows that the calibration pipeline is robust and can be applied to future much larger field theory simulations.

preprint2020arXiv

Revisiting Proxima with ESPRESSO

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

preprint2020arXiv

WASP-127b: A misaligned planet with a partly cloudy atmosphere and tenuous sodium signature seen by ESPRESSO

The study of exoplanet atmospheres is essential to understand the formation, evolution and composition of exoplanets. The transmission spectroscopy technique is playing a significant role in this domain. In particular, the combination of state-of-the-art spectrographs at low- and high-spectral resolution is key to our understanding of atmospheric structure and composition. Two transits of the close-in sub Saturn-mass planet,WASP-127b, have been observed with ESPRESSO in the frame of the Guaranteed Time Observations Consortium. Transit observations allow us to study simultaneously the system architecture and the exoplanet atmosphere. We found that this planet is orbiting its slowly rotating host star (veq sin(i)=0.53+/-0.07 km/s) on a retrograde misaligned orbit (lambda=-128.41+/-5.60 deg). We detected the sodium line core at the 9-sigma confidence level with an excess absorption of 0.3+/-0.04%, a blueshift of 2.7+/-0.79 km/s and a FWHM of 15.18+/-1.75 km/s. However, we did not detect the presence of other atomic species but set upper-limits of only few scale heights. Finally, we put a 3-sigma upper limit, to the average depth of the 1600 strongest water lines at equilibrium temperature in the visible band, of 38 ppm. This constrains the cloud-deck pressure between 0.3 and 0.5 mbar by combining our data with low-resolution data in the near-infrared and models computed for this planet. To conclude, WASP-127b, with an age of about 10 Gyr, is an unexpected exoplanet by its orbital architecture but also by the small extension of its sodium atmosphere (~7 scale heights). ESPRESSO allows us to take a step forward in the detection of weak signals, thus bringing strong constraints on the presence of clouds in exoplanet atmospheres. The framework proposed in this work can be applied to search for molecular species and study cloud-decks in other exoplanets.

preprint2019arXiv

Astrophysical and local constraints on string theory: runaway dilaton models

One of the clear predictions of string theory is the presence of a dynamical scalar partner of the spin-2 graviton, known as the dilaton. This will violate the Einstein Equivalence Principle, leading to a plethora of possibly observable consequences which is a cosmological context include dynamical dark energy and spacetime variations of nature&#39;s fundamental constants. The runaway dilaton scenario of Damour, Piazza and Veneziano is a particularly interesting class of string theory inspired models which can in principle reconcile a massless dilaton with experimental data. Here we use the latest background cosmology observations, astrophysical and laboratory tests of the stability of the fine-structure constant and local tests of the Weak Equivalence Principle to provide updated constraints on this scenario, under various simplifying assumptions. Overall we find consistency with the standard $Λ$CDM paradigm, and we improve the existing constraints on the coupling of the dilaton to baryonic matter by a factor of six, and to the dark sector by a factor of two. At the one sigma level the current data already excludes dark sector couplings of order unity, which would be their natural value.

preprint2019arXiv

Cosmological Evolution of Semilocal String Networks

Semilocal strings -- a particular limit of electroweak strings -- are an interesting example of a stable non-topological defect whose properties resemble those of their topological cousins, the Abrikosov-Nielsen-Olesen vortices. There is, however, one important difference: a network of semilocal strings will contain segments. These are &#39;dumbbells&#39; whose ends behave almost like global monopoles that are strongly attracted to one another. While closed loops of string will eventually shrink and disappear, the segments can either shrink or grow, and a cosmological network of semilocal strings will reach a scaling regime. We discuss attempts to find a &#34;thermodynamic&#34; description of the cosmological evolution and scaling of a network of semilocal strings, by analogy with well-known descriptions for cosmic strings and for monopoles. We propose a model for the time evolution of an overall lengthscale and typical velocity for the network as well as for its segments, and some supporting (preliminary) numerical evidence.

preprint2019arXiv

Extending and Calibrating the Velocity dependent One-Scale model for Cosmic Strings with One Thousand Field Theory Simulations

Understanding the evolution and cosmological consequences of topological defect networks requires a combination of analytic modeling and numerical simulations. The canonical analytic model for defect network evolution is the Velocity-dependent One-Scale (VOS) model. For the case of cosmic strings, this has so far been calibrated using small numbers of Goto-Nambu and field theory simulations, in the radiation and matter eras, as well as in Minkowski spacetime. But the model is only as good as the available simulations, and it should be extended as further simulations become available. In previous work we presented a General Purpose Graphics Processing Unit implementation of the evolution of cosmological domain wall networks, and used it to obtain an improved VOS model for domain walls. Here we continue this effort, exploiting a more recent analogous code for local Abelian-Higgs string networks. The significant gains in speed afforded by this code enabled us to carry out 1032 field theory simulations of $512^3$ size, with 43 different expansion rates. This detailed exploration of the effects of the expansion rate on the network properties in turn enables a statistical separation of various dynamical processes affecting the evolution of the network. We thus extend and accurately calibrate the VOS model for cosmic strings, including separate terms for energy losses due to loop production and scalar/gauge radiation. By comparing this newly calibrated VOS model with the analogous one for domain walls we quantitatively show that energy loss mechanisms are different for the two types of defects.