Showing 151–200 of 502 results for author: Morales, J

Weighing stars from birth to death: mass determination methods across the HRD

Authors: Aldo Serenelli, Achim Weiss, Conny Aerts, George C. Angelou, David Baroch, Nate Bastian, Paul G. Beck, Maria Bergemann, Joachim M. Bestenlehner, Ian Czekala, Nancy Elias-Rosa, Ana Escorza, Vincent Van Eylen, Diane K. Feuillet, Davide Gandolfi, Mark Gieles, Leo Girardi, Yveline Lebreton, Nicolas Lodieu, Marie Martig, Marcelo M. Miller Bertolami, Joey S. G. Mombarg, Juan Carlos Morales, Andres Moya, Benard Nsamba , et al. (9 additional authors not shown)

Abstract: The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approac… ▽ More The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,\msun$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars. △ Less

Submitted 9 April, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Invited review article for The Astronomy and Astrophysics Review. 146 pages, 16 figures, 11 tables. Accepted version by the Journal. It includes summary figure of accuracy/precision of methods for mass ranges and summary table for individual methods

The CARMENES search for exoplanets around M dwarfs. Measuring precise radial velocities in the near infrared: the example of the super-Earth CD Cet b

Authors: F. F. Bauer, M. Zechmeister, A. Kaminski, C. Rodríguez López, J. A. Caballero, M. Azzaro, O. Stahl, D. Kossakowski, A. Quirrenbach, S. Becerril Jarque, E. Rodríguez, P. J. Amado, W. Seifert, A. Reiners, S. Schäfer, I. Ribas, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, A. Hatzes, T. Henning, S. V. Jeffers, M. Kürster, M. Lafarga, D. Montes , et al. (4 additional authors not shown)

Abstract: The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and… ▽ More The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and the work flow used to derive precise Doppler measurements from the spectra. The capability of both CARMENES channels to detect small exoplanets is demonstrated with the example of the nearby M5.0 V star CD Cet (GJ 1057), around which we announce a super-Earth ($4.0\pm0.4\,M_\oplus$) companion on a 2.29 d orbit. △ Less

Submitted 3 June, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

User behavior and token adoption on ERC20

Authors: Alfredo J. Morales, Shahar Somin, Yaniv Altshuler, Alex 'Sandy' Pentland

Abstract: Cryptocurrencies and Blockchain-based technologies are disrupting all markets. While the potential of such technologies remains to be seen, there is a current need to understand emergent patterns of user behavior and token adoption in order to design future products. In this paper we analyze the social dynamics taking place during one arbitrary day on the ERC20 platform. We characterize the networ… ▽ More Cryptocurrencies and Blockchain-based technologies are disrupting all markets. While the potential of such technologies remains to be seen, there is a current need to understand emergent patterns of user behavior and token adoption in order to design future products. In this paper we analyze the social dynamics taking place during one arbitrary day on the ERC20 platform. We characterize the network of token transactions among agents. We show heterogeneous profiles of user behavior, portfolio diversity, and token adoption. While most users are specialized in transacting with a few tokens, those that have diverse portfolios are bridging across large parts of the network and may jeopardize the system stability. We believe this work to be a foundation for unveiling the usage dynamics of crypto-currencies networks. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 12 pages, 4 figures

The CARMENES search for exoplanets around M dwarfs. The He I infrared triplet lines in PHOENIX models of M2-3 V stars

Authors: D. Hintz, B. Fuhrmeister, S. Czesla, J. H. M. M. Schmitt, A. Schweitzer, E. Nagel, E. N. Johnson, J. A. Caballero, M. Zechmeister, S. V. Jeffers, A. Reiners, I. Ribas, P. J. Amado, A. Quirrenbach, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, D. Galadí-Enríquez, E. W. Guenther, P. H. Hauschildt, A. Kaminski, M. Kürster, M. Lafarga , et al. (3 additional authors not shown)

Abstract: The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believe… ▽ More The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believed to be dominant in solar-type stars. For this purpose, we use published PHOENIX models for stars of spectral types M2 V and M3 V and also compute new series of models with different levels of activity following an ansatz developed for the case of the Sun. We perform a detailed analysis of the behavior of the He I IR line within these models. We evaluate the line in relation to other chromospheric lines and also the influence of the extreme ultraviolet (EUV) radiation field. The analysis of the He I IR line strengths as a function of the respective EUV radiation field strengths suggests that the mechanism of photoionization and recombination is necessary to form the line for inactive models, while collisions start to play a role in our most active models. Moreover, the published model set, which is optimized in the ranges of the Na I D2, H$α$, and the bluest Ca II IR triplet line, gives an adequate prediction of the He I IR line for most stars of the stellar sample. Because especially the most inactive stars with weak He I IR lines are fit worst by our models, it seems that our assumption of a 100% filling factor of a single inactive component no longer holds for these stars. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A, 13 pages, 9 figures, 4 tables

Journal ref: A&A 638, A115 (2020)

About apparent superluminal drives in generic gravity theories

Authors: J. Osorio Morales, O. Santillán

Abstract: As is well known, there exists warp drives in GR, such as the Alcubierre bubbles, which achieve an apparent faster than light travel \cite{alcubierre}. A result due to Gao and Wald \cite{gaowald} suggests that such a travel is unlikely for GR with matter satisfying both the Null Energy and the Null Generic Conditions. There exists a generalization of this statement due to Galloway, that ensures th… ▽ More As is well known, there exists warp drives in GR, such as the Alcubierre bubbles, which achieve an apparent faster than light travel \cite{alcubierre}. A result due to Gao and Wald \cite{gaowald} suggests that such a travel is unlikely for GR with matter satisfying both the Null Energy and the Null Generic Conditions. There exists a generalization of this statement due to Galloway, that ensures that the Gao-Wald result is true regardless the underlying gravity model, unless there exists at least one inextendible null geodesic with achronal image in the space time (a null line). The proof of this proposition is based on techniques of causal theories, and has never been released. In the present work an independent proof of this result is presented by use of the Raychaudhuri equation, and avoiding several technical complications described along the text. Some consequences of these affirmations are discussed at last section, in particular their potential use in problems of causality. △ Less

Submitted 14 April, 2022; v1 submitted 26 April, 2020; originally announced April 2020.

Comments: Some typos corrected

A novel embedded min-max approach for feature selection in nonlinear support vector machine classification

Authors: Asunción Jiménez-Cordero, Juan Miguel Morales, Salvador Pineda

Abstract: In recent years, feature selection has become a challenging problem in several machine learning fields, such as classification problems. Support Vector Machine (SVM) is a well-known technique applied in classification tasks. Various methodologies have been proposed in the literature to select the most relevant features in SVM. Unfortunately, all of them either deal with the feature selection probl… ▽ More In recent years, feature selection has become a challenging problem in several machine learning fields, such as classification problems. Support Vector Machine (SVM) is a well-known technique applied in classification tasks. Various methodologies have been proposed in the literature to select the most relevant features in SVM. Unfortunately, all of them either deal with the feature selection problem in the linear classification setting or propose ad-hoc approaches that are difficult to implement in practice. In contrast, we propose an embedded feature selection method based on a min-max optimization problem, where a trade-off between model complexity and classification accuracy is sought. By leveraging duality theory, we equivalently reformulate the min-max problem and solve it without further ado using off-the-shelf software for nonlinear optimization. The efficiency and usefulness of our approach are tested on several benchmark data sets in terms of accuracy, number of selected features and interpretability. △ Less

Submitted 15 January, 2021; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: Published at European Journal of Operational Research

Forecasting the Price-Response of a Pool of Buildings via Homothetic Inverse Optimization

Authors: Ricardo Fernández-Blanco, Juan Miguel Morales, Salvador Pineda

Abstract: This paper focuses on the day-ahead forecasting of the aggregate power of a pool of smart buildings equipped with thermostatically-controlled loads. We first propose the modeling of the aggregate behavior of its power trajectory by using a geometric approach. Specifically, we assume that the aggregate power is a homothet of a prototype building, whose physical and technical parameters are chosen t… ▽ More This paper focuses on the day-ahead forecasting of the aggregate power of a pool of smart buildings equipped with thermostatically-controlled loads. We first propose the modeling of the aggregate behavior of its power trajectory by using a geometric approach. Specifically, we assume that the aggregate power is a homothet of a prototype building, whose physical and technical parameters are chosen to be the mean of those in the pool. This allows us to preserve the building thermal dynamics of the pool. We then apply inverse optimization to estimate the homothetic parameters with bilevel programming. The lower level characterizes the price-response of the ensemble by a set of marginal utility curves and a homothet of the prototype building, which, in turn, are inferred in the upper-level problem. The upper level minimizes the mean absolute error over a training sample. This bilevel program is transformed into a regularized nonlinear problem that is initialized with the solution given by an efficient heuristic procedure. This heuristic consists in solving two linear programs and its solution is deemed a suitable proxy for the original bilevel problem. The results have been compared to state-of-the-art methodologies. △ Less

Submitted 24 March, 2021; v1 submitted 21 April, 2020; originally announced April 2020.

Precise mass and radius of a transiting super-Earth planet orbiting the M dwarf TOI-1235: a planet in the radius gap?

Authors: P. Bluhm, R. Luque, N. Espinoza, E. Palle, J. A. Caballero, S. Dreizler, J. H. Livingston, S. Mathur, A. Quirrenbach, S. Stock, V. Van Eylen, G. Nowak, E. Lopez, Sz. Csizmadia, M. R. Zapatero Osorio, P. Schoefer, J. Lillo-Box, M. Oshagh, P. J. Amado, D. Barrado, V. J. S. Bejar, B. Cale, P. Chaturvedi, C. Cifuentes, W. D. Cochran , et al. (56 additional authors not shown)

