Date: Thu, 9 May 2024 | Total: 14
The aim of this project is to recover the CMB anisotropies maps in temperature and polarized intensity by means of a deep convolutional neural network (CNN) which, after appropiate training, can remove the foregrounds from Planck and QUIJOTE data. The results are then compared with those obtained by COMMANDER, based on Bayesian parametric component separation. The CNN successfully recovered the CMB signal for both All Sky and Partial Sky maps showing frequency dependant results, being optimum for central frequencies where there is less contamination by foregrounds emissions such as galactic synchrotron and thermal dust emissions. Recovered maps in temperature are consistent with those obtained by Planck Collaboration, while polarized intensity has been recovered as a new observable. The polarized intensity maps recovered from QUIJOTE experiment are novel and of potential interest to the scientific community for the detection of primordial gravitational waves. The way forward will be to recover the maps at higher NSIDE and make them available to the scientific community.
Artificial broadcasts from extraterrestrial intelligences (ETIs) are a hypothetical class of celestial phenomena. Unlike known astrophysical objects, the societies that generate them may be able to replicate on galactic scales through interstellar travel. Different galaxies could thus have drastically different populations, with abundance variations of many orders of magnitude. I present a probabilistic formalism to treat this shared history, in which societies and their broadcasts are described by distributions over basic properties like lifespan and energy released. The framework contains a hierarchy of objects related by a tree structure. Discrete societies, the sources of broadcasts, are organized into potentially interstellar "metasocieties." The population of each type of object is represented by a random point process in an abstract parameter hyperspace, a "haystack." When a selection like an observation draws a sample, the point process is thinned. Given assumptions of interchangeability and independence, observables are modeled with compound Poisson random variables. I present an example of how selection bias can favor sampling longer-lived objects. I rederive the Drake Equation for societies in the limit of no expansion. When interstellar replication is present, however, the mean number of detected broadcasts can depend quadratically on stellar mass, suggesting a search strategy favoring large galaxies.
The search for extraterrestrial intelligence includes efforts to constrain populations of artificial broadcasts in other galaxies. Previous efforts use individualist methods, searching for single broadcasts with high signal-to-noise ratio. These would be detected as observables with extreme values. This approach is limited to very bright broadcasts and also is subject to confusion, where a large number of broadcasts blend together to form a noise continuum. The mean value of the total emission provides an additional collective bound: the luminosity of the transmitters is no higher than the galaxy's observed luminosity. Using the framework developed in Paper I, I evaluate how confusion affects individualist searches. I then compare individualist and collective approaches for radio broadcasts from the Milky Way, M31, and three Virgo Cluster elliptical galaxies. For current observations, confusion blurs narrowband radio broadcasts together in the Virgo ellipticals when there is one broadcast per gigahertz per 1000 stars. The collective bound implies fewer than $\sim 10^6 (\overline{\ell}/10^{13} W)^{-1}$ L-band broadcasts per star gigahertz GHz in the Milky Way and is about 10 and 400 times stronger in M31 and M59, respectively. Applying the collective bound to the far-infrared--radio correlation yields constraints on radio broadcast populations in star-forming galaxies throughout the Universe. The collective bound allows us to rule out large regions of broadcast population parameter space even for distant galaxies. It also imposes constraints on gamma-ray, neutrino, and gravitational-wave broadcasts in the nearest galaxies.
The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called "Erebor," designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog $\sim12000$ Galactic Binaries ($\sim8000$ as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.
Trans-Neptunian objects (TNOs) are small, icy bodies in the outer solar system. They are observed to have a complex orbital distribution that was shaped by the early dynamical history and migration of the giant planets. Comparisons between the different dynamical classes of modeled and observed TNOs can help constrain the history of the outer solar system. Because of the complex dynamics of TNOs, particularly those in and near mean motion resonances with Neptune, classification has traditionally been done by human inspection of plots of the time evolution of orbital parameters. This is very inefficient. The Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) is expected to increase the number of known TNOs by a factor of $\sim$10, necessitating a much more automated process. In this chapter we present an improved supervised machine learning classifier for TNOs. Using a large and diverse training set as well as carefully chosen, dynamically motivated data features calculated from numerical integrations of TNO orbits, our classifier returns results that match those of a human classifier 98% of the time, and dynamically relevant classifications 99.7% of the time. This classifier is dramatically more efficient than human classification, and it will improve classification of both observed and modeled TNO data.
In this work, we assess the potential detectability of solar panels made of silicon on an Earth-like exoplanet as a potential technosignature. Silicon-based photovoltaic cells have high reflectance in the UV-VIS and in the near-IR, within the wavelength range of a space-based flagship mission concept like the Habitable Worlds Observatory (HWO). Assuming that only solar energy is used to provide the 2022 human energy needs with a land cover of ~2.4%, and projecting the future energy demand assuming various growth-rate scenarios, we assess the detectability with an 8 m HWO-like telescope. Assuming the most favorable viewing orientation, and focusing on the strong absorption edge in the ultraviolet-to-visible (0.34 - 0.52 um), we find that several 100s of hours of observation time is needed to reach a SNR of 5 for an Earth-like planet around a Sun-like star at 10pc, even with a solar panel coverage of ~23% land coverage of a future Earth. We discuss the necessity of concepts like Kardeshev Type I/II civilizations and Dyson spheres, which would aim to harness vast amounts of energy. Even with much larger populations than today, the total energy use of human civilization would be orders of magnitude below the threshold for causing direct thermal heating or reaching the scale of a Kardashev Type I civilization. Any extraterrrestrial civilization that likewise achieves sustainable population levels may also find a limit on its need to expand, which suggests that a galaxy-spanning civilization as imagined in the Fermi paradox may not exist.
Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of "true" redshift posteriors that can be used for photo-z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo-z estimation. We discuss how PZFlow will be used to validate the photo-z estimation pipeline of the Dark Energy Science Collaboration (DESC), and the wider applicability of PZFlow for statistical modeling of any tabular data.
The digital frequency domain multiplexing (DfMux) technique is widely used for astrophysical instruments with large detector arrays. Detailed detector characterization is required for instrument calibration and systematics control. We conduct the TES complex electrothermal-feedback (ETF) response measurement with the DfMux readout system as follows. By injecting a single sideband signal, we induce modulation in TES power dissipation over a frequency range encompassing the detector response. The modulated current signal induced by TES heating effect is measured, allowing for the ETF response characterization of the detector. With the injection of an upper sideband, the TES readout current shows both an upper and a lower sideband. We model the upper and lower sideband complex ETF response and verify the model by fitting to experimental data. The model not only can fit for certain physical parameters of the detector, such as loop gain, temperature sensitivity, current sensitivity, and time constant, but also enables us to estimate the systematic effect introduced by the multiplexed readout. The method is therefore useful for in-situ detector calibration and for estimating systematic effects during astronomical telescope observations, such as those performed by the upcoming LiteBIRD satellite.
We developed the X-ray, Gamma-ray and Relativistic Electron detector (XGRE) onboard the TARANIS satellite, to investigate high-energy phenomena associated with lightning discharges such as terrestrial gamma-ray flashes and terrestrial electron beams. XGRE consisted of three sensors. Each sensor has one layer of LaBr$_{3}$ crystals for X-ray/gamma-ray detections, and two layers of plastic scintillators for electron and charged-particle discrimination. Since 2018, the flight model of XGRE was developed, and validation and calibration tests, such as a thermal cycle test and a calibration test with the sensors onboard the satellite were performed before the launch of TARANIS on 17 November 2020. The energy range of the LaBr$_{3}$ crystals sensitive to X-rays and gamma rays was determined to be 0.04-11.6 MeV, 0.08-11.0 MeV, and 0.08-11.3 MeV for XGRE1, 2, and 3, respectively. The energy resolution at 0.662 MeV (full width at half maximum) was to be 20.5%, 25.9%, and 28.6%, respectively. Results from the calibration test were then used to validate a simulation model of XGRE and TARANIS. By performing Monte Carlo simulations with the verified model, we calculated effective areas of XGRE to X-rays, gamma rays, electrons, and detector responses to incident photons and electrons coming from various elevation and azimuth angles.
The Cherenkov Telescope Array (CTA) is a future ground-based observatory for gamma-ray astronomy providing unparalleled sensitivity in the energy range from 20 GeV up to 300 TeV. CTA will consist of telescopes with three different sizes. The Medium-Sized Telescopes (MSTs) will have 12 m reflectors with a tessellated mirror design of 86 mirror facets each. Each mirror facet is mounted on the mirror support structure with two actuators that are adjustable in length to align the mirrors, and a freely rotating fixpoint. Image resolution and pointing accuracy constraints impose limits on the backlash and deformation of the actuators and the fixpoint under various weight and wind loads. In this contribution, the test stand to measure the backlash and deformation behaviour of actuators and fixpoints is described and the measurement procedure is explained.
We present the ContEvol (continuous evolution) formalism, a family of implicit numerical methods which only need to solve linear equations and are almost symplectic. Combining values and derivatives of functions, ContEvol outputs allow users to recover full history and render full distributions. Using classic harmonic oscillator as a prototype case, we show that ContEvol methods lead to lower-order errors than two commonly used Runge--Kutta methods. Applying first-order ContEvol to simple celestial mechanics problems, we demonstrate that deviation from equation(s) of motion of ContEvol tracks is still $\mathcal{O}(h^5)$ ($h$ is the step length) by our definition. Numerical experiments with an eccentric elliptical orbit indicate that first-order ContEvol is a viable alternative to classic Runge--Kutta or the symplectic leapfrog integrator. Solving stationary Schr\"odinger equation in quantum mechanics, we manifest ability of ContEvol to handle boundary value or eigenvalue problems. Important directions for future work, including mathematical foundation, higher dimensions, and technical improvements, are discussed at the end of this article.
Microwave Kinetic Inductance Detectors (MKIDs) are photon detectors comprised of superconducting LC resonators with unique resonant frequencies corresponding to their geometrical structure. As each pixel has its own geometry, electromagnetic simulations by hand of every pixel in a kilo-pixel array are impractical. Simulating fewer pixels and interpolating in between risks reduced pixel yield in arrays due to overlapping resonant frequencies. We introduce a new software called AEM (Automated Electromagnetic MKID simulations) that automates the constructions and simulations of every simulated MKID pixel in an array according to specified resonant frequencies and a Qc range. We show automated designs to have an increased pixel yield (avoiding loses due to interpolation completely), increased accuracy in resonance frequency and Qc values when compared to interpolated structures. We also demonstrate a simulated trial of AEM for 100 MKIDs between 4 & 8 GHz to produce MKIDs with accuracies of +-0.2 MHz with a runtime of 10 hrs 45 mins.
We present an analysis of the LISA Pathfinder differential acceleration performance over the entire mission . We show that the Brownian noise level, detected for frequencies $f\gtrsim \SI{1}{mHz}$, has been evolving consistently with the outgassing of a single gaseous species, with an activation temperature of $(7.0\pm 0.2)\,\text{kK}$. In excess to the Brownian noise, the acceleration amplitude spectral density (ASD) always shows a sub-mHz tail which is reasonably well fit, between $f=\SI{36}{\micro\hertz}$ and $\SI{1}{\milli\hertz}$, to $\widetilde{S}_{\Delta g}^{1/2}(1\, \text{mHz}/f)$. A Bayesian estimate of $\widetilde{S}_{\Delta g}^{1/2}$ on a partition of the entire set of measurements in 27 data stretches, each 2.75\,d long, gives $\widetilde{S}_{\Delta g}^{1/2}=(1.1\pm0.3)\,\si{\femto\meter\,\second^{-2}/\rtHz}$, with no particular time pattern over the course of the mission. The width the posterior contains, in excess of the statistical uncertainty, a true physical fluctuation of $\widetilde{S}_{\Delta g}^{1/2}$ from run to run, of about $\SI{0.2}{\femto\meter\,\second^{-2}/\rtHz}$, with no correlation with specific operating conditions. At the lowest considered frequency of $f=\SI{18}{\micro\hertz}$, the ASD significantly deviates from the $1/f$ behavior, because of temperature fluctuations that appear to modulate a quasi-static pressure gradient, sustained by the asymmetries of outgassing . We also present a projection of acceleration noise on the sources for which we had either a correlation measurement, or an estimate from dedicated experiments.These sources account for about 40\% of the noise power the $1/f$ tail. We discuss the possible sources of the unaccounted-for fraction, present a series of analyses that rule many of them out, and identify the possible measures that may be taken to keep the remaining ones under control in LISA.
We introduce a novel technique to mitigate the adverse effects of atmospheric turbulence on astronomical imaging. Utilizing a video-to-image neural network trained on simulated data, our method processes a sliding sequence of short-exposure ($\sim$0.2s) stellar field images to reconstruct an image devoid of both turbulence and noise. We demonstrate the method with simulated and observed stellar fields, and show that the brief exposure sequence allows the network to accurately associate speckles to their originating stars and effectively disentangle light from adjacent sources across a range of seeing conditions, all while preserving flux to a lower signal-to-noise ratio than an average stack. This approach results in a marked improvement in angular resolution without compromising the astrometric stability of the final image.