A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

Classifying GRB 170817A/GW170817 in a Fermi duration–hardness plane

GRB 170817A, associated with the LIGO-Virgo GW170817 neutron-star merger event, lacks the short duration and hard spectrum of a Short gamma-ray burst (GRB) expected from long-standing classification models. Correctly identifying the class to which this burst belongs requires comparison with other GRBs detected by the Fermi GBM. The aim of our analysis is to classify Fermi GRBs and to test whether or not GRB 170817A belongs—as suggested—to the Short GRB class. The Fermi GBM catalog provides a large database with many measured variables that can be used to explore gamma-ray burst classification. We use statistical techniques to look for clustering in a sample of 1298 gamma-ray bursts described by duration and spectral hardness. Classification of the detected bursts shows that GRB 170817A most likely belongs to the Intermediate, rather than the Short GRB class. We discuss this result in light of theoretical neutron-star merger models and existing GRB classification schemes. It appears that GRB classification schemes may not yet be linked to appropriate theoretical models, and that theoretical models may not yet adequately account for known GRB class properties. We conclude that GRB 170817A may not fit into a simple phenomenological classification scheme.     

Propagation of Leaky MHD Waves at Discontinuities with Tilted Magnetic Field

We investigate the characteristics of magneto-acoustic surface waves propagating at a single density interface, in the presence of an inclined magnetic field. For linear wave propagation, the dispersion relation is obtained and analytical solutions are derived for small inclination angle. The inclination of the field renders the frequency of the waves complex, where the imaginary part describes wave attenuation, due to lateral energy leakage.     

A recent history of science cases for optical interferometry

Optical long-baseline interferometry is a unique and powerful technique for astronomical research. Since the 1980’s (with I2T, GI2T, Mark I to III, SUSI, ...), optical interferometers have produced an increasing number of scientific papers covering various fields of astrophysics. As current interferometric facilities are reaching their maturity, we take the opportunity in this paper to summarize the conclusions of a few key meetings, workshops, and conferences dedicated to interferometry. We present the most persistent recommendations related to science cases and discuss some key technological developments required to address them. In the era of extremely large telescopes, optical long-baseline interferometers will remain crucial to probe the smallest spatial scales and make breakthrough discoveries.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

Fuel optimization for low-thrust Earth–Moon transfer via indirect optimal control

The problem of designing low-energy transfers between the Earth and the Moon has attracted recently a major interest from the scientific community. In this paper, an indirect optimal control approach is used to determine minimum-fuel low-thrust transfers between a low Earth orbit and a Lunar orbit in the Sun–Earth–Moon Bicircular Restricted Four-Body Problem. First, the optimal control problem is formulated and its necessary optimality conditions are derived from Pontryagin’s Maximum Principle. Then, two different solution methods are proposed to overcome the numerical difficulties arising from the huge sensitivity of the problem’s state and costate equations. The first one consists in the use of continuation techniques. The second one is based on a massive exploration of the set of unknown variables appearing in the optimality conditions. The dimension of the search space is reduced by considering adapted variables leading to a reduction of the computational time. The trajectories found are classified in several families according to their shape, transfer duration and fuel expenditure. Finally, an analysis based on the dynamical structure provided by the invariant manifolds of the two underlying Circular Restricted Three-Body Problems, Earth–Moon and Sun–Earth is presented leading to a physical interpretation of the different families of trajectories.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Eight luminous early-type galaxies in nearby pairs and sparse groups. I. Stellar populations spatially analysed