Abstract: We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and H… ▽ More We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and HARPS-N spectrographs. A comparison of the properties derived for TOI-1235 b's with theoretical models reveals that the planet has a rocky composition, with a bulk density slightly higher than Earth's. In particular, we measure a mass of M_p = 5.9+/-0.6 M_Earth and a radius of R_p = 1.69+/-0.08 R_Earth, which together result in a density of rho_p = 6.7+1.3-1.1 g/cm3. When compared with other well-characterized exoplanetary systems, the particular combination of planetary radius and mass puts our discovery in the radius gap, a transition region between rocky planets and planets with significant atmospheric envelopes, with few known members. While the exact location of the radius gap for M dwarfs is still a matter of debate, our results constrain it to be located at around 1.7 R_Earth or larger at the insolation levels received by TOI-1235 b (~60 S_Earth), which makes it an extremely interesting object for further studies of planet formation and atmospheric evolution. △ Less

Submitted 20 June, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: A&A, in press

Journal ref: A&A 639, A132 (2020)

The CARMENES search for exoplanets around M dwarfs. A super-Earth planet orbiting HD 79211 (GJ 338 B)

Authors: E. González-Álvarez, M. R. Zapatero Osorio, J. A. Caballero, J. Sanz-Forcada, V. J. S. Béjar, L. González-Cuesta, S. Dreizler, F. F. Bauer, E. Rodríguez, L. Tal-Or, M. Zechmeister, D. Montes, M. J. López-González, I. Ribas, A. Reiners, A. Quirrenbach, P. J. Amado, G. Anglada-Escudé, M. Azzaro, M. Cortés-Contreras, A. P. Hatzes, T. Henning, S. V. Jeffers, A. Kaminski, M. Kürster , et al. (5 additional authors not shown)

Abstract: We report on radial velocity time series for two M0.0V stars, GJ338B and GJ338A, using the CARMENES spectrograph, complemented by ground-telescope photometry from Las Cumbres and Sierra Nevada observatories. We aim to explore the presence of small planets in tight orbits using the spectroscopic radial velocity technique. We obtained 159 and 70 radial velocity measurements of GJ338B and A, respecti… ▽ More We report on radial velocity time series for two M0.0V stars, GJ338B and GJ338A, using the CARMENES spectrograph, complemented by ground-telescope photometry from Las Cumbres and Sierra Nevada observatories. We aim to explore the presence of small planets in tight orbits using the spectroscopic radial velocity technique. We obtained 159 and 70 radial velocity measurements of GJ338B and A, respectively, with the CARMENES visible channel. We also compiled additional relative radial velocity measurements from the literature and a collection of astrometric data that cover 200 a of observations to solve for the binary orbit. We found dynamical masses of 0.64$\pm$0.07M$_\odot$ for GJ338B and 0.69$\pm$0.07M$_\odot$ for GJ338A. The CARMENES radial velocity periodograms show significant peaks at 16.61$\pm$0.04 d (GJ338B) and 16.3$^{+3.5}_{-1.3}$ d (GJ338A), which have counterparts at the same frequencies in CARMENES activity indicators and photometric light curves. We attribute these to stellar rotation. GJ338B shows two additional, significant signals at 8.27$\pm$0.01 and 24.45$\pm$0.02 d, with no obvious counterparts in the stellar activity indices. The former is likely the first harmonic of the star's rotation, while we ascribe the latter to the existence of a super-Earth planet with a minimum mass of 10.27$^{+1.47}_{-1.38}$$M_{\oplus}$ orbiting GJ338B. GJ338B b lies inside the inner boundary of the habitable zone around its parent star. It is one of the least massive planets ever found around any member of stellar binaries. The masses, spectral types, brightnesses, and even the rotational periods are very similar for both stars, which are likely coeval and formed from the same molecular cloud, yet they differ in the architecture of their planetary systems. △ Less

Submitted 29 March, 2020; originally announced March 2020.

Comments: All the figures were size reduced for arXiv preprint archive. astro-ph.EP astro-ph.SR

Journal ref: A&A 637, A93 (2020)

The CARMENES search for exoplanets around M dwarfs. Radial velocities and activity indicators from cross-correlation functions with weighted binary masks

Authors: M. Lafarga, I. Ribas, C. Lovis, M. Perger, M. Zechmeister, F. F. Bauer, M. Kürster, M. Cortés-Contreras, J. C. Morales, E. Herrero, A. Rosich, D. Baroch, A. Reiners, J. A. Caballer, A. Quirrenbach, P. J. Amado, J. M. Alacid, V. J. S. Béjar, S. Dreizler, A. P. Hatzes, T. Henning, S. V. Jeffers, A. Kaminski, D. Montes, S. Pedraz , et al. (2 additional authors not shown)

Abstract: For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology… ▽ More For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology to first build weighted binary masks and, second, to compute the CCF of spectral observations with these masks from which we derive radial velocities and activity indicators. These methods are implemented in a python code that is publicly available. To build the masks, we selected a large number of sharp absorption lines based on the profile of the minima present in high signal-to-noise ratio (S/N) spectrum templates built from observations of reference stars. We computed the CCFs of observed spectra and derived RVs and the following three standard activity indicators: full-width-at-half-maximum as well as contrast and bisector inverse slope.We applied our methodology to CARMENES high-resolution spectra and obtain RV and activity indicator time series of more than 300 M dwarf stars observed for the main CARMENES survey. Compared with the standard CARMENES template matching pipeline, in general we obtain more precise RVs in the cases where the template used in the standard pipeline did not have enough S/N. We also show the behaviour of the three activity indicators for the active star YZ CMi and estimate the absolute RV of the M dwarfs analysed using the CCF RVs. △ Less

Submitted 16 March, 2020; originally announced March 2020.

Comments: 24 pages, 14 figures. Accepted for publication in A&A. Code available at https://github.com/mlafarga/raccoon

Journal ref: A&A 636, A36 (2020)

The CARMENES search for exoplanets around M dwarfs. Two planets on the opposite sides of the radius gap transiting the nearby M dwarf LTT 3780

Authors: G. Nowak, R. Luque, H. Parviainen, E. Pallé, K. Molaverdikhani, V. J. S. Béjar, J. Lillo-Box, C. Rodríguez-López, J. A. Caballero, M. Zechmeister, V. M. Passegger, C. Cifuentes, A. Schweitzer, N. Narita, B. Cale, N. Espinoza, F. Murgas, D. Hidalgo, M. R. Zapatero Osorio, F. J. Pozuelos, F. J. Aceituno, P. J. Amado, K. Barkaoui, D. Barrado, F. F. Bauer , et al. (75 additional authors not shown)

Abstract: We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LTT 3780 (TOI-732). We confirm both planets and their association with LTT 3780 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrogra… ▽ More We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LTT 3780 (TOI-732). We confirm both planets and their association with LTT 3780 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LTT 3780 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short-period planet LTT 3780 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LTT 3780 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LTT 3780 c is an ideal object for atmospheric studies with the James Webb Space Telescope. △ Less

Submitted 8 October, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A

Journal ref: A&A 642, A173 (2020)

An Efficient Robust Approach to the Day-ahead Operation of an Aggregator of Electric Vehicles

Authors: Álvaro Porras, Ricardo Fernández-Blanco, Juan M. Morales, Salvador Pineda

Abstract: The growing use of electric vehicles (EVs) may hinder their integration into the electricity system as well as their efficient operation due to the intrinsic stochasticity associated with their driving patterns. In this work, we assume a profit-maximizer EV-aggregator who participates in the day-ahead electricity market. The aggregator accounts for the technical aspects of each individual EV and t… ▽ More The growing use of electric vehicles (EVs) may hinder their integration into the electricity system as well as their efficient operation due to the intrinsic stochasticity associated with their driving patterns. In this work, we assume a profit-maximizer EV-aggregator who participates in the day-ahead electricity market. The aggregator accounts for the technical aspects of each individual EV and the uncertainty in its driving patterns. We propose a hierarchical optimization approach to represent the decision-making of this aggregator. The upper level models the profit-maximizer aggregator's decisions on the EV-fleet operation, while a series of lower-level problems computes the worst-case EV availability profiles in terms of battery draining and energy exchange with the market. Then, this problem can be equivalently transformed into a mixed-integer linear single-level equivalent given the totally unimodular character of the constraint matrices of the lower-level problems and their convexity. Finally, we thoroughly analyze the benefits of the hierarchical model compared to the results from stochastic and deterministic models. △ Less

Submitted 19 June, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 11 pages, 6 figures

MSC Class: 90C90

Chaos at the rim of black hole and fuzzball shadows

Authors: Massimo Bianchi, Alfredo Grillo, Jose Francisco Morales

Abstract: We study the scattering of massless probes in the vicinity of the {\it photon-sphere} of asymptotically AdS black holes and horizon-free microstate geometries (fuzzballs). We find that these exhibit a chaotic behaviour characterised by exponentially large deviations of nearby trajectories. We compute the Lyapunov exponent $λ$ governing the exponential growth in $d$ dimensions and show that it is b… ▽ More We study the scattering of massless probes in the vicinity of the {\it photon-sphere} of asymptotically AdS black holes and horizon-free microstate geometries (fuzzballs). We find that these exhibit a chaotic behaviour characterised by exponentially large deviations of nearby trajectories. We compute the Lyapunov exponent $λ$ governing the exponential growth in $d$ dimensions and show that it is bounded from above by $λ_b = \sqrt{d{-}3}/2b_{\rm min}$ where $b_{\rm min}$ is the minimal impact parameter under which a massless particle is swallowed by the black hole or gets trapped in the fuzzball for a very long time. Moreover we observe that $λ$ is typically below the advocated bound on chaos $λ_H=2πκ_B T/\hbar$, that in turn characterises the radial fall into the horizon, but the bound is violated in a narrow window near extremality, where the photon-sphere coalesces with the horizon. Finally, we find that fuzzballs are characterised by Lyapunov exponents smaller than those of the corresponding BH's suggesting the possibility of discriminating the existence of micro-structures at horizon scales via the detection of ring-down modes with time scales $λ^{-1}$ longer than those expected for a BH of the given mass and spin. △ Less

