Showing 1–50 of 245 results for author: Mathis, S

Stochastic excitation of waves in magnetic stars -- I. Scaling laws for the modes amplitudes

Authors: Leïla Bessila, Stéphane Mathis

Abstract: Stellar oscillations are key to unravelling stars' properties, such as their mass, radius and age. Amplitudes of acoustic modes in solar-like stars are intrinsically linked to their convective turbulent excitation source, which in turn is influenced by magnetism. In the observations of the Sun and stars, the amplitude of the modes is modulated following their magnetic activity cycles: the higher t… ▽ More Stellar oscillations are key to unravelling stars' properties, such as their mass, radius and age. Amplitudes of acoustic modes in solar-like stars are intrinsically linked to their convective turbulent excitation source, which in turn is influenced by magnetism. In the observations of the Sun and stars, the amplitude of the modes is modulated following their magnetic activity cycles: the higher the magnetic field, the lower the modes' amplitudes. When the magnetic field is strong, it can even inhibit the acoustic modes, which are not detected in a majority of solar-like stars presenting a strong magnetic activity. Magnetic fields are known to freeze convection when stronger than a critical value: an "on-off" approach is used in the literature. In this work, we investigate the impact of magnetic fields on the stochastic excitation of acoustic modes. First, we generalise the forced wave equation formalism, including the effects of magnetic fields. Second, we assess how convection is affected by magnetic fields using results from Magnetic Mixing-Length Theory. We provide the source terms of stochastic excitation, including a new magnetic source term and the Reynolds stresses. We provide scaling laws for the amplitudes of the modes, taking into account both the driving and the damping. Those scalings are based on the inverse Alfvén dimensionless parameter: the damping increases with the magnetic field and reaches a saturation threshold when the magnetic field is strong. The driving of the modes diminishes when the magnetic field becomes stronger, the turbulent convection being weaker. As expected from the observations, we find that a higher magnetic field diminishes the resulting modes amplitudes. Evaluating the inverse Alfvén number in stellar models provides a means to estimate the expected amplitudes of acoustic modes in magnetic active solar-type stars. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Astronomy & Astrophysics, 13 pages, 6 figures

Evaluating representation learning on the protein structure universe

Authors: Arian R. Jamasb, Alex Morehead, Chaitanya K. Joshi, Zuobai Zhang, Kieran Didi, Simon V. Mathis, Charles Harris, Jian Tang, Jianlin Cheng, Pietro Lio, Tom L. Blundell

Abstract: We introduce ProteinWorkshop, a comprehensive benchmark suite for representation learning on protein structures with Geometric Graph Neural Networks. We consider large-scale pre-training and downstream tasks on both experimental and predicted structures to enable the systematic evaluation of the quality of the learned structural representation and their usefulness in capturing functional relations… ▽ More We introduce ProteinWorkshop, a comprehensive benchmark suite for representation learning on protein structures with Geometric Graph Neural Networks. We consider large-scale pre-training and downstream tasks on both experimental and predicted structures to enable the systematic evaluation of the quality of the learned structural representation and their usefulness in capturing functional relationships for downstream tasks. We find that: (1) large-scale pretraining on AlphaFold structures and auxiliary tasks consistently improve the performance of both rotation-invariant and equivariant GNNs, and (2) more expressive equivariant GNNs benefit from pretraining to a greater extent compared to invariant models. We aim to establish a common ground for the machine learning and computational biology communities to rigorously compare and advance protein structure representation learning. Our open-source codebase reduces the barrier to entry for working with large protein structure datasets by providing: (1) storage-efficient dataloaders for large-scale structural databases including AlphaFoldDB and ESM Atlas, as well as (2) utilities for constructing new tasks from the entire PDB. ProteinWorkshop is available at: github.com/a-r-j/ProteinWorkshop. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: ICLR 2024

RNA-FrameFlow: Flow Matching for de novo 3D RNA Backbone Design

Authors: Rishabh Anand, Chaitanya K. Joshi, Alex Morehead, Arian R. Jamasb, Charles Harris, Simon V. Mathis, Kieran Didi, Bryan Hooi, Pietro Liò

Abstract: We introduce RNA-FrameFlow, the first generative model for 3D RNA backbone design. We build upon SE(3) flow matching for protein backbone generation and establish protocols for data preparation and evaluation to address unique challenges posed by RNA modeling. We formulate RNA structures as a set of rigid-body frames and associated loss functions which account for larger, more conformationally fle… ▽ More We introduce RNA-FrameFlow, the first generative model for 3D RNA backbone design. We build upon SE(3) flow matching for protein backbone generation and establish protocols for data preparation and evaluation to address unique challenges posed by RNA modeling. We formulate RNA structures as a set of rigid-body frames and associated loss functions which account for larger, more conformationally flexible RNA backbones (13 atoms per nucleotide) vs. proteins (4 atoms per residue). Toward tackling the lack of diversity in 3D RNA datasets, we explore training with structural clustering and cropping augmentations. Additionally, we define a suite of evaluation metrics to measure whether the generated RNA structures are globally self-consistent (via inverse folding followed by forward folding) and locally recover RNA-specific structural descriptors. The most performant version of RNA-FrameFlow generates locally realistic RNA backbones of 40-150 nucleotides, over 40% of which pass our validity criteria as measured by a self-consistency TM-score >= 0.45, at which two RNAs have the same global fold. Open-source code: https://github.com/rish-16/rna-backbone-design △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: To be presented as an Oral at ICML 2024 Structured Probabilistic Inference & Generative Modeling Workshop, and a Spotlight at ICML 2024 AI4Science Workshop

Improving Antibody Design with Force-Guided Sampling in Diffusion Models

Authors: Paulina Kulytė, Francisco Vargas, Simon Valentin Mathis, Yu Guang Wang, José Miguel Hernández-Lobato, Pietro Liò

Abstract: Antibodies, crucial for immune defense, primarily rely on complementarity-determining regions (CDRs) to bind and neutralize antigens, such as viruses. The design of these CDRs determines the antibody's affinity and specificity towards its target. Generative models, particularly denoising diffusion probabilistic models (DDPMs), have shown potential to advance the structure-based design of CDR regio… ▽ More Antibodies, crucial for immune defense, primarily rely on complementarity-determining regions (CDRs) to bind and neutralize antigens, such as viruses. The design of these CDRs determines the antibody's affinity and specificity towards its target. Generative models, particularly denoising diffusion probabilistic models (DDPMs), have shown potential to advance the structure-based design of CDR regions. However, only a limited dataset of bound antibody-antigen structures is available, and generalization to out-of-distribution interfaces remains a challenge. Physics based force-fields, which approximate atomic interactions, offer a coarse but universal source of information to better mold designs to target interfaces. Integrating this foundational information into diffusion models is, therefore, highly desirable. Here, we propose a novel approach to enhance the sampling process of diffusion models by integrating force field energy-based feedback. Our model, DiffForce, employs forces to guide the diffusion sampling process, effectively blending the two distributions. Through extensive experiments, we demonstrate that our method guides the model to sample CDRs with lower energy, enhancing both the structure and sequence of the generated antibodies. △ Less

Submitted 9 June, 2024; originally announced June 2024.

The PLATO Mission

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

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

Submitted 8 June, 2024; originally announced June 2024.

Effect of Rotation on Wave Mixing in Intermediate Mass Stars

Authors: A. Varghese, R. P. Ratnasingam, R. Vanon, P. V. F. Edelmann, S. Mathis, T. M. Rogers

Abstract: Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$,… ▽ More Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$, $2\times 10^{-5}$, $3\times 10^{-5}$, $4\times 10^{-5}$ and $1\times 10^{-4}$ rad.s$^{-1}$) and observe that the mixing decreases with decreasing Rossby number. This can be attributed to the effect of rotation on convection which influences the amplitude with which the waves are excited near the convective-radiative interface. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 13 pages, 8 figures, accepted at ApJ

DEFT: Efficient Finetuning of Conditional Diffusion Models by Learning the Generalised $h$-transform

Authors: Alexander Denker, Francisco Vargas, Shreyas Padhy, Kieran Didi, Simon Mathis, Vincent Dutordoir, Riccardo Barbano, Emile Mathieu, Urszula Julia Komorowska, Pietro Lio

Abstract: Generative modelling paradigms based on denoising diffusion processes have emerged as a leading candidate for conditional sampling in inverse problems. In many real-world applications, we often have access to large, expensively trained unconditional diffusion models, which we aim to exploit for improving conditional sampling. Most recent approaches are motivated heuristically and lack a unifying f… ▽ More Generative modelling paradigms based on denoising diffusion processes have emerged as a leading candidate for conditional sampling in inverse problems. In many real-world applications, we often have access to large, expensively trained unconditional diffusion models, which we aim to exploit for improving conditional sampling. Most recent approaches are motivated heuristically and lack a unifying framework, obscuring connections between them. Further, they often suffer from issues such as being very sensitive to hyperparameters, being expensive to train or needing access to weights hidden behind a closed API. In this work, we unify conditional training and sampling using the mathematically well-understood Doob's h-transform. This new perspective allows us to unify many existing methods under a common umbrella. Under this framework, we propose DEFT (Doob's h-transform Efficient FineTuning), a new approach for conditional generation that simply fine-tunes a very small network to quickly learn the conditional $h$-transform, while keeping the larger unconditional network unchanged. DEFT is much faster than existing baselines while achieving state-of-the-art performance across a variety of linear and non-linear benchmarks. On image reconstruction tasks, we achieve speedups of up to 1.6$\times$, while having the best perceptual quality on natural images and reconstruction performance on medical images. △ Less

Submitted 3 June, 2024; originally announced June 2024.

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

TimelinePTC: Development of a unified interface for pathways to care collection, visualization, and collaboration in first episode psychosis

Authors: Walter S. Mathis, Maria Ferrara, John Cahill, Sneha Karmani, Sümeyra N. Tayfur, Vinod Srihari

Abstract: This paper presents TimelinePTC, a web-based tool developed to improve the collection and analysis of Pathways to Care (PTC) data in first episode psychosis (FEP) research. Accurately measuring the duration of untreated psychosis (DUP) is essential for effective FEP treatment, requiring detailed understanding of the patient's journey to care. However, traditional PTC data collection methods, mainl… ▽ More This paper presents TimelinePTC, a web-based tool developed to improve the collection and analysis of Pathways to Care (PTC) data in first episode psychosis (FEP) research. Accurately measuring the duration of untreated psychosis (DUP) is essential for effective FEP treatment, requiring detailed understanding of the patient's journey to care. However, traditional PTC data collection methods, mainly manual and paper-based, are time-consuming and often fail to capture the full complexity of care pathways. TimelinePTC addresses these limitations by providing a digital platform for collaborative, real-time data entry and visualization, thereby enhancing data accuracy and collection efficiency. Initially created for the Specialized Treatment Early in Psychosis (STEP) program in New Haven, Connecticut, its design allows for straightforward adaptation to other healthcare contexts, facilitated by its open-source codebase. The tool significantly simplifies the data collection process, making it more efficient and user-friendly. It automates the conversion of collected data into a format ready for analysis, reducing manual transcription errors and saving time. By enabling more detailed and consistent data collection, TimelinePTC has the potential to improve healthcare access research, supporting the development of targeted interventions to reduce DUP and improve patient outcomes. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Tidal Dissipation in Giant Planets