We present a detailed spatial analysis of stellar populations based on long-slit optical spectra in a sample of eight luminous early-type galaxies selected from nearby sparse groups and pairs, three of them may have interaction with another galaxy of similar mass. We have spatially measured luminosity-weighted averages of age, [M/H], [Fe/H], and [\(\alpha /\mbox{Fe}\)] in the sample galaxies to add empirical data relative to the influence of galaxy mass, environment, interaction, and AGN feedback in their formation and evolution. The stellar population of the individual galaxies were determined through the well-established stellar population synthesis code starlight using semi-empirical simple stellar population models. Radial variations of luminosity- weighted means of age, [M/H], [Fe/H], and [\(\alpha /\mbox{Fe}\)] were quantified up to half of the effective radius of each galaxy. We found trends between representative values of age, [M/H], [\(\alpha /\mbox{Fe}\)], and the nuclear stellar velocity dispersion. There are also relations between the metallicity/age gradients and the velocity dispersion. Contributions of 1–4 Gyr old stellar populations were quantified in IC?5328 and NGC?6758 as well as 4–8 Gyr old ones in NGC?5812. Extended gas is present in IC?5328, NGC?1052, NGC?1209, and NGC?6758, and the presence of a LINER is identified in all these galaxies. The regions up to one effective radius of all galaxies are basically dominated by \(\alpha \)-enhanced metal-rich old stellar populations likely due to rapid star formation episodes that induced efficient chemical enrichment. On average, the age and [\(\alpha /\mbox{Fe}\)] gradients are null and the [M/H] gradients are negative, although discordant cases were found. We found no correlation between the stellar population properties and the LINER presence as well as between the stellar properties and environment or gravitational interaction, suggesting that the influence of progenitor mass cannot be discarded in the formation and evolution of early-type galaxies.     

Implications for behavior of volatile elements during impacts—Zinc and copper systematics in sediments from the Ries impact structure and central European tektites

Moldavites are tektites genetically related to the Ries impact structure, located in Central Europe, but the source materials and the processes related to the chemical fractionation of moldavites are not fully constrained. To further understand moldavite genesis, the Cu and Zn abundances and isotope compositions were measured in a suite of tektites from four different substrewn fields (South Bohemia, Moravia, Cheb Basin, Lusatia) and chemically diverse sediments from the surroundings of the Ries impact structure. Moldavites are slightly depleted in Zn (~10–20%) and distinctly depleted in Cu (>90%) relative to supposed sedimentary precursors. Moreover, the moldavites show a wide range in δ⁶⁶Zn values between 1.7 and 3.7‰ (relative to JMC 3‐0749 Lyon) and δ⁶⁵Cu values between 1.6 and 12.5‰ (relative to NIST SRM 976) and are thus enriched in heavy isotopes relative to their possible parent sedimentary sources (δ⁶⁶Zn = ?0.07 to +0.64‰; δ⁶⁵Cu = ?0.4 to +0.7‰). In particular, the Cheb Basin moldavites show some of the highest δ⁶⁵Cu values (up to 12.5‰) ever observed in natural samples. The relative magnitude of isotope fractionation for Cu and Zn seen here is opposite to oxygen‐poor environments such as the Moon where Zn is significantly more isotopically fractionated than Cu. One possibility is that monovalent Cu diffuses faster than divalent Zn in the reduced melt and diffusion will not affect the extent of Zn isotope fractionation. These observations imply that the capability of forming a redox environment may aid in volatilizing some elements, accompanied by isotope fractionation, during the impact process. The greater extent of elemental depletion, coupled with isotope fractionation of more refractory Cu relative to Zn, may also hinge on the presence of carbonyl species of transition metals and electromagnetic charge, which could exist in the impact‐induced high‐velocity jet of vapor and melts.     

Porous,S‐bearing silica in metal‐sulfide nodules and in the interchondrule clastic matrix in two EH3 chondrites

Two new occurrences of porous, S‐bearing, amorphous silica are described within metal‐sulfide nodules (MSN) and as interchondrule patches in EH3 chondrites SAH 97072 and ALH 84170. This porous amorphous material, which was first reported from sulfide‐bearing chondrules, consists of sinewy SiO₂‐rich areas containing S with minor Na or Ca as well as Fe, Mg, and Al. Some pores contain minerals including pyrite, pyrrhotite, and anhydrite. Most pores appear vacant or contain unidentified material that is unstable under analytical conditions. Niningerite, olivine, enstatite, albite, and kumdykolite occur enclosed within porous silica patches. Porous silica is commonly interfingered with cristobalite suggesting its amorphous structure resulted from high‐temperature quenching. We interpret the S‐bearing porous silica to be a product of silicate sulfidation, and the Na, Ca, Fe, Mg, and Al detectable within this material are chemical residues of sulfidized silicates and metal. The occurrence of porous silica in the cores of MSN, which are considered to be pre‐accretionary objects, suggests the sulfidizing conditions occurred prior to final parent‐body solidification. Ubiquitous S‐bearing porous silica among sulfide‐bearing chondrules, MSN, and in the interchondrule clastic matrix, suggests that similar sulfidizing conditions affected all the constituents of these EH3 chondrites.