Submitted 3 June, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: 26 pages, 9 figures

The CARMENES search for exoplanets around M dwarfs. Characterization of the nearby ultra-compact multiplanetary system YZ Ceti

Authors: S. Stock, J. Kemmer, S. Reffert, T. Trifonov, A. Kaminski, S. Dreizler, A. Quirrenbach, J. A. Caballero, A. Reiners, S. V. Jeffers, G. Anglada-Escudé, I. Ribas, P. J. Amado, D. Barrado, J. R. Barnes, F. F. Bauer, Z. M. Berdiñas, V. J. S. Béjar, G. A. L. Coleman, M. Cortés-Contreras, E. Díez-Alonso, A. J. Domínguez-Fernández, N. Espinoza, C. A. Haswell, A. Hatzes , et al. (16 additional authors not shown)

Abstract: The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since t… ▽ More The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 636, A119 (2020)

Three planets transiting the evolved star EPIC 249893012: a hot 8.8-M$_\oplus$ super-Earth and two warm 14.7 and 10.2-M$_\oplus$ sub-Neptunes

Authors: D. Hidalgo, E. Pallé, R. Alonso, D. Gandolfi, M. Fridlund, G. Nowak, R. Luque, T. Hirano, A. B. Justesen, W. D. Cochran, O. Barragan, L. Spina, F. Rodler, S. Albrecht, D. Anderson, P. Amado, E. Bryant, J. A. Caballero, J. Cabrera, Sz. Csizmadia, F. Dai, J. De Leon, H. J. Deeg, Ph. Eigmuller, M. Endl , et al. (35 additional authors not shown)

Abstract: We report the discovery of a new planetary system with three transiting planets, one super-Earth and two sub-Neptunes, that orbit EPIC\,249893012, a G8\,IV-V evolved star ($M_\star$\,=\,1.05\,$\pm$\,0.05\,$M_\odot$, $R_\star$\,=\,1.71\,$\pm$\,0.04\,$R_\odot$, $T_\mathrm{eff}$\,=5430\,$\pm$\,85\,K). The star is just leaving the main sequence. We combined \ktwo \ photometry with IRCS adaptive-optics… ▽ More We report the discovery of a new planetary system with three transiting planets, one super-Earth and two sub-Neptunes, that orbit EPIC\,249893012, a G8\,IV-V evolved star ($M_\star$\,=\,1.05\,$\pm$\,0.05\,$M_\odot$, $R_\star$\,=\,1.71\,$\pm$\,0.04\,$R_\odot$, $T_\mathrm{eff}$\,=5430\,$\pm$\,85\,K). The star is just leaving the main sequence. We combined \ktwo \ photometry with IRCS adaptive-optics imaging and HARPS, HARPS-N, and CARMENES high-precision radial velocity measurements to confirm the planetary system, determine the stellar parameters, and measure radii, masses, and densities of the three planets. With an orbital period of $3.5949^{+0.0007}_{-0.0007}$ days, a mass of $8.75^{+1.09}_{-1.08}\ M_{\oplus}$ , and a radius of $1.95^{+0.09}_{-0.08}\ R_{\oplus}$, the inner planet b is compatible with nickel-iron core and a silicate mantle ($ρ_b= 6.39^{+1.19}_{-1.04}$ g cm$^{-3}$). Planets c and d with orbital periods of $15.624^{+0.001}_{-0.001}$ and $35.747^{+0.005}_{-0.005}$ days, respectively, have masses and radii of $14.67^{+1,84}_{-1.89}\ M_{\oplus}$ and $3.67^{+0.17}_{-0.14}\ R_{\oplus}$ and $10.18^{+2.46}_{-2.42}\ M_{\oplus}$ and $3.94^{+0.13}_{-0.12}\ R_{\oplus}$, respectively, yielding a mean density of $1.62^{+0.30}_{-0.29}$ and $0.91^{+0.25}_{-0.23}$ g cm$^{-3}$, respectively. The radius of planet b lies in the transition region between rocky and gaseous planets, but its density is consistent with a rocky composition. Its semimajor axis and the corresponding photoevaporation levels to which the planet has been exposed might explain its measured density today. In contrast, the densities and semimajor axes of planets c and d suggest a very thick atmosphere. The singularity of this system, which orbits a slightly evolved star that is just leaving the main sequence, makes it a good candidate for a deeper study from a dynamical point of view. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A\& A

Journal ref: A&A 636, A89 (2020)

The CARMENES search for exoplanets around M dwarfs. Dynamical characterization of the multiple planet system GJ 1148 and prospects of habitable exomoons around GJ 1148 b

Authors: T. Trifonov, M. H. Lee, M. Kürster, Th. Henning, E. Grishin, S. Stock, J. Tjoa, J. A. Caballero, K. H. Wong, F. F. Bauer, A. Quirrenbach, M. Zechmeister, I. Ribas, S. Reffert, A. Reiners, P. J. Amado, D. Kossakowski, M. Azzaro, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, A. P. Hatzes, S. V. Jeffers, A. Kaminski, M. Lafarga , et al. (7 additional authors not shown)

Abstract: Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements fro… ▽ More Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A

Design and construction of a new detector to measure ultra-low radioactive-isotope contamination of argon

Authors: The DarkSide Collaboration, C. E. Aalseth, S. Abdelhakim, F. Acerbi, P. Agnes, R. Ajaj, I. F. M. Albuquerque, T. Alexander, A. Alici, A. K. Alton, P. Amaudruz, F. Ameli, J. Anstey, P. Antonioli, M. Arba, S. Arcelli, R. Ardito, I. J. Arnquist, P. Arpaia, D. M. Asner, A. Asunskis, M. Ave, H. O. Back, A. Barrado Olmedo, G. Batignani , et al. (306 additional authors not shown)

Abstract: Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioa… ▽ More Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $β$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $γ$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Comments: 13 pages, 8 figures. Corresponding author: E. Sánchez García

How connected is too connected? Impact of network topology on systemic risk and collapse of complex economic systems

Authors: Aymeric Vié, Alfredo J. Morales

Abstract: Economic interdependencies have become increasingly present in globalized production, financial and trade systems. While establishing interdependencies among economic agents is crucial for the production of complex products, they may also increase systemic risks due to failure propagation. It is crucial to identify how network connectivity impacts both the emergent production and risk of collapse… ▽ More Economic interdependencies have become increasingly present in globalized production, financial and trade systems. While establishing interdependencies among economic agents is crucial for the production of complex products, they may also increase systemic risks due to failure propagation. It is crucial to identify how network connectivity impacts both the emergent production and risk of collapse of economic systems. In this paper we propose a model to study the effects of network structure on the behavior of economic systems by varying the density and centralization of connections among agents. The complexity of production increases with connectivity given the combinatorial explosion of parts and products. Emergent systemic risks arise when interconnections increase vulnerabilities. Our results suggest a universal description of economic collapse given in the emergence of tipping points and phase transitions in the relationship between network structure and risk of individual failure. This relationship seems to follow a sigmoidal form in the case of increasingly denser or centralized networks. The model sheds new light on the relevance of policies for the growth of economic complexity, and highlights the trade-off between increasing the potential production of the system and its robustness to collapse. We discuss the policy implications of intervening in the organization of interconnections and system features, and stress how different network structures and node characteristics suggest different directions in order to promote complex and robust economic systems. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Study of the effects of radiation on the CMS Drift Tubes Muon Detector for the HL-LHC

Authors: G. Abbiendi, J. Alcaraz Maestre, A. Álvarez Fernández, B. Álvarez González, N. Amapane, I. Bachiller, J. M. Barcala, L. Barcellan, C. Battilana, M. Bellato, G. Bencze, M. Benettoni, N. Beni, A. Benvenuti, L. C. Blanco Ramos, A. Boletti, A. Bragagnolo, J. A. Brochero Cifuentes, V. Cafaro, A. Calderon, E. Calvo, A. Cappati, R. Carlin, C. A. Carrillo Montoya, F. R. Cavallo , et al. (118 additional authors not shown)

Abstract: The CMS drift tubes (DT) muon detector, built for withstanding the LHC expected integrated and instantaneous luminosities, will be used also in the High Luminosity LHC (HL-LHC) at a 5 times larger instantaneous luminosity and, consequently, much higher levels of radiation, reaching about 10 times the LHC integrated luminosity. Initial irradiation tests of a spare DT chamber at the CERN gamma irrad… ▽ More The CMS drift tubes (DT) muon detector, built for withstanding the LHC expected integrated and instantaneous luminosities, will be used also in the High Luminosity LHC (HL-LHC) at a 5 times larger instantaneous luminosity and, consequently, much higher levels of radiation, reaching about 10 times the LHC integrated luminosity. Initial irradiation tests of a spare DT chamber at the CERN gamma irradiation facility (GIF++), at large ($\sim$O(100)) acceleration factor, showed ageing effects resulting in a degradation of the DT cell performance. However, full CMS simulations have shown almost no impact in the muon reconstruction efficiency over the full barrel acceptance and for the full integrated luminosity. A second spare DT chamber was moved inside the GIF++ bunker in October 2017. The chamber was being irradiated at lower acceleration factors, and only 2 out of the 12 layers of the chamber were switched at working voltage when the radioactive source was active, being the other layers in standby. In this way the other non-aged layers are used as reference and as a precise and unbiased telescope of muon tracks for the efficiency computation of the aged layers of the chamber, when set at working voltage for measurements. An integrated dose equivalent to two times the expected integrated luminosity of the HL-LHC run has been absorbed by this second spare DT chamber and the final impact on the muon reconstruction efficiency is under study. Direct inspection of some extracted aged anode wires presented a melted resistive deposition of materials. Investigation on the outgassing of cell materials and of the gas components used at the GIF++ are underway. Strategies to mitigate the ageing effects are also being developed. From the long irradiation measurements of the second spare DT chamber, the effects of radiation in the performance of the DTs expected during the HL-LHC run will be presented. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 10 pages, 7 figures, to be published in JINST, editor I. González Caballero