Authors: Jim Fuller, Tristan Guillot, Stephane Mathis, Carl Murray

Abstract: Tidal interactions between moons and planets can have major effects on the orbits, spins, and thermal evolution of the moons. In the Saturn system, tidal dissipation in the planet transfers angular momentum from Saturn to the moons, causing them to migrate outwards. The rate of migration is determined by the mechanism of dissipation within the planet, which is closely tied to the planet's uncertai… ▽ More Tidal interactions between moons and planets can have major effects on the orbits, spins, and thermal evolution of the moons. In the Saturn system, tidal dissipation in the planet transfers angular momentum from Saturn to the moons, causing them to migrate outwards. The rate of migration is determined by the mechanism of dissipation within the planet, which is closely tied to the planet's uncertain structure. We review current knowledge of giant planet internal structure and evolution, which has improved thanks to data from the \textit{Juno} and \textit{Cassini} missions. We discuss general principles of tidal dissipation, describing both equilibrium and dynamical tides, and how dissipation can occur in a solid core or a fluid envelope. Finally, we discuss the possibility of resonance locking, whereby a moon can lock into resonance with a planetary oscillation mode, producing enhanced tidal migration relative to classical theories, and possibly explaining recent measurements of moon migration rates. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for Space Science Reviews. Chapter in the book based on the ISSI workshop "New Vision of the Saturnian System in the Context of a Highly Dissipative Saturn" (9-13 May 2022)

Reflections on Inductive Thematic Saturation as a potential metric for measuring the validity of an inductive Thematic Analysis with LLMs

Authors: Stefano De Paoli, Walter Stan Mathis

Abstract: This paper presents a set of reflections on saturation and the use of Large Language Models (LLMs) for performing Thematic Analysis (TA). The paper suggests that initial thematic saturation (ITS) could be used as a metric to assess part of the transactional validity of TA with LLM, focusing on the initial coding. The paper presents the initial coding of two datasets of different sizes, and it refl… ▽ More This paper presents a set of reflections on saturation and the use of Large Language Models (LLMs) for performing Thematic Analysis (TA). The paper suggests that initial thematic saturation (ITS) could be used as a metric to assess part of the transactional validity of TA with LLM, focusing on the initial coding. The paper presents the initial coding of two datasets of different sizes, and it reflects on how the LLM reaches some form of analytical saturation during the coding. The procedure proposed in this work leads to the creation of two codebooks, one comprising the total cumulative initial codes and the other the total unique codes. The paper proposes a metric to synthetically measure ITS using a simple mathematical calculation employing the ratio between slopes of cumulative codes and unique codes. The paper contributes to the initial body of work exploring how to perform qualitative analysis with LLMs. △ Less

Submitted 6 January, 2024; originally announced January 2024.

A framework for conditional diffusion modelling with applications in motif scaffolding for protein design

Authors: Kieran Didi, Francisco Vargas, Simon V Mathis, Vincent Dutordoir, Emile Mathieu, Urszula J Komorowska, Pietro Lio

Abstract: Many protein design applications, such as binder or enzyme design, require scaffolding a structural motif with high precision. Generative modelling paradigms based on denoising diffusion processes emerged as a leading candidate to address this motif scaffolding problem and have shown early experimental success in some cases. In the diffusion paradigm, motif scaffolding is treated as a conditional… ▽ More Many protein design applications, such as binder or enzyme design, require scaffolding a structural motif with high precision. Generative modelling paradigms based on denoising diffusion processes emerged as a leading candidate to address this motif scaffolding problem and have shown early experimental success in some cases. In the diffusion paradigm, motif scaffolding is treated as a conditional generation task, and several conditional generation protocols were proposed or imported from the Computer Vision literature. However, most of these protocols are motivated heuristically, e.g. via analogies to Langevin dynamics, and lack a unifying framework, obscuring connections between the different approaches. In this work, we unify conditional training and conditional sampling procedures under one common framework based on the mathematically well-understood Doob's h-transform. This new perspective allows us to draw connections between existing methods and propose a new variation on existing conditional training protocols. We illustrate the effectiveness of this new protocol in both, image outpainting and motif scaffolding and find that it outperforms standard methods. △ Less

Submitted 13 March, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 9 pages

A Hitchhiker's Guide to Geometric GNNs for 3D Atomic Systems

Authors: Alexandre Duval, Simon V. Mathis, Chaitanya K. Joshi, Victor Schmidt, Santiago Miret, Fragkiskos D. Malliaros, Taco Cohen, Pietro Liò, Yoshua Bengio, Michael Bronstein

Abstract: Recent advances in computational modelling of atomic systems, spanning molecules, proteins, and materials, represent them as geometric graphs with atoms embedded as nodes in 3D Euclidean space. In these graphs, the geometric attributes transform according to the inherent physical symmetries of 3D atomic systems, including rotations and translations in Euclidean space, as well as node permutations.… ▽ More Recent advances in computational modelling of atomic systems, spanning molecules, proteins, and materials, represent them as geometric graphs with atoms embedded as nodes in 3D Euclidean space. In these graphs, the geometric attributes transform according to the inherent physical symmetries of 3D atomic systems, including rotations and translations in Euclidean space, as well as node permutations. In recent years, Geometric Graph Neural Networks have emerged as the preferred machine learning architecture powering applications ranging from protein structure prediction to molecular simulations and material generation. Their specificity lies in the inductive biases they leverage - such as physical symmetries and chemical properties - to learn informative representations of these geometric graphs. In this opinionated paper, we provide a comprehensive and self-contained overview of the field of Geometric GNNs for 3D atomic systems. We cover fundamental background material and introduce a pedagogical taxonomy of Geometric GNN architectures: (1) invariant networks, (2) equivariant networks in Cartesian basis, (3) equivariant networks in spherical basis, and (4) unconstrained networks. Additionally, we outline key datasets and application areas and suggest future research directions. The objective of this work is to present a structured perspective on the field, making it accessible to newcomers and aiding practitioners in gaining an intuition for its mathematical abstractions. △ Less

Submitted 13 March, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Evaluating Zero-Shot Scoring for In Vitro Antibody Binding Prediction with Experimental Validation

Authors: Divya Nori, Simon V. Mathis, Amir Shanehsazzadeh

Abstract: The success of therapeutic antibodies relies on their ability to selectively bind antigens. AI-based antibody design protocols have shown promise in generating epitope-specific designs. Many of these protocols use an inverse folding step to generate diverse sequences given a backbone structure. Due to prohibitive screening costs, it is key to identify candidate sequences likely to bind in vitro. H… ▽ More The success of therapeutic antibodies relies on their ability to selectively bind antigens. AI-based antibody design protocols have shown promise in generating epitope-specific designs. Many of these protocols use an inverse folding step to generate diverse sequences given a backbone structure. Due to prohibitive screening costs, it is key to identify candidate sequences likely to bind in vitro. Here, we compare the efficacy of 8 common scoring paradigms based on open-source models to classify antibody designs as binders or non-binders. We evaluate these approaches on a novel surface plasmon resonance (SPR) dataset, spanning 5 antigens. Our results show that existing methods struggle to detect binders, and performance is highly variable across antigens. We find that metrics computed on flexibly docked antibody-antigen complexes are more robust, and ensembles scores are more consistent than individual metrics. We provide experimental insight to analyze current scoring techniques, highlighting that the development of robust, zero-shot filters is an important research gap. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Direct Fabrication of Atomically Defined Pores in MXenes

Authors: Matthew G. Boebinger, Dundar E. Yilmaz, Ayana Ghosh, Sudhajit Misra, Tyler S. Mathis, Sergei V. Kalinin, Stephen Jesse, Yury Gogotsi, Adri C. T. van Duin, Raymond R. Unocic

Abstract: Controlled fabrication of nanopores in atomically thin two-dimensional material offers the means to create robust membranes needed for ion transport, nanofiltration, and DNA sensing. Techniques for creating nanopores have relied upon either plasma etching or direct irradiation using electrons or ions; however, aberration-corrected scanning transmission electron microscopy (STEM) offers the advanta… ▽ More Controlled fabrication of nanopores in atomically thin two-dimensional material offers the means to create robust membranes needed for ion transport, nanofiltration, and DNA sensing. Techniques for creating nanopores have relied upon either plasma etching or direct irradiation using electrons or ions; however, aberration-corrected scanning transmission electron microscopy (STEM) offers the advantage of combining a highly energetic, sub-angstrom sized electron beam for atomic manipulation along with atomic resolution imaging. Here, we utilize a method for automated nanopore fabrication with real-time atomic visualization to enhance our mechanistic understanding of beam-induced transformations. Additionally, an electron beam simulation technique, Electron-Beam Simulator (E-BeamSim) was developed to observe the atomic movements and interactions resulting from electron beam irradiation. Using the 2D MXene Ti3C2Tx, we explore the influence of temperature on nanopore fabrication by tracking atomic transformation pathways and find that at room temperature, electron beam irradiation induces random displacement of atoms and results in a pileup of titanium atoms at the nanopore edge. This pileup was confirmed and demonstrated in E-BeamSim simulations around the small, milled area in the MXene monolayer. At elevated temperatures, the surface functional groups on MXene are effectively removed, and the mobility of atoms increases, which results in atomic transformations that lead to the selective removal of atoms layer by layer. Through controllable manufacture using e-beam milling fabrication, the production and then characterization of the fabricated defects can be better understood for future work. This work can lead to the development of defect engineering techniques within functionalized MXene layers. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Experimental and simulations on the electron beam interactions with MXene monolayers to form nanopores as a function of temperature

Hydrodynamic modelling of dynamical tides dissipation in Jupiter's interior as revealed by Juno

Authors: Hachem Dhouib, Clément Baruteau, Stéphane Mathis, Florian Debras, Aurélie Astoul, Michel Rieutord

