Search

The Orion Nebula Cluster (ONC) is the closest site of very young ($\sim$ 1 Myrs) massive star formation. The ONC hosts more than 1600 young and X-ray bright stars with masses ranging from $\sim$ 0.1 to 35 $M_\odot$. The Chandra HETGS Orion Legacy Project observed the ONC with the Chandra high energy transmission grating spectrometer (HETGS) for $2.1\,$Ms. We describe the spectral extraction and cleaning processes necessary to separate overlapping spectra. We obtained 36 high resolution spectra which includes a high brilliance X-ray spectrum of $\theta^1$ Ori C with over 100 highly significant X-ray lines. The lines show Doppler broadening between 300 and $400\;\mathrm{km}\;\mathrm{s}^{-1}$. Higher spectral diffraction orders allow us to resolve line components of high Z He-like triplets in $\theta^1$ Ori C with unprecedented spectral resolution. Long term light curves spanning $\sim$20 years show all stars to be highly variable, including the massive stars. Spectral fitting with thermal coronal emission line models reveals that most sources show column densities of up to a few times $10^{22}\,$cm$^{-2}$ and high coronal temperatures of 10 to 90 MK. We observe a bifurcation of the high temperature component where some stars show a high component of 40 MK, while others show above 60 MK indicating heavy flaring activity. Some lines are resolved with Doppler broadening above our threshold of $\sim200\;\mathrm{km}\;\mathrm{s}^{-1}$, up to $500\;\mathrm{km}\;\mathrm{s}^{-1}$. This data set represents the largest collection of HETGS high resolution X-ray spectra from young pre-MS stars in a single star-forming region to date., Comment: Accepted by ApJ, fixed typo in the xlabel of Fig 12 (left panel)

High-resolution X-ray spectra of $\pi\,$Aqr, a $\gamma\,$Cas-type star, obtained with the Chandra/HETG grating spectrometer, revealed emission lines of H-like ions of Mg, Si, S, and Fe, a strong, hard continuum, and a lack of He-like ions, indicating the presence of very hot thermal plasma. The X-ray light curve showed significant fluctuations, with coherent variability at period of about 3400 seconds in one observation. The hardness ratio was relatively constant except for one observation in which the spectrum was much harder and more absorbed. We interpret the X-ray emission as arising from accretion onto the secondary, which is likely a magnetic white dwarf, an intermediate polar system., Comment: 7 pages, 3 figures, 3 tables

Evidence suggests that when compact objects such as black holes and neutron stars form, they may receive a ``natal kick,'' where the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial, yet limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks. Here, we report the discovery that the canonical black hole low-mass X-ray binary V404 Cygni is part of a wide hierarchical triple with a tertiary companion at least 3500 astronomical units away from the inner binary. Given the orbital configuration, the black hole likely received a sub-5 kilometer per second kick to have avoided unbinding the tertiary. This discovery reveals that at least some black holes form with nearly no natal kick. Furthermore, the tertiary in this system lends credence to evolutionary models of low-mass X-ray binaries involving a hierarchical triple structure. Remarkably, the tertiary is evolved, indicating that the system formed 3-5 billion years ago, and that the black hole has removed at least half a solar mass of matter from its evolved secondary companion. During the event in which the black hole formed, it is likely that at least half of the mass of the black hole progenitor collapsed into the black hole; it may even have undergone a complete implosion, enabling the tertiary to remain loosely bound., Comment: 36 pages, 10 figures, 1 table