Report number: CMS CR-2019/159

Journal ref: JINST 14 C12010 (2019)

Segregated interactions in urban and online space

Authors: Xiaowen Dong, Alfredo J. Morales, Eaman Jahani, Esteban Moro, Bruno Lepri, Burcin Bozkaya, Carlos Sarraute, Yaneer Bar-Yam, Alex Pentland

Abstract: Urban income segregation is a widespread phenomenon that challenges societies across the globe. Classical studies on segregation have largely focused on the geographic distribution of residential neighborhoods rather than on patterns of social behaviors and interactions. In this study, we analyze segregation in economic and social interactions by observing credit card transactions and Twitter ment… ▽ More Urban income segregation is a widespread phenomenon that challenges societies across the globe. Classical studies on segregation have largely focused on the geographic distribution of residential neighborhoods rather than on patterns of social behaviors and interactions. In this study, we analyze segregation in economic and social interactions by observing credit card transactions and Twitter mentions among thousands of individuals in three culturally different metropolitan areas. We show that segregated interaction is amplified relative to the expected effects of geographic segregation in terms of both purchase activity and online communication. Furthermore, we find that segregation increases with difference in socio-economic status but is asymmetric for purchase activity, i.e., the amount of interaction from poorer to wealthier neighborhoods is larger than vice versa. Our results provide novel insights into the understanding of behavioral segregation in human interactions with significant socio-political and economic implications. △ Less

Submitted 19 April, 2020; v1 submitted 10 November, 2019; originally announced November 2019.

A giant exoplanet orbiting a very low-mass star challenges planet formation models

Authors: J. C. Morales, A. J. Mustill, I. Ribas, M. B. Davies, A. Reiners, F. F. Bauer, D. Kossakowski, E. Herrero, E. Rodríguez, M. J. López-González, C. Rodríguez-López, V. J. S. Béjar, L. González-Cuesta, R. Luque, E. Pallé, M. Perger, D. Baroch, A. Johansen, H. Klahr, C. Mordasini, G. Anglada-Escudé, J. A. Caballero, M. Cortés-Contreras, S. Dreizler, M. Lafarga , et al. (157 additional authors not shown)

Abstract: Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an… ▽ More Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: Manuscript author version. 41 pages, 11 figures

Segregation Dynamics with Reinforcement Learning and Agent Based Modeling

Authors: Egemen Sert, Yaneer Bar-Yam, Alfredo J. Morales

Abstract: Societies are complex. Properties of social systems can be explained by the interplay and weaving of individual actions. Incentives are key to understand people's choices and decisions. For instance, individual preferences of where to live may lead to the emergence of social segregation. In this paper, we combine Reinforcement Learning (RL) with Agent Based Models (ABM) in order to address the sel… ▽ More Societies are complex. Properties of social systems can be explained by the interplay and weaving of individual actions. Incentives are key to understand people's choices and decisions. For instance, individual preferences of where to live may lead to the emergence of social segregation. In this paper, we combine Reinforcement Learning (RL) with Agent Based Models (ABM) in order to address the self-organizing dynamics of social segregation and explore the space of possibilities that emerge from considering different types of incentives. Our model promotes the creation of interdependencies and interactions among multiple agents of two different kinds that want to segregate from each other. For this purpose, agents use Deep Q-Networks to make decisions based on the rules of the Schelling Segregation model and the Predator-Prey model. Despite the segregation incentive, our experiments show that spatial integration can be achieved by establishing interdependencies among agents of different kinds. They also reveal that segregated areas are more probable to host older people than diverse areas, which attract younger ones. Through this work, we show that the combination of RL and ABMs can create an artificial environment for policy makers to observe potential and existing behaviors associated to incentives. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 14 pages, 4 figures + supplemental material, in review

Water vapor detection in the transmission spectra of HD 209458 b with the CARMENES NIR channel

Authors: A. Sánchez-López, F. J. Alonso-Floriano, M. López-Puertas, I. A. G. Snellen, B. Funke, E. Nagel, F. F. Bauer, P. J. Amado, J. A. Caballero, S. Czesla, L. Nortmann, E. Pallé, M. Salz, A. Reiners, I. Ribas, A. Quirrenbach, G. Anglada-Escudé, V. J. S. Béjar, N. Casasayas-Barris, D. Galadí-Enríquez, E. W. Guenther, Th. Henning, A. Kaminski, M. Kürster, M. Lampón , et al. (9 additional authors not shown)

Abstract: Aims: We aim at detecting H$_2$O in the atmosphere of the hot Jupiter HD 209458 b and perform a multi-band study in the near infrared with CARMENES. Methods: The H$_2$O absorption lines from the planet's atmosphere are Doppler-shifted due to the large change in its radial velocity during transit. This shift is of the order of tens of km s$^{-1}$, whilst the Earth's telluric and the stellar lines… ▽ More Aims: We aim at detecting H$_2$O in the atmosphere of the hot Jupiter HD 209458 b and perform a multi-band study in the near infrared with CARMENES. Methods: The H$_2$O absorption lines from the planet's atmosphere are Doppler-shifted due to the large change in its radial velocity during transit. This shift is of the order of tens of km s$^{-1}$, whilst the Earth's telluric and the stellar lines can be considered quasi-static. We took advantage of this to remove the telluric and stellar lines using SYSREM, a principal component analysis algorithm. The residual spectra contain the signal from thousands of planetary molecular lines well below the noise level. We retrieve this information by cross-correlating the spectra with models of the atmospheric absorption. Results: We find evidence of H$_2$O in HD 209458 b with a signal-to-noise ratio (S/N) of 6.4. The signal is blueshifted by --5.2 $^{+2.6}_{-1.3}$ km s$^{-1}$, which, despite the error bars, is a firm indication of day-to-night winds at the terminator of this hot Jupiter. Additionally, we performed a multi-band study for the detection of H$_2$O individually from the three NIR bands covered by CARMENES. We detect H$_2$O from its 1.0 $μ$m band with a S/N of 5.8, and also find hints from the 1.15 $μ$m band, with a low S/N of 2.8. No clear planetary signal is found from the 1.4 $μ$m band. Conclusions: Our significant signal from the 1.0 $μ$m band in HD 209458 b represents the first detection of H$_2$O from this band, the bluest one to date. The unfavorable observational conditions might be the reason for the inconclusive detection from the stronger 1.15 and 1.4 $μ$m bands. H$_2$O is detected from the 1.0 $μ$m band in HD 209458 b, but hardly in HD 189733 b, which supports a stronger aerosol extinction in the latter. △ Less

Submitted 23 August, 2019; originally announced August 2019.

Comments: 11 pages, 10 figures; accepted for publication in A&A

Journal ref: A&A 630, A53 (2019)

On the trend to global equilibrium for Kuramoto Oscillators

Authors: Javier Morales, David Poyato

Abstract: In this paper, we study the convergence to the stable equilibrium for Kuramoto oscillators. Specifically, we derive estimates on the rate of convergence to the global equilibrium for solutions of the Kuramoto-Sakaguchi equation in a large coupling strength regime from generic initial data. As a by-product, using the stability of the equation in the Wasserstein distance, we quantify the rate at whi… ▽ More In this paper, we study the convergence to the stable equilibrium for Kuramoto oscillators. Specifically, we derive estimates on the rate of convergence to the global equilibrium for solutions of the Kuramoto-Sakaguchi equation in a large coupling strength regime from generic initial data. As a by-product, using the stability of the equation in the Wasserstein distance, we quantify the rate at which discrete Kuramoto oscillators concentrate around the global equilibrium. In doing this, we achieve a quantitative estimate in which the probability that the oscillators will concentrate at the given rate tends to one as the number of oscillators increases. Among the essential steps in our proof are: 1) An entropy production estimate inspired by the formal Riemannian structure of the space of probability measures, first introduced by F. Otto in [35]; 2) A new quantitative estimate on the instability of equilibria with antipodal oscillators based on the dynamics of norms of the solution in sets evolving by the continuity equation; 3) The use of generalized local logarithmic Sobolev and Talagrand type inequalities, similar to the ones derived by F. Otto and C. Villani in [36]; 4) The study of a system of coupled differential inequalities, by a treatment inspired by the work of L. Desvillettes and C. Villani [13]. Since the Kuramoto-Sakaguchi equation is not a gradient flow with respect to the Wasserstein distance, we derive such inequalities under a suitable fibered transportation distance. △ Less

Submitted 20 August, 2019; originally announced August 2019.

MSC Class: 34C15; 34D06; 35B40; 35Q70; 35Q83; 70F99; 92B20; 92B25

Journal ref: Ann. Inst. H. Poincaré Anal. Non Linéaire 40 (2023), no. 3, pp. 631-716

Visualization of Abel's Impossibility Theorem

Authors: J. Morales, V. Kalicki, R. Ostrander

Abstract: In this paper we construct a visualization of the Abel's Impossibility Theorem also known as the Abel-Ruffini Theorem. Using the canvas object in JavaScript along with the p5.js library, and given any expression that uses analytic functions and radicals one can always construct closed paths such that the expression evaluated at the coefficients of a general polynomial returns to it's initial posit… ▽ More In this paper we construct a visualization of the Abel's Impossibility Theorem also known as the Abel-Ruffini Theorem. Using the canvas object in JavaScript along with the p5.js library, and given any expression that uses analytic functions and radicals one can always construct closed paths such that the expression evaluated at the coefficients of a general polynomial returns to it's initial position, while the roots of the polynomial undergo a non-trivial permutation. Hence, such expression does not reconstruct the roots from the coefficients. Using the visualization we begin by considering the necessity of radicals to solve second degree polynomial equations and build towards degree five polynomial equations. In eventuality our program shows that there is no formula for an arbitrary fifth degree polynomial equation that uses analytic functions, finite field operations, and radicals that reconstructs the roots of the polynomial from it's coefficients. This theorem was partially completed by Paolo Ruffini in 1799 and completed by Niels Abel in 1824. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 11 pages, 4 figures