Abstract: The Juno spacecraft has acquired exceptionally precise data on Jupiter's gravity field, offering invaluable insights into Jupiter's tidal response, interior structure, and dynamics, establishing crucial constraints. We develop a new model for calculating Jupiter's tidal response based on its latest interior model, while also examining the significance of different dissipation processes for the evo… ▽ More The Juno spacecraft has acquired exceptionally precise data on Jupiter's gravity field, offering invaluable insights into Jupiter's tidal response, interior structure, and dynamics, establishing crucial constraints. We develop a new model for calculating Jupiter's tidal response based on its latest interior model, while also examining the significance of different dissipation processes for the evolution of its system. We study the dissipation of dynamical tides in Jupiter by thermal, viscous and molecular diffusivities acting on gravito-inertial waves in stably stratified zones and inertial waves in convection ones. We solve the linearised equations for the equilibrium tide. Next, we compute the dynamical tides using linear hydrodynamical simulations based on a spectral method. The Coriolis force is fully taken into account, but the centrifugal effect is neglected. We study the dynamical tides occurring in Jupiter using internal structure models that respect Juno's constraints. We study specifically the dominant quadrupolar tidal components and our focus is on the frequency range that corresponds to the tidal frequencies associated with Jupiter's Galilean satellites. By incorporating the different dissipation mechanisms, we calculate the total dissipation and determine the imaginary part of the tidal Love number. We find a significant frequency dependence in dissipation spectra, indicating a strong relationship between dissipation and forcing frequency. Furthermore, our analysis reveals that, in the chosen parameter regime in which kinematic viscosity, thermal and molecular diffusivities are equal, the dominant mechanism contributing to dissipation is viscosity, exceeding in magnitude both thermal and chemical dissipation. We find that the presence of stably stratified zones plays an important role in explaining the high dissipation observed in Jupiter. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 21 pages, 15 figures. Accepted for publication in A&A

Stellar spectral-type (mass) dependence of the dearth of close-in planets around fast-rotating stars. Architecture of Kepler confirmed single-exoplanet systems compared to star-planet evolution models

Authors: R. A. García, C. Gourvès, A. R. G. Santos, A. Strugarek, D. Godoy-Rivera, S. Mathur, V. Delsanti, S. N. Breton, P. G. Beck, A. S. Brun, S. Mathis

Abstract: In 2013 a dearth of close-in planets around fast-rotating host stars was found using statistical tests on Kepler data. The addition of more Kepler and Transiting Exoplanet Survey Satellite (TESS) systems in 2022 filled this region of the diagram of stellar rotation period (Prot) versus the planet orbital period (Porb). We revisited the Prot extraction of Kepler planet-host stars, we classify the s… ▽ More In 2013 a dearth of close-in planets around fast-rotating host stars was found using statistical tests on Kepler data. The addition of more Kepler and Transiting Exoplanet Survey Satellite (TESS) systems in 2022 filled this region of the diagram of stellar rotation period (Prot) versus the planet orbital period (Porb). We revisited the Prot extraction of Kepler planet-host stars, we classify the stars by their spectral type, and we studied their Prot-Porb relations. We only used confirmed exoplanet systems to minimize biases. In order to learn about the physical processes at work, we used the star-planet evolution code ESPEM (French acronym for Evolution of Planetary Systems and Magnetism) to compute a realistic population synthesis of exoplanet systems and compared them with observations. Because ESPEM works with a single planet orbiting around a single main-sequence star, we limit our study to this population of Kepler observed systems filtering out binaries, evolved stars, and multi-planets. We find in both, observations and simulations, the existence of a dearth in close-in planets orbiting around fast-rotating stars, with a dependence on the stellar spectral type (F, G, and K), which is a proxy of the mass in our sample of stars. There is a change in the edge of the dearth as a function of the spectral type (and mass). It moves towards shorter Prot as temperature (and mass) increases, making the dearth look smaller. Realistic formation hypotheses included in the model and the proper treatment of tidal and magnetic migration are enough to qualitatively explain the dearth of hot planets around fast-rotating stars and the uncovered trend with spectral type. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted in A&A. 13 pages, 8 figures

Journal ref: A&A 679, L12 (2023)

Asteroseismic g-mode period spacings in strongly magnetic rotating stars

Authors: Nicholas Z. Rui, J. M. Joel Ong, Stéphane Mathis

Abstract: Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively.… ▽ More Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively. We provide a compact asymptotic formula for the g-mode period spacing, and universally find that strong magnetism decreases this period spacing substantially more than is predicted by perturbation theory. These results are validated with explicit numerical mode calculations for realistic stellar models. The approach we present is highly versatile: once the eigenvalues $λ$ of a certain differential operator are precomputed as a function of the magnetogravity and rotational frequencies (in units of the mode frequency), the non-perturbative impact of the Coriolis and Lorentz forces is understood under a broad domain of validity, and is readily incorporated into asteroseismic modeling. △ Less

Submitted 7 November, 2023; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 18 pages, 8 figures; accepted to MNRAS

Hydrodynamical modelling of tidal dissipation in gas giant planets at the time of space missions

Authors: Hachem Dhouib, Clément Baruteau, Stéphane Mathis, Florian Debras, Aurélie Astoul, Michel Rieutord

Abstract: Gas giant planets are differentially rotating magnetic objects that have strong and complex interactions with their environment. In our Solar system, they interact with their numerous moons while exoplanets with very short orbital periods (hot Jupiters), interact with their host star. The dissipation of waves excited by tidal forces in their interiors shapes the orbital architecture and the rotati… ▽ More Gas giant planets are differentially rotating magnetic objects that have strong and complex interactions with their environment. In our Solar system, they interact with their numerous moons while exoplanets with very short orbital periods (hot Jupiters), interact with their host star. The dissipation of waves excited by tidal forces in their interiors shapes the orbital architecture and the rotational dynamics of these systems. Recently, astrometric observations of Jupiter and Saturn systems have challenged our understanding of their formation and evolution, with stronger tidal dissipation in these planets than previously predicted, in contrast to what appears to be weaker in gas giant exoplanets. These new constraints are motivating the development of realistic models of tidal dissipation inside these planets. At the same time, the Juno and Cassini space missions have revolutionised our knowledge of the interiors of Jupiter and Saturn, whose structure is a combination of stably stratified zones and convective regions. In this work, we present results of hydrodynamical calculations modelling tidal waves and their dissipation in Jupiter, taking for the first time the latest, state-of-the-art interior model of the planet. We performed 2D numerical simulations of linear tidal gravito-inertial waves that propagate and dissipate within Jupiter interior by taking into account viscous, thermal and chemical diffusions. This new model allows us to explore the properties of the dissipation and the associated tidal torque as a function of all the key hydrodynamical and structural parameters. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 4 pages, 2 figures. Proceeding of the Annual meeting of the French Society of Astronomy and Astrophysics (SF2A 2023)

In search of gravity mode signatures in main sequence solar-type stars observed by Kepler

Authors: Sylvain N. Breton, Hachem Dhouib, Rafael A. García, Allan Sacha Brun, Stéphane Mathis, Fernando Pérez Hernández, Savita Mathur, Achrène Dyrek, Angela R. G. Santos, Pere L. Pallé

Abstract: Gravity modes (g modes), mixed gravito-acoustic modes (mixed modes), and gravito-inertial modes (gi modes) possess unmatched properties as probes for stars with radiative interiors. The structural and dynamical constraints that they are able to provide cannot be accessed by other means. While they provide precious insights into the internal dynamics of evolved stars as well as massive and intermed… ▽ More Gravity modes (g modes), mixed gravito-acoustic modes (mixed modes), and gravito-inertial modes (gi modes) possess unmatched properties as probes for stars with radiative interiors. The structural and dynamical constraints that they are able to provide cannot be accessed by other means. While they provide precious insights into the internal dynamics of evolved stars as well as massive and intermediate-mass stars, their non-detection in main sequence (MS) solar-type stars make them a crucial missing piece in our understanding of angular momentum transport in radiative zones and stellar rotational evolution. In this work, we aim to apply certain analysis tools originally developed for helioseismology in order to look for g-mode signatures in MS solar-type stars. We select a sample of the 34 most promising MS solar-type stars with Kepler four-year long photometric time series. All these stars are well-characterised late F-type stars with thin convective envelopes, fast convective flows, and stochastically excited acoustic modes (p modes). For each star, we compute the background noise level of the Fourier power spectrum to identify significant peaks at low frequency. After successfully detecting individual peaks in 12 targets, we further analyse four of them and observe distinct patterns of surrounding peaks with a low probability of being noise artifacts. Comparisons with the predictions from reference models suggest that these patterns are compatible with the presence of non-asymptotic low-order pure g modes, pure p modes, and mixed modes. Given their sensitivity to both the convective core interface stratification and the coupling between p- and g-mode resonant cavities, such modes are able to provide strong constraints on the structure and evolutionary states of the related targets. [abridged] △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 19 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 679, A104 (2023)

DiffHopp: A Graph Diffusion Model for Novel Drug Design via Scaffold Hopping

Authors: Jos Torge, Charles Harris, Simon V. Mathis, Pietro Lio

Abstract: Scaffold hopping is a drug discovery strategy to generate new chemical entities by modifying the core structure, the \emph{scaffold}, of a known active compound. This approach preserves the essential molecular features of the original scaffold while introducing novel chemical elements or structural features to enhance potency, selectivity, or bioavailability. However, there is currently a lack of… ▽ More Scaffold hopping is a drug discovery strategy to generate new chemical entities by modifying the core structure, the \emph{scaffold}, of a known active compound. This approach preserves the essential molecular features of the original scaffold while introducing novel chemical elements or structural features to enhance potency, selectivity, or bioavailability. However, there is currently a lack of generative models specifically tailored for this task, especially in the pocket-conditioned context. In this work, we present DiffHopp, a conditional E(3)-equivariant graph diffusion model tailored for scaffold hopping given a known protein-ligand complex. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Benchmarking Generated Poses: How Rational is Structure-based Drug Design with Generative Models?

Authors: Charles Harris, Kieran Didi, Arian R. Jamasb, Chaitanya K. Joshi, Simon V. Mathis, Pietro Lio, Tom Blundell

Abstract: Deep generative models for structure-based drug design (SBDD), where molecule generation is conditioned on a 3D protein pocket, have received considerable interest in recent years. These methods offer the promise of higher-quality molecule generation by explicitly modelling the 3D interaction between a potential drug and a protein receptor. However, previous work has primarily focused on the quali… ▽ More Deep generative models for structure-based drug design (SBDD), where molecule generation is conditioned on a 3D protein pocket, have received considerable interest in recent years. These methods offer the promise of higher-quality molecule generation by explicitly modelling the 3D interaction between a potential drug and a protein receptor. However, previous work has primarily focused on the quality of the generated molecules themselves, with limited evaluation of the 3D molecule \emph{poses} that these methods produce, with most work simply discarding the generated pose and only reporting a "corrected" pose after redocking with traditional methods. Little is known about whether generated molecules satisfy known physical constraints for binding and the extent to which redocking alters the generated interactions. We introduce PoseCheck, an extensive analysis of multiple state-of-the-art methods and find that generated molecules have significantly more physical violations and fewer key interactions compared to baselines, calling into question the implicit assumption that providing rich 3D structure information improves molecule complementarity. We make recommendations for future research tackling identified failure modes and hope our benchmark can serve as a springboard for future SBDD generative modelling work to have a real-world impact. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Mode coupling coefficients between the convective core and radiative envelope of $γ\,$Doradus and slowly pulsating B stars