Solar radio bursts (SRBs), are emitted by electrons propagating through the corona and interplanetary space. Tracking such bursts is key to understanding the properties of accelerated electrons and radio wave propagation as well as the local plasma environment that they propagate through. Here, we present a novel multilateration algorithm called BayEsian LocaLisation Algorithm (BELLA). In addition, apparent SRB positions from BELLA are compared with comparable localisation methods and the predictions of solar wind models. BELLA uses Bayesian inference to create probabilistic distributions of source positions and their uncertainties. This facilitates the estimation of algorithmic, instrumental, and physical uncertainties in a quantitative manner. We validated BELLA using simulations and a Type III SRB observed by STEREO A/B and Wind. BELLA tracked the Type III source from $\sim$ 10--150 $R_{sun}$ (2-0.15 MHz) along a spiral trajectory. This allowed for an estimate of an apparent solar wind speed of $v_{sw} \sim$ 400 km s$^{-1}$ and a source longitude of $\phi_0 \sim$ 30deg. We compared these results with well-established methods of positioning: Goniopolarimetric (GP), analytical time-difference-of-arrival (TDOA), and Solar radio burst Electron Motion Tracker (SEMP). We found them to be in agreement with the results obtained by BELLA. Additionally, the results aligned with solar wind properties assimilated by the Heliospheric Upwind Extrapolation with time dependence (HUXt) model. We have validated BELLA and used it to identify apparent source positions as well as velocities and densities of the solar wind. Furthermore, we identified higher than expected electron densities, suggesting that the true emission sources were at lower altitudes than those identified by BELLA, an effect that may be due to appreciable scattering of electromagnetic waves by electrons in interplanetary space., Comment: 13 pages, 12 figures, in press, accepted by Astronomy and Astrophysics for publication

X-ray binaries are known to launch powerful accretion disk winds that can have significant impact on the binary systems and their surroundings. To quantify the impact and determine the launching mechanisms of these outflows, we need to measure the wind plasma number density, an important ingredient in the theoretical disk wind models. While X-ray spectroscopy is a crucial tool to understanding the wind properties, such as their velocity and ionization, in nearly all cases, we lack the signal-to-noise to constrain the plasma number density, weakening the constraints on outflow location and mass outflow rate. We present a new approach to determine this number density in the X-ray binary Hercules X-1 by measuring the speed of the wind ionization response to time-variable illuminating continuum. Hercules X-1 is powered by a highly magnetized neutron star, pulsating with a period of 1.24 s. We show that the wind number density in Hercules X-1 is sufficiently high to respond to these pulsations by modeling the ionization response with the time-dependent photoionization model TPHO. We then perform a pulse-resolved analysis of the best-quality XMM-Newton observation of Hercules X-1 and directly detect the wind response, confirming that the wind density is at least $10^{12}$ cm$^{-3}$. Finally, we simulate XRISM observations of Hercules X-1 and show that they will allow us to accurately measure the number density at different locations within the outflow. With XRISM we will rule out $\sim3$ orders of magnitude in density parameter space, constraining the wind mass outflow rate, energetics, and its launching mechanism., Comment: Accepted in ApJ. 24 pages, 12 figures, 3 tables

The focus of this paper is the study of the inverse point-source scattering problem, specifically in relation to a certain class of electric potentials. Our research provides a novel uniqueness result for the inverse problem with local data, obtained from the near field pattern. Our work improves the work of Caro and Garcia, who investigated both the direct problem and the inverse problem with global near field data for critically singular and $\delta$-shell potentials. The primary contribution of our research is the introduction of a Runge approximation result for the near field data on the scattering problem which, in combination with an interior regularity argument, enables us to establish a uniqueness result for the inverse problem with local data. Additionaly, we manage to consider a slightly wider class of potentials., Comment: 33 pages. Preprint

This manuscript presents novel techniques for identifying the switch states, phase identification, and estimation of equipment parameters in multi-phase low voltage electrical grids, which is a major challenge in long-standing German low voltage grids that lack observability and are heavily impacted by modelling errors. The proposed methods are tailored for systems with a limited number of spatially distributed measuring devices, which measure voltage magnitudes at specific nodes and some line current magnitudes. The overall approach employs a problem decomposition strategy to divide the problem into smaller subproblems, which are addressed independently. The techniques for identifying switch states and system phases are based on heuristics and a binary optimization problem using correlation analysis of the measured time series. The estimation of equipment parameters is achieved through a data-driven regression approach and by an optimization problem, and the identification of cable types is solved using a Mixed-Integer Quadratic Programming solver. To validate the presented methods, a realistic grid is used and the presented techniques are evaluated for their resilience to data quality and time resolution, discussing the limitations of the proposed methods.