Day-ahead Operation of an Aggregator of Electric Vehicles via Optimization under Uncertainty

Authors: Álvaro Porras, Ricardo Fernández-Blanco, Juan Miguel Morales, Salvador Pineda

Abstract: We pose the aggregator's problem as a bilevel model, where the upper level minimizes the total operation costs of the fleet of EVs, while each lower level minimizes the energy available to each vehicle for transportation given a certain charging plan. Thanks to the totally unimodular character of the constraint matrix in the lower-level problems, the model can be mathematically recast as a computa… ▽ More We pose the aggregator's problem as a bilevel model, where the upper level minimizes the total operation costs of the fleet of EVs, while each lower level minimizes the energy available to each vehicle for transportation given a certain charging plan. Thanks to the totally unimodular character of the constraint matrix in the lower-level problems, the model can be mathematically recast as a computationally efficient mixed-integer program that delivers charging schedules that are robust against the uncertain availability of the EVs. Finally, we use synthetic data from the National Household Travel Survey 2017 to analyze the behavior of the EV aggregator from both economic and technical viewpoints and compare it with the results from a deterministic approach. △ Less

Submitted 25 September, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

How Can Smart Buildings Be Price-Responsive?

Authors: Ricardo Fernández-Blanco, Juan Miguel Morales, Salvador Pineda

Abstract: The prospective participation of smart buildings in the electricity system is strongly related to the increasing active role of demand-side resources in the electrical grid. In addition, the growing penetration of smart meters and recent advances on home automation technologies will spur the development of new mathematical tools to help optimize the local resources of these buildings. Within this… ▽ More The prospective participation of smart buildings in the electricity system is strongly related to the increasing active role of demand-side resources in the electrical grid. In addition, the growing penetration of smart meters and recent advances on home automation technologies will spur the development of new mathematical tools to help optimize the local resources of these buildings. Within this context, this paper first provides a comprehensive model to determine the electrical consumption of a single-zone household based on economic model predictive control. The goal of this problem is to minimize the electricity consumption cost while accounting for the heating dynamics of the building, smart home appliances, and comfort constraints. This paper then identifies and analyzes the key parameters responsible for the price-responsive behaviour of smart households. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Inverse Optimization with Kernel Regression: Application to the Power Forecasting and Bidding of a Fleet of Electric Vehicles

Authors: Ricardo Fernández-Blanco, Juan Miguel Morales, Salvador Pineda, Álvaro Porras

Abstract: This paper considers an aggregator of Electric Vehicles (EVs) who aims to learn the aggregate power of his/her fleet while also participating in the electricity market. The proposed approach is based on a data-driven inverse optimization (IO) method, which is highly nonlinear. To overcome such a caveat, we use a two-step estimation procedure which requires solving two convex programs. Both program… ▽ More This paper considers an aggregator of Electric Vehicles (EVs) who aims to learn the aggregate power of his/her fleet while also participating in the electricity market. The proposed approach is based on a data-driven inverse optimization (IO) method, which is highly nonlinear. To overcome such a caveat, we use a two-step estimation procedure which requires solving two convex programs. Both programs depend on penalty parameters that can be adjusted by using grid search. In addition, we propose the use of kernel regression to account for the nonlinear relationship between the behaviour of the pool of EVs and the explanatory variables, i.e., the past electricity prices and EV fleet's driving patterns. Unlike any other forecasting method, the proposed IO framework also allows the aggregator to derive a bid/offer curve, i.e. the tuple of price-quantity to be submitted to the electricity market, according to the market rules. We show the benefits of the proposed method against the machine-learning techniques that are reported to exhibit the best forecasting performance for this application in the technical literature. △ Less

Submitted 5 March, 2021; v1 submitted 1 August, 2019; originally announced August 2019.

Towards a General Framework for Static Cost Analysis of Parallel Logic Programs

Authors: Maximiliano Klemen, Pedro Lopez-Garcia, John P. Gallagher, Jose F. Morales, Manuel V. Hermenegildo

Abstract: The estimation and control of resource usage is now an important challenge in an increasing number of computing systems. In particular, requirements on timing and energy arise in a wide variety of applications such as internet of things, cloud computing, health, transportation, and robots. At the same time, parallel computing, with (heterogeneous) multi-core platforms in particular, has become the… ▽ More The estimation and control of resource usage is now an important challenge in an increasing number of computing systems. In particular, requirements on timing and energy arise in a wide variety of applications such as internet of things, cloud computing, health, transportation, and robots. At the same time, parallel computing, with (heterogeneous) multi-core platforms in particular, has become the dominant paradigm in computer architecture. Predicting resource usage on such platforms poses a difficult challenge. Most work on static resource analysis has focused on sequential programs, and relatively little progress has been made on the analysis of parallel programs, or more specifically on parallel logic programs. We propose a novel, general, and flexible framework for setting up cost equations/relations which can be instantiated for performing resource usage analysis of parallel logic programs for a wide range of resources, platforms and execution models. The analysis estimates both lower and upper bounds on the resource usage of a parallel program (without executing it) as functions on input data sizes. In addition, it also infers other meaningful information to better exploit and assess the potential and actual parallelism of a system. We develop a method for solving cost relations involving the max function that arise in the analysis of parallel programs. Finally, we instantiate our general framework for the analysis of logic programs with Independent And-Parallelism, report on an implementation within the CiaoPP system, and provide some experimental results. To our knowledge, this is the first approach to the cost analysis of parallel logic programs. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 19 pages, 3 tables; submitted to ICLP'19, accepted as technical communication

Report number: CLIP-1/2019.0

Computing Abstract Distances in Logic Programs

Authors: Ignacio Casso, Jose F. Morales, Pedro Lopez-Garcia, Manuel V. Hermenegildo

Abstract: Abstract interpretation is a well-established technique for performing static analyses of logic programs. However, choosing the abstract domain, widening, fixpoint, etc. that provides the best precision-cost trade-off remains an open problem. This is in a good part because of the challenges involved in measuring and comparing the precision of different analyses. We propose a new approach for measu… ▽ More Abstract interpretation is a well-established technique for performing static analyses of logic programs. However, choosing the abstract domain, widening, fixpoint, etc. that provides the best precision-cost trade-off remains an open problem. This is in a good part because of the challenges involved in measuring and comparing the precision of different analyses. We propose a new approach for measuring such precision, based on defining distances in abstract domains and extending them to distances between whole analyses of a given program, thus allowing comparing precision across different analyses. We survey and extend existing proposals for distances and metrics in lattices or abstract domains, and we propose metrics for some common domains used in logic program analysis, as well as extensions of those metrics to the space of whole program analysis. We implement those metrics within the CiaoPP framework and apply them to measure the precision of different analyses over both benchmarks and a realistic program. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 21 pages, 8 figures; submitted to ICLP'19, accepted as technical communication

Report number: CLIP-2/2019.0

Approximate Bayesian inference for a "steps and turns" continuous-time random walk observed at regular time intervals

Authors: Sofia Ruiz-Suarez, Vianey Leos-Barajas, Ignacio Alvarez-Castro, Juan M. Morales

Abstract: The study of animal movement is challenging because it is a process modulated by many factors acting at different spatial and temporal scales. Several models have been proposed which differ primarily in the temporal conceptualization, namely continuous and discrete time formulations. Naturally, animal movement occurs in continuous time but we tend to observe it at fixed time intervals. To account… ▽ More The study of animal movement is challenging because it is a process modulated by many factors acting at different spatial and temporal scales. Several models have been proposed which differ primarily in the temporal conceptualization, namely continuous and discrete time formulations. Naturally, animal movement occurs in continuous time but we tend to observe it at fixed time intervals. To account for the temporal mismatch between observations and movement decisions, we used a state-space model where movement decisions (steps and turns) are made in continuous time. The movement process is then observed at regular time intervals. As the likelihood function of this state-space model turned out to be complex to calculate yet simulating data is straightforward, we conduct inference using a few variations of Approximate Bayesian Computation (ABC). We explore the applicability of these methods as a function of the discrepancy between the temporal scale of the observations and that of the movement process in a simulation study. We demonstrate the application of this model to a real trajectory of a sheep that was reconstructed in high resolution using information from magnetometer and GPS devices. Our results suggest that accurate estimates can be obtained when the observations are less than 5 times the average time between changes in movement direction. The state-space model used here allowed us to connect the scales of the observations and movement decisions in an intuitive and easy to interpret way. Our findings underscore the idea that the time scale at which animal movement decisions are made needs to be considered when designing data collection protocols, and that sometimes high-frequency data may not be necessary to have good estimates of certain movement processes. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Jo: The Smart Journal

Authors: Vivian Li, Alon Halevy, Adi Zief-Balteriski Ph. D, Wang-Chiew Tan, George Mihaila, John Morales, Natalie Nuno, Huining Liu, Chen Chen, Xiaojuan Ma, Shani Robins Ph. D., Jessica Johnson

Abstract: We introduce Jo, a mobile application that attempts to improve user's well-being. Jo is a journaling application--users log their important moments via short texts and optionally an attached photo. Unlike a static journal, Jo analyzes these moments and helps users take action towards increased well-being. For example, Jo annotates each moment with a set of values (e.g., family, socialization, mind… ▽ More We introduce Jo, a mobile application that attempts to improve user's well-being. Jo is a journaling application--users log their important moments via short texts and optionally an attached photo. Unlike a static journal, Jo analyzes these moments and helps users take action towards increased well-being. For example, Jo annotates each moment with a set of values (e.g., family, socialization, mindfulness), thereby giving the user insights about the balance in their lives. In addition, Jo helps the user create reminders that enable them to create additional happy moments. We describe the results of fielding Jo in a study of 39 participants. The results illustrate the promise of a journaling application that provides personalized feedback, and points at further research. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: 7 pages, 3 figures