Authors: Conny Aerts, Stéphane Mathis

Abstract: Signatures of coupling between an inertial mode in the convective core and a gravito-inertial mode in the envelope have been found in four-year Kepler light curves of 16 rapidly rotating $γ\,$Doradus ($γ\,$Dor) stars. This makes it possible to obtain a measurement of the rotation frequency in their convective core. Despite their similar internal structure and available data, inertial modes have no… ▽ More Signatures of coupling between an inertial mode in the convective core and a gravito-inertial mode in the envelope have been found in four-year Kepler light curves of 16 rapidly rotating $γ\,$Doradus ($γ\,$Dor) stars. This makes it possible to obtain a measurement of the rotation frequency in their convective core. Despite their similar internal structure and available data, inertial modes have not yet been reported for slowly pulsating B (SPB) stars. We aim to provide a numerical counterpart of the recently published theoretical expressions for the mode-coupling coefficients, $\varepsilon$ and $\tilde{\varepsilon}$. These coefficients represent the two cases of a continuous and a discontinuous Brunt-Väisälä frequency profile at the core-envelope interface, respectively. We used asteroseismic forward models of two samples consisting of 26 SPB stars and 37 $γ\,$Dor stars to infer their numerical values of $\varepsilon$. The asteroseismically inferred values of $\varepsilon$ for the two samples are between 0.0 and 0.34. While $\varepsilon$ is most strongly correlated with the near-core rotation frequency for $γ\,$Dor stars, the fractional radius of the convective core instead provides the tightest correlation for SPB stars. We find $\varepsilon$ to decrease mildly as the stars evolve. Our asteroseismic results for the mode coupling support the theoretical interpretation and reveal that young, fast-rotating $γ\,$Dor stars are most suitable for undergoing couplings between inertial modes in the rotating convective core and gravito-inertial modes in the radiative envelope. The phenomenon has been found in 2.4\% of such pulsators with detected period spacing patterns, whereas it has not been seen in any of the SPB stars so far. (shortened abstract to meet the arXiv limits) △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: Manuscript in press at A&A, 9 pages, 8 figures

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

Constraining stellar and orbital co-evolution through ensemble seismology of solar-like oscillators in binary systems -- A census of oscillating red-giants and main-sequence stars in Gaia DR3 binaries

Authors: P. G. Beck, D. H. Grossmann, L. Steinwender, L. S. Schimak, N. Muntean, M. Vrard, R. A. Patton, J. Merc, S. Mathur, R. A. Garcia, M. H. Pinsonneault, D. M. Rowan, P. Gaulme, C. Allende Prieto, K. Z. Arellano-Córdova, L. Cao, E. Corsaro, O. Creevey, K. M. Hambleton, A. Hanslmeier, B. Holl, J. Johnson, S. Mathis, D. Godoy-Rivera, S. Símon-Díaz , et al. (1 additional authors not shown)

Abstract: Binary systems constitute a valuable astrophysics tool for testing our understanding of stellar structure and evolution. Systems containing a oscillating component are interesting as asteroseismology offers independent parameters for the oscillating component that aid the analysis. About 150 of such systems are known in the literature. To enlarge the sample of these benchmark objects, we crossmatc… ▽ More Binary systems constitute a valuable astrophysics tool for testing our understanding of stellar structure and evolution. Systems containing a oscillating component are interesting as asteroseismology offers independent parameters for the oscillating component that aid the analysis. About 150 of such systems are known in the literature. To enlarge the sample of these benchmark objects, we crossmatch the Two-Body-Orbit Catalogue (TBO) of Gaia DR3, with catalogs of confirmed solar-like oscillators on the main-sequence and red-giant phase from NASA Kepler and TESS. We obtain 954 new binary system candidates hosting solar-like oscillators, of which 45 and 909 stars are on the main sequence and red-giant, resp., including 2 new red giants in eclipsing systems. 918 oscillators in potentially long-periodic systems are reported. We increase the sample size of known solar-like oscillators in binary systems by an order of magnitude. We present the seismic properties of the full sample and conclude that the grand majority of the orbital elements in the TBO is physically reasonable. 82% of all TBO binary candidates with multiple times with APOGEE are confirmed from radial-velocity measurement. However, we suggest that due to instrumental noise of the TESS satellite the seismically inferred masses and radii of stars with $ν_\textrm{max}$$\lesssim$30$μ$Hz could be significantly overestimated. For 146 giants the seismically inferred evolutionary state has been determined and shows clear differences in their distribution in the orbital parameters, which are accounted the accumulative effect of the equilibrium tide acting in these evolved binary systems. For other 146 systems hosting oscillating stars values for the orbital inclination were found in the TBO. From testing the TBO on the SB9 catalogue, we obtain a completeness factor of 1/3. △ Less

Submitted 6 November, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Comments: accepted for publication in Astronomy&Astrophysics (23 pages + 4 pages of appendix, 21 figures, 33 pages of tables in the Appendix)

Artificial Intelligence for Science in Quantum, Atomistic, and Continuum Systems

Authors: Xuan Zhang, Limei Wang, Jacob Helwig, Youzhi Luo, Cong Fu, Yaochen Xie, Meng Liu, Yuchao Lin, Zhao Xu, Keqiang Yan, Keir Adams, Maurice Weiler, Xiner Li, Tianfan Fu, Yucheng Wang, Haiyang Yu, YuQing Xie, Xiang Fu, Alex Strasser, Shenglong Xu, Yi Liu, Yuanqi Du, Alexandra Saxton, Hongyi Ling, Hannah Lawrence , et al. (38 additional authors not shown)

Abstract: Advances in artificial intelligence (AI) are fueling a new paradigm of discoveries in natural sciences. Today, AI has started to advance natural sciences by improving, accelerating, and enabling our understanding of natural phenomena at a wide range of spatial and temporal scales, giving rise to a new area of research known as AI for science (AI4Science). Being an emerging research paradigm, AI4Sc… ▽ More Advances in artificial intelligence (AI) are fueling a new paradigm of discoveries in natural sciences. Today, AI has started to advance natural sciences by improving, accelerating, and enabling our understanding of natural phenomena at a wide range of spatial and temporal scales, giving rise to a new area of research known as AI for science (AI4Science). Being an emerging research paradigm, AI4Science is unique in that it is an enormous and highly interdisciplinary area. Thus, a unified and technical treatment of this field is needed yet challenging. This work aims to provide a technically thorough account of a subarea of AI4Science; namely, AI for quantum, atomistic, and continuum systems. These areas aim at understanding the physical world from the subatomic (wavefunctions and electron density), atomic (molecules, proteins, materials, and interactions), to macro (fluids, climate, and subsurface) scales and form an important subarea of AI4Science. A unique advantage of focusing on these areas is that they largely share a common set of challenges, thereby allowing a unified and foundational treatment. A key common challenge is how to capture physics first principles, especially symmetries, in natural systems by deep learning methods. We provide an in-depth yet intuitive account of techniques to achieve equivariance to symmetry transformations. We also discuss other common technical challenges, including explainability, out-of-distribution generalization, knowledge transfer with foundation and large language models, and uncertainty quantification. To facilitate learning and education, we provide categorized lists of resources that we found to be useful. We strive to be thorough and unified and hope this initial effort may trigger more community interests and efforts to further advance AI4Science. △ Less

Submitted 15 November, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Asteroseismology with the Roman Galactic Bulge Time-Domain Survey

Authors: Daniel Huber, Marc Pinsonneault, Paul Beck, Timothy R. Bedding, Joss Bland-Hawthorn, Sylvain N. Breton, Lisa Bugnet, William J. Chaplin, Rafael A. Garcia, Samuel K. Grunblatt, Joyce A. Guzik, Saskia Hekker, Steven D. Kawaler, Stephane Mathis, Savita Mathur, Travis Metcalfe, Benoit Mosser, Melissa K. Ness, Anthony L. Piro, Aldo Serenelli, Sanjib Sharma, David R. Soderblom, Keivan G. Stassun, Dennis Stello, Jamie Tayar , et al. (2 additional authors not shown)

Abstract: Asteroseismology has transformed stellar astrophysics. Red giant asteroseismology is a prime example, with oscillation periods and amplitudes that are readily detectable with time-domain space-based telescopes. These oscillations can be used to infer masses, ages and radii for large numbers of stars, providing unique constraints on stellar populations in our galaxy. The cadence, duration, and spat… ▽ More Asteroseismology has transformed stellar astrophysics. Red giant asteroseismology is a prime example, with oscillation periods and amplitudes that are readily detectable with time-domain space-based telescopes. These oscillations can be used to infer masses, ages and radii for large numbers of stars, providing unique constraints on stellar populations in our galaxy. The cadence, duration, and spatial resolution of the Roman galactic bulge time-domain survey (GBTDS) are well-suited for asteroseismology and will probe an important population not studied by prior missions. We identify photometric precision as a key requirement for realizing the potential of asteroseismology with Roman. A precision of 1 mmag per 15-min cadence or better for saturated stars will enable detections of the populous red clump star population in the Galactic bulge. If the survey efficiency is better than expected, we argue for repeat observations of the same fields to improve photometric precision, or covering additional fields to expand the stellar population reach if the photometric precision for saturated stars is better than 1 mmag. Asteroseismology is relatively insensitive to the timing of the observations during the mission, and the prime red clump targets can be observed in a single 70 day campaign in any given field. Complementary stellar characterization, particularly astrometry tied to the Gaia system, will also dramatically expand the diagnostic power of asteroseismology. We also highlight synergies to Roman GBTDS exoplanet science using transits and microlensing. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Roman Core Community Survey White Paper, 3 pages, 4 figures

Predicting protein variants with equivariant graph neural networks

Authors: Antonia Boca, Simon Mathis

Abstract: Pre-trained models have been successful in many protein engineering tasks. Most notably, sequence-based models have achieved state-of-the-art performance on protein fitness prediction while structure-based models have been used experimentally to develop proteins with enhanced functions. However, there is a research gap in comparing structure- and sequence-based methods for predicting protein varia… ▽ More Pre-trained models have been successful in many protein engineering tasks. Most notably, sequence-based models have achieved state-of-the-art performance on protein fitness prediction while structure-based models have been used experimentally to develop proteins with enhanced functions. However, there is a research gap in comparing structure- and sequence-based methods for predicting protein variants that are better than the wildtype protein. This paper aims to address this gap by conducting a comparative study between the abilities of equivariant graph neural networks (EGNNs) and sequence-based approaches to identify promising amino-acid mutations. The results show that our proposed structural approach achieves a competitive performance to sequence-based methods while being trained on significantly fewer molecules. Additionally, we find that combining assay labelled data with structure pre-trained models yields similar trends as with sequence pre-trained models. Our code and trained models can be found at: https://github.com/semiluna/partIII-amino-acid-prediction. △ Less

