The ionization fraction in neutral interstellar clouds is a key physical parameter controlling multiple physical and chemical processes, and varying by orders of magnitude from the UV irradiated ...surface of the cloud to its cosmic-ray dominated central regions. Traditional observational tracers of the ionization fraction, which mostly rely on deuteration ratios of molecules like HCO
+
, suffer from the fact that the deuterated molecules are only detected in a tiny fraction of a given Giant Molecular Cloud (GMC). In 1, we propose a
machine learning
-based, semi-automatic method to search in a large dataset of astrochemical model results for new tracers of the ionization fraction, and propose several new tracers relevant in different ranges of physical conditions.
In this paper, we propose a modification of an acoustic-transport operator splitting Lagrange-projection method for simulating compressible flows with gravity. The original method involves two steps ...that respectively account for acoustic and transport effects. Our work proposes a simple modification of the transport step, and the resulting modified scheme turns out to be a flux-splitting method. This new numerical method is less computationally expensive in the low-Mach regime, more memory efficient, and easier to implement than the original one. We prove stability properties for this new scheme by showing that under classical CFL conditions, the method is positivity preserving for mass, energy and entropy satisfying. The flexible flux-splitting structure of the method enables straightforward extensions of the method to multi-dimensional problems (with respect to space) and high-order discretizations that are presented in this work. We also propose an interpretation of the flux-splitting solver as a relaxation approximation. Both the stability and the accuracy of the new method are tested against one-dimensional and two-dimensional numerical experiments that involve highly compressible flows and low-Mach regimes.
•The all-regime acoustic/transport splitting solver proposed in 16 can be modified into a flux-splitting, one-step version.•The resulting method is less computationally expensive, more memory efficient, and easier to implement than the original one.•Stability properties are derived by interpreting the scheme as a convex combination of the advection/pressure steps.•This recasting allows an easy high-order extension as the new method is a standard flux-based solver.•Following the lines of 32, we equip our scheme with the well-balanced property for gravity source terms.
In this paper, we propose a modification of an acoustic-transport operator splitting Lagrange-projection method for simulating compressible flows with gravity. The original method involves two steps ...that respectively account for acoustic and transport effects. Our work proposes a simple modification of the transport step, and the resulting modified scheme turns out to be a flux-splitting method. This new numerical method is less computationally expensive, more memory efficient, and easier to implement than the original one. We prove stability properties for this new scheme by showing that under classical CFL conditions, the method is positivity preserving for mass, energy and entropy satisfying. The flexible flux-splitting structure of the method enables straightforward extensions of the method to multi-dimensional problems (with respect to space) and high-order discretizations that are presented in this work. We also propose an interpretation of the flux-splitting solver as a relaxation approximation. Both the stability and the accuracy of the new method are tested against one-dimensional and two-dimensional numerical experiments that involve highly compressible flows and low-Mach regimes.
The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer (WISE) led to the discovery of extremely cold, low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to ...lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST, ESO, Spitzer, and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put onto the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that nonequilibrium atmospheric chemistry is important for these objects. Nonequilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 ≤ Teff (K) ≤ 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 ≤ Teff (K) ≤ 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the Teff = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5-8 Jupiter masses and an age of 0.3-6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6-8.5 Gyr.
Atoms and molecules have long been thought to be versatile tracers of the cold neutral gas in the universe, from high-redshift galaxies to star forming regions and proto-planetary disks, because ...their internal degrees of freedom bear the signature of the physical conditions where these species reside. However, the promise that molecular emission has a strong diagnostic power of the underlying physical and chemical state is still hampered by the difficulty to combine sophisticated chemical codes with gas dynamics. It is therefore important 1) to acquire self-consistent data sets that can be used as templates for this theoretical work, and 2) to reveal the diagnostic capabilities of molecular lines accurately. The advent of sensitive wideband spectrometers in the (sub)- millimeter domain (e.g., IRAM-30m/EMIR, NOEMA, …) during the 2010s has allowed us to image a significant fraction of a Giant Molecular Cloud with enough sensitivity to detect tens of molecular lines in the 70 – 116 GHz frequency range. Machine learning techniques applied to these data start to deliver the next generation of molecular line diagnostics of mass, density, temperature, and radiation field.
Abstract
James Webb Space Telescope (JWST) mid-infrared spectra of cold brown dwarfs (BDs) are becoming available. In Leggett & Tremblin we demonstrated that an ATMO2020++ synthetic spectrum provided ...an excellent fit to the first published spectrum. Here we show that these models can reproduce the 1–14
μ
m spectra of BDs with a range of effective temperature: 350 ≤
T
eff
K ≤ 500. A grid of these models, which include disequilibrium chemistry and non-adiabatic pressure–temperature profiles, is available. In addition, a number of models without phosphine were calculated, as PH
3
is absent in the JWST data. Two of the six BDs studied are particularly interesting; one appears to be a young low-mass BD, the other is best fit as a binary. More complete analyses will be possible with an extended model grid, and with additional JWST data. The ATMO2020++ models are recommended for analyses of all BD data.
We present a geometric multigrid solver for the M1 model of radiative transfer without source terms. In radiative hydrodynamics applications, the radiative transfer needs to be solved implicitly ...because of the fast propagation speed of photons relative to the fluid velocity. The M1 model is hyperbolic and can be discretized with an HLL solver, whose time implicit integration can be done using a nonlinear Jacobi method. One can show that this iterative method always preserves the admissible states, such as positive radiative energy and reduced flux less than 1. To decrease the number of iterations required for the solver to converge, and therefore to decrease the computational cost, we propose a geometric multigrid algorithm. Unfortunately, this method is not able to preserve the admissible states. In order to preserve the admissible state states, we introduce a pseudo-time such that the solution of the problem on the coarse grid is the steady state of a differential equation in pseudo-time. We present preliminary results showing the decrease of the number of iterations and computational cost as a function of the number of multigrid levels used in the method. These results suggest that nonlinear multigrid methods can be used as a robust implicit solver for hyperbolic systems such as the M1 model.
•Jacobi method to solve the nonlinear M1 model for radiative transfer.•Geometric multigrid method to improve the convergence rate.•Pseudo-time to preserve the admissible states with large time steps.•Decrease of computational time compared to Newton-Raphson method.
Context.
The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection ...of potential temperature leads to a hotter adiabatic temperature profile in the deep atmosphere than the profile obtained with standard 1D models.
Aims.
In order to confirm the viability of that scenario, we extend this investigation to 3D, time-dependent models.
Methods.
We use a 3D general circulation model DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both the upper and deep atmospheric forcing.
Results.
In agreement with the previous, 2D study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of the order of 1500 yr are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling timescale is longer than ~3000 yr at 200 bar.
Conclusions.
We conclude that steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain the inflated radii of hot Jupiters. We suggest that future models of hot Jupiters be evolved for a longer time than currently done, and when possible that models initialised with a hot deep adiabat be included. We stress that this mechanism stems from the advection of entropy by irradiation-induced mass flows and does not require a (finely tuned) dissipative process, in contrast with most previously suggested scenarios.
We present here a reanalysis of the Spitzer Space Telescope phase curves of the hot Jupiter WASP43 b, using the wavelet pixel-independent component analysis, a blind signal-source separation method. ...The data analyzed were recorded with the Infrared Array Camera and consisted of two visits at 3.6 m, and one visit at 4.5 m, each visit covering one transit and two eclipse events. To test the robustness of our technique we repeated the analysis on smaller portions of the phase curves, and by employing different instrument ramp models. Our reanalysis presents significant updates of the planetary parameters compared to those reported in the original phase curve study of WASP43 b. In particular, we found (1) higher nightside temperatures, (2) smaller hotspot offsets, (3) a greater consistency (∼1 ) between the two 3.6 m visits, and (4) a greater similarity with the predictions of the atmospheric circulation models. Our parameter results are consistent within 1 with those reported by a recent reanalysis of the same data sets. For each visit we studied the variation of the retrieved transit parameters as a function of various sets of stellar limb-darkening coefficients, finding significant degeneracy between the limb-darkening models and the analysis output. Furthermore, we performed the analysis of the single transit and eclipse events, and we examined the differences between these results with the ones obtained with the whole phase curve. Finally we provide a formula useful to optimize the trade-off between precision and duration of observations of transiting exoplanets.