Feature-driven Improvement of Renewable Energy Forecasting and Trading

Authors: Miguel Á. Muñoz, Juan M. Morales, Salvador Pineda

Abstract: Inspired from recent insights into the common ground of machine learning, optimization and decision-making, this paper proposes an easy-to-implement, but effective procedure to enhance both the quality of renewable energy forecasts and the competitive edge of renewable energy producers in electricity markets with a dual-price settlement of imbalances. The quality and economic gains brought by the… ▽ More Inspired from recent insights into the common ground of machine learning, optimization and decision-making, this paper proposes an easy-to-implement, but effective procedure to enhance both the quality of renewable energy forecasts and the competitive edge of renewable energy producers in electricity markets with a dual-price settlement of imbalances. The quality and economic gains brought by the proposed procedure essentially stem from the utilization of valuable predictors (also known as features) in a data-driven newsvendor model that renders a computationally inexpensive linear program. We illustrate the proposed procedure and numerically assess its benefits on a realistic case study that considers the aggregate wind power production in the Danish DK1 bidding zone as the variable to be predicted and traded. Within this context, our procedure leverages, among others, spatial information in the form of wind power forecasts issued by transmission system operators (TSO) in surrounding bidding zones and publicly available in online platforms. We show that our method is able to improve the quality of the wind power forecast issued by the Danish TSO by several percentage points (when measured in terms of the mean absolute or the root mean square error) and to significantly reduce the balancing costs incurred by the wind power producer. △ Less

Submitted 16 January, 2020; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 10 pages, 6 figures

Data-Driven Screening of Network Constraints for Unit Commitment

Authors: Salvador Pineda, Juan Miguel Morales, Asunción Jiménez-Cordero

Abstract: The transmission-constrained unit commitment (TC-UC) problem is one of the most relevant problems solved by independent system operators for the daily operation of power systems. Given its computational complexity, this problem is usually not solved to global optimality for real-size power systems. In this paper, we propose a data-driven method that leverages historical information to screen out n… ▽ More The transmission-constrained unit commitment (TC-UC) problem is one of the most relevant problems solved by independent system operators for the daily operation of power systems. Given its computational complexity, this problem is usually not solved to global optimality for real-size power systems. In this paper, we propose a data-driven method that leverages historical information to screen out network constraints in the TC-UC problem. First, past data on demand and renewable generation throughout the network are used to learn the congestion status of transmission lines. Then, we infer the lines that will not become congested for upcoming operating conditions based on such learning and disregard their capacity constraints. This way, we formulate a reduced TC-UC problem that is easier to solve. Numerical results on a medium- and a large-size power system show that the proposed approach outperforms existing ones by significantly reducing the computational time while obtaining solutions that are equal or close to the one obtained with the original TC-UC problem. Furthermore, the purely data-driven method we propose can be seamlessly complemented with a constraint generation procedure to guarantee that the optimal solution to the original TC-UC problem is eventually recovered. △ Less

Submitted 24 January, 2020; v1 submitted 10 July, 2019; originally announced July 2019.

The CARMENES search for exoplanets around M dwarfs -- Photospheric parameters of target stars from high-resolution spectroscopy. II. Simultaneous multiwavelength range modeling of activity insensitive lines

Authors: V. M. Passegger, A. Schweitzer, D. Shulyak, E. Nagel, P. H. Hauschildt, A. Reiners, P. J. Amado, J. A. Caballero, M. Cortés-Contreras, A. J. Domínguez-Fernández, A. Quirrenbach, I. Ribas, M. Azzaro, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, S. Dreizler, E. W. Guenther, T. Henning, S. V. Jeffers, A. Kaminski, M. Kürster, M. Lafarga, E. L. Martín, D. Montes , et al. (3 additional authors not shown)