Submitted 24 July, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 4 pages, 2 figures, accepted to the 2023 ICML Workshop on Computational Biology

Asymmetries of frequency splittings of dipolar mixed modes: a window on the topology of deep magnetic fields

Authors: Stéphane Mathis, Lisa Bugnet

Abstract: Space asteroseismology is revolutionizing our knowledge of the internal structure and dynamics of stars. A breakthrough is ongoing with the recent discoveries of signatures of strong magnetic fields in the core of red giant stars. The key signature for such a detection is the asymmetry these fields induce in the frequency splittings of observed dipolar mixed gravito-acoustic modes. We investigate… ▽ More Space asteroseismology is revolutionizing our knowledge of the internal structure and dynamics of stars. A breakthrough is ongoing with the recent discoveries of signatures of strong magnetic fields in the core of red giant stars. The key signature for such a detection is the asymmetry these fields induce in the frequency splittings of observed dipolar mixed gravito-acoustic modes. We investigate the ability of the observed asymmetries of the frequency splittings of dipolar mixed modes to constrain the geometrical properties of deep magnetic fields. We use the powerful analytical Racah-Wigner algebra used in Quantum Mechanics to characterize the geometrical couplings of dipolar mixed oscillation modes with various possible realistic fossil magnetic fields' topologies and compute the induced perturbation of their frequencies. First, in the case of an oblique magnetic dipole, we provide the exact analytical expression of the asymmetry as a function of the angle between the rotation and magnetic axes. Its value provides a direct measure of this angle. Second, considering a combination of axisymmetric dipolar and quadrupolar fields, we show how the asymmetry is blind to unravel the relative strength and sign of each component. Finally, in the case of a given multipole, we show that a negative asymmetry is a signature of non-axisymmetric topologies. Therefore, asymmetries of dipolar mixed modes provide key but only partial information on the geometrical topology of deep fossil magnetic fields. Asteroseismic constraints should therefore be combined with spectropolarimetric observations and numerical simulations, which aim to predict the more probable stable large-scale geometries. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 10 pages, 3 figures, Letter accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 676, L9 (2023)

TOI-1416: A system with a super-Earth planet with a 1.07d period

Authors: H. J. Deeg, I. Y. Georgieva, G. Nowak, C. M. Persson, B. L. Cale, F. Murgas, E. Pallé, D. Godoy Rivera, F. Dai, D. R. Ciardi, J. M. Akana Murphy, P. G. Beck, C. J. Burke, J. Cabrera, I. Carleo, W. D. Cochran, K. A. Collins, Sz. Csizmadia, M. El Mufti, M. Fridlund, A. Fukui, D. Gandolfi, R. A. García, E. W. Guenther, P. Guerra , et al. (27 additional authors not shown)

Abstract: TOI 1416 (BD+42 2504, HIP 70705) is a V=10 late G or early K-type dwarf star with transits detected by TESS. Radial velocities verify the presence of the transiting planet TOI-1416 b, with a period of 1.07d, a mass of $3.48 M_{Earth}$ and a radius of $1.62 R_{Earth}$, implying a slightly sub-Earth density of $4.50$ g cm$^{-3}$. The RV data also further indicate a tentative planet c with a period o… ▽ More TOI 1416 (BD+42 2504, HIP 70705) is a V=10 late G or early K-type dwarf star with transits detected by TESS. Radial velocities verify the presence of the transiting planet TOI-1416 b, with a period of 1.07d, a mass of $3.48 M_{Earth}$ and a radius of $1.62 R_{Earth}$, implying a slightly sub-Earth density of $4.50$ g cm$^{-3}$. The RV data also further indicate a tentative planet c with a period of 27.4 or 29.5 days, whose nature cannot be verified due to strong suspicions about contamination by a signal related to the Moon's synodic period of 29.53 days. The near-USP (Ultra Short Period) planet TOI-1416 b is a typical representative of a short-period and hot ($T_{eq} \approx$ 1570 K) super-Earth like planet. A planet model of an interior of molten magma containing a significant fraction of dissolved water provides a plausible explanation for its composition, and its atmosphere could be suitable for transmission spectroscopy with JWST. The position of TOI-1416 b within the radius-period distribution corroborates that USPs with periods of less than one day do not form any special group of planets. Rather, this implies that USPs belong to a continuous distribution of super-Earth like planets with periods ranging from the shortest known ones up to ~ 30 days, whose period-radius distribution is delimitated against larger radii by the Neptune desert and by the period-radius valley that separates super-Earths from sub-Neptune planets. In the abundance of small-short periodic planets against period, a plateau between periods of 0.6 to 1.4 days has however become notable that is compatible with the low-eccentricity formation channel. For the Neptune desert, its lower limits required a revision due to the increasing population of short period planets and new limits are provided. These limits are also given in terms of the planets' insolation and effective temperatures. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Comments: 31 pages, 31 figures, 8 tables, accepted for publication in A&A

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

Secular dipole-dipole stability of magnetic binaries

Authors: Christopher Aykroyd, Adrien Bourgoin, Christophe Le Poncin-Lafitte, Stéphane Mathis, Marie-Christine Angonin

Abstract: The presence of strong large-scale stable magnetic fields in a significant portion of early-type stars, white dwarfs, and neutron stars is well established. Despite this, the origins of these fields remain unresolved, with leading propositions advocating fossil fields, mergers, and shear-driven dynamos as the main mechanism. A potential key for further insight could lie in the connection with bina… ▽ More The presence of strong large-scale stable magnetic fields in a significant portion of early-type stars, white dwarfs, and neutron stars is well established. Despite this, the origins of these fields remain unresolved, with leading propositions advocating fossil fields, mergers, and shear-driven dynamos as the main mechanism. A potential key for further insight could lie in the connection with binarity: notably, magnetism can play a role in the long-term orbital and rotational dynamics of binaries. In gravitational wave astronomy, the advanced sensitivity of upcoming detectors such as LISA and the Einstein Telescope will enable the characterisation of the orbital inspirals of compact systems, including their magnetic properties. A comprehensive understanding of the dynamics of magnetism in these systems is required for the interpretation of the gravitational wave signals and to avoid calibration biases. Furthermore, this knowledge can be used to create new magnetic population models and to provide insight into the nature of their internal fields. The aim of this study is to investigate the secular spin precession dynamics of binary systems under pure magnetic dipole interactions, focusing on stars with strong, stable, dipolar fields. We employ an orbit-averaging procedure for the spin equations and obtain an effective secular description. By minimising the magnetic energy, we derive the configurations of equilibrium. We show that among the four states, only one is stable, consisting of the spin and magnetic axes of one star reversed with respect to the companions', and orthogonal to the orbital plane. Our long-term stability results disagree with usual methods, which tend to neglect orbital motion. Finally, we provide analytical solutions for the system out of equilibrium, which can be used to derive secular orbital evolution in the context of gravitational wave astronomy. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Comments: 16 pages, 5 figures

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

gRNAde: Geometric Deep Learning for 3D RNA inverse design

Authors: Chaitanya K. Joshi, Arian R. Jamasb, Ramon Viñas, Charles Harris, Simon V. Mathis, Alex Morehead, Rishabh Anand, Pietro Liò

Abstract: Computational RNA design tasks are often posed as inverse problems, where sequences are designed based on adopting a single desired secondary structure without considering 3D geometry and conformational diversity. We introduce gRNAde, a geometric RNA design pipeline operating on 3D RNA backbones to design sequences that explicitly account for structure and dynamics. Under the hood, gRNAde is a mul… ▽ More Computational RNA design tasks are often posed as inverse problems, where sequences are designed based on adopting a single desired secondary structure without considering 3D geometry and conformational diversity. We introduce gRNAde, a geometric RNA design pipeline operating on 3D RNA backbones to design sequences that explicitly account for structure and dynamics. Under the hood, gRNAde is a multi-state Graph Neural Network that generates candidate RNA sequences conditioned on one or more 3D backbone structures where the identities of the bases are unknown. On a single-state fixed backbone re-design benchmark of 14 RNA structures from the PDB identified by Das et al. [2010], gRNAde obtains higher native sequence recovery rates (56% on average) compared to Rosetta (45% on average), taking under a second to produce designs compared to the reported hours for Rosetta. We further demonstrate the utility of gRNAde on a new benchmark of multi-state design for structurally flexible RNAs, as well as zero-shot ranking of mutational fitness landscapes in a retrospective analysis of a recent RNA polymerase ribozyme structure. Open source code: https://github.com/chaitjo/geometric-rna-design △ Less

Submitted 25 May, 2024; v1 submitted 24 May, 2023; originally announced May 2023.

Comments: Previously titled 'Multi-State RNA Design with Geometric Multi-Graph Neural Networks', presented at ICML 2023 Computational Biology Workshop

Detection of magnetic galactic binaries in quasi-circular orbit with LISA

Authors: Etienne Savalle, Adrien Bourgoin, Christophe Le Poncin-Lafitte, StÉphane Mathis, Marie-Christine Angonin, Christopher Aykroyd

Abstract: Laser Interferometer Space Antenna (LISA) will observe gravitational waves from galactic binaries (GBs) of white dwarfs or neutron stars. Some of these objects are among the most magnetic astrophysical objects in the Universe. Magnetism, by secularly disrupting the orbit, can eventually affect the gravitational waves emission and could then be potentially detected and characterized after several y… ▽ More Laser Interferometer Space Antenna (LISA) will observe gravitational waves from galactic binaries (GBs) of white dwarfs or neutron stars. Some of these objects are among the most magnetic astrophysical objects in the Universe. Magnetism, by secularly disrupting the orbit, can eventually affect the gravitational waves emission and could then be potentially detected and characterized after several years of observations by LISA. Currently, the data processing pipeline of the LISA Data Challenge (LDC) for GBs does not consider either magnetism or eccentricity. Recently, it was shown [Bourgoin et al. PRD 105, 124042 (2022)] that magnetism induces a shift on the gravitational wave frequencies. Additionally, it was argued that, if the binary's orbit is eccentric, the presence of magnetism could be detected by LISA. In this work, we explore the consequences of a future data analysis conducted on quasi-circular and magnetic GB systems using the current LDC tools. We first show that a single eccentric GB can be interpreted as several GBs and this can eventually bias population studies deduced from LISA's future catalog. Then, we confirm that for quasi-circular orbits, the secular magnetic energy of the system can be inferred if the signal-to-noise ratio of the second harmonic is high enough to be detected by traditional quasi-monochromatic source searching algorithms. LISA observations could therefore bring new insights on the nature and origin of magnetic fields in white dwarfs or neutron stars. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 18 pages, 6 figures