One of the greatest challenge facing current space weather monitoring operations is forecasting the arrival of coronal mass ejections (CMEs) and Solar Energetic Particles (SEPs) within their Earth-Sun propagation timescales. Current campaigns mainly rely on extreme ultra-violet and white light observations to create forecasts, missing out many potential events that may be hazardous to Earth's infrastructure undetectable at these wavelengths. Here we introduce the SURROUND mission, a constellation of CubeSats each with identical radio spectrometers, and the results of the initial Phase-0 study for the concept. The main goal of SURROUND is to monitor and track solar radio bursts (SRBs), widely utilised as a useful diagnostic for space weather activity, and revolutionise current forecasting capabilities. The Phase-0 study concludes that SURROUND can achieve its mission objectives using 3 - 5 spacecraft using current technologies with feasible SEP and CME forecasting potential: a first for heliospheric monitors., Comment: 12 pages, 3 figures, 2 tables, submitted to Planetary, Solar and Heliospheric Radio Emissions IX, held 26-28 September, 2022 [ACCEPTED]; Edited by: C. K. Louis, C. M. Jackman, G. Fischer, A. H. Sulaiman, P. Zucca, Publishers: Dublin Institute for Advanced Studies and Trinity College Dublin. (these authors contributed equally to this work: D. M. Weigt, L. A. Ca\~nizares)

This paper presents novel methods for parameter identification in electrical grids with small numbers of spatially distributed measuring devices, which is an issue for distribution system operators managing aged and not properly mapped underground Low Voltage (LV) grids, especially in Germany. For this purpose, the total impedance of individual branches of the overall system is estimated by measuring currents and voltages at a subset of all system nodes over time. It is shown that, under common assumptions for electrical distsribution systems, an estimate of the total impedance can be made using readily computable proxies. Different regression methods are then used and compared to estimate the total impedance of the respective branches, with varying weights of the input data. The results on realistic LV feeders with different branch lengths and number of unmeasured segments are discussed and multiple influencing factors are investigated through simulations. It is shown that estimates of the total impedances can be obtained with acceptable quality under realistic assumptions.

The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether by radiation pressure, magnetic forces, thermal irradiation, or a combination thereof. An important means of studying disk winds produced near the central compact object is through X-ray absorption line spectroscopy, which allows us to probe outflow properties along a single line of sight, but usually provides little information about the global 3D disk wind structure that is vital for understanding the launching mechanism and total wind energy budget. Here, we study Hercules X-1, a unique, nearly edge-on X-ray binary with a warped accretion disk precessing with a period of about 35 days. This disk precession results in changing sightlines towards the neutron star, through the ionized outflow. We perform time-resolved X-ray spectroscopy over the precession phase and detect a strong decrease in the wind column density by three orders of magnitude as our sightline progressively samples the wind at greater heights above the accretion disk. The wind becomes clumpier as it rises upwards and expands away from the neutron star. Modelling the warped disk shape, we create a 2D map of wind properties. This unique measurement of the vertical structure of an accretion disk wind allows direct comparisons to 3D global simulations to reveal the outflow launching mechanism., Comment: Published in Nature Astronomy on April 10, available at https://www.nature.com/articles/s41550-023-01929-7

Integrated information theory (IIT) is a theoretical framework that provides a quantitative measure to estimate when a physical system is conscious, its degree of consciousness, and the complexity of the qualia space that the system is experiencing. Formally, IIT rests on the assumption that if a surrogate physical system can fully embed the phenomenological properties of consciousness, then the system properties must be constrained by the properties of the qualia being experienced. Following this assumption, IIT represents the physical system as a network of interconnected elements that can be thought of as a probabilistic causal graph, $\mathcal{G}$, where each node has an input-output function and all the graph is encoded in a transition probability matrix. Consequently, IIT's quantitative measure of consciousness, $\Phi$, is computed with respect to the transition probability matrix and the present state of the graph. In this paper, we provide a random search algorithm that is able to optimize $\Phi$ in order to investigate, as the number of nodes increases, the structure of the graphs that have higher $\Phi$. We also provide arguments that show the difficulties of applying more complex black-box search algorithms, such as Bayesian optimization or metaheuristics, in this particular problem. Additionally, we suggest specific research lines for these techniques to enhance the search algorithm that guarantees maximal $\Phi$.

