-
The Most Common Habitable Planets III -- Modeling Temperature Forcing and Surface Conditions on Rocky Exoplanets and Exomoons
Authors:
Beatriz B. Siffert,
Raquel G. Gonçalves Farias,
Matias Garcia,
Luiz Felipe Melo de Menezes,
Gustavo F. Porto de Mello,
Marcelo Borges Fernandes,
Rafael Pinotti
Abstract:
Small rocky planets, as well as larger planets that suffered extensive volatile loss, tend to be drier and have thinner atmospheres as compared to Earth. Such planets probably outnumber worlds better endowed with volatiles, being the most common habitable planets. For the subgroup of fast rotators following eccentric orbits, atmospheres suffer radiative forcing and their heat capacity provides a m…
▽ More
Small rocky planets, as well as larger planets that suffered extensive volatile loss, tend to be drier and have thinner atmospheres as compared to Earth. Such planets probably outnumber worlds better endowed with volatiles, being the most common habitable planets. For the subgroup of fast rotators following eccentric orbits, atmospheres suffer radiative forcing and their heat capacity provides a method for gauging atmospheric thickness and surface conditions. We further explore the model presented in a previous paper and apply it to real and hypothetical exoplanets in the habitable zone of various classes of stars, simulating atmospheric and orbital characteristics. For planetary eccentricities e ~0.3, the forcing-induced hypothetical temperature variation would reach ~80 K for airless planets and ~10 K for planets with substantial atmospheres. For Kepler-186 f and Kepler-442 b, assuming e ~0.1, temperature variations can reach ~24 K. We also consider habitable exomoons in circular orbits around gas giants within the habitable zone, which suffer radiative forcing due to their epicyclic motion. We study several combinations of parameters for the characterization of planets (mass, eccentricity and semi-major axis) and exomoons (mass, orbital radius, albedo and atmospheric characteristics) for different stellar types. For e ~0.3, exomoon temperature varies up to ~90 K, while for ~0.6 variations can reach ~200 K. Such exomoons may plausibly retain their volatiles by continued volcanic activity fueled by tidal dissipation. Although currently undetectable, such effects might be within reach of future Extremely Large Telescope-class telescopes and space missions with mid-infrared and coronagraphic capabilities.
△ Less
Submitted 9 May, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: Solar and stellar observations
Authors:
Sven Wedemeyer,
Miroslav Barta,
Roman Brajsa,
Yi Chai,
Joaquim Costa,
Dale Gary,
Guillermo Gimenez de Castro,
Stanislav Gunar,
Gregory Fleishman,
Antonio Hales,
Hugh Hudson,
Mats Kirkaune,
Atul Mohan,
Galina Motorina,
Alberto Pellizzoni,
Maryam Saberi,
Caius L. Selhorst,
Paulo J. A. Simoes,
Masumi Shimojo,
Ivica Skokic,
Davor Sudar,
Fabian Menezes,
Stephen White,
Mark Booth,
Pamela Klaassen
, et al. (13 additional authors not shown)
Abstract:
Observations at (sub-)millimeter wavelengths offer a complementary perspective on our Sun and other stars, offering significant insights into both the thermal and magnetic composition of their chromospheres. Despite the fundamental progress in (sub-)millimeter observations of the Sun, some important aspects require diagnostic capabilities that are not offered by existing observatories. In particul…
▽ More
Observations at (sub-)millimeter wavelengths offer a complementary perspective on our Sun and other stars, offering significant insights into both the thermal and magnetic composition of their chromospheres. Despite the fundamental progress in (sub-)millimeter observations of the Sun, some important aspects require diagnostic capabilities that are not offered by existing observatories. In particular, simultaneously observations of the radiation continuum across an extended frequency range would facilitate the mapping of different layers and thus ultimately the 3D structure of the solar atmosphere. Mapping large regions on the Sun or even the whole solar disk at a very high temporal cadence would be crucial for systematically detecting and following the temporal evolution of flares, while synoptic observations, i.e., daily maps, over periods of years would provide an unprecedented view of the solar activity cycle in this wavelength regime. As our Sun is a fundamental reference for studying the atmospheres of active main sequence stars, observing the Sun and other stars with the same instrument would unlock the enormous diagnostic potential for understanding stellar activity and its impact on exoplanets. The Atacama Large Aperture Submillimeter Telescope (AtLAST), a single-dish telescope with 50\,m aperture proposed to be built in the Atacama desert in Chile, would be able to provide these observational capabilities. Equipped with a large number of detector elements for probing the radiation continuum across a wide frequency range, AtLAST would address a wide range of scientific topics including the thermal structure and heating of the solar chromosphere, flares and prominences, and the solar activity cycle. In this white paper, the key science cases and their technical requirements for AtLAST are discussed.
△ Less
Submitted 6 March, 2024; v1 submitted 1 March, 2024;
originally announced March 2024.
-
Trajectories of Coronal Mass Ejection from Solar-type Stars
Authors:
Fabian Menezes,
Adriana Valio,
Yuri Netto,
Alexandre Araújo,
Christina Kay,
Merav Opher
Abstract:
The Sun and other solar-type stars have magnetic fields that permeate their interior and surface, extends through the interplanetary medium, and is the main driver of stellar activity. Stellar magnetic activity affects physical processes and conditions of the interplanetary medium and orbiting planets. Coronal mass ejections (CMEs) are the most impacting of these phenomena in near-Earth space weat…
▽ More
The Sun and other solar-type stars have magnetic fields that permeate their interior and surface, extends through the interplanetary medium, and is the main driver of stellar activity. Stellar magnetic activity affects physical processes and conditions of the interplanetary medium and orbiting planets. Coronal mass ejections (CMEs) are the most impacting of these phenomena in near-Earth space weather, and consist of plasma clouds, with magnetic field, ejected from the solar corona. Precisely predicting the trajectory of CMEs is crucial in determining whether a CME will hit a planet and impact its magnetosphere and atmosphere. Despite the rapid developments in the search for stellar CMEs, their detection is still very incipient. In this work we aim to better understand the propagation of CMEs by analysing the influence of initial parameters on CME trajectories, such as position, velocities, and stellar magnetic field's configuration. We reconstruct magnetograms for Kepler-63 (KIC 11554435) and Kepler-411 (KIC 11551692) from spot transit mapping, and use a CME deflection model, ForeCAT, to simulate trajectories of hypothetical CMEs launched into the interplanetary medium from Kepler-63 and Kepler-411. We apply the same methodology to the Sun, for comparison. Our results show that in general, deflections and rotations of CMEs decrease with their radial velocity, and increase with ejection latitude. Moreover, magnetic fields stronger than the Sun's, such as Kepler-63's, tend to cause greater CME deflections.
△ Less
Submitted 11 May, 2023;
originally announced May 2023.
-
Subterahertz Radius and Limb Brightening of the Sun Derived from SST and ALMA
Authors:
Fabian Menezes,
Caius L. Selhorst,
Carlos Guillermo Giménez de Castro,
Adriana Valio
Abstract:
Measurements of the radius and limb brightening of the Sun provide important information about the solar atmosphere structure and temperature. The solar radius increases as the observation at radio frequency decreases, indicating that each emission originates higher in the atmosphere. Thus, different layers of the solar atmosphere can be probed by observing at multiple wavelengths. In this work, w…
▽ More
Measurements of the radius and limb brightening of the Sun provide important information about the solar atmosphere structure and temperature. The solar radius increases as the observation at radio frequency decreases, indicating that each emission originates higher in the atmosphere. Thus, different layers of the solar atmosphere can be probed by observing at multiple wavelengths. In this work, we determined the average radius and limb brightening at 100, 212, 230, and 405 GHz, using data from the Solar Submillimeter Telescope and ALMA's single-dish observations. For the first time, limb brightening values for frequencies of 212 and 405 GHz were estimated. At sub-THz frequencies, the observed limb brightening may affect the solar radius measurements. We use two different and well known approaches to determine the radius: the half-power method and the inflection-point method. We investigate how the antenna beam size and the limb brightening level, LB, can affect the radius measurements using both methods. Our results showed that the inflection-point method is the least affected by these parameters, and should thus be used for solar radius estimates at radio wavelengths. The measured average radii are 968"~$\pm$~3" (100 GHz), 963"~$\pm$~3" (212 GHz), 963"~$\pm$~2" (230 GHz), and 963"~$\pm$~5" (405 GHz). Finally, we used forward modeling to estimate the ranges of LB of the solar disk resulting in 5%-19% (100 GHz), 2%-12% (212 GHz), 6%-18% (230 GHz), and 3%-17% (405 GHz). Both radius and limb brightening estimates agree with previous measurements reported in the literature.
△ Less
Submitted 30 November, 2021;
originally announced November 2021.
-
A Model for Bimodal Rates and Proportions
Authors:
Roberto Vila,
Lucas Alfaia,
André F. B. Menezes,
Mehmet N. Çankaya,
Marcelo Bourguignon
Abstract:
The beta model is the most important distribution for fitting data with the unit interval. However, the beta distribution is not suitable to model bimodal unit interval data. In this paper, we propose a bimodal beta distribution constructed by using an approach based on the alpha-skew-normal model. We discuss several properties of this distribution such as bimodality, real moments, entropy measure…
▽ More
The beta model is the most important distribution for fitting data with the unit interval. However, the beta distribution is not suitable to model bimodal unit interval data. In this paper, we propose a bimodal beta distribution constructed by using an approach based on the alpha-skew-normal model. We discuss several properties of this distribution such as bimodality, real moments, entropy measures and identifiability. Furthermore, we propose a new regression model based on the proposed model and discuss residuals. Estimation is performed by maximum likelihood. A Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples with a discussion of the results. An application is provided to show the modelling competence of the proposed distribution when the data sets show bimodality.
△ Less
Submitted 17 August, 2021;
originally announced August 2021.
-
Non-intrusive polynomial chaos expansion for topology optimization using polygonal meshes
Authors:
Nilton Cuellar,
Anderson Pereira,
Ivan F. M. Menezes,
Americo Cunha Jr
Abstract:
This paper deals with the applications of stochastic spectral methods for structural topology optimization in the presence of uncertainties. A non-intrusive polynomial chaos expansion is integrated into a topology optimization algorithm to calculate low-order statistical moments of the mechanical-mathematical model response. This procedure, known as robust topology optimization, can optimize the m…
▽ More
This paper deals with the applications of stochastic spectral methods for structural topology optimization in the presence of uncertainties. A non-intrusive polynomial chaos expansion is integrated into a topology optimization algorithm to calculate low-order statistical moments of the mechanical-mathematical model response. This procedure, known as robust topology optimization, can optimize the mean of the compliance while simultaneously minimizing its standard deviation. In order to address possible variabilities in the loads applied to the mechanical system of interest, magnitude and direction of the external forces are assumed to be uncertain. In this probabilistic framework, forces are described as a random field or a set of random variables. Representation of the random objects and propagation of load uncertainties through the model are efficiently done through Karhunen-Loève and polynomial chaos expansions. We take advantage of using polygonal elements, which have been shown to be effective in suppressing checkerboard patterns and reducing mesh dependency in the solution of topology optimization problems. Accuracy and applicability of the proposed methodology are demonstrated by means of several topology optimization examples. The obtained results, which are in excellent agreement with reference solutions computed via Monte Carlo method, show that load uncertainties play an important role in optimal design of structural systems, so that they must be taken into account to ensure a reliable optimization process.
△ Less
Submitted 14 July, 2021;
originally announced August 2021.
-
The subterahertz solar cycle: Polar and equatorial radii derived from SST and ALMA
Authors:
Fabian Menezes,
Caius L. Selhorst,
Carlos Guillermo Giménez de Castro,
Adriana Valio
Abstract:
At subterahertz frequencies -- i.e., millimeter and submillimeter wavelengths -- there is a gap of measurements of the solar radius as well as other parameters of the solar atmosphere. As the observational wavelength changes, the radius varies because the altitude of the dominant electromagnetic radiation is produced at different heights in the solar atmosphere. Moreover, radius variations through…
▽ More
At subterahertz frequencies -- i.e., millimeter and submillimeter wavelengths -- there is a gap of measurements of the solar radius as well as other parameters of the solar atmosphere. As the observational wavelength changes, the radius varies because the altitude of the dominant electromagnetic radiation is produced at different heights in the solar atmosphere. Moreover, radius variations throughout long time series are indicative of changes in the solar atmosphere that may be related to the solar cycle. Therefore, the solar radius is an important parameter for the calibration of solar atmospheric models enabling a better understanding of the atmospheric structure. In this work we use data from the Solar Submillimeter-wave Telescope (SST) and from the Atacama Large Millimeter/submillimeter Array (ALMA), at the frequencies of 100, 212, 230, and 405 GHz, to measure the equatorial and polar radii of the Sun. The radii measured with extensive data from the SST agree with the radius-vs-frequency trend present in the literature. The radii derived from ALMA maps at 230 GHz also agree with the radius-vs-frequency trend, whereas the 100-GHz radii are slightly above the values reported by other authors. In addition, we analyze the equatorial and polar radius behavior over the years, by determining the correlation coefficient between solar activity and subterahertz radii time series at 212 and 405 GHz (SST). The variation of the SST-derived radii over 13 years are correlated to the solar activity when considering equatorial regions of the solar atmosphere, and anticorrelated when considering polar regions. The ALMA derived radii time series for 100 and 230 GHz show very similar behaviors with those of SST.
△ Less
Submitted 12 February, 2021; v1 submitted 8 February, 2021;
originally announced February 2021.
-
Accelerating engineering design by automatic selection of simulation cases through Pool-Based Active Learning
Authors:
J. H. Gaspar Elsas,
N. A. G. Casaprima,
I. F. M. Menezes
Abstract:
A common workflow for many engineering design problems requires the evaluation of the design system to be investigated under a range of conditions. These conditions usually involve a combination of several parameters. To perform a complete evaluation of a single candidate configuration, it may be necessary to perform hundreds to thousands of simulations. This can be computationally very expensive,…
▽ More
A common workflow for many engineering design problems requires the evaluation of the design system to be investigated under a range of conditions. These conditions usually involve a combination of several parameters. To perform a complete evaluation of a single candidate configuration, it may be necessary to perform hundreds to thousands of simulations. This can be computationally very expensive, particularly if several configurations need to be evaluated, as in the case of the mathematical optimization of a design problem. Although the simulations are extremely complex, generally, there is a high degree of redundancy in them, as many of the cases vary only slightly from one another. This redundancy can be exploited by omitting some simulations that are uninformative, thereby reducing the number of simulations required to obtain a reasonable approximation of the complete system. The decision of which simulations are useful is made through the use of machine learning techniques, which allow us to estimate the results of "yet-to-be-performed" simulations from the ones that are already performed. In this study, we present the results of one such technique, namely active learning, to provide an approximate result of an entire offshore riser design simulation portfolio from a subset that is 80% smaller than the original one. These results are expected to facilitate a significant speed-up in the offshore riser design.
△ Less
Submitted 16 September, 2020; v1 submitted 2 September, 2020;
originally announced September 2020.
-
Solar polar brightening and radius at 100 and 230 GHz observed by ALMA
Authors:
Caius L. Selhorst,
Paulo J. A. Simões,
Roman Brajša,
Adriana Valio,
C. G. Giménez de Castro,
Joaquim E. R. Costa,
Fabian Menezes,
Jean Pierre Rozelot,
Antonio S. Hales,
Kazumasa Iwai,
Stephen White
Abstract:
Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models…
▽ More
Polar brightening of the Sun at radio frequencies has been studied for almost fifty years and yet a disagreement persists between solar atmospheric models and observations. Some observations reported brightening values much smaller than the expected values obtained from the models, with discrepancies being particularly large at millimeter wavelengths. New clues to calibrate the atmospheric models can be obtained with the advent of the Atacama Large Millimeter/submillimeter Array (ALMA) radio interferometer. In this work, we analyzed the lower limit of the polar brightening observed at 100 and 230 GHz by ALMA, during its Science Verification period, 2015 December 16-20. We find that the average polar intensity is higher than the disk intensity at 100 and 230 GHz, with larger brightness intensities at the South pole in eight of the nine maps analyzed. The observational results were compared with calculations of the millimetric limb brightnening emission for two semi-empirical atmospheric models, FAL- C (Fontenla et al. 1993) and SSC (Selhorst et al. 2005a). Both models presented larger limb intensities than the averaged observed values. The intensities obtained with the SSC model were closer to the observations, with polar brightenings of 10.5% and 17.8% at 100 and 230 GHz, respectively. This discrepancy may be due to the presence of chromospheric features (like spicules) at regions close to the limb.
△ Less
Submitted 29 November, 2018;
originally announced November 2018.
-
On the one parameter unit-Lindley distribution and its associated regression model for proportion data
Authors:
J. Mazucheli,
A. F. B. Menezes,
S. Chakraborty
Abstract:
In this paper considering the transformation $X=\frac{Y}{1+Y}$, where $Y \sim\text{Lindley}(θ)$, we propose the unit-Lindley distribution and investigate some of its mathematical properties. A important fact associated with this new distribution is that is possible to obtain the analytical expression for bias correction of the maximum likelihood estimator. Moreover, it belongs to the exponential f…
▽ More
In this paper considering the transformation $X=\frac{Y}{1+Y}$, where $Y \sim\text{Lindley}(θ)$, we propose the unit-Lindley distribution and investigate some of its mathematical properties. A important fact associated with this new distribution is that is possible to obtain the analytical expression for bias correction of the maximum likelihood estimator. Moreover, it belongs to the exponential family. This distribution allows us to incorporate covariates directly in the mean and consequently to quantify the influence on the average of the response variable. Finally, a practical application is present and it is shown that our model fits much better than the Beta regression.
△ Less
Submitted 8 January, 2018;
originally announced January 2018.
-
Solar Radius at Subterahertz Frequencies and its Relation to Solar Activity
Authors:
Fabian Menezes,
Adriana Valio
Abstract:
The Sun emits radiation at several wavelengths of the electromagnetic spectrum. In the optical band, the solar radius is 695,700 km and this is what defines the photosphere, the visible surface of the Sun. However, as the altitude increases, the electromagnetic radiation is produced at other frequencies, causing the solar radius to change as function of wavelength. These measurements enable a bett…
▽ More
The Sun emits radiation at several wavelengths of the electromagnetic spectrum. In the optical band, the solar radius is 695,700 km and this is what defines the photosphere, the visible surface of the Sun. However, as the altitude increases, the electromagnetic radiation is produced at other frequencies, causing the solar radius to change as function of wavelength. These measurements enable a better understanding of the solar atmosphere and the radius dependence on the solar cycle is a good indicator of the changes that occur in the atmospheric structure. We measure the solar radius at the subterahertz frequencies of 0.212 and 0.405 THz -- i.e., the altitude where these emissions are primarily generated -- and also analyse the radius variation over the 11-year solar activity cycle. For this, we used radio maps of the solar disk for the period between 1999 and 2017, reconstructed from daily scans made by the Solar Submillimeter-wave Telescope (SST), installed at El Leoncito Astronomical Complex (CASLEO) in the Argentinean Andes. Our measurements yield radii of 966.5" $\pm$ 2.8" for 0.2 THz and 966.5" $\pm$ 2.7" for 0.4 THz. This implies a height of $(5.0 \pm 2.0 \times 10^6)$ m above the photosphere. Furthermore, we also observed strong anti-correlation between the radius variation and the solar activity at both frequencies.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
Topology Optimization Using Polytopes
Authors:
Arun L. Gain,
Glaucio H. Paulino,
Leonardo Duarte,
Ivan F. M. Menezes
Abstract:
Meshing complex engineering domains is a challenging task. Arbitrary polyhedral meshes can provide the much needed flexibility in automated discretization of such domains. The geometric property of the polyhedral meshes such as the unstructured nature and the facial connectivity between elements makes them specially attractive for topology optimization applications. Numerical anomalies in designs…
▽ More
Meshing complex engineering domains is a challenging task. Arbitrary polyhedral meshes can provide the much needed flexibility in automated discretization of such domains. The geometric property of the polyhedral meshes such as the unstructured nature and the facial connectivity between elements makes them specially attractive for topology optimization applications. Numerical anomalies in designs such as the single node connections and checkerboard pattern, which are difficult to manufacture physically, are naturally alleviated with polyhedrons. Special interpolants such as Wachspress, mean value coordinates, maximum entropy shape functions are available to handle arbitrary shaped elements. But the finite elements approaches based on these shape functions face some challenges such as accurate and efficient computation of the shape functions and their derivatives for the numerical evaluation of the weak form integrals. In the current work, we solve the governing three-dimensional elasticity state equation using a Virtual Element Method (VEM) approach. The main characteristic difference between VEM and standard finite element methods (FEM) is that in VEM the canonical basis functions are not constructed explicitly. Rather the stiffness matrix is computed directly utilizing a projection map which extracts the linear component of the deformation. Such a construction guarantees the satisfaction of the patch test (used by engineers as an indicator of optimal convergence of numerical solutions under mesh refinement). Finally, the computations reduce to the evaluation of matrices which contain purely geometric surface facet quantities. The present work focuses on the first-order VEM in which the degrees of freedom associated with the vertices. Utilizing polyhedral elements for topology optimization, we show that the mesh bias in the member orientation is alleviated.
△ Less
Submitted 25 December, 2013;
originally announced December 2013.
-
Nonuniversality in the pair contact process with diffusion
Authors:
Ronald Dickman,
Marcio Argollo Ferreira de Menezes
Abstract:
We study the static and dynamic behavior of the one dimensional pair contact process with diffusion. Several critical exponents are found to vary with the diffusion rate, while the order-parameter moment ratio m=\bar{rho^2} /\bar{rho}^2 grows logarithmically with the system size. The anomalous behavior of m is traced to a violation of scaling in the order parameter probability density, which in…
▽ More
We study the static and dynamic behavior of the one dimensional pair contact process with diffusion. Several critical exponents are found to vary with the diffusion rate, while the order-parameter moment ratio m=\bar{rho^2} /\bar{rho}^2 grows logarithmically with the system size. The anomalous behavior of m is traced to a violation of scaling in the order parameter probability density, which in turn reflects the presence of two distinct sectors, one purely diffusive, the other reactive, within the active phase. Studies restricted to the reactive sector yield precise estimates for exponents beta and nu_perp, and confirm finite size scaling of the order parameter. In the course of our study we determine, for the first time, the universal value m_c = 1.334 associated with the parity-conserving universality class in one dimension.
△ Less
Submitted 30 July, 2002;
originally announced July 2002.
-
The Onset of Mackey-Glass Leukemia at the Edge of Chaos
Authors:
M. F. Menezes,
R. M. Zorzenon dos Santos
Abstract:
In this paper we revisit the Mackey-Glass model for blood-forming process, which was proposed to describe the spontaneous fluctuations of the blood cell counts in normal individuals and the first stage of chronic myelocytic (or granylocytic) leukemia (CML). We obtain the bifurcation diagram as a function of the time delay parameter and show that the onset of leukemia is related to instabilities…
▽ More
In this paper we revisit the Mackey-Glass model for blood-forming process, which was proposed to describe the spontaneous fluctuations of the blood cell counts in normal individuals and the first stage of chronic myelocytic (or granylocytic) leukemia (CML). We obtain the bifurcation diagram as a function of the time delay parameter and show that the onset of leukemia is related to instabilities associated to the presence of periodic windows in the midst of a chaotic regime. We also introduce a very simple modification in the death rate parameter in order to simulate the accumulation of cells and the progressive increase of the minima counts experimentally observed in the final stage of the disease in CML patients. The bifurcation diagram as a function of the death rate parameter is also obtained and we discuss the effects of treatments like leukapheresis.
△ Less
Submitted 12 July, 2001;
originally announced July 2001.