Spin evolution of Venus-like planets subjected to gravitational and thermal tides

Authors: Alexandre Revol, Émeline Bolmont, Gabriel Tobie, Caroline Dumoulin, Yann Musseau, Stéphane Mathis, Antoine Strugarek, Allan-Sacha Brun

Abstract: The arrival of powerful instruments will provide valuable data for the characterization of rocky exoplanets. It is then crucial to accurately model the dynamical state of exoplanets. Rocky planets with sufficiently large orbits should have non-zero eccentricities and/or obliquities. Realistic models of tides for rocky planets can allow for higher spin states than the synchronization state in the p… ▽ More The arrival of powerful instruments will provide valuable data for the characterization of rocky exoplanets. It is then crucial to accurately model the dynamical state of exoplanets. Rocky planets with sufficiently large orbits should have non-zero eccentricities and/or obliquities. Realistic models of tides for rocky planets can allow for higher spin states than the synchronization state in the presence of eccentricities or obliquities. This work explores the secular evolution of a star-planet system under tidal interactions, both gravitational and thermal, induced respectively by the quadrupolar component of the gravitational potential and the irradiation of the planet's surface. We use the formalism of Kaula associated with an Andrade rheology to model a relevant response of a rocky planet to gravitational tides and a prescription of thermal tides fitted for Venus to model the response of the atmosphere to the thermal tides. We implemented the general secular evolution equations of tidal interactions in the secular code ESPEM (French acronym for Evolution of Planetary System and Magnetism). We show the possible spin-orbit evolution and resonances for eccentric orbits and explore the possible spin orbit resonances raised by the obliquity of the planet. Our simulations have shown that the secular evolution of the spin and obliquity can lead to the retrograde spin of the Venus-like planet if the system starts from a high spin obliquity, in agreement to previous studies. Taking into account the luminosity evolution of the Sun changes the picture. We find that the planet never reaches the equilibrium: the timescale of rotation evolution is longer than the luminosity variation timescale, which suggests that Venus may never reach a spin equilibrium state but may still evolve. △ Less

Submitted 20 March, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 19 pages, 16 figues

Linking the Interiors and Surfaces of Magnetic Stars

Authors: Jim Fuller, Stephane Mathis

Abstract: Strong magnetic fields are observed in a substantial fraction of upper main sequence stars and white dwarfs. Many such stars are observed to exhibit photometric modulations as the magnetic poles rotate in and out of view, which could be a consequence of magnetic perturbations to the star's thermal structure. The magnetic pressure is typically larger than the gas pressure at the star's photosphere,… ▽ More Strong magnetic fields are observed in a substantial fraction of upper main sequence stars and white dwarfs. Many such stars are observed to exhibit photometric modulations as the magnetic poles rotate in and out of view, which could be a consequence of magnetic perturbations to the star's thermal structure. The magnetic pressure is typically larger than the gas pressure at the star's photosphere, but much smaller than the gas pressure in the star's interior, so the expected surface flux perturbations are not clear. We compute magnetically perturbed stellar structures of young $3 \, M_\odot$ stars that are in both hydrostatic and thermal equilibrium, and which contain both poloidal and toroidal components of a dipolar magnetic field as expected for stable fossil fields. This provides semi-analytical models of such fields in baroclinic stably stratified regions. The star's internal pressure, temperature, and flux perturbations can have a range of magnitudes, though we argue the most likely configurations exhibit flux perturbations much smaller than the ratio of surface magnetic pressure to surface gas pressure, but much larger than the ratio of surface magnetic pressure to central gas pressure. The magnetic pole is hotter than the equator in our models, but a cooler magnetic pole is possible depending on the magnetic field configuration. The expected flux variations for observed field strengths are $δL/L \! \lesssim \! 10^{-6}$, much smaller than those observed in magnetic stars, suggesting that observed perturbations stem from changes to the emergent spectrum rather than changes to the bolometric flux. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: Submitted to MNRAS, comments welcome!

On the Expressive Power of Geometric Graph Neural Networks

Authors: Chaitanya K. Joshi, Cristian Bodnar, Simon V. Mathis, Taco Cohen, Pietro Liò

Abstract: The expressive power of Graph Neural Networks (GNNs) has been studied extensively through the Weisfeiler-Leman (WL) graph isomorphism test. However, standard GNNs and the WL framework are inapplicable for geometric graphs embedded in Euclidean space, such as biomolecules, materials, and other physical systems. In this work, we propose a geometric version of the WL test (GWL) for discriminating geo… ▽ More The expressive power of Graph Neural Networks (GNNs) has been studied extensively through the Weisfeiler-Leman (WL) graph isomorphism test. However, standard GNNs and the WL framework are inapplicable for geometric graphs embedded in Euclidean space, such as biomolecules, materials, and other physical systems. In this work, we propose a geometric version of the WL test (GWL) for discriminating geometric graphs while respecting the underlying physical symmetries: permutations, rotation, reflection, and translation. We use GWL to characterise the expressive power of geometric GNNs that are invariant or equivariant to physical symmetries in terms of distinguishing geometric graphs. GWL unpacks how key design choices influence geometric GNN expressivity: (1) Invariant layers have limited expressivity as they cannot distinguish one-hop identical geometric graphs; (2) Equivariant layers distinguish a larger class of graphs by propagating geometric information beyond local neighbourhoods; (3) Higher order tensors and scalarisation enable maximally powerful geometric GNNs; and (4) GWL's discrimination-based perspective is equivalent to universal approximation. Synthetic experiments supplementing our results are available at \url{https://github.com/chaitjo/geometric-gnn-dojo} △ Less

Submitted 3 March, 2024; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: ICML 2023

Journal ref: Proceedings of the 40th International Conference on Machine Learning, PMLR 202:15330-15355, 2023

Gravitational waves radiated by magnetic galactic binaries and detection by LISA

Authors: Adrien Bourgoin, Christophe Le Poncin-Lafitte, Stéphane Mathis, Marie-Christine Angonin

Abstract: In the context of the future Laser Interferometer Space Antenna (LISA) mission, galactic binary systems of white dwarfs and neutron stars will represent the dominant source of Gravitational Waves (GWs) within the $10^{-4}-10^{-1}\,\mathrm{Hz}$ frequency band. It is expected that LISA will measure simultaneously, the GWs from more than ten thousands of these compact galactic binaries. The analysis… ▽ More In the context of the future Laser Interferometer Space Antenna (LISA) mission, galactic binary systems of white dwarfs and neutron stars will represent the dominant source of Gravitational Waves (GWs) within the $10^{-4}-10^{-1}\,\mathrm{Hz}$ frequency band. It is expected that LISA will measure simultaneously, the GWs from more than ten thousands of these compact galactic binaries. The analysis of such a superposition of signals will represent one of the greatest challenge for the mission. Currently, in the LISA Datacode Challenge, each galactic binary is modeled as a quasi-monochromatic source of GWs. This corresponds to the circular motion of two point-masses at the 2.5 post-Newtonian approximation. If this picture is expected to be an accurate description for most of the galactic binaries that LSIA will detect, we nevertheless expect to observe eccentric systems with complex physical properties beyond the point-mass approximation. In this work, we investigate how a binary system of highly magnetic objects in quasi-circular orbit could affect the quasi-monochromatic picture of the GW signal detected by LISA. We demonstrate that the eccentricity generates additional frequency peaks at harmonics of the mean motion and that magnetism is responsible for shifting each frequency peak with respect to the case without magnetism. We provide analytical estimates and argue that LISA will be able to detect magnetism if it can measure the main peaks at two and three times the mean motion with a sufficient accuracy. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 4 pages, 1 figure, proceedings Les rencontres de Moriond

Impact of magnetism on gravitational waves emitted by compact galactic binaries in quasi-circular orbits

Authors: Adrien Bourgoin, Etienne Savalle, Christophe Le Poncin-Lafitte, Stéphane Mathis, Marie-Christine Angonin, Antoine Strugarek

Abstract: The LISA (Laser Interferometer Space Antenna) mission will observe in the low frequency band from 0.1 mHz to 1 Hz. In this regime, we expect the galactic binaries to be the dominant (by number) sources of gravitational waves signal. Considering that galactic binaries are composed of the most magnetized astrophysical objects in the universe (i.e., the white dwarfs and the neutron stars), LISA is ex… ▽ More The LISA (Laser Interferometer Space Antenna) mission will observe in the low frequency band from 0.1 mHz to 1 Hz. In this regime, we expect the galactic binaries to be the dominant (by number) sources of gravitational waves signal. Considering that galactic binaries are composed of the most magnetized astrophysical objects in the universe (i.e., the white dwarfs and the neutron stars), LISA is expected to bring new informations about the origin and the nature of magnetism inside degenerated stars. Currently, the data processing assumes that the galactic binary systems are non-magnetic and in circular orbits which can potentially biased the determination of the parameters of the sources and also the calibration of the detector. In this work, we investigate the impact of magnetism on gravitational waves emitted by compact galactic binaries assuming quasi-circular orbits. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 4 pages, 2 figures, proceedings SF2A

How do tidal waves interact with convective vortices in rapidly-rotating planets and stars?

Authors: Virgile Dandoy, Junho Park, Kyle Augustson, Aurélie Astoul, Stéphane Mathis

Abstract: The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal wave… ▽ More The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal waves is modelled most of the time by an effective eddy-viscosity. This approach may be valid when the characteristic length scales of convective eddies are smaller than those of tidal waves. However, it becomes highly questionable in the case where tidal waves interact with potentially stable large-scale vortices, as those observed at the pole of Jupiter and Saturn. They are potentially triggered by convection in rapidly-rotating bodies in which the Coriolis acceleration forms the flow in columnar vortical structures along the direction of the rotation axis. In this paper, we investigate the complex interactions between a tidal inertial wave and a columnar convective vortex. We use a quasi-geostrophic semi-analytical model of a convective columnar vortex. We perform linear stability analysis to identify the unstable regime and conduct linear numerical simulations for the interactions between the convective vortex and an incoming tidal inertial wave. We verify that in the unstable regime, an incoming tidal inertial wave triggers the most unstable mode of the vortex leading to turbulent dissipation. For stable vortices, the wave-vortex interaction leads to the momentum mixing while it creates a low-velocity region around the vortex core and a new wave-like perturbation in the form of a progressive wave radiating in the far field. The emission of this secondary wave is the strongest when the wavelength of the incoming wave is close to the characteristic size of the vortex. △ Less

Submitted 8 March, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: 20 pages, 15 figures, accepted in Astronomy & Astrophysics

Journal ref: A&A 673, A6 (2023)

Can we detect deep axisymmetric toroidal magnetic fields in stars?

Authors: Hachem Dhouib, Stéphane Mathis, Lisa Bugnet, Timothy Van Reeth, Conny Aerts