Although electric power networks and district heating networks are physically coupled, they are not operated in a coordinated manner. With increasing penetration of renewable energy sources, a coordinated market-based operation of the two networks can yield significant advantages, as reduced need for grid reinforcements, by optimizing the power flows in the coupled systems. Transactive control has been developed as a promising approach based on market and control mechanisms to coordinate supply and demand in energy systems, which when applied to power systems is being referred to as transactive energy. However, this approach has not been fully investigated in the context of market-based operation of coupled electric power and district heating networks. Therefore, this paper proposes a transactive control approach to coordinate flexible producers and consumers while taking into account the operational aspects of both networks, for the benefit of all participants and considering their privacy. A nonlinear model predictive control approach is applied in this work to maximize the social welfare of both networks, taking into account system operational limits, while reducing losses and considering system dynamics and forecasted power supply and demand of inflexible producers and consumers. A subtle approximation of the operational optimization problem is used to enable the practical application of the proposed approach in real time. The presented technique is implemented, tested, and demonstrated in a realistic test system, illustrating its benefits., Comment: 35 pages, 16 Figures

Hercules X-1 is a nearly edge-on accreting X-ray pulsar with a warped accretion disk, precessing with a period of about 35 days. The disk precession allows for unique and changing sightlines towards the X-ray source. To investigate the accretion flow at a variety of sightlines, we obtained a large observational campaign on Her X-1 with XMM-Newton (380 ks exposure) and Chandra (50 ks exposure) for a significant fraction of a single disk precession cycle, resulting in one of the best datasets taken to date on a neutron star X-ray binary. Here we present the spectral analysis of the High State high-resolution grating and CCD datasets, including the extensive archival data available for this famous system. The observations reveal a complex Fe K region structure, with three emission line components of different velocity widths. Similarly, the high-resolution soft X-ray spectra reveal a number of emission lines of various widths. We correct for the uncertain gain of the EPIC-pn Timing mode spectra, and track the evolution of these spectral components with Her X-1 precession phase and observed luminosity. We find evidence for three groups of emission lines: one originates in the outer accretion disk (10^5 RG from the neutron star). The second line group plausibly originates at the boundary between the inner disk and the pulsar magnetosphere (10^3 RG). The last group is too broad to arise in the magnetically-truncated disk and instead must originate very close to the neutron star surface, likely from X-ray reflection from the accretion curtain (~10^2 RG)., Comment: Accepted for publication in ApJ. 26 pages, 17 figures, 4 tables

In recent years, the electrification of Canadian Remote Communities (RCs) has received significant attention, as their current electric energy systems are not only expensive, but are also highly polluting due to the prevalence of diesel generators. In addition, RCs' inherent geographic characteristics impose a series of challenges that must be considered when planning their electricity supply. Thus, in this paper, an optimization model for the long-term planning of RC Microgrids (MGs) including Renewable Energy Sources (RESs) and Energy Storage Systems (ESSs) is proposed, with the objective of reducing costs and emissions. The proposed model considers lithium-ion batteries and hydrogen systems as part of ESSs technologies. The model is used to investigate the feasibility of integrating RESs and ESSs in an MG in Sanikiluaq, an RC in the Nunavut territory in Northern Canada. The results show that wind resources along with solar and storage technologies can play a key role in satisfying RC electricity demand, while significantly reducing costs and Greenhouse Gas Emissions (GHG). In addition, insights on sustainable and affordable policies for RC MGs are provided., Comment: In Proceedings of the 11th Bulk Power Systems Dynamics and Control Symposium (IREP 2022), July 25-30, 2022, Banff, Canada

This paper presents an Energy Management System (EMS) that considers power exchanges between a set of interconnected microgrids (MGs) and the main grid, in the context of Multi-MG (MMG) systems. The model is first formulated as a centralized optimization problem, which is then decomposed into subproblems corresponding to each MG, using Lagrangian relaxation, and solved through a distributed approach using a subgradient method. The proposed model determines the power exchanges minimizing the operation cost of each MG, considering grid constraints and preserving the privacy of each MG by not revealing their generation cost and demand information. The distributed approach is validated with respect to the centralized problem, and various case studies are presented to demonstrate the performance of the proposed approach, comparing the costs of the MGs operating individually and cooperatively. The results show that all MGs in the MMG system improve their cost as consequence of the power exchanges, thus demonstrating the advantages of interconnecting MGs., Comment: In proceedings of the 11th Bulk Power Systems Dynamics and Control Symposium (IREP 2022), July 25-30, 2022, Banff, Canada

