FOM: A Framework for Metaheuristic Optimization

Most metaheuristic approaches for discrete optimization are usually implemented from scratch. In this paper, we introduce and discuss FOM, an object-oriented framework for metaheuristic optimization to be used as a general tool for the development and the implementation of metaheuristic algorithms. The basic idea behind the framework is to separate the problem side from the metaheuristic algorithms, allowing this to reuse different metaheuristic components in different problems. In addition to describing the design and functionality of the framework, we apply it to illustrative examples. Finally, we present our conclusions and discuss futures development

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.Comment: 45 pages, 15 figure

Soil-landscape and climatic relationships in the middle Miocene of the Madrid Basin

The Miocene alluvial-lacustrine sequences of the Madrid Basin, Spain, formed in highly varied landscapes. The presence of various types of palaeosols allows assessment of the effects of local and external factors onsedimentation, pedogenesis and geomorphological development. In the northern, more arid, tectonicallyactive arca, soils were weakly developed in aggrading alluvial fans, dominated by mass flows. reflecting high sedimentation rates. In more distal parts of the fans and in playa lakes calcretes and dolocretes developed: the former were associated with Mg-poor fan sediments whitc: the latter formed on Mg-rich lake clays exposed during minar lake lowstands. The nonh-east part of the basin had a less arid climate. Alluvial fans in this area were dominated by stream Aood deposits, sourced by carbonate terrains. Floodplain and freshwater lakc deposits formed in distal areas. The high local supply of calcium carbonate may have contributed to the preferential developmenl on calcretes on the fans. Both the fan and floodplain palaeosols exhibit pedofacies relationships and more mature soils developed in settings more distant from the sediment sources. Palaeosols also developed on pond and lake margin carbonates, and led to the formation of palustrine limestones. The spatial distributions and stratigraphies of palaeosols in the Madrid Basin alluvial fans suggest that soil formation was controlled by local factors. These palaeosols differ from those seen in Quatemary fans. Which are characterized by climatically induced periods of stability and instability

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further.Comment: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysic

Gaia Early Data Release 3 Acceleration of the Solar System from Gaia astrometry

Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar systembarycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. Theeffect of the acceleration was obtained as a part of the general expansion of the vector field of proper motions in vector spherical harmonics (VSH). Various versions of the VSH fit and various subsets of the sources were tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution were used to get a better idea of the possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is (2.32 +/- 0.16) x 10(-10) m s(-2) (or 7.33 +/- 0.51 km s(-1) Myr-1) towards alpha = 269.1 degrees +/- 5.4 degrees, delta = -31.6 degrees +/- 4.1 degrees, corresponding to a proper motion amplitude of 5.05 +/- 0.35 mu as yr(-1). This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 mu as yr(-1).Peer reviewe

Gaia Data Release 2 Mapping the Milky Way disc kinematics

Context. The second Gaia data release (Gaia DR2) contains high-precision positions, parallaxes, and proper motions for 1.3 billion sources as well as line-of-sight velocities for 7.2 million stars brighter than G(RVS) = 12 mag. Both samples provide a full sky coverage. Aims. To illustrate the potential of Gaia DR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. Methods. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes (sigma((omega) over bar)/(omega) over bar Results. Gaia DR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the U - V plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. Conclusions. Gaia DR2 provides the largest existing full 6D phase-space coordinates catalogue. It also vastly increases the number of available distances and transverse velocities with respect to Gaia DR1. Gaia DR2 offers a great wealth of information on the Milky Way and reveals clear non-axisymmetric kinematic signatures within the Galactic disc, for instance. It is now up to the astronomical community to explore its full potential.Peer reviewe

Gaia Focused Product Release: Radial velocity time series of long-period variables

The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged)Comment: 36 pages, 38 figure

Gaia Early Data Release 3: The celestial reference frame (Gaia-CRF3)

Context. Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. Aims. We describe the construction of Gaia-CRF3 and its properties in terms of the distributions in magnitude, colour, and astrometric quality. Methods. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasi-stellar objects (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio-optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). Results. Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13-21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 μas yr-1 on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but it exceeds 4 mas in either coordinate for 127 sources. We outline the future of Gaia-CRF in the next Gaia data releases. Appendices give further details on the external catalogues used, how to extract information about the Gaia-CRF3 sources, potential (Galactic) confusion sources, and the estimation of the spin and orientation of an astrometric solution

Pulsations in main sequence OBAF-type stars

CONTEXT: The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators (M ≥  1.3 M⊙) of spectral types O, B, A, or F, known as β Cep, slowly pulsating B (SPB), δ Sct, and γ Dor stars. These stars are often multi-periodic and display low amplitudes, making them challenging targets to analyse with sparse time series. AIMS: We investigate the extent to which the sparse Gaia DR3 data can be used to detect OBAF-type pulsators and discriminate them from other types of variables. We aim to probe the empirical instability strips and compare them with theoretical predictions. The most populated variability class is that of the δ Sct variables. For these stars, we aim to confirm their empirical period-luminosity (PL) relation, and verify the relation between their oscillation amplitude and rotation. METHODS: All datasets used in this analysis are part of the Gaia DR3 data release. The photometric time series were used to perform a Fourier analysis, while the global astrophysical parameters necessary for the empirical instability strips were taken from the Gaia DR3 gspphot tables, and the v sin i data were taken from the Gaia DR3 esphs tables. The δ Sct PL relation was derived using the same photometric parallax method as the one recently used to establish the PL relation for classical Cepheids using Gaia data. RESULTS: We show that for nearby OBAF-type pulsators, the Gaia DR3 data are precise and accurate enough to pinpoint them in the Hertzsprung-Russell (HR) diagram. We find empirical instability strips covering broader regions than theoretically predicted. In particular, our study reveals the presence of fast rotating gravity-mode pulsators outside the strips, as well as the co-existence of rotationally modulated variables inside the strips as reported before in the literature. We derive an extensive period–luminosity relation for δ Sct stars and provide evidence that the relation features different regimes depending on the oscillation period. We demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of δ Sct stars. CONCLUSIONS: The Gaia DR3 time-series photometry already allows for the detection of the dominant (non-)radial oscillation mode in about 100 000 intermediate- and high-mass dwarfs across the entire sky. This detection capability will increase as the time series becomes longer, allowing the additional delivery of frequencies and amplitudes of secondary pulsation modes