Abstract: One of the major discoveries of asteroseismology is the signature of a strong extraction of angular momentum (AM) in the radiative zones of stars across the entire Hertzsprung-Russell diagram, resulting in weak core-to-surface rotation contrasts. Despite all efforts, a consistent AM transport theory, which reproduces both the internal rotation and mixing probed thanks to the seismology of stars, r… ▽ More One of the major discoveries of asteroseismology is the signature of a strong extraction of angular momentum (AM) in the radiative zones of stars across the entire Hertzsprung-Russell diagram, resulting in weak core-to-surface rotation contrasts. Despite all efforts, a consistent AM transport theory, which reproduces both the internal rotation and mixing probed thanks to the seismology of stars, remains one of the major open problems in modern stellar astrophysics. A possible key ingredient to figure out this puzzle is magnetic field with its various possible topologies. Among them, strong axisymmetric toroidal fields, which are subject to the so-called Tayler MHD instability, could play a major role. They could trigger a dynamo action in radiative layers while the resulting magnetic torque allows an efficient transport of AM. But is it possible to detect signatures of these deep toroidal magnetic fields? The only way to answer this question is asteroseismology and the best laboratories of study are intermediate-mass and massive stars because of their external radiative envelope. Since most of these are rapid rotators during their main-sequence, we have to study stellar pulsations propagating in stably stratified, rotating, and potentially strongly magnetised radiative zones. For that, we generalise the traditional approximation of rotation, which provides in its classic version a flexible treatment of the adiabatic propagation of gravito-inertial modes, by taking simultaneously general axisymmetric differential rotation and toroidal magnetic fields into account. Using this new non-perturbative formalism, we derive the asymptotic properties of magneto-gravito-inertial modes and we explore the different possible field configurations. We found that the magnetic effects should be detectable for equatorial fields using high-precision asteroseismic data. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 4 pages, 2 figures. Proceeding of the Annual meeting of the French Society of Astronomy and Astrophysics (SF2A 2022)

Thought experiments in a quantum computer

Authors: Nuriya Nurgalieva, Simon Mathis, Lídia del Rio, Renato Renner