Type III radio bursts are the result of plasma emission from mildly relativistic electron beams propagating from the low solar corona into the heliosphere where they can eventually be detected in situ if they align with the location of a heliospheric spacecraft. Here we observe a type III radio burst from 0.1-16 MHz using the Parker Solar Probe (PSP) FIELDS Radio Frequency Spectrometer (RFS), and from 20-80 MHz using the Low Frequency Array (LOFAR). This event was not associated with any detectable flare activity but was part of an ongoing type III and noise storm that occurred during PSP encounter 2. A deprojection of the LOFAR radio sources into 3D space shows that the type III radio burst sources were located on open magnetic field from 1.6-3 R$_\odot$ and originated from a near-equatorial active region around longitude E48$^o$. Combining PSP/RFS observations with WIND/WAVES and STEREO/WAVES, we reconstruct the type III radio source trajectory in the heliosphere interior to PSP's position, assuming ecliptic confinement. An energetic electron enhancement is subsequently detected in situ at the STEREO-A spacecraft at compatible times although the onset and duration suggests the individual burst contributes a subset of the enhancement. This work shows relatively small-scale flux emergence in the corona can cause the injection of electron beams from the low corona into the heliosphere, without needing a strong solar flare. The complementary nature of combined ground and space-based radio observations, especially in the era of PSP, is also clearly highlighted by this study., Comment: 16 pages, 11 figures, Submitted to ApJ: April 15 2022, Accepted: September 7 2022

Modern radio interferometers such as the LOw Frequency ARray (LOFAR) are capable of producing data at hundreds of gigabits to terabits per second. This high data rate makes the analysis of radio data cumbersome and computationally expensive. While high performance computing facilities exist for large national and international facilities, that may not be the case for instruments operated by a single institution or a small consortium. Data rates for next generation radio telescopes are set to eclipse those currently in operation, hence local processing of data will become all the more important. Here, we introduce the REAL-time Transient Acquisition backend (REALTA), a computing backend at the Irish LOFAR station (I-LOFAR) which facilitates the recording of data in near real-time and post-processing. We also present first searches and scientific results of a number of radio phenomena observed by I-LOFAR and REALTA, including pulsars, fast radio bursts (FRBs), rotating radio transients (RRATs), the search for extraterrestrial intelligence (SETI), Jupiter, and the Sun., Comment: 12 pages, 10 figures, accepted for publication in Astronomical instrumentation section of Astronomy & Astrophysics 24/08/2021

We report the temperature dependent THz to mid-infrared response of hexagonal-ErMnO3 using absorption, reflectivity, and emissivity techniques from 2 K to 1700 K. At low temperatures, lowest frequency vibrational modes coexist with paramagnon excitations associated with well-defined crystal field Rare Earth pure magnetic replicas in an intriguing phonon-magnetic convergence. Increasing the temperature, a number of vibrational bands close to the space group predicted undergo profile broadening and softening. In particular, a distinctive set of bands in the 288-329 cm-1 (300 K) range has a component whose profile is carrier screened becoming nearly fully blurred in the intermediate phase between ~830 K and ~1500 K. Below TC ~830 K this asymmetric band further splits as spin phonon interaction and the tripling of the unit cell takes place revealing a delicate balance of long- and short-range interactions. Ambient Raman scattering brings up evidence of a Fr\"ohlich resonance due to Coulomb interactions between carriers and the macroscopic field of the longitudinal optical phonon mode. We found it is dynamically correlated to the hexagonal c-axis negative thermal expansion. Mid-infrared (MIR) optical conductivity show that also plays a role in small polarons and mediates in high temperature bipolarones.At higher temperatures a low frequency Drude contribution is triggered by electron hopping signaling an insulator-metal phase transition at ~1600 K while the MIR response suggests coexistence between single small polarons and bipolarons. We draw a parallel with improper ferroelectrics sustaining a lattice incommensurate intermediate phase and unit cell tripling. We argue that in the h-RMnO3 (R=Rare Earth, Y) family of compounds the intermediate phase be considered incommensurate with onset at TINC ~1500 K and ferroelectric lock-in at TC ~830 K delimiting this regime in h-ErMnO3., Comment: Manuscript and Supplemental Material (66 pages)