Abstract: We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $μ$m) and near-infrared wavelength range (0.96 - 1.71 $μ$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the targ… ▽ More We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $μ$m) and near-infrared wavelength range (0.96 - 1.71 $μ$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $χ^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Journal ref: A&A 627, A161 (2019)

The CARMENES search for exoplanets around M dwarfs. Two temperate Earth-mass planet candidates around Teegarden's Star

Authors: M. Zechmeister, S. Dreizler, I. Ribas, A. Reiners, J. A. Caballero, F. F. Bauer, V. J. S. Béjar, L. González-Cuesta, E. Herrero, S. Lalitha, M. J. López-González, R. Luque, J. C. Morales, E. Pallé, E. Rodríguez, C. Rodríguez López, L. Tal-Or, G. Anglada-Escudé, A. Quirrenbach, P. J. Amado, M. Abril, F. J. Aceituno, J. Aceituno, F. J. Alonso-Floriano, M. Ammler-von Eiff , et al. (160 additional authors not shown)

Abstract: Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of… ▽ More Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities. △ Less

Submitted 13 September, 2019; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: A&A 627, A49. 26 pages, 17 figures, 6 tables. Press release available at http://www.astro.physik.uni-goettingen.de/~zechmeister/teegarden/teegarden.html. v2: two authors and one reference added

Journal ref: A&A 627, A49 (2019)

The CARMENES search for exoplanets around M dwarfs Detection of a mini-Neptune around LSPM J2116+0234 and refinement of orbital parameters of a super-Earth around GJ 686 (BD+18 3421)

Authors: S. Lalitha, D. Baroch, J. C. Morales, V. M. Passegger, F. F. Bauer, C. Cardona Guillén, S. Dreizler, M. Oshagh, A. Reiners, I. Ribas, J. A. Caballero, A. Quirrenbach, P. J. Amado, V. J. S. Béjar, J. Colomé, M. Cortés-Contreras, D. Galadí-Enríquez, L. González-Cuesta, E. W. Guenther, H. -J. Hagen, T. Henning, E. Herrero, T. -O. Husser, S. V. Jeffers, A. Kaminski , et al. (12 additional authors not shown)

Abstract: Although M dwarfs are known for high levels of stellar activity, they are ideal targets for the search of low-mass exoplanets with the radial velocity (RV) method. We report the discovery of a planetary-mass companion around LSPM J2116+0234 (M3.0 V) and confirm the existence of a planet orbiting GJ 686 (BD+18 3421; M1.0 V). The discovery of the planet around LSPM J2116+0234 is based on CARMENES RV… ▽ More Although M dwarfs are known for high levels of stellar activity, they are ideal targets for the search of low-mass exoplanets with the radial velocity (RV) method. We report the discovery of a planetary-mass companion around LSPM J2116+0234 (M3.0 V) and confirm the existence of a planet orbiting GJ 686 (BD+18 3421; M1.0 V). The discovery of the planet around LSPM J2116+0234 is based on CARMENES RV observations in the visual and near-infrared channels. We confirm the planet orbiting around GJ 686 by analyzing the RV data spanning over two decades of observations from CARMENES VIS, HARPS-N, HARPS, and HIRES. We find planetary signals at 14.44 and 15.53 d in the RV data for LSPM J2116+0234 and GJ 686, respectively. Additionally, the RV, photometric time series, and various spectroscopic indicators show hints of variations of 42 d for LSPM J2116+0234 and 37 d for GJ 686, which we attribute to the stellar rotation periods. The orbital parameters of the planets are modeled with Keplerian fits together with correlated noise from the stellar activity. A mini-Neptune with a minimum mass of 11.8 Me orbits LSPM J2116+0234 producing an RV semi-amplitude of 6.19 m/s, while a super-Earth of mass 6.6 Me orbits GJ 686 and produces an RV semi-amplitude of 3.0 m/s. Both LSPM J2116+0234 and GJ 686 have planetary companions populating the regime of exoplanets with masses lower than 15 Me and orbital periods <20 d. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 23 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A 627, A116 (2019)

Planetary system around the nearby M dwarf GJ 357 including a transiting, hot, Earth-sized planet optimal for atmospheric characterization

Authors: R. Luque, E. Pallé, D. Kossakowski, S. Dreizler, J. Kemmer, N. Espinoza, J. Burt, G. Anglada-Escudé, V. J. S. Béjar, J. A. Caballero, K. A. Collins, K. I. Collins, M. Cortés-Contreras, E. Díez-Alonso, F. Feng, A. Hatzes, C. Hellier, T. Henning, S. V. Jeffers, L. Kaltenegger, M. Kürster, J. Madden, K. Molaverdikhani, D. Montes, N. Narita , et al. (51 additional authors not shown)

Abstract: We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HI… ▽ More We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HIRES, UVES, and HARPS also display a 3.93-day periodicity, confirming the planetary nature and leading to a planetary mass of Mb=1.84+-0.31 Me. In addition to the radial velocity signal for GJ 357 b, more periodicities are present in the data indicating the presence of two further planets in the system: GJ 357 c, with a minimum mass of Mc=3.40+-0.46 Me in a 9.12 d orbit, and GJ 357 d, with a minimum mass of Md=6.1+-1.0 Me in a 55.7 d orbit inside the habitable zone. The host is relatively inactive and exhibits a photometric rotation period of Prot=78+-2 d. GJ 357 b is to date the second closest transiting planet to the Sun, making it a prime target for further investigations such as transmission spectroscopy. Therefore, GJ 357 b represents one of the best terrestrial planets suitable for atmospheric characterization with the upcoming JWST and ground-based ELTs. △ Less

Submitted 30 July, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Journal ref: A&A 628, A39 (2019)

Magnetic fields in M dwarfs from the CARMENES survey

Authors: D. Shulyak, A. Reiners, E. Nagel, L. Tal-Or, J. A. Caballero, M. Zechmeister, V. J. S. Béjar, M. Cortés-Contreras, E. L. Martin, A. Kaminski, I. Ribas, A. Quirrenbach, P. J. Amado, G. Anglada-Escudé, F. F. Bauer, S. Dreizler, E. W. Guenther, T. Henning, S. V. Jeffers, M. Kürster, M. Lafarga, J. C. Morales, S. Pedraz

Abstract: M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity still lacks a consistent picture. In this work we measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample o… ▽ More M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity still lacks a consistent picture. In this work we measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample of 29 active M dwarfs and compare our results against stellar parameters. We use the state-of-the-art radiative transfer code to measure total magnetic flux densities from the Zeeman broadening of spectral lines and filling factors. We detect strong kG magnetic fields in all our targets. In 16 stars the magnetic fields were measured for the first time. Our measurements are consistent with the magnetic field saturation in stars with rotation periods P<4d. The analysis of the magnetic filling factors reveal two different patterns of either very smooth distribution or a more patchy one, which can be connected to the dynamo state of the stars and/or stellar mass. Our measurements extend the list of M dwarfs with strong surface magnetic fields. They also allow us to better constrain the interplay between the magnetic energy, stellar rotation, and underlying dynamo action. The high spectral resolution and observations at near-infrared wavelengths are the beneficial capabilities of the CARMENES instrument that allow us to address important questions about the stellar magnetism. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: 13 pages of main text, 14 pages of online material, 2 tables

Journal ref: A&A 626, A86 (2019)

Mobility profiles and calendars for food security and livelihoods analysis

Authors: Pedro J. Zufiria, David Pastor-Escuredo, Luis Ubeda Medina, Miguel A. Hernandez Medina, Iker Barriales Valbuena, Alfredo J. Morales, Wilfred Nkwambi, John Quinn, Paula Hidalgo Sanchis, Miguel Luengo-Oroz

Abstract: Social vulnerability is defined as the capacity of individuals and social groups to respond to any external stress placed on their livelihoods and wellbeing. Mobility and migrations are relevant when assessing vulnerability since the movements of a population reflect on their livelihoods, coping strategies and social safety nets. Although in general migration characterization is complex and open t… ▽ More Social vulnerability is defined as the capacity of individuals and social groups to respond to any external stress placed on their livelihoods and wellbeing. Mobility and migrations are relevant when assessing vulnerability since the movements of a population reflect on their livelihoods, coping strategies and social safety nets. Although in general migration characterization is complex and open to controversy, changes in mobility patterns for vulnerable population groups are likely to indicate a change in livelihoods or coping strategies. These changes can also indicate that the population groups may be exposed to new shocks; hence, monitoring of changes in mobility patterns can be a powerful early warning mechanism. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: Submitted and accepted at D4D'15 Senegal NetMob. Final version now published in PLOS ONE https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0195714

Journal ref: PloS one. 2018;13(4)

F-algebra--Rinehart Pairs and Super F-algebroids

Authors: John Alexander Cruz Morales, Javier A. Gutierrez, Alexander Torres-Gomez

Abstract: We define F-algebra--Rinehart pairs and super F-algebroids and study the connection between them. We define F-algebra--Rinehart pairs and super F-algebroids and study the connection between them. △ Less

Submitted 21 May, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: 14 pages. Version 2, the main result has not been changed; a discussion about tensor product and direct product of F-algebras has been added and the theorem about monoidal structure of the category of F-algebras has been removed

The CARMENES search for exoplanets around M dwarfs: Different roads to radii and masses of the target stars

Authors: Andreas Schweitzer, V. M. Passegger, C. Cifuentes, V. J. S. Bejar, M. Cortes-Contreras, J. A. Caballero, C. del Burgo, S. Czesla, M. Kuerster, D. Montes, M. R. Zapatero Osorio, I. Ribas, A. Reiners, A. Quirrenbach, P. J. Amado, J. Aceituno, G. Anglada-Escud, F. F. Bauer, S. Dreizler, S. V. Jeffers, E. W. Guenther, T. Henning, A. Kaminski, M. Lafarga, E. Marfil , et al. (6 additional authors not shown)

Abstract: We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we… ▽ More We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we obtained the required luminosities L from integrated broadband photometry together with the Gaia DR2 parallaxes. The mass was then determined using a mass-radius relation that we derived from eclipsing binaries known in the literature. We compared this method with three other methods: (1) We calculated the mass from the radius and the surface gravity log g, which was obtained from the same spectral analysis as $T_{\rm eff}$. (2) We used a widely used infrared mass-magnitude relation. (3) We used a Bayesian approach to infer stellar parameters from the comparison of the absolute magnitudes and colors of our targets with evolutionary models. Between spectral types M0V and M7V our radii cover the range $0.1\,R_{\normalsize\odot}<R<0.6\,R_{\normalsize\odot}$ with an error of 2-3% and our masses cover $0.09\,{\mathcal M}_{\normalsize\odot}<{\mathcal M}<0.6\,{\mathcal M}_{\normalsize\odot}$ with an error of 3-5%. We find good agreement between the masses determined with these different methods for most of our targets. Only the masses of very young objects show discrepancies. This can be well explained with the assumptions that we used for our methods. △ Less

Submitted 5 April, 2019; originally announced April 2019.

Comments: Accepted by A&A, 27 pages, 20 figures

Journal ref: A&A 625, A68 (2019)

Gliese 49: Activity evolution and detection of a super-Earth

Authors: M Perger, G Scandariato, I Ribas, J C Morales, L Affer, M Azzaro, P J Amado, G Anglada-Escudé, D Baroch, D Barrado, F F Bauer, V J S Béjar, J A Caballero, M Cortés-Contreras, M Damasso, S Dreizler, L González-Cuesta, J I González Hernández, E W Guenther, T Henning, E Herrero, S. V Jeffers, A Kaminski, M Kürster, M Lafarga , et al. (18 additional authors not shown)

Abstract: Small planets around low-mass stars often show orbital periods in a range that corresponds to the temperate zones of their host stars which are therefore of prime interest for planet searches. Surface phenomena such as spots and faculae create periodic signals in radial velocities and in observational activity tracers in the same range, so they can mimic or hide true planetary signals. We aim to d… ▽ More Small planets around low-mass stars often show orbital periods in a range that corresponds to the temperate zones of their host stars which are therefore of prime interest for planet searches. Surface phenomena such as spots and faculae create periodic signals in radial velocities and in observational activity tracers in the same range, so they can mimic or hide true planetary signals. We aim to detect Doppler signals corresponding to planetary companions, determine their most probable orbital configurations, and understand the stellar activity and its impact on different datasets. We analyze 22 years of data of the M1.5V-type star Gl49 (BD+61 195) including HARPS-N and CARMENES spectrographs, complemented by APT2 and SNO photometry. Activity indices are calculated from the observed spectra, and all datasets are analyzed with periodograms and noise models. We investigate how the variation of stellar activity imprints on our datasets. We further test the origin of the signals and investigate phase shifts between the different sets. To search for the best-fit model we maximize the likelihood function in a Markov Chain Monte Carlo approach. As a result of this study, we are able to detect the super-Earth Gl49b with a minimum mass of 5.6 Ms. It orbits its host star with a period of 13.85d at a semi-major axis of 0.090 au and we calculate an equilibrium temperature of 350 K and a transit probability of 2.0%. The contribution from the spot-dominated host star to the different datasets is complex, and includes signals from the stellar rotation at 18.86d, evolutionary time-scales of activity phenomena at 40-80d, and a long-term variation of at least four years. △ Less

Submitted 3 December, 2020; v1 submitted 12 March, 2019; originally announced March 2019.

Journal ref: A&A 624, A123 (2019)

Partition-based Distributionally Robust Optimization via Optimal Transport with Order Cone Constraints

Authors: Adrián Esteban-Pérez, Juan M. Morales

Abstract: In this paper we wish to tackle stochastic programs affected by ambiguity about the probability law that governs their uncertain parameters. Using optimal transport theory, we construct an ambiguity set that exploits the knowledge about the distribution of the uncertain parameters, which is provided by: (1) sample data and (2) a-priori information on the order among the probabilities that the true… ▽ More In this paper we wish to tackle stochastic programs affected by ambiguity about the probability law that governs their uncertain parameters. Using optimal transport theory, we construct an ambiguity set that exploits the knowledge about the distribution of the uncertain parameters, which is provided by: (1) sample data and (2) a-priori information on the order among the probabilities that the true data-generating distribution assigns to some regions of its support set. This type of order is enforced by means of order cone constraints and can encode a wide range of information on the shape of the probability distribution of the uncertain parameters such as information related to monotonicity or multi-modality. We seek decisions that are distributionally robust. In a number of practical cases, the resulting distributionally robust optimization (DRO) problem can be reformulated as a finite convex problem where the a-priori information translates into linear constraints. In addition, our method inherits the finite-sample performance guarantees of the Wasserstein-metric-based DRO approach proposed by Mohajerin Esfahani and Kuhn (Math Program 171(1-2):115-166. https://doi.org/10.1007/s10107-017-1172-1, 2018), while generalizing this and other popular DRO approaches. Finally, we have designed numerical experiments to analyze the performance of our approach with the newsvendor problem and the problem of a strategic firm competing à la Cournot in a market. △ Less

Submitted 13 June, 2021; v1 submitted 5 March, 2019; originally announced March 2019.

MSC Class: 90C15 Stochastic programming; 90C47 Minimax problems

Voyageur: An Experiential Travel Search Engine

Authors: Sara Evensen, Aaron Feng, Alon Halevy, Jinfeng Li, Vivian Li, Yuliang Li, Huining Liu, George Mihaila, John Morales, Natalie Nuno, Ekaterina Pavlovic, Wang-Chiew Tan, Xiaolan Wang

Abstract: We describe Voyageur, which is an application of experiential search to the domain of travel. Unlike traditional search engines for online services, experiential search focuses on the experiential aspects of the service under consideration. In particular, Voyageur needs to handle queries for subjective aspects of the service (e.g., quiet hotel, friendly staff) and combine these with objective attr… ▽ More We describe Voyageur, which is an application of experiential search to the domain of travel. Unlike traditional search engines for online services, experiential search focuses on the experiential aspects of the service under consideration. In particular, Voyageur needs to handle queries for subjective aspects of the service (e.g., quiet hotel, friendly staff) and combine these with objective attributes, such as price and location. Voyageur also highlights interesting facts and tips about the services the user is considering to provide them with further insights into their choices. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: Demo paper accepted to the Web Conference

First heat flux estimation in the lower divertor of WEST with embedded thermal measurements

Authors: J. Gaspar, Y. Corre, M. Firdaouss, J. Gardarein, J. Gerardin, J. Gunn, M. Houry, G. Laffont, T. Loarer, M. Missirlian, J. Morales, P. Moreau, C. Pocheau, E. Tsitrone, West Team

Abstract: The present paper deals with the surface heat flux estimation with thermocouples (TC) and fiber Bragg grating (FBG) embedded in the plasma facing components (PFC) of the WEST tokamak. A 2D heat transfer model combined with the conjugate gradient method (CGM) and the adjoint state is used to estimate the plasma heat flux deposited on the PFC. The plasma heat flux is characterized by the time evolut… ▽ More The present paper deals with the surface heat flux estimation with thermocouples (TC) and fiber Bragg grating (FBG) embedded in the plasma facing components (PFC) of the WEST tokamak. A 2D heat transfer model combined with the conjugate gradient method (CGM) and the adjoint state is used to estimate the plasma heat flux deposited on the PFC. The plasma heat flux is characterized by the time evolution of its amplitude and spatial shape on the target (heat flux decay length $λ^t_q$, power spreading in the private flux region $S^t$ and the strike point location $x_0$). As a first step, five ohmic pulses have been investigated with different magnetic configuration and divertor X-point height varying from 44 to 68 mm from the surface. Despite an outboard shift, the relative displacements of the outer strike point as well as the heat flux decay length derived from the TC/FBG systems are consistent with the magnetic equilibrium reconstruction. △ Less

Submitted 28 February, 2019; originally announced February 2019.

Comments: Fusion Engineering and Design, Elsevier, 2019

The CARMENES search for exoplanets around M dwarfs. Chromospheric modeling of M2-3 V stars with PHOENIX

Authors: D. Hintz, B. Fuhrmeister, S. Czesla, J. H. M. M. Schmitt, E. N. Johnson, A. Schweitzer, J. A. Caballero, M. Zechmeister, S. V. Jeffers, A. Reiners, I. Ribas, P. J. Amado, A. Quirrenbach, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, D. Galadí-Enríquez, E. W. Guenther, P. H. Hauschildt, A. Kaminski, M. Kürster, M. Lafarga, M. López del Fresno , et al. (4 additional authors not shown)

Abstract: Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framewo… ▽ More Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framework of the CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, exoplanet survey. The spectral comparison between observed data and models is performed in the chromospheric lines of Na I D2, H$α$, and the bluest Ca II infrared triplet line to obtain best-fit models for each star in the sample. We find that for inactive stars a single model with a VAL C-like temperature structure is sufficient to describe simultaneously all three lines adequately. Active stars are rather modeled by a combination of an inactive and an active model, also giving the filling factors of inactive and active regions. Moreover, the fitting of linear combinations on variable stars yields relationships between filling factors and activity states, indicating that more active phases are coupled to a larger portion of active regions on the surface of the star. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 16 pages + 6 pages appendix, 13 figures, 2 + 4 tables

Journal ref: A&A 623, A136 (2019)

Two-point correlators in non-conformal $\mathcal{N}=2$ gauge theories

Authors: M. Billo, F. Fucito, G. P. Korchemsky, A. Lerda, J. F. Morales

Abstract: We study the two-point correlation functions of chiral/anti-chiral operators in $N=2$ supersymmetric Yang-Mills theories on $R^4$ with gauge group SU(N) and $N_f$ massless hypermultiplets in the fundamental representation. We compute them in perturbation theory, using dimensional regularization up to two loops, and show that field-theory observables built out of dimensionless ratios of two-point r… ▽ More We study the two-point correlation functions of chiral/anti-chiral operators in $N=2$ supersymmetric Yang-Mills theories on $R^4$ with gauge group SU(N) and $N_f$ massless hypermultiplets in the fundamental representation. We compute them in perturbation theory, using dimensional regularization up to two loops, and show that field-theory observables built out of dimensionless ratios of two-point renormalized correlators on $R^4$ are in perfect agreement with the same quantities computed using localization on the four-sphere, even in the non-conformal case $N_f\not=2N$. △ Less

Submitted 28 May, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 57 pages, many figures. v2: version to appear on JHEP. Title changed, a few comments added in sections 2,3,4 to clarify some points

Report number: ROM2F/2019/02

The CARMENES search for exoplanets around M dwarfs. Activity indicators at visible and near-infrared wavelengths

Authors: P. Schöfer, S. V. Jeffers, A. Reiners, D. Shulyak, B. Fuhrmeister, E. N. Johnson, M. Zechmeister, I. Ribas, A. Quirrenbach, P. J. Amado, J. A. Caballero, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, E. W. Guenther, A. Kaminski, M. Kürster, M. Lafarga, D. Montes, J. C. Morales, S. Pedraz, L. Tal-Or

Abstract: The CARMENES survey is searching for Earth-like planets orbiting M dwarfs using the radial velocity method. Studying the stellar activity of the target stars is important to avoid false planet detections and to improve our understanding of the atmospheres of late-type stars. In this work we present measurements of activity indicators at visible and near-infrared wavelengths for 331 M dwarfs observ… ▽ More The CARMENES survey is searching for Earth-like planets orbiting M dwarfs using the radial velocity method. Studying the stellar activity of the target stars is important to avoid false planet detections and to improve our understanding of the atmospheres of late-type stars. In this work we present measurements of activity indicators at visible and near-infrared wavelengths for 331 M dwarfs observed with CARMENES. Our aim is to identify the activity indicators that are most sensitive and easiest to measure, and the correlations among these indicators. We also wish to characterise their variability. Using a spectral subtraction technique, we measured pseudo-equivalent widths of the He I D3, H$α$, He I $λ$10833 Å, and Pa$β$ lines, the Na I D doublet, and the Ca II infrared triplet, which have a chromospheric component in active M dwarfs. In addition, we measured an index of the strength of two TiO and two VO bands, which are formed in the photosphere. We also searched for periodicities in these activity indicators for all sample stars using generalised Lomb-Scargle periodograms. We find that the most slowly rotating stars of each spectral subtype have the strongest H$α$ absorption. H$α$ is correlated most strongly with He I D3, whereas Na I D and the Ca II infrared triplet are also correlated with H$α$. He I $λ$10833 Å and Pa$β$ show no clear correlations with the other indicators. The TiO bands show an activity effect that does not appear in the VO bands. We find that the relative variations of H$α$ and He I D3 are smaller for stars with higher activity levels, while this anti-correlation is weaker for Na I D and the Ca II infrared triplet, and is absent for He I $λ$10833 Å and Pa$β$. Periodic variation with the rotation period most commonly appears in the TiO bands, H$α$, and in the Ca II infrared triplet. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: 16 pages + 2 pages appendix, 11 figures, 3+1 tables, accepted for publication by A&A

Journal ref: A&A 623, A44 (2019)

The CARMENES search for exoplanets around M dwarfs - The enigmatic planetary system GJ 4276: One eccentric planet or two planets in a 2:1 resonance?

Authors: E. Nagel, S. Czesla, J. H. M. M. Schmitt, S. Dreizler, G. Anglada-Escudé, E. Rodríguez, I. Ribas, A. Reiners, A. Quirrenbach, P. J. Amado, J. A. Caballero, J. Aceituno, V. J. S. Béjar, M. Cortés-Contreras, L. González-Cuesta, E. W. Guenther, T. Henning, S. V. Jeffers, A. Kaminski, M. Kürster, M. Lafarga, M. J. López-González, D. Montes, J. C. Morales, V. M. Passegger , et al. (3 additional authors not shown)

Abstract: We report the detection of a Neptune-mass exoplanet around the M4.0 dwarf GJ 4276 (G 232-070) based on radial velocity (RV) observations obtained with the CARMENES spectrograph. The RV variations of GJ 4276 are best explained by the presence of a planetary companion that has a minimum mass of $m_{\rm b}\sin i \approx 16\, M_\oplus$ on a $P_{\rm b}=13.35$ day orbit. The analysis of the activity ind… ▽ More We report the detection of a Neptune-mass exoplanet around the M4.0 dwarf GJ 4276 (G 232-070) based on radial velocity (RV) observations obtained with the CARMENES spectrograph. The RV variations of GJ 4276 are best explained by the presence of a planetary companion that has a minimum mass of $m_{\rm b}\sin i \approx 16\, M_\oplus$ on a $P_{\rm b}=13.35$ day orbit. The analysis of the activity indicators and spectral diagnostics exclude stellar induced RV perturbations and prove the planetary interpretation of the RV signal. We show that a circular single-planet solution can be excluded by means of a likelihood ratio test. Instead, we find that the RV variations can be explained either by an eccentric orbit or interpreted as a pair of planets on circular orbits near a period ratio of 2:1. Although the eccentric single-planet solution is slightly preferred, our statistical analysis indicates that none of these two scenarios can be rejected with high confidence using the RV time series obtained so far. Based on the eccentric interpretation, we find that GJ 4276 b is the most eccentric ($e_{\rm b} = 0.37$) exoplanet around an M dwarf with such a short orbital period known today. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Comments: 15 pages, 12 figures, 3 tables, accepted for publication by A&A

Journal ref: A&A 622, A153 (2019)