Abstract: We introduce a software package that allows users to design and run simulations of thought experiments in quantum theory. In particular, it covers cases where several reasoning agents are modelled as quantum systems, such as Wigner's friend experiment. Users can customize the protocol of the experiment, the inner workings of agents (including a quantum circuit that models their reasoning process),… ▽ More We introduce a software package that allows users to design and run simulations of thought experiments in quantum theory. In particular, it covers cases where several reasoning agents are modelled as quantum systems, such as Wigner's friend experiment. Users can customize the protocol of the experiment, the inner workings of agents (including a quantum circuit that models their reasoning process), the abstract logical system used (which may or not allow agents to combine premises and make inferences about each other's reasoning), and the interpretation of quantum theory used by different agents. Our open-source software is written in a quantum programming language, ProjectQ, and runs on classical or quantum hardware. As an example, we model the Frauchiger-Renner extended Wigner's friend thought experiment, where agents are allowed to measure each other's physical memories, and make inferences about each other's reasoning. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 18 + 10 pages

Spinning up the Surface: Evidence for Planetary Engulfment or Unexpected Angular Momentum Transport?

Authors: Jamie Tayar, Facundo D. Moyano, Melinda Soares-Furtado, Ana Escorza, Meridith Joyce, Sarah L. Martell, Rafael A. García, Sylvain N. Breton, Stéphane Mathis, Savita Mathur, Vincent Delsanti, Sven Kiefer, Sabine Reffert, Dominic M. Bowman, Timothy Van Reeth, Shreeya Shetye, Charlotte Gehan, Samuel K. Grunblatt

Abstract: In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower enve… ▽ More In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower envelope and surface. KIC 9267654, however, seems to show a surface rotation rate that is faster than its bulk envelope rotation rate, in conflict with this simple angular momentum conservation argument. We improve the spectroscopic surface constraint, show that the pulsation frequencies are consistent with the previously published core and envelope rotation rates, and demonstrate that the star does not show strong chemical peculiarities. We discuss the evidence against any tidally interacting stellar companion. Finally, we discuss the possible origin of this unusual rotation profile, including the potential ingestion of a giant planet or unusual angular momentum transport by tidal inertial waves triggered by a close substellar companion, and encourage further observational and theoretical efforts. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 18 pages, 9 figures, submitted to AAS Journals

Energy self-sufficient systems for monitoring sewer networks

Authors: Simon Mathis, Juan-Mario Gruber, Christian Ebi, Simon Bloem, Jörg Rieckermann, Frank Blumensaat

Abstract: Underground infrastructure networks form the backbone of vital supply and disposal systems. However, they are under-monitored in comparison to their value. This is due, in large part, to the lack of energy supply for monitoring and data transmission. In this paper, we investigate a novel, energy harvesting system used to power underground sewer infrastructure monitoring networks. The system collec… ▽ More Underground infrastructure networks form the backbone of vital supply and disposal systems. However, they are under-monitored in comparison to their value. This is due, in large part, to the lack of energy supply for monitoring and data transmission. In this paper, we investigate a novel, energy harvesting system used to power underground sewer infrastructure monitoring networks. The system collects the required energy from ambient sources, such as temperature differences or residual light in sewer networks. A prototype was developed that could use either a thermoelectric generator (TEG) or a solar cell to capture the energy needed to acquire and transmit ultrasonic water level data via LoRaWAN. Real-world field trials were satisfactory and showed the potential power output, as well as, possibilities to improve the system. Using an extrapolation model, we proved that the developed solution could work reliably throughout the year. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: To be published in proceedings of the conference "21. ITG/GMA- Fachtagung Sensoren und Messsysteme 2022", 10.-11. Mai 2022, Nürnberger CongressCenter, Nuremberg, Germany, or IEEE explore

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Detecting deep axisymmetric toroidal magnetic fields in stars. The traditional approximation of rotation for differentially rotating deep spherical shells with a general azimuthal magnetic field

Authors: Hachem Dhouib, Stéphane Mathis, Lisa Bugnet, Timothy Van Reeth, Conny Aerts

Abstract: Asteroseismology has revealed small core-to-surface rotation contrasts in stars in the whole HR diagram. This is the signature of strong transport of angular momentum (AM) in stellar interiors. One of the plausible candidates to efficiently carry AM is magnetic fields with various topologies that could be present in stellar radiative zones. Among them, strong axisymmetric azimuthal magnetic fields… ▽ More Asteroseismology has revealed small core-to-surface rotation contrasts in stars in the whole HR diagram. This is the signature of strong transport of angular momentum (AM) in stellar interiors. One of the plausible candidates to efficiently carry AM is magnetic fields with various topologies that could be present in stellar radiative zones. Among them, strong axisymmetric azimuthal magnetic fields have received a lot of interest. Indeed, if they are subject to the so-called Tayler instability, the accompanying triggered Maxwell stresses can transport AM efficiently. In addition, the electromotive force induced by the fluctuations of magnetic and velocity fields could potentially sustain a dynamo action that leads to the regeneration of the initial strong axisymmetric azimuthal magnetic field. The key question we aim to answer is: can we detect signatures of these deep strong azimuthal magnetic fields? The only way to answer this question is asteroseismology and the best laboratories of study are intermediate-mass and massive stars. Most of these are rapid rotators during their main-sequence. Therefore, we have to study stellar pulsations propagating in stably stratified, rotating, and potentially strongly magnetised radiative zones. We generalise the traditional approximation of rotation by simultaneously taking general axisymmetric differential rotation and azimuthal magnetic fields into account in a non-perturbative way. Using this new formalism, we derive the asymptotic properties of magneto-gravito-inertial (MGI) waves and their period spacings. We find that toroidal magnetic fields induce a shift in the period spacings of MGI modes. An equatorial azimuthal magnetic field with an amplitude of the order of $10^5\,\rm G$ leads to signatures that can be detectable thanks to modern space photometry. More complex hemispheric configurations are more difficult to observe. △ Less

Submitted 22 May, 2022; v1 submitted 21 February, 2022; originally announced February 2022.

Comments: 21 pages, 15 figures, 1 table, abstract shortened for arXiv. Published in A&A

Journal ref: A&A 661, A133 (2022)

99 oscillating red-giant stars in binary systems with NASA TESS and NASA Kepler identified from the SB9-Catalogue

Authors: P. G. Beck, S. Mathur, K. Hambleton, R. A. García, L. Steinwender, N. L. Eisner, J. -D. do Nascimento, P. Gaulme, S. Mathis

Abstract: Oscillating red-giant stars in binary systems are an ideal testbed for investigating the structure and evolution of stars in the advanced phases of evolution. With 83 known red giants in binary systems, of which only ~40 have determined global seismic parameters and orbital parameters, the sample is small compared to the numerous known oscillating stars. The detection of red-giant binary systems i… ▽ More Oscillating red-giant stars in binary systems are an ideal testbed for investigating the structure and evolution of stars in the advanced phases of evolution. With 83 known red giants in binary systems, of which only ~40 have determined global seismic parameters and orbital parameters, the sample is small compared to the numerous known oscillating stars. The detection of red-giant binary systems is typically obtained from the signature of stellar binarity in space photometry. The time base of such data biases the detection towards systems with shorter periods and orbits of insufficient size to allow a red giant to fully extend as it evolves up the red-giant branch. Consequently, the sample shows an excess of H-shell burning giants while containing very few stars in the He-core burning phase. From the ninth catalogue of spectroscopic binary orbits (SB9), we identified candidate systems hosting a red-giant primary component. Searching space photometry from the NASA missions Kepler, K2, and TESS (Transiting Exoplanet Survey Satellite) as well as the BRITE (BRIght Target Explorer) constellation mission, we find 99 systems, which were previously unknown to host an oscillating giant component. The revised search strategy allowed us to extend the range of orbital periods of systems hosting oscillating giants up to 26000days. Such wide orbits allow a rich population of He-core burning primaries, which are required for a complete view of stellar evolution from binary studies. Tripling the size of the sample of known oscillating red-giant stars in binary systems is an important step towards an ensemble approach for seismology and tidal studies. While for non-eclipsing binaries the inclination is unknown, such a seismically well-characterized sample will be a treasure trove in combination with Gaia astrometric orbits for binary~systems. △ Less

Submitted 27 April, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: A&A (accepted for publication)

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

Impact of dipolar magnetic fields on gravitational wave strain by galactic binaries

Authors: Adrien Bourgoin, Christophe Le Poncin-Lafitte, Stéphane Mathis, Marie-Christine Angonin

Abstract: White dwarfs (WDs) and neutron stars (NSs) are among the most magnetized astrophysical objects in the universe, with magnetic fields that can reach up to $10^9\,\mathrm{G}$ for WDs and up to $10^{15}\,\mathrm{G}$ for NSs. The galaxy is expected to be populated with approximately one hundred million of double WD and millions of NS-WD binaries. Throughout the duration of the mission, the Laser Inter… ▽ More White dwarfs (WDs) and neutron stars (NSs) are among the most magnetized astrophysical objects in the universe, with magnetic fields that can reach up to $10^9\,\mathrm{G}$ for WDs and up to $10^{15}\,\mathrm{G}$ for NSs. The galaxy is expected to be populated with approximately one hundred million of double WD and millions of NS-WD binaries. Throughout the duration of the mission, the Laser Interferometer Space Antenna (LISA) will observe gravitational waves (GWs) emitted simultaneously by more than ten thousand of such galactic binaries. In this paper, we investigate the effect of the magnetic dipole-dipole interaction on the GW signal emitted by magnetic galactic binaries. We derive the secular equations governing the orbital and rotational motion of these objects. Then, we integrate these equations both numerically and analytically. We conclude that the overall visible effect is an additional secular drift of the mean longitude. This drift is proportional to the product of the magnetic moments and is inversely proportional to the $7/2$ power of the semi-major axis. Finally, we show that, at zeroth-order in eccentricity, the magnetic dipole-dipole interaction shifts the main frequency of the gravitational strain measured by LISA. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 21 pages, 8 figures

Rossby numbers and stiffness values inferred from gravity-mode asteroseismology of rotating F- and B-type dwarfs: consequences for mixing, transport, magnetism, and convective penetration

Authors: C. Aerts, K. Augustson, S. Mathis, M. G. Pedersen, J. S. G. Mombarg, V. Vanlaer, J. Van Beeck, T. Van Reeth

Abstract: Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of inte… ▽ More Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of intermediate-mass stars. We deduce these quantities for rotating dwarfs from the observed properties of their identified gravity and gravito-inertial modes. We rely on near-core rotation rates and asteroseismic models of 26 B- and 37 F-type dwarf pulsators derived from 4-year Kepler space photometry, high-resolution spectroscopy and Gaia astrometry in the literature to deduce their convective and wave Rossby numbers. We compute the stiffness at the convection/radiation interface from the inferred maximum buoyancy frequency at the interface and the convective turnover frequency in the core. We use those asteroseismically inferred quantities to make predictions of convective penetration levels, local flux levels of gravito-inertial waves triggered by the convective core, and of the cores' potential rotational and magnetic states. Our sample of 63 gravito-inertial mode pulsators covers near-core rotation rates from almost zero up to the critical rate. The frequencies of their identified modes lead to models with stiffness values between $10^{2.69}$ and $10^{3.60}$ for the B-type pulsators, while those of F-type stars cover the range from $10^{3.47}$ to $10^{4.52}$. The convective Rossby numbers derived from the maximum convective diffusion coefficient in the convective core, based on mixing length theory and a value of the mixing length coefficient relevant for these pulsators, vary between $10^{-2.3}$ and $10^{-0.8}$ for B-type stars and $10^{-3}$ and $10^{-1.5}$ for F-type stars. (abridged) △ Less

Submitted 12 December, 2021; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 10 pages, 8 figures, accepted for publication in Astronomy & Astrophysics, version adapted to proofs

Journal ref: A&A 656, A121 (2021)

The traditional approximation of rotation for rapidly rotating stars and planets. II. Deformation and differential rotation

Authors: Hachem Dhouib, Vincent Prat, Timothy Van Reeth, Stéphane Mathis

Abstract: We examine the dynamics of low-frequency gravito-inertial waves (GIWs) in differentially rotating deformed radiation zones in stars and planets by generalising the traditional approximation of rotation (TAR). The TAR treatment was built on the assumptions that the star is spherical and uniformly rotating. However, it has been generalised in our previous work by including the effects of the centrif… ▽ More We examine the dynamics of low-frequency gravito-inertial waves (GIWs) in differentially rotating deformed radiation zones in stars and planets by generalising the traditional approximation of rotation (TAR). The TAR treatment was built on the assumptions that the star is spherical and uniformly rotating. However, it has been generalised in our previous work by including the effects of the centrifugal deformation using a non-perturbative approach in the uniformly rotating case. We aim to carry out a new generalisation of the TAR treatment to account for the differential rotation and the strong centrifugal deformation simultaneously. We generalise our previous work by taking into account the differential rotation in the derivation of our complete analytical formalism that allows the study of the dynamics of GIWs in differentially and rapidly rotating stars. We derived the complete set of equations that generalises the TAR, simultaneously taking the full centrifugal acceleration and the differential rotation into account. Within the validity domain of the TAR, we derived a generalised Laplace tidal equation for the horizontal eigenfunctions and asymptotic wave periods of the GIWs, which can be used to probe the structure and dynamics of differentially rotating deformed stars with asteroseismology. A new generalisation of the TAR, which simultaneously takes into account the differential rotation and the centrifugal acceleration in a non-perturbative way, was derived. This generalisation allowed us to study the detectability and the signature of the differential rotation on GIWs in rapidly rotating deformed stars and planets. We found that the effects of the differential rotation in early-type deformed stars on GIWs is theoretically largely detectable in modern space photometry using observations from $\textit{Kepler}$ and TESS. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 13 pages, 11 figures, 1 table, abstract shortened for arXiv. Accepted for publication in A&A. arXiv admin note: substantial text overlap with arXiv:2104.09302

Journal ref: A&A 656, A122 (2021)

Dipolar magnetic fields in binaries and gravitational waves

Authors: Adrien Bourgoin, Christophe Le Poncin-Lafitte, Stéphane Mathis, Marie-Christine Angonin

Abstract: The LISA mission will observe gravitational waves emitted from tens of thousands of galactic binaries, in particular white dwarf binary systems. These objects are known to have intense magnetic fields. However, these fields are usually not considered as their influence on the orbital and rotational motion of the binary is assumed for being too weak. It turns out that magnetic fields modify the orb… ▽ More The LISA mission will observe gravitational waves emitted from tens of thousands of galactic binaries, in particular white dwarf binary systems. These objects are known to have intense magnetic fields. However, these fields are usually not considered as their influence on the orbital and rotational motion of the binary is assumed for being too weak. It turns out that magnetic fields modify the orbits, in particular their geometry with respect to the observer. In this work, we revisit the issue, assuming magnetostatic approximation, and we show how the magnetic fields within a binary system generate a secular drift in the argument of the periastron, leading then, to modifications of the gravitational waveforms that are potentially detectable by LISA. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 4 pages, 2 figures, conference paper

A 20-Second Cadence View of Solar-Type Stars and Their Planets with TESS: Asteroseismology of Solar Analogs and a Re-characterization of pi Men c

Authors: Daniel Huber, Timothy R. White, Travis S. Metcalfe, Ashley Chontos, Michael M. Fausnaugh, Cynthia S. K. Ho, Vincent Van Eylen, Warrick Ball, Sarbani Basu, Timothy R. Bedding, Othman Benomar, Diego Bossini, Sylvain Breton, Derek L. Buzasi, Tiago L. Campante, William J. Chaplin, Joergen Christensen-Dalsgaard, Margarida S. Cunha, Morgan Deal, Rafael A. Garcia, Antonio Garcia Munoz, Charlotte Gehan, Lucia Gonzalez-Cuesta, Chen Jiang, Cenk Kayhan , et al. (28 additional authors not shown)

Abstract: We present an analysis of the first 20-second cadence light curves obtained by the TESS space telescope during its extended mission. We find a precision improvement of 20-second data compared to 2-minute data for bright stars when binned to the same cadence (~10-25% better for T<~8 mag, reaching equal precision at T~13 mag), consistent with pre-flight expectations based on differences in cosmic ra… ▽ More We present an analysis of the first 20-second cadence light curves obtained by the TESS space telescope during its extended mission. We find a precision improvement of 20-second data compared to 2-minute data for bright stars when binned to the same cadence (~10-25% better for T<~8 mag, reaching equal precision at T~13 mag), consistent with pre-flight expectations based on differences in cosmic ray mitigation algorithms. We present two results enabled by this improvement. First, we use 20-second data to detect oscillations in three solar analogs (gamma Pav, zeta Tuc and pi Men) and use asteroseismology to measure their radii, masses, densities and ages to ~1%, ~3%, ~1% and ~20% respectively, including systematic errors. Combining our asteroseismic ages with chromospheric activity measurements we find evidence that the spread in the activity-age relation is linked to stellar mass and thus convection-zone depth. Second, we combine 20-second data and published radial velocities to re-characterize pi Men c, which is now the closest transiting exoplanet for which detailed asteroseismology of the host star is possible. We show that pi Men c is located at the upper edge of the planet radius valley for its orbital period, confirming that it has likely retained a volatile atmosphere and that the "asteroseismic radius valley" remains devoid of planets. Our analysis favors a low eccentricity for pi Men c (<0.1 at 68% confidence), suggesting efficient tidal dissipation (Q/k <~ 2400) if it formed via high-eccentricity migration. Combined, these early results demonstrate the strong potential of TESS 20-second cadence data for stellar astrophysics and exoplanet science. △ Less

Submitted 13 October, 2021; v1 submitted 20 August, 2021; originally announced August 2021.

Comments: 17 pages (excluding references), 13 figures, 6 tables; accepted for publication in AJ. Data and scripts to reproduce results are archived at https://zenodo.org/record/5555456

Gauge invariant quantum circuits for $U(1)$ and Yang-Mills lattice gauge theories

Authors: Giulia Mazzola, Simon V. Mathis, Guglielmo Mazzola, Ivano Tavernelli

Abstract: Quantum computation represents an emerging framework to solve lattice gauge theories (LGT) with arbitrary gauge groups, a general and long-standing problem in computational physics. While quantum computers may encode LGT using only polynomially increasing resources, a major openissue concerns the violation of gauge-invariance during the dynamics and the search for groundstates. Here, we propose a… ▽ More Quantum computation represents an emerging framework to solve lattice gauge theories (LGT) with arbitrary gauge groups, a general and long-standing problem in computational physics. While quantum computers may encode LGT using only polynomially increasing resources, a major openissue concerns the violation of gauge-invariance during the dynamics and the search for groundstates. Here, we propose a new class of parametrized quantum circuits that can represent states belonging only to the physical sector of the total Hilbert space. This class of circuits is compact yet flexible enough to be used as a variational ansatz to study ground state properties, as well as representing states originating from a real-time dynamics. Concerning the first application, the structure of the wavefunction ansatz guarantees the preservation of physical constraints such as the Gauss law along the entire optimization process, enabling reliable variational calculations. As for the second application, this class of quantum circuits can be used in combination with timedependent variational quantum algorithms, thus drastically reducing the resource requirements to access dynamical properties. △ Less

Submitted 20 January, 2022; v1 submitted 12 May, 2021; originally announced May 2021.

Comments: 16 pages, 13 figures, published version

Journal ref: Phys. Rev. Research, 2021