Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 11, August
Previous Issue
Volume 11, April
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 

Galaxies, Volume 11, Issue 3 (June 2023) – 19 articles

Cover Story (view full-size image): Z CMa has intrigued astronomers for decades. It is an active early-type binary with a Herbig Be primary and an FU Orionis-type secondary. Both of the stars exhibit sub-arcsecond jet-like ejecta. In addition, the primary is associated with an extended jet as well as with a large-scale outflow. In this paper, we further investigate the nature of the large-scale outflow, which has not been studied since its discovery almost three and a half decades ago. We present proper motion measurements of individual features of the large-scale outflow and determine their kinematical ages. Furthermore, with our newly acquired deep images, we have discovered additional faint arc-shaped features that can be associated with the central binary. View this paper
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Select all
Export citation of selected articles as:
21 pages, 1509 KiB  
Article
Discovering New B[e] Supergiants and Candidate Luminous Blue Variables in Nearby Galaxies
by Grigoris Maravelias, Stephan de Wit, Alceste Z. Bonanos, Frank Tramper, Gonzalo Munoz-Sanchez and Evangelia Christodoulou
Galaxies 2023, 11(3), 79; https://doi.org/10.3390/galaxies11030079 - 19 Jun 2023
Cited by 2 | Viewed by 1551
Abstract
Mass loss is one of the key parameters that determine stellar evolution. Despite the progress we have achieved over the last decades we still cannot match the observational derived values with theoretical predictions. Even worse, there are certain phases, such as the B[e] [...] Read more.
Mass loss is one of the key parameters that determine stellar evolution. Despite the progress we have achieved over the last decades we still cannot match the observational derived values with theoretical predictions. Even worse, there are certain phases, such as the B[e] supergiants (B[e]SGs) and the Luminous Blue Variables (LBVs), where significant mass is lost through episodic or outburst activity. This leads to various structures forming around them that permit dust formation, making these objects bright IR sources. The ASSESS project aims to determine the role of episodic mass in the evolution of massive stars, by examining large numbers of cool and hot objects (such as B[e]SGs/LBVs). For this purpose, we initiated a large observation campaign to obtain spectroscopic data for ∼1000 IR-selected sources in 27 nearby galaxies. Within this project we successfully identified seven B[e] supergiants (one candidate) and four Luminous Blue Variables of which six and two, respectively, are new discoveries. We used spectroscopic, photometric, and light curve information to better constrain the nature of the reported objects. We particularly noted the presence of B[e]SGs at metallicity environments as low as 0.14 Z. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>Spectra of objects classified as B[e]SGs (including the B[e]SG candidate NGC7793-1). (Left) The full spectra for all stars with small offsets for better illustration purposes. The most prominent emission features are indicated. (Right) The region around H<math display="inline"><semantics> <mi>α</mi> </semantics></math> is highlighted to emphasize the relative strength of the emission compared to the continuum.</p>
Full article ">Figure 2
<p>Similar to <a href="#galaxies-11-00079-f001" class="html-fig">Figure 1</a>, but for LBVc. We note the lack of forbidden emission lines.</p>
Full article ">Figure 3
<p>The region between the [O <span class="html-small-caps">i</span>] and H<math display="inline"><semantics> <mi>α</mi> </semantics></math> line, that showcases multiple Fe <span class="html-small-caps">ii</span> emission lines. We note the clear presence of [O <span class="html-small-caps">i</span>] <math display="inline"><semantics> <mi>λ</mi> </semantics></math>6300 line for the B[e]SGs (left and middle panels, with the exception of the candidate NGC7793-1, due to the problematic spectrum; see text for more) and its absence from the LBVc spectra (right panel).</p>
Full article ">Figure 4
<p>The region around the [Ca <span class="html-small-caps">ii</span>] emission doublet. Its presence is evident in some B[e]SGs (left and middle panels, including NGC7793-1 candidate source, that suffers from data reduction artifacts due to slit overlaps), while LBVc (right panel) do not typically exhibit these lines (except for NGC55-3).</p>
Full article ">Figure 5
<p>The light curves from the Pan-STARRS survey for the B[e]SGs WLM-1 and NGC247-1. Each panel (per filter) shows the difference of each epoch from the mean value (noted on the y-axis label). See text for more.</p>
Full article ">Figure 6
<p>Same as <a href="#galaxies-11-00079-f005" class="html-fig">Figure 5</a>, but for the candidate LBVs NGC3109-1 and NGC247-2.</p>
Full article ">Figure 7
<p>(Left) The mid-IR <span class="html-italic">WISE</span> CCD for B[e]SGs and LBVs, including sources from the MCs (after [<a href="#B7-galaxies-11-00079" class="html-bibr">7</a>]) and our sample (for 5 out of 11 sources with <span class="html-italic">WISE</span> data). In general, the separation also holds for the new sources, with the exception of NGC55-1 (see text for more). (Right) IR CMD combining near-IR <span class="html-italic">J</span>-band (available for only five of our sources) with <span class="html-italic">Spitzer</span> [3.6]. We notice that, in this case, the newly found sources are consistent with the positions of the MC sources.</p>
Full article ">Figure 8
<p>(Left) The optical (<span class="html-italic">Gaia</span>) CMD, plotting BP–RP vs. M<math display="inline"><semantics> <msub> <mrow/> <mi>G</mi> </msub> </semantics></math> magnitude. We included all our sample and the MC sources from [<a href="#B7-galaxies-11-00079" class="html-bibr">7</a>] (except for two sources without a complete dataset in both <span class="html-italic">Gaia</span> and <span class="html-italic">Spitzer</span> surveys). (Right) The mid-IR (<span class="html-italic">Spitzer</span>) CMD using the IR color [3.6]–[4.5] vs. M<math display="inline"><semantics> <msub> <mrow/> <mrow> <mo>[</mo> <mn>4.5</mn> <mo>]</mo> </mrow> </msub> </semantics></math>. In this case, there is a significant improvement in the separation between the two classes. The position of NGC7793-1 favors a B[e]SG nature (see text for more).</p>
Full article ">Figure 9
<p>Similar to <a href="#galaxies-11-00079-f008" class="html-fig">Figure 8</a> but plotting the IR color [3.6]–[4.5] vs. the optical M<math display="inline"><semantics> <msub> <mrow/> <mi>G</mi> </msub> </semantics></math> magnitude. Similar to the IR CMD we saw relatively good separation between the two classes, with LBVs being brighter in the optical and less dusty compared to the B[e]SGs.</p>
Full article ">Figure 10
<p>The cumulative distribution function of the B[e]SGs and LBVs (including candidates) from this work and the literature. We notice (for the first time) the presence of B[e]SGs in lower metallicity environments and the fact that the two populations are not totally different (see text for more).</p>
Full article ">
22 pages, 33967 KiB  
Article
On the Evolution of, and Hot Gas in, Wind-Blown Bubbles around Massive Stars - Wind Bubbles Are Not Energy-Conserving
by Vikram V. Dwarkadas
Galaxies 2023, 11(3), 78; https://doi.org/10.3390/galaxies11030078 - 19 Jun 2023
Cited by 4 | Viewed by 1590
Abstract
The structure and evolution of wind-blown bubbles (WBBs) around massive stars has primarily been investigated using an energy-conserving model of wind-blown bubbles. While this model is useful in explaining the general properties of the evolution, several problems remain, including inconsistencies between observed wind [...] Read more.
The structure and evolution of wind-blown bubbles (WBBs) around massive stars has primarily been investigated using an energy-conserving model of wind-blown bubbles. While this model is useful in explaining the general properties of the evolution, several problems remain, including inconsistencies between observed wind luminosities and those derived using this formulation. Major difficulties include the low X-ray temperature and X-ray luminosity, compared to the model. In this paper, we re-examine the evolution, dynamics, and kinematics of WBBs around massive stars, using published ionization gasdynamic simulations of wind-blown bubbles. We show that WBBs can cool efficiently due to the presence of various instabilities and turbulence within the bubble. The expansion of WBBs is more consistent with a momentum-conserving solution, rather than an energy-conserving solution. This compares well with the dynamics and kinematics of observed wind bubbles. Despite the cooling of the bubble, the shocked wind temperature is not reduced to the observed values. We argue that the X-ray emission arise mainly from clumps and filaments within the hot shocked wind region, with temperatures just above 106 K. The remainder of the plasma can contribute to a lesser extent. Full article
Show Figures

Figure 1

Figure 1
<p>Snapshots of the number density (in cm<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </semantics></math>) with time. From the calculation of the evolution of a wind bubble around a 40 <math display="inline"><semantics> <msub> <mi>M</mi> <mo>⊙</mo> </msub> </semantics></math> star (600 × 400 zones), computed using an expanding grid ([<a href="#B30-galaxies-11-00078" class="html-bibr">30</a>,<a href="#B31-galaxies-11-00078" class="html-bibr">31</a>]). Time increases from left to right and top to bottom. The time in years is listed at the top of each panel. The scale shows the log of the number density. The dense swept-up shell, the ionized HII region, and the shocked wind region are all marked in the figure. The top two panels depict the MS phase. In the bottom left one, the RSG wind can be seen expanding near the star. It does not go too far given its low velocity. The bottom right panel depicts the clumps and filaments in the shocked wind region during the W-R phase.</p>
Full article ">Figure 2
<p>The evolution of the outer shell radius over time (shown in brown) from the simulation described herein. Overplotted is a cyan curve with Radius ∝ Time<math display="inline"><semantics> <msup> <mrow/> <mrow> <mn>0.48</mn> </mrow> </msup> </semantics></math>. The good fit clearly shows that the radius is not consistent with that predicted by the Weaver et al. [<a href="#B10-galaxies-11-00078" class="html-bibr">10</a>] model, but is consistent with a momentum-conserving bubble.</p>
Full article ">Figure 3
<p>Snapshots showing the velocity vectors (in blue) plotted over the density contours (in dark red). The formation of vortices in the hot shocked wind region is clearly seen. Time increases from left to right and is listed in years at the top of each panel. The vector in the top right corner of each panel denotes a velocity of 500 km s<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math>. The panels all correspond to the MS phase; the last one is at the transition between MS and RSG phases.</p>
Full article ">Figure 4
<p>Snapshots showing the pressure within the wind-blown bubble in various phases. Note the variations in pressure within the shocked wind region, as opposed to the uniform pressure expected in the Weaver et al. [<a href="#B10-galaxies-11-00078" class="html-bibr">10</a>] model. Time (in years) increases from left to right and is listed at the top of each panel. The first panel corresponds to the MS phase, the second to the RSG phase, and the third to the W-R phase.</p>
Full article ">Figure 5
<p>Snapshots showing the temperature of the wind-blown bubble in various phases. Note the variations in temperature in the hot shocked wind, with low temperatures particularly at the sites of instabilities, and mixing of cool ionized material with the hot interior. Time (in years) increases from left to right and is listed at the top of each panel. The first panel corresponds to the MS phase, the second to the RSG phase, and the third to the W-R phase.</p>
Full article ">Figure 6
<p>Panels showing the ionization fraction within the wind-blown bubble over time. Time (in years) increases from left to right, and is listed at the top of each panel. The first panel corresponds to the MS phase, the second to the RSG phase, and the third to the W-R phase. The ionization fraction decreases in the RSG phase, then steadily increases in the W-R phase.</p>
Full article ">Figure 7
<p>Evolution of the X-ray spectra within the W-R phase of the bubble. Time increases from left to right. The absorption column due to the nebular material decreases with time. Adapted from Figure 2 in [<a href="#B30-galaxies-11-00078" class="html-bibr">30</a>].</p>
Full article ">
11 pages, 1186 KiB  
Article
Vertical Structure of the Milky Way Disk with Gaia DR3
by Katherine Vieira, Vladimir Korchagin, Giovanni Carraro and Artem Lutsenko
Galaxies 2023, 11(3), 77; https://doi.org/10.3390/galaxies11030077 - 16 Jun 2023
Cited by 4 | Viewed by 1988
Abstract
Using a complete sample of about 330,000 dwarf stars, well measured by Gaia DR3, limited to the galactic north and south solid angles |b|<75° and up to a vertical distance of 2 kpc, we analyze the vertical structure [...] Read more.
Using a complete sample of about 330,000 dwarf stars, well measured by Gaia DR3, limited to the galactic north and south solid angles |b|<75° and up to a vertical distance of 2 kpc, we analyze the vertical structure of the Milky Way stellar disks, based on projected tangential velocities. From selected subsamples dominated by their corresponding population, we obtain the thin and thick disk scale heights as hZ=279.76±12.49 pc and HZ=797.23±12.34 pc, respectively. Then from the simultaneous fitting of the sum of two populations over the whole sample, assuming these scale heights, we estimate the thick-to-thin disk number density ratio at the galactic plane to be ρT/ρt=0.750±0.049, which is consistent with a previous result by the authors: in the galactic plane there is a significant number of thick disk stars, possibly as many as thin disk ones, which also points to the existence of more thick disk stars than generally thought. The overall fit does not closely follow the data for |Z|>700 pc and points to the presence of more stars beyond the thin disk that cannot be accounted for by the two-disk model. Full article
Show Figures

Figure 1

Figure 1
<p><math display="inline"><semantics> <mrow> <mo>(</mo> <mi>U</mi> <mo>,</mo> <mi>V</mi> <mo>)</mo> </mrow> </semantics></math> velocities of the sample extracted from Gaia. The disk sample is limited to the stars within the circle, having <math display="inline"><semantics> <mrow> <msqrt> <mrow> <msup> <mi>U</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>V</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&lt;</mo> <mn>180</mn> </mrow> </semantics></math> km s<sup>−1</sup>.</p>
Full article ">Figure 2
<p>Color-magnitude diagram of the disk sample, corrected by extinction and reddening. Disk stars are shown in grey while disk dwarfs stars are in black.</p>
Full article ">Figure 3
<p>Completeness evaluation of the 339,876 disk dwarfs stars sample. The black line shows the incompleteness percentage estimate (<b>left</b> axis), while the grey dashed line counts the number of stars (<b>right</b> axis), both evaluated at 100-pc height <span class="html-italic">coins</span> centered at each <span class="html-italic">Z</span>.</p>
Full article ">Figure 4
<p>KDE plot of the vertical distribution at various <math display="inline"><semantics> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> </semantics></math> intervals. Each curve corresponds to a <math display="inline"><semantics> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> </semantics></math> bin of 20 km s<sup>−1</sup> width, centered at −290 to −50 km s<sup>−1</sup>. The thick grey lines are subsamples dominated by the thin disk and the thick black lines are dominated by the thick disk. The thin black lines are a mix of both and the thin dashed line is the fastest rotating sample that exhibits significant asymmetry. The filled light and dark grey plots correspond to subsamples as defined in the legend: thin (<math display="inline"><semantics> <mrow> <mo>−</mo> <mn>280</mn> <mo>&lt;</mo> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> <mo>&lt;</mo> <mo>−</mo> <mn>200</mn> </mrow> </semantics></math> km s<sup>−1</sup>) and thick disk (<math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> <mo>&gt;</mo> <mo>−</mo> <mn>140</mn> </mrow> </semantics></math> km s<sup>−1</sup>) dominated samples.</p>
Full article ">Figure 5
<p>Thick disk sample (1000 &lt; <math display="inline"><semantics> <mrow> <mo>|</mo> <mi>Z</mi> <mo>|</mo> <mo>&lt;</mo> <mn>2000</mn> </mrow> </semantics></math> pc and <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> <mo>&gt;</mo> <mo>−</mo> <mn>140</mn> </mrow> </semantics></math> km s<sup>−1</sup>) vertical number density fit. Obtained scale height is <math display="inline"><semantics> <mrow> <mn>797.23</mn> <mo>±</mo> <mn>12.34</mn> </mrow> </semantics></math> pc.</p>
Full article ">Figure 6
<p>Thin disk sample (<math display="inline"><semantics> <mrow> <mn>100</mn> <mo>&lt;</mo> <mo>|</mo> <mi>Z</mi> <mo>|</mo> <mo>&lt;</mo> <mn>200</mn> </mrow> </semantics></math> pc and <math display="inline"><semantics> <mrow> <mo>−</mo> <mn>280</mn> <mo>&lt;</mo> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> <mo>&lt;</mo> <mo>−</mo> <mn>200</mn> </mrow> </semantics></math> km s<sup>−1</sup>) vertical number density fit. Obtained scale height is <math display="inline"><semantics> <mrow> <mn>279.76</mn> <mo>±</mo> <mn>12.49</mn> </mrow> </semantics></math> pc.</p>
Full article ">Figure 7
<p>Disk sample vertical number density simultaneous fit of thin and thick disk using Equation (<a href="#FD1-galaxies-11-00077" class="html-disp-formula">1</a>). Black points are the computed densities at each truncated cone centered at <span class="html-italic">Z</span> and of 100 pc height. All dwarfs stars with <math display="inline"><semantics> <mrow> <mn>100</mn> <mo>&lt;</mo> <mo>|</mo> <mi>Z</mi> <mo>|</mo> <mo>&lt;</mo> <mn>2000</mn> </mrow> </semantics></math> pc and no cut in <math display="inline"><semantics> <msub> <mi>V</mi> <mi>ϕ</mi> </msub> </semantics></math> were considered.</p>
Full article ">
23 pages, 968 KiB  
Article
Dense Molecular Environments of B[e] Supergiants and Yellow Hypergiants
by Michaela Kraus, Michalis Kourniotis, María Laura Arias, Andrea F. Torres and Dieter H. Nickeler
Galaxies 2023, 11(3), 76; https://doi.org/10.3390/galaxies11030076 - 16 Jun 2023
Cited by 3 | Viewed by 1422
Abstract
Massive stars expel large amounts of mass during their late evolutionary phases. We aim to unveil the physical conditions within the warm molecular environments of B[e] supergiants (B[e]SGs) and yellow hypergiants (YHGs), which are known to be embedded in circumstellar shells and disks. [...] Read more.
Massive stars expel large amounts of mass during their late evolutionary phases. We aim to unveil the physical conditions within the warm molecular environments of B[e] supergiants (B[e]SGs) and yellow hypergiants (YHGs), which are known to be embedded in circumstellar shells and disks. We present K-band spectra of two B[e]SGs from the Large Magellanic Cloud and four Galactic YHGs. The CO band emission detected from the B[e]SGs LHA 120-S 12 and LHA 120-S 134 suggests that these stars are surrounded by stable rotating molecular rings. The spectra of the YHGs display a rather diverse appearance. The objects 6 Cas and V509 Cas lack any molecular features. The star [FMR2006] 15 displays blue-shifted CO bands in emission, which might be explained by a possible close to pole-on oriented bipolar outflow. In contrast, HD 179821 shows blue-shifted CO bands in absorption. While the star itself is too hot to form molecules in its outer atmosphere, we propose that it might have experienced a recent outburst. We speculate that we currently can only see the approaching part of the expelled matter because the star itself might still block the receding parts of a (possibly) expanding gas shell. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>HR diagram with evolutionary tracks for LMC (<math display="inline"><semantics> <mrow> <mi>Z</mi> <mo>=</mo> <mn>0.006</mn> </mrow> </semantics></math>, [<a href="#B76-galaxies-11-00076" class="html-bibr">76</a>], left panel) and solar (<math display="inline"><semantics> <mrow> <mi>Z</mi> <mo>=</mo> <mn>0.014</mn> </mrow> </semantics></math>, [<a href="#B42-galaxies-11-00076" class="html-bibr">42</a>], right panel) metallicity for models of rotating stars (<math display="inline"><semantics> <mrow> <mi>v</mi> <mo>/</mo> <msub> <mi>v</mi> <mi>crit</mi> </msub> <mo>=</mo> <mn>0.4</mn> </mrow> </semantics></math>) with initial masses from <math display="inline"><semantics> <mrow> <mn>20</mn> <mrow> <mo>−</mo> </mrow> <mn>40</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math>. Shown are the positions of LMC (blue symbols) and Galactic (red symbols) objects. Only B[e]SGs with hot circumstellar molecular CO gas are shown. These populate similar evolutionary tracks as the YHGs. The minimum and maximum temperature values (where known) of the YHGs are connected by dashed lines. YHGs with reported (at least once) CO band emission are shown with triangles. The stars of the current study are labeled.</p>
Full article ">Figure 2
<p>Best fitting model (red) to the normalized Phoenix spectra (black) of LHA 120-S 12 (<b>top</b>) and LHA 120-S 134 (<b>bottom</b>). For each star we display the entire fit to the total spectrum (top panels) and the zoom to the band heads (bottom panels). Emission from the Pfund series (blue dashed line), detected in the spectrum of LHA 120-S 134, is included in the total fit.</p>
Full article ">Figure 3
<p>Normalized medium-resolution <span class="html-italic">K</span>-band spectra of the YHGs taken with GNIRS. For better visualization, the spectra have been offset along the flux axis. Positions of the CO band heads and of the lines from the Pfund series are marked by ticks. The lines of Br<math display="inline"><semantics> <mi>γ</mi> </semantics></math> and of the Na <span class="html-small-caps">i</span> doublet are labeled as well.</p>
Full article ">Figure 4
<p>Best fitting model (red) to the observed (black) <span class="html-italic">K</span>-band spectrum of [FMR2006] 15 for the model including a rotational component (<b>top</b>) and a pure Gaussian broadening (<b>bottom</b>).</p>
Full article ">Figure 5
<p>Residual spectrum of [FMR2006] 15 after subtraction of the blue-shifted CO band emission. For illustration purposes, a synthetic model spectrum of a cool supergiant is shown as well (shifted down along the flux axis for better visualization) with parameters similar to those determined for [FMR2006] 15.</p>
Full article ">Figure 6
<p>Comparison of the <span class="html-italic">K</span>-band spectrum of HD 179821 (<b>top</b>) with synthetic spectra for effective temperatures of 5400 K (<b>middle</b>) and 6600 K (<b>bottom</b>). For illustration purposes, the synthetic model spectra are included in this figure and shifted down along the flux axis for better visualization.</p>
Full article ">
31 pages, 2008 KiB  
Review
A Review of the Mixing Length Theory of Convection in 1D Stellar Modeling
by Meridith Joyce and Jamie Tayar
Galaxies 2023, 11(3), 75; https://doi.org/10.3390/galaxies11030075 - 16 Jun 2023
Cited by 36 | Viewed by 3070
Abstract
We review the application of the one-dimensional Mixing Length Theory (MLT) model of convection in stellar interiors and low-mass stellar evolution. We summarize the history of MLT, present a derivation of MLT in the context of 1D stellar structure equations, and discuss the [...] Read more.
We review the application of the one-dimensional Mixing Length Theory (MLT) model of convection in stellar interiors and low-mass stellar evolution. We summarize the history of MLT, present a derivation of MLT in the context of 1D stellar structure equations, and discuss the physical regimes in which MLT is relevant. We review attempts to improve and extend the formalism, including to higher dimensions. We discuss the interactions of MLT with other modeling physics, and demonstrate the impact of introducing variations in the convective mixing length, αMLT, on stellar tracks and isochrones. We summarize the process of performing a solar calibration of αMLT and state-of-the-art on calibrations to non-solar targets. We discuss the scientific implications of changing the mixing length, using recent analyses for demonstration. We review the most prominent successes of MLT, and the remaining challenges, and we conclude by speculating on the future of this treatment of convection. Full article
(This article belongs to the Special Issue The Structure and Evolution of Stars)
Show Figures

Figure 1

Figure 1
<p>Classification of stars according to convective structure in the main sequence. While all regions of the stellar model are treated using the MLT formalism, it is only fully convective stars and stars with convective envelopes (including those with both convective envelopes and convective cores) for which the <span class="html-italic">choice</span> of <math display="inline"><semantics> <msub> <mi>α</mi> <mi>MLT</mi> </msub> </semantics></math> is relevant. This is because the temperature gradient is approximately adiabatic in, e.g., convective cores (see <a href="#sec6-galaxies-11-00075" class="html-sec">Section 6</a> for more detail).</p>
Full article ">Figure 2
<p>The distance a “convective blob” can travel is measured in multiples of the pressure scale height, <math display="inline"><semantics> <msub> <mi>H</mi> <mi>P</mi> </msub> </semantics></math>. The upward motion of the parcel is driven by the thermal excess of the fluid parcel, compared to its surroundings. A larger mixing length implies that the parcel travels over a larger pressure differential before denaturing, which corresponds to more efficient transport of the flux by convection.</p>
Full article ">Figure 3
<p>Stellar tracks computed using the Dartmouth Stellar Evolution Program (DSEP) for a range of masses and mixing lengths.</p>
Full article ">Figure 4
<p>MESA: all one solar mass, range of MLT values. Other assumptions: photosphere tables for atmospheric boundary conditions and the Asplund et al. [<a href="#B113-galaxies-11-00075" class="html-bibr">113</a>] solar abundance scale. The pre-ZAMS evolutionary tracks are printed in fainter colors to avoid crowding and overlap with the ZAMS-RGB tracks. Note that in the case of <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mi>MLT</mi> </msub> <mo>=</mo> <mn>0.40</mn> </mrow> </semantics></math>, the model failed before reaching the ZAMS.</p>
Full article ">Figure 5
<p>A set of isochrones generated using MESA. All have identical compositions and ages, but the mixing length is varied. The colors of the curves correspond to the mixing length value assigned (uniformly) in the constituent tracks.</p>
Full article ">Figure 6
<p>(<b>Left</b>): Solar-normalized <math display="inline"><semantics> <msub> <mi>α</mi> <mi>MLT</mi> </msub> </semantics></math> fits to Alpha Cen A and B. This panel appears as Figure 3 in Joyce and Chaboyer [<a href="#B86-galaxies-11-00075" class="html-bibr">86</a>], and is reproduced with permission. (<b>Right</b>): Optimal mixing length as a function of mass for Alpha Cen A, B, and the Sun.</p>
Full article ">Figure 7
<p>(<b>Left</b>): Isochrones produced using the Dartmouth Stellar Evolution Program (DSEP) assume a range of mixing lengths. HST data for the metal-poor ([Fe/H]∼−2.4) globular cluster M92 are shown in gray. This figure appears as Figure 5 from Joyce and Chaboyer [<a href="#B85-galaxies-11-00075" class="html-bibr">85</a>], and is reproduced with permission. (<b>Right</b>): Stellar tracks from DSEP designed to fit tight constraints for the metal-poor subdwarf HD 140283 are computed with a range of assumptions about composition and <math display="inline"><semantics> <msub> <mi>α</mi> <mi>MLT</mi> </msub> </semantics></math>, demonstrating that variations in one can mimic the other. This figure appears as Figure 2 in Joyce and Chaboyer [<a href="#B85-galaxies-11-00075" class="html-bibr">85</a>], and is reproduced with permission.</p>
Full article ">Figure 8
<p>Custom isochrones based on MESA stellar tracks assume a range of mixing length values. In the upper right and lower two panels, spectroscopic parameters for micro-lensed subdwarfs determined by Bensby et al. [<a href="#B123-galaxies-11-00075" class="html-bibr">123</a>] are shown as black points, with their <math display="inline"><semantics> <mrow> <mn>1</mn> <mi>σ</mi> </mrow> </semantics></math> uncertainties shown as horizontal and vertical black lines. This figure appears as Figure 11 in Joyce et al. [<a href="#B4-galaxies-11-00075" class="html-bibr">4</a>], and is reproduced with permission.</p>
Full article ">Figure 9
<p>(<b>Left</b>): The effective mixing length required to match observations in a set of 1D stellar evolution models (black points) does not agree with the predictions of 3D simulations (gray bands) as a function of metallicity. (<b>Right</b>): The effective mixing length required to match observations in a set of 1D stellar evolution models (black points) does not agree with the predictions of 3D simulations (teal and red points) as a function of surface gravity. Reproduced with permission from Magic et al. [<a href="#B74-galaxies-11-00075" class="html-bibr">74</a>] and Tayar et al. [<a href="#B84-galaxies-11-00075" class="html-bibr">84</a>].</p>
Full article ">
43 pages, 7503 KiB  
Review
The Dynamics and Energetics of Remnant and Restarting RLAGN
by Vijay H. Mahatma
Galaxies 2023, 11(3), 74; https://doi.org/10.3390/galaxies11030074 - 13 Jun 2023
Cited by 3 | Viewed by 1606
Abstract
In this article, I review past, current, and future advances on the study of radio-loud AGN (RLAGN; radio-loud quasars and radio galaxies) lifecycles exclusively in the remnant and restarting phases. I focus on their dynamics and energetics as inferred from radio observations while [...] Read more.
In this article, I review past, current, and future advances on the study of radio-loud AGN (RLAGN; radio-loud quasars and radio galaxies) lifecycles exclusively in the remnant and restarting phases. I focus on their dynamics and energetics as inferred from radio observations while discussing their radiative lifetimes, population statistics, and trends in their physical characteristics. I briefly summarise multi-wavelength observations, particularly X-rays, that have enabled studies of the large-scale environments of RLAGN in order to understand their role in feedback. Furthermore, I discuss analytic and numerical simulations that predict key properties of remnant and restarting sources as found in wide-area surveys, and discuss the prospects of future surveys that may shed further light on these elusive subpopulations of RLAGN. Full article
(This article belongs to the Special Issue New Perspectives on Radio Galaxy Dynamics)
Show Figures

Figure 1

Figure 1
<p>1.4–5 GHz spectral index map of 3C 388, presented by [<a href="#B57-galaxies-11-00074" class="html-bibr">57</a>]. The sharp transition from flat (~<math display="inline"><semantics> <mrow> <mn>0.7</mn> </mrow> </semantics></math>) to steep (~<math display="inline"><semantics> <mrow> <mn>1.5</mn> </mrow> </semantics></math>) spectral indices can be seen between the regions coloured in blue and orange. © AAS. Reproduced with permission.</p>
Full article ">Figure 2
<p>Radio maps of the restarting source B1144 + 352 at 1.4 GHz from [<a href="#B75-galaxies-11-00074" class="html-bibr">75</a>]. The inset shows a high resolution image of the central GPS source, showing evidence of an edge-brightened restarting jet. Reproduced with permission © ESO.</p>
Full article ">Figure 3
<p>The cluster-centre remnant source of WNB1829 + 6911. <b>Left</b>: X-ray ROSAT All Sky Survey image with WENSS 327 MHz radio contours overlaid. <b>Right</b>: Digitized Sky Survey optical image with VLA 1.4 GHz radio contours overlaid, showing the outer steep-spectrum lobe material and the compact inner flat-spectrum source at the location of the host, with clear evidence of restarting activity. Image credit: [<a href="#B144-galaxies-11-00074" class="html-bibr">144</a>]. Reproduced with permission © ESO.</p>
Full article ">Figure 4
<p>The well-studied remnant source B2 0924 + 30. <b>Top</b>: 144 MHz total intensity map observed in the LoTSS survey at 6 arcsec resolution, from the publicly available archive of DR2, accessed in January 2023. <a href="https://lofar-surveys.org/dr2_release.html" target="_blank">https://lofar-surveys.org/dr2_release.html</a>, accessed on 25 May 2023. <b>Bottom</b>: radiative age map by [<a href="#B95-galaxies-11-00074" class="html-bibr">95</a>] at 1 arcmin angular resolution.</p>
Full article ">Figure 5
<p><b>Left panel</b>: distribution of lobe spectral indices (<math display="inline"><semantics> <mrow> <mi>S</mi> <mo>∝</mo> <msup> <mi>ν</mi> <mi>α</mi> </msup> </mrow> </semantics></math> here, not <math display="inline"><semantics> <mrow> <mi>S</mi> <mo>∝</mo> <msup> <mi>ν</mi> <mrow> <mo>−</mo> <mi>α</mi> </mrow> </msup> </mrow> </semantics></math> as for the right panel) between 150 MHz (LOFAR) and 1400 MHz (NVSS) for the candidate remnant sample selected by [<a href="#B162-galaxies-11-00074" class="html-bibr">162</a>]. Blue: candidate remnants with no high frequency core [<a href="#B163-galaxies-11-00074" class="html-bibr">163</a>]. Red: candidate remnants with a high frequency core that was not detected by FIRST. Black: initial flux-complete sample with FIRST cores (see text). <b>Right panel</b>: distribution of lobe spectral indices between 150 MHz and 1400 MHz for the mock radio galaxy catalogue using radiative and dynamical evolution for active (green), remnant (red), and remnant with ultra-steep spectra (blue) from [<a href="#B167-galaxies-11-00074" class="html-bibr">167</a>].</p>
Full article ">Figure 6
<p>WISE colour–colour plots for the host galaxies of remnants and restarting sources (<b>top</b> [<a href="#B118-galaxies-11-00074" class="html-bibr">118</a>]; <b>bottom</b> [<a href="#B120-galaxies-11-00074" class="html-bibr">120</a>]) and DDRGs (<b>middle</b> [<a href="#B117-galaxies-11-00074" class="html-bibr">117</a>]) compared with those for the parent samples of RLAGN. Both studies that select samples heterogeneously from the same survey but with different selection methods (top and middle) replicate the conclusion that interrupted radio jet activity is not a consequence of being hosted by different types of galaxies, a conclusion based on optical and infrared photometric properties.</p>
Full article ">Figure 7
<p>The source ILTJ131403.17 + 543939.6 from the DDRG sample of [<a href="#B117-galaxies-11-00074" class="html-bibr">117</a>], showing evidence of a third episode of activity based on the diffuse extended plasma beyond the outer double. Greyscale image: 6 GHz VLA data showing the inner restarting double. Orange contours: 144 MHz LOFAR data.</p>
Full article ">Figure 8
<p><b>Top</b>: 1.5 GHz radio image (red contours) of 3C 293 superposed on a near-infrared 2 <math display="inline"><semantics> <mi mathvariant="sans-serif">μ</mi> </semantics></math>m image (from [<a href="#B192-galaxies-11-00074" class="html-bibr">192</a>]). The inset (a) shows 5 GHz MERLIN contours (red), with the inner restarting jet and CO distribution showing the molecular gas (blue) overlaid on a <span class="html-italic">HST</span> image with large dust lanes (© AAS. Reproduced with permission). <b>Bottom</b>: Spectral index map between 144 MHz and 1360 MHz at 0.28 × 0.23 arcsec resolution (from [<a href="#B106-galaxies-11-00074" class="html-bibr">106</a>]), clearly showing the distinct diffuse older electron population beyond the flat spectrum inner lobes.</p>
Full article ">Figure 9
<p>Composite optical (background), X-ray (blue), and radio (red) image of the galaxy group NGC 507 (from [<a href="#B195-galaxies-11-00074" class="html-bibr">195</a>]). Remnant emission newly detected by LOFAR, interpreted as remnant emission from a previous AGN episode, is suggested to be transported by gas sloshing induced by merger events. Reproduced with permission © ESO.</p>
Full article ">Figure 10
<p>Radio source size function from the dynamical models of intermittent radio sources by [<a href="#B243-galaxies-11-00074" class="html-bibr">243</a>] (dotted line and binned as black filled squares). Open squares show the data from the CSS/GPS sample by [<a href="#B56-galaxies-11-00074" class="html-bibr">56</a>], showing the plateau in sizes at ~1 kpc, approximately consistent with theoretical models of intermittent radio sources showing fine structure related to the number of intermittency cycles.</p>
Full article ">Figure 11
<p>Time-stamps (<b>left</b> to <b>right</b>, <b>top</b> to <b>bottom</b>) of the evolution of intermittent RLAGN plasma (taken from [<a href="#B249-galaxies-11-00074" class="html-bibr">249</a>]). The colour-coding is the fraction of material that originated in the jet relative to the environment (dark red = jet plasma, white = ICM). Vectors denote gas velocities. The jet is off in the top left panels, restarting in the top right, turned off in the second row right panel, and remaining as a remnant for the rest. It can be seen that vortices remain intact for the entire remnant phase, and mixing between the jet material and the ICM happens slowly, smoothing out the sporadic nature of intermittency. © AAS, reproduced with permission.</p>
Full article ">Figure 12
<p>Anticipated sensitivity of the SKA1 and SKA2 telescopes as a function of frequency, compared with existing telescopes. Figure taken from <a href="https://www.skao.int/en/science-users/" target="_blank">https://www.skao.int/en/science-users/</a>, accessed on 1 February 2023. Credit: SKA Observatory.</p>
Full article ">Figure 13
<p>Modelled spectral energy distribution of a typical remnant at <math display="inline"><semantics> <mrow> <mi>z</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> and with a jet power of <math display="inline"><semantics> <mrow> <msub> <mi>L</mi> <mi>j</mi> </msub> <mo>=</mo> <msup> <mn>10</mn> <mn>45</mn> </msup> </mrow> </semantics></math> erg s<sup>−1</sup>, taken from [<a href="#B239-galaxies-11-00074" class="html-bibr">239</a>]. They predict that while the total flux from the remnant lobes (thick black line) may not be visible at GHz frequencies with existing telescopes, the bow shock (thick dark red line) may be marginally visible with both SKA1−LOW and SKA1−MID for an integration time of 10 h.</p>
Full article ">
46 pages, 5709 KiB  
Review
Recent Progress in Modeling the Macro- and Micro-Physics of Radio Jet Feedback in Galaxy Clusters
by Martin A. Bourne and Hsiang-Yi Karen Yang
Galaxies 2023, 11(3), 73; https://doi.org/10.3390/galaxies11030073 - 13 Jun 2023
Cited by 15 | Viewed by 2048
Abstract
Radio jets and the lobes they inflate are common in cool-core clusters and are known to play a critical role in regulating the heating and cooling of the intracluster medium (ICM). This is an inherently multi-scale problem, and much effort has been made [...] Read more.
Radio jets and the lobes they inflate are common in cool-core clusters and are known to play a critical role in regulating the heating and cooling of the intracluster medium (ICM). This is an inherently multi-scale problem, and much effort has been made to understand the processes governing the inflation of lobes and their impact on the cluster, as well as the impact of the environment on the jet–ICM interaction, on both macro- and microphysical scales. The developments of new numerical techniques and improving computational resources have seen simulations of jet feedback in galaxy clusters become ever more sophisticated. This ranges from modeling ICM plasma physics processes such as the effects of magnetic fields, cosmic rays, and viscosity to including jet feedback in cosmologically evolved cluster environments in which the ICM thermal and dynamic properties are shaped by large-scale structure formation. In this review, we discuss the progress made over the last ∼decade in capturing both the macro- and microphysical processes in numerical simulations, highlighting both the current state of the field, as well as the open questions and potential ways in which these questions can be addressed in the future. Full article
(This article belongs to the Special Issue New Perspectives on Radio Galaxy Dynamics)
Show Figures

Figure 1

Figure 1
<p>Diagram illustrating the general processes that occur during and after lobe inflation. (<b>A</b>): A fast jet drives into the ambient medium, forms a bow shock, and inflates a hot lobe that expands into the ICM. The lobe morphology can depend sensitively on the injected jet properties (e.g., content, velocity, geometry). The expanding shock wave results in a layer of shocked ICM material surrounding the jet lobe. (<b>B</b>): As the lobe expansion slows, which may or may not be accompanied by the jet switching off, the shock driven into the ICM broadens into a sound wave that can detach from the lobes. (<b>C</b>): Once the jet has ceased and as the lobe buoyantly rises through the ICM, dense, low entropy material can be entrained and pulled up in the wake. Moreover, instabilities can lead to mixing of the lobe and ICM material. Sound waves generated by the lobe expansion can continue to propagate to large distances depending on the ICM viscosity. (<b>D</b>): This process continues at late times, with mixing continuing to dilute the lobe material, although the rate at which this occurs can sensitively depend on the ICM physical processes including magnetic fields, viscosity, and cluster weather (i.e., the ICM and sub-structure motions, see <a href="#sec2dot2-galaxies-11-00073" class="html-sec">Section 2.2</a>).</p>
Full article ">Figure 2
<p>Thin temperature projections illustrate how jet injection parameters impact jet and lobe morphologies. All jets are 100 Myr old, with each row illustrating the effect of changing one parameter: from top to bottom, these are jet power, half opening angle, velocity, and resolution, respectively. These quantities take fiducial values of <math display="inline"><semantics> <msup> <mn>10</mn> <mn>46</mn> </msup> </semantics></math> erg s<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math>, <math display="inline"><semantics> <mrow> <mn>10</mn> <mo>°</mo> </mrow> </semantics></math>, 15,000 km s<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math>, and <math display="inline"><semantics> <mrow> <mn>1.81</mn> <mo>×</mo> <msup> <mn>10</mn> <mn>5</mn> </msup> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math>, unless being varied. Color bars shown to the right of each row extend down to <math display="inline"><semantics> <msup> <mn>10</mn> <mn>7</mn> </msup> </semantics></math> K, and the gas below this is shown in black. (Figure 6 from Huško and Lacey [<a href="#B71-galaxies-11-00073" class="html-bibr">71</a>], CC BY 4.0.)</p>
Full article ">Figure 3
<p>An illustration of the jet-driven shocks within the central 100 kpc of a simulated cluster. The simulation includes self-regulated AGN feedback that assumes the jet power is coupled to the SMBH accretion of cold gas with an efficiency parameter of <math display="inline"><semantics> <mrow> <mi>ϵ</mi> <mo>=</mo> <mn>1</mn> <mo>%</mo> </mrow> </semantics></math> (as defined in Equation (<a href="#FD1-galaxies-11-00073" class="html-disp-formula">1</a>)). The jets are assumed to have a small angle precession (<math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>0.15</mn> </mrow> </semantics></math>), with a 10 Myr period around the <span class="html-italic">z</span>-axis in the simulation. The left panel shows the energy dissipation rate, while the right-hand panel shows shock Mach numbers, illustrating that stronger shocks and hence higher dissipation rates are seen closer to the jet; however, overall the shocks are typically quite weak. (Figure 1 from Li et al. [<a href="#B96-galaxies-11-00073" class="html-bibr">96</a>], ©AAS, reproduced with permission).</p>
Full article ">Figure 4
<p>Volume rendering of jet lobes (green) and cold material (red) within a cosmologically evolved galaxy cluster. The top, middle, and bottom rows show low, medium (labeled “Fid”), and high-power jets, respectively. The small panels show the evolution of the jet lobes for two different viewing angles (rotated by <math display="inline"><semantics> <mrow> <mn>90</mn> <mo>°</mo> </mrow> </semantics></math> about the <span class="html-italic">z</span>-axis with respect to each other), while the large panels also overlay the gas velocity field. Overall, the jet lobes can be displaced, disrupted, and mixed by cluster weather and cold structures, with lower-power jets being more susceptible. (Figure 2 from Bourne and Sijacki [<a href="#B228-galaxies-11-00073" class="html-bibr">228</a>], CC BY 4.0.)</p>
Full article ">Figure 5
<p>Comparisons of the simulation results with varied jet composition and assumptions for modeling CR transport. The simulations model a single outburst of jet activity (jet power of <math display="inline"><semantics> <mrow> <mn>5</mn> <mspace width="3.33333pt"/> <mo>×</mo> <mspace width="3.33333pt"/> <msup> <mn>10</mn> <mn>45</mn> </msup> </mrow> </semantics></math> erg s<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math> and jet duration of 10 Myr) in a Perseus-like hydrostatic atmosphere. Rows from top to bottom show the results from kinetic-energy-dominated jets (KIN), CR-dominated jets (CR), and CR-dominated jets with diffusion and heating (CRdh). The morphology of jet-inflated lobes tend to be elongated for the KIN case, whereas bubbles inflated by CR-dominated jets are wider. Comparisons between the bottom two panels show that CR heating is more efficient and the amount of cold gas is less for the CRdh case as the CRs can diffuse outside the bubbles and heat the ICM due to hadronic and Coulomb interactions. (Figure 3 from Yang et al. [<a href="#B90-galaxies-11-00073" class="html-bibr">90</a>], ©AAS, reproduced with permission).</p>
Full article ">Figure 6
<p>Impact of assumptions about ICM viscosity on the evolution of AGN-jet-inflated bubbles. Cases (<b>A</b>–<b>D</b>) show the simulations with no viscosity, isotropic viscosity with full Spitzer values, anisotropic viscosity with full Braginskii values, and anisotropic viscosity limited by microinstabilities, respectively. While the hydrodynamic instabilities are suppressed by viscosity in cases (<b>B</b>,<b>C</b>), when the parallel viscosity along the magnetic field lines is suppressed by microinstabilities (case (<b>D</b>)), the viscosity is strongly limited and the bubbles are deformed as in the inviscid case (<b>A</b>). This illustrates the importance of modeling the ICM microphysics in AGN simulations. (Figure 1 from Kingsland et al. [<a href="#B117-galaxies-11-00073" class="html-bibr">117</a>], ©AAS, reproduced with permission).</p>
Full article ">Figure 7
<p>Illustration of the sound waves generated by jets. The left-hand panel shows the jet and cocoon temperature and density structure, with key features labeled. The right-hand panel shows jet entropy and the acoustic flux density, with the structure in the latter illustrating the production of sound waves within the shocked ICM material. (Figure 1 from Bambic and Reynolds [<a href="#B239-galaxies-11-00073" class="html-bibr">239</a>], ©AAS, reproduced with permission).</p>
Full article ">
20 pages, 936 KiB  
Article
New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy
by Andrea Fabiana Torres, María Laura Arias, Michaela Kraus, Lorena Verónica Mercanti and Tõnis Eenmäe
Galaxies 2023, 11(3), 72; https://doi.org/10.3390/galaxies11030072 - 10 Jun 2023
Viewed by 1427
Abstract
The B[e] phenomenon is manifested by a heterogeneous group of stars surrounded by gaseous and dusty circumstellar envelopes with similar physical conditions. Among these stars, the FS CMa-type objects are suspected to be binary systems, which could be experiencing or have undergone a [...] Read more.
The B[e] phenomenon is manifested by a heterogeneous group of stars surrounded by gaseous and dusty circumstellar envelopes with similar physical conditions. Among these stars, the FS CMa-type objects are suspected to be binary systems, which could be experiencing or have undergone a mass-transfer process that could explain the large amount of material surrounding them. We aim to contribute to the knowledge of a recently confirmed binary, MWC 645, which could be undergoing an active mass-transfer process. We present near-infrared and optical spectra, identify atomic and molecular spectral features, and derive different quantitative properties of line profiles. Based on publicly available photometric data, we search for periodicity in the light curve and model the spectral energy distribution. We have detected molecular bands of CO in absorption at 1.62 μm and 2.3 μm for the first time. We derive an upper limit for the effective temperature of the cool binary component. We found a correlation between the enhancement of the Hα emission and the decrease in optical brightness that could be associated with mass-ejection events or an increase in mass loss. We outline the global properties of the envelope, possibly responsible for brightness variations due to a variable extinction, and briefly speculate on different possible scenarios. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS on July 2018 from 8500 Å to 13,600 Å. The normalized Ondřejov spectrum from 8400 Å to 8870 Å, acquired on September 2018, is shown in cyan. Main spectral lines are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, calcium in pink, nitrogen in violet, sulfur in green, and helium in cyan. Wavelengths are given in angstroms.</p>
Full article ">Figure 2
<p>Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS in 2018, covering the <span class="html-italic">H</span>-(upper panel) and <span class="html-italic">K</span>-bands (lower panel). Main spectral lines and molecular bands are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) or molecule (of different isotopes) are joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, and carbon monoxide in pink. Wavelengths are given in angstroms.</p>
Full article ">Figure 3
<p>Normalized <span class="html-italic">L</span>-band spectrum of MWC 645 obtained in 2022. The emission lines of H I and Fe II are marked in red and blue, respectively. The “bump” longward of the Br<math display="inline"><semantics> <mi>α</mi> </semantics></math> line is a remnant from telluric correction. Wavelengths are in angstroms.</p>
Full article ">Figure 4
<p>Strongest H I lines detected in our IR spectra of MWC 645. They display an asymmetric profile, where the red flank is steeper than the blue one. Wavelengths are in angstroms.</p>
Full article ">Figure 5
<p>CO molecular bands of MWC 645. Some emission and absorption lines are also identified. The spectral features of a given element or molecule (of different isotopes) are indicated by colored markings and joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, calcium in cyan, and carbon monoxide in pink. Upper panel: <math display="inline"><semantics> <msup> <mrow/> <mn>12</mn> </msup> </semantics></math>CO second-overtone band heads seen in the <span class="html-italic">H</span>-band spectrum taken in 2018. Lower panel: <math display="inline"><semantics> <msup> <mrow/> <mn>12</mn> </msup> </semantics></math>CO and <math display="inline"><semantics> <msup> <mrow/> <mn>13</mn> </msup> </semantics></math>CO band heads in absorption detected in the <span class="html-italic">K</span>-band. The spectra obtained in 2017 (in red) and 2018 (in black) revealed the variability in the strength of the observed bands. Wavelengths are given in angstroms.</p>
Full article ">Figure 6
<p>Comparison between the degraded spectrum of MWC 645 to R ∼ 2000 (solid black line) and a G0 Ib-II star, HD 185018 (dashed red line), where the CO(2-0) band head fits well. The other <math display="inline"><semantics> <msup> <mrow/> <mn>12</mn> </msup> </semantics></math>CO absorption bands from MWC 645 are shallower than those of the template; perhaps they are filled by emission. Wavelengths are given in angstroms.</p>
Full article ">Figure 7
<p>Light curve of MWC 645 from 16 December 2014 up to 11 November 2022 taken from ASAS-SN. The purple squares indicate the <span class="html-italic">V</span>-band magnitudes, and the green circles represent the <span class="html-italic">g</span>-band measurements converted to the <span class="html-italic">V</span>-band magnitudes. The conversion has been carried out with the relationship found by Nodyarov et al. [<a href="#B20-galaxies-11-00072" class="html-bibr">20</a>]. Their optical photometric observations are also included (blue circles). Vertical solid blue lines mark the dates of our IR observations (2017, 2018, and 2022, respectively); dotted blue lines mark the dates of our optical spectra (2018 and 2021, respectively); and those dotted in green, gray, and red correspond to the spectra downloaded from the BeSS database taken in 2019, 2020, and 2022, respectively. Time is given in heliocentric Julian dates (HJD) minus 2.45 × 10<math display="inline"><semantics> <msup> <mrow/> <mn>6</mn> </msup> </semantics></math> days.</p>
Full article ">Figure 8
<p>Fourier power spectrum of the ASAS-SN light curve of MWC 645. The periods in days are on a logarithmic scale. The green line shows the confidence level.</p>
Full article ">Figure 9
<p>H<math display="inline"><semantics> <mi>α</mi> </semantics></math> line variation of MWC 645. Spectra taken in 2018, 2019, 2020, and 2022 are displayed in blue, green, gray, and red, respectively. The spectrum with the lowest resolution is plotted with a dashed-line. The heliocentric radial velocity scale is shown in km s<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math>.</p>
Full article ">Figure 10
<p>Example of the emission lines in the surroundings of the H<math display="inline"><semantics> <mi>α</mi> </semantics></math> line of MWC 645. Ondřejov normalized spectra taken in 2018 on 12 September (in red line) and 13 September (in black line) and the 2019 low-resolution spectrum (in green) are shown with the main lines identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, nitrogen in violet, and sulfur in green. Wavelengths are given in angstroms.</p>
Full article ">Figure 11
<p>Spectral energy distribution of MWC 645. The open triangles represent the observed photometric data: optical bands (yellow), 2MASS (green) [<a href="#B42-galaxies-11-00072" class="html-bibr">42</a>], WISE (violet) [<a href="#B43-galaxies-11-00072" class="html-bibr">43</a>], MSX (light blue) [<a href="#B44-galaxies-11-00072" class="html-bibr">44</a>], IRAS (red) [<a href="#B45-galaxies-11-00072" class="html-bibr">45</a>], and AKARI (black) [<a href="#B46-galaxies-11-00072" class="html-bibr">46</a>]. The error bars of the photometric data are included (in most photometric bands, they fall inside the symbols). The low-resolution spectrum acquired in 2021 over the H<math display="inline"><semantics> <mi>α</mi> </semantics></math> region is also displayed (with a dashed black line). The solid red line shows the SED modeled considering the contribution of the photospheric fluxes from both stars, the thermal emission from a gaseous shell close to the system and the effect of the interstellar medium extinction. The solid blue line shows our best-fitting theoretical SED, obtained by adding to the SED plotted in red the contribution of three dusty shells surrounding the stellar system. The flux is normalized to that at <math display="inline"><semantics> <msub> <mi>λ</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> </semantics></math> = 0.55 <math display="inline"><semantics> <mi mathvariant="sans-serif">μ</mi> </semantics></math>m and displayed on a logarithmic scale. The wavelengths are in microns.</p>
Full article ">
12 pages, 338 KiB  
Article
Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope
by Lindita Hamolli, Mimoza Hafizi, Francesco De Paolis and Esmeralda Guliqani
Galaxies 2023, 11(3), 71; https://doi.org/10.3390/galaxies11030071 - 1 Jun 2023
Viewed by 1399
Abstract
In this work, we investigate the possibility of observing quasars, particularly lensed quasars, by the Nancy Grace Roman Space Telescope (Roman). To this aim, based on the capabilities of the Roman Space Telescope and the results from the quasar luminosity function (QLF) in [...] Read more.
In this work, we investigate the possibility of observing quasars, particularly lensed quasars, by the Nancy Grace Roman Space Telescope (Roman). To this aim, based on the capabilities of the Roman Space Telescope and the results from the quasar luminosity function (QLF) in the infrared band of the Spitzer Space Telescope imaging survey, we calculated the number of quasars expected to be in its field of view. In order to estimate the number of lensed quasars, we develop a Monte Carlo simulation to estimate the probability that a quasar is lensed once or more times by foreground galaxies. Using the mass–luminosity distribution function of galaxies and the redshift distributions of galaxies and quasars, we find that 1 per 180 observed quasars will be lensed by foreground galaxies. Further on, adopting a singular isothermal sphere (SIS) model for lens galaxies, we calculate the time delay between lensed images for single and multiple lensing systems and present their distributions. We emphasize that detailed studies of these lensing systems will provide a powerful probe of the physical properties of quasars and may allow testing the mass distribution models of galaxies in addition to being extremely helpful for constraining the cosmological parameters. Full article
Show Figures

Figure 1

Figure 1
<p>Redshift distributions in single lensing events expected to be observed by Roman Space Telescope. The red line shows the redshift distribution of the galaxy and the gray shadow shows the redshift distribution of quasars.</p>
Full article ">Figure 2
<p>Redshift distributions in multiple lensing events expected to be observed by the Roman Space Telescope. The red line shows the redshift distribution of the first galaxy, the dashed black line shows the redshift distribution of the second galaxy and gray shadow shows the redshift distribution of quasars.</p>
Full article ">Figure 3
<p>Distribution of number of images for strong lensing events caused by multiple lenses. As one can see, up to six images can form (see text for details).</p>
Full article ">Figure 4
<p>Time delay distribution between two lensed images by the Einstein radius.</p>
Full article ">
13 pages, 627 KiB  
Article
Hawking Radiation and Lifetime of Primordial Black Holes in Braneworld
by Bobur Turimov, Akhror Mamadjanov and Ozodbek Rahimov
Galaxies 2023, 11(3), 70; https://doi.org/10.3390/galaxies11030070 - 31 May 2023
Cited by 3 | Viewed by 1639
Abstract
The paper explores the thermodynamic properties of primordial black holes (PBHs) in the braneworld. Specifically, the researchers examined Hawking radiation and the lifetime of PBHs. Through their analysis, an exact analytical expression for the Bekenstein–Hawking entropy, temperature, and heat capacity was derived. Their [...] Read more.
The paper explores the thermodynamic properties of primordial black holes (PBHs) in the braneworld. Specifically, the researchers examined Hawking radiation and the lifetime of PBHs. Through their analysis, an exact analytical expression for the Bekenstein–Hawking entropy, temperature, and heat capacity was derived. Their findings suggest that the lifetime of PBHs in the early universe is reduced by at least one order of magnitude, ultimately leading to their evaporation. This could explain why we have not observed the final rapid evaporation of PBHs in the recent epoch of the universe. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>Left panel</b>) Dependence of radius of horizon of PBH form its mass for the different values of brane tension. (<b>Right panel</b>) Dependence of radius of horizon of PBH form brane tension parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>(<b>Left panel</b>) Dependence of temperature of PBH from its mass for the different values of brane parameter. (<b>Right panel</b>) Dependence of temperature of PBH from brane parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>(<b>Left panel</b>) Fractional lifetime of PBH is a function of brane charge parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>. (<b>Right panel</b>) Fractional lifetime of PBH is a function of its mass for different values of brane charge parameter.</p>
Full article ">Figure 4
<p>(<b>Left panel</b>) Dependence of entropy of PBH (in Equation (<a href="#FD14-galaxies-11-00070" class="html-disp-formula">14</a>)) from its mass for the different values of brane tension. (<b>Right panel</b>) Dependence of entropy of PBH from brane tension parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Dependence of free energy from brane tension <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>(<b>Left panel</b>) Dependence of specific heat of PBH from its mass for the different values of brane tension. (<b>Right panel</b>) Dependence of specific heat of PBH from brane tension parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>(Color online) Radial dependence of the effective potential for the different values of brane tension parameter <math display="inline"><semantics> <mrow> <mi>b</mi> <mo>/</mo> <mi>M</mi> </mrow> </semantics></math> at <math display="inline"><semantics> <mrow> <mo>ℓ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>. Solid lines represent scalar field (<math display="inline"><semantics> <mrow> <mi>s</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>) while dashed lines vector field (<math display="inline"><semantics> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>) one.</p>
Full article ">Figure 8
<p>(<b>Left panel</b>) The transmission and reflection coefficients of PBH for the different values of brane tension at <math display="inline"><semantics> <mrow> <mo>ℓ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>. (<b>Right panel</b>) Dependence of transmission and reflection coefficients of PBH from the brane tension for the scalar and vector fields at <math display="inline"><semantics> <mrow> <mo>ℓ</mo> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">
40 pages, 90706 KiB  
Article
Synthetic Light Curve Design for Pulsating Binary Stars to Compare the Efficiency in the Detection of Periodicities
by Aldana Alberici Adam, Gunther F. Avila Marín, Alejandra Christen and Lydia Sonia Cidale
Galaxies 2023, 11(3), 69; https://doi.org/10.3390/galaxies11030069 - 31 May 2023
Cited by 1 | Viewed by 1275
Abstract
B supergiant stars pulsate in regular and quasi-regular oscillations resulting in intricate light variations that might conceal their binary nature. To discuss possible observational bias in a light curve, we performed a simulation design of a binary star affected by sinusoidal functions emulating [...] Read more.
B supergiant stars pulsate in regular and quasi-regular oscillations resulting in intricate light variations that might conceal their binary nature. To discuss possible observational bias in a light curve, we performed a simulation design of a binary star affected by sinusoidal functions emulating pulsation phenomena. The Period04 tool and the WaveletComp package of R were used for this purpose. Thirty-two models were analysed based on a combination of two values on each of the k = 6 variables, such as multiple pulsations, the amplitude of the pulsation, the pulsation frequency, the beating phenomenon, the light-time effect, and regular or quasi-regular periods. These synthetic models, unlike others, consider an ARMA (1, 1) statistical noise, irregular sampling, and a gap of about 4 days. Comparing Morlet wavelet with Fourier methods, we observed that the orbital period and its harmonics were well detected in most cases. Although the Fourier method provided more accurate period detection, the wavelet analysis found it more times. Periods seen with the wavelet method have a shift due to the slightly irregular time scale used. The pulsation period hitting rate depends on the wave amplitude and frequency with respect to eclipse depth and orbital period. None of the methods was able to distinguish accurate periods leading to a beating phenomenon when they were longer than the orbital period, resulting, in both cases, in an intermediate value. When the beating period was shorter, the Fourier analysis found it in all cases except for unsolved quasi-regular periods. Overall, the Morlet wavelet analysis performance was lower than the Fourier analysis. Considering the strengths and disadvantages found in these methods, we recommend using at least two diagnosis tools for a detailed time series data analysis to obtain confident results. Moreover, a fine-tuning of trial periods by applying phase diagrams would be helpful for recovering accurate values. The combined analysis could reduce observational bias in searching binaries using photometric techniques. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>Synthetic light curve of an eclipsing binary affected by ARMA (1, 1) noise.</p>
Full article ">Figure 2
<p>Example of applying wavelet analysis for a multi-periodic time series. Left panel: sinusoidal series with periods 30, 38, and 80. Central panel: Scalogram or wavelet power spectrum of the series. The top in the color code indicates a higher significance for periods. Periods with a higher power are detected with solid black lines. Right panel: average wavelet power of the series in logarithm scale. In this diagram, the beating periods 30 and 38 are more difficult to identify individually, unlike period 80, but they show a striking feature in the scalogram. Although all the periods represented in red have a significance level less than or equal to 0.05, the most relevant ones, about 30, 38, and 80, are the three highest average wavelet power. They were obtained using the WaveletComp Monte Carlo sampling-based test.</p>
Full article ">Figure 3
<p>Synthetic light curves obtained from the simulation design. The model number is placed according to the list in <a href="#galaxies-11-00069-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 4
<p>Results of the wavelet analysis for Model No. 0, 1, 5, 17, and 21. From left to right, the wavelet power spectrum and the average wavelet power for the second part of each synthetic light curve.</p>
Full article ">Figure 5
<p>Periodograms obtained for the models No. 0, 1, 5, 17, and 21.</p>
Full article ">Figure 6
<p>Comparison between the proportion of correctness of the original period obtained by wavelet (cyan) and Fourier (grey) analysis. The thick solid black line represents the median of the data, and the red symbol stands for an outlier or extreme value.</p>
Full article ">Figure 7
<p>Boxplots of hit ratio of the original periods with the wavelet (cyan) and Fourier (grey) analysis for each variable involved. The thick solid black line represents the median of the data, and the red symbol stands for an outlier or extreme value. The hit rate is on the <span class="html-italic">y</span>-axis and the two categories (indicated by <math display="inline"><semantics> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </semantics></math> and 1) for each variable are shown on <span class="html-italic">x</span>-axis: beating phenomenon (absence/presence), period of the pulsation (minor/major), amplitude of the pulsation (<math display="inline"><semantics> <mrow> <mn>4</mn> <mo>%</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mn>20</mn> <mo>%</mo> </mrow> </semantics></math>), number of pulsations periods (1/2), and quasi-regular pulsations (absence/presence).</p>
Full article ">Figure 8
<p>Light curves obtained by TESS of pulsating binary stars.</p>
Full article ">Figure A1
<p>Periodograms obtained from the Fourier analysis for models No. 0 to No. 23.</p>
Full article ">Figure A1 Cont.
<p>Periodograms obtained from the Fourier analysis for models No. 0 to No. 23.</p>
Full article ">Figure A2
<p>Periodograms obtained from the Fourier analysis for models No. 24 to No. 32.</p>
Full article ">Figure A3
<p>Results of the wavelet analysis for models No. 0 to No. 5. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">Figure A4
<p>Results of the wavelet analysis for models No. 6 to No. 11. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">Figure A5
<p>Results of the wavelet analysis for models No. 12 to No. 17. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">Figure A6
<p>Results of the wavelet analysis for models No. 18 to No. 23. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">Figure A7
<p>Results of the wavelet analysis for models No. 24 to No. 29. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">Figure A8
<p>Results of the wavelet analysis for models No. 30 to No. 32. From left to right the wavelet power spectrum and the average wavelet power for the first and second part of each synthetic curve.</p>
Full article ">
33 pages, 2909 KiB  
Review
Radiation-Driven Wind Hydrodynamics of Massive Stars: A Review
by Michel Curé and Ignacio Araya
Galaxies 2023, 11(3), 68; https://doi.org/10.3390/galaxies11030068 - 12 May 2023
Cited by 6 | Viewed by 2190
Abstract
Mass loss from massive stars plays a determining role in their evolution through the upper Hertzsprung–Russell diagram. The hydrodynamic theory that describes their steady-state winds is the line-driven wind theory (m-CAK). From this theory, the mass loss rate and the velocity profile of [...] Read more.
Mass loss from massive stars plays a determining role in their evolution through the upper Hertzsprung–Russell diagram. The hydrodynamic theory that describes their steady-state winds is the line-driven wind theory (m-CAK). From this theory, the mass loss rate and the velocity profile of the wind can be derived, and estimating these properly will have a profound impact on quantitative spectroscopy analyses from the spectra of these objects. Currently, the so-called β law, which is an approximation for the fast solution, is widely used instead of m-CAK hydrodynamics, and when the derived value is β1.2, there is no hydrodynamic justification for these values. This review focuses on (1) a detailed topological analysis of the equation of motion (EoM), (2) solving the EoM numerically for all three different (fast and two slow) wind solutions, (3) deriving analytical approximations for the velocity profile via the LambertW function and (4) presenting a discussion of the applicability of the slow solutions. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>Function <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math> for a typical O5 V star without rotation. The (<b>left panel</b>) shows only the function <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math>, while the (<b>right panel</b>) is similar to the (<b>left panel</b>), but the plane <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> is also plotted in light grey. Furthermore, the intersection of both curves (black solid lines) shows two loci of singular points.</p>
Full article ">Figure 2
<p>Velocity profile for a typical O5 V star without rotation. The velocity profile is plotted as a function of <math display="inline"><semantics> <mrow> <mo form="prefix">log</mo> <mo>(</mo> <mi>r</mi> <mo>/</mo> <msub> <mi>R</mi> <mo>*</mo> </msub> <mo>−</mo> <mn>1</mn> <mo>)</mo> </mrow> </semantics></math> (<b>left panel</b>) and as a function of <span class="html-italic">u</span> (<b>right panel</b>). The location of the singular point is shown with a red dot, while the sonic point is in black.</p>
Full article ">Figure 3
<p>The radiative acceleration, <math display="inline"><semantics> <msup> <mi>g</mi> <mi>line</mi> </msup> </semantics></math>, for a typical O5 V star without rotation as function of <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>/</mo> <msub> <mi>R</mi> <mo>*</mo> </msub> </mrow> </semantics></math> for <math display="inline"><semantics> <mrow> <mi>r</mi> <mo>&lt;</mo> <mn>10</mn> <mspace width="0.166667em"/> <msub> <mi>R</mi> <mo>*</mo> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Dependence of the wind parameters as a function of the line force parameter <math display="inline"><semantics> <mi>α</mi> </semantics></math>. Terminal velocity (<b>left panel</b>) and mass loss rate (<b>right panel</b>). The values obtained for our typical O5 V star without rotation are shown in red.</p>
Full article ">Figure 5
<p>Dependence of the wind parameters as a function of the line force parameter <span class="html-italic">k</span>. Terminal velocity (<b>left panel</b>) and mass loss rate (<b>right panel</b>). The values obtained for our typical O5 V star without rotation are shown in red.</p>
Full article ">Figure 6
<p>Dependence of the wind parameters as a function of the line force parameter <math display="inline"><semantics> <mi>δ</mi> </semantics></math>. Upper panels are for <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mn>0.124</mn> </mrow> </semantics></math> and lower panels are for <math display="inline"><semantics> <mrow> <mi>k</mi> <mo>=</mo> <mn>0.0124</mn> </mrow> </semantics></math>. The values obtained for our typical O5 V star without rotation are shown in red.</p>
Full article ">Figure 7
<p>The function <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math> (topology) of the m-CAK theory as function of <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>: <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math> (<b>upper left panel</b>), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math> (<b>upper right panel</b>), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.7</mn> </mrow> </semantics></math> (<b>lower left panel</b>) and <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.9</mn> </mrow> </semantics></math> (<b>lower right panel</b>). The plane <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> is shown in light grey, and its intersection with the surface <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math> (locus of singular points) is plotted with a black line.</p>
Full article ">Figure 8
<p>Velocity profiles, <math display="inline"><semantics> <mrow> <mi>v</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math>, as function of the rotational rate <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>. (<b>Left panel</b>): fast solutions (<math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>≲</mo> <mn>0.74</mn> </mrow> </semantics></math>) are plotted in grey lines, while the <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>-slow solutions are in coloured lines: <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.74</mn> </mrow> </semantics></math> (red line), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.76</mn> </mrow> </semantics></math> (blue line), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.80</mn> </mrow> </semantics></math> (cyan line), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.82</mn> </mrow> </semantics></math> (magenta line), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.85</mn> </mrow> </semantics></math> (green line), <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.90</mn> </mrow> </semantics></math> (black line) and <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">Ω</mi> <mo>=</mo> <mn>0.95</mn> </mrow> </semantics></math> (orange line). (<b>Right panel</b>): The same <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>-slow solutions, but zoomed and including the location of the singular points (red dots).</p>
Full article ">Figure 9
<p>Wind density profiles, <math display="inline"><semantics> <mrow> <mi>ρ</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math> (in gr/cm<math display="inline"><semantics> <msup> <mrow/> <mn>3</mn> </msup> </semantics></math>) versus <span class="html-italic">u</span>, for fast and <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>-slow solutions. The colour scheme is the same as the one used in <a href="#galaxies-11-00068-f008" class="html-fig">Figure 8</a>.</p>
Full article ">Figure 10
<p>The topological function <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math> of the m-CAK theory as function of <math display="inline"><semantics> <mi>δ</mi> </semantics></math>. (<b>Upper left panel</b>): <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>. (<b>Upper right panel</b>): <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.12</mn> </mrow> </semantics></math>. (<b>Lower left panel</b>): <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.24</mn> </mrow> </semantics></math>. (<b>Lower right panel</b>): <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math>. The plane <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math> is shown in light grey, and its intersection with the surface <math display="inline"><semantics> <mrow> <mi>R</mi> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>Z</mi> <mo>)</mo> </mrow> </semantics></math> (locus of singular points) is plotted with black lines.</p>
Full article ">Figure 11
<p>Velocity profiles, <math display="inline"><semantics> <mrow> <mi>v</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math>, for different values of the line force parameter <math display="inline"><semantics> <mi>δ</mi> </semantics></math>. (<b>Left panel</b>): fast solutions are plotted in grey lines for <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.0</mn> <mo>,</mo> <mn>0.1</mn> <mo>,</mo> <mn>0.2</mn> <mo>,</mo> <mn>0.24</mn> </mrow> </semantics></math>, while <math display="inline"><semantics> <mi>δ</mi> </semantics></math>-slow solutions are in coloured lines. The red line corresponds to <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math>, the blue line to <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.31</mn> </mrow> </semantics></math>, the cyan line to <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.32</mn> </mrow> </semantics></math> and the magenta line to <math display="inline"><semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0.33</mn> </mrow> </semantics></math>. The (<b>right panel</b>) shows only <math display="inline"><semantics> <mi>δ</mi> </semantics></math>-slow solutions, and the location of the singular points for each solution is shown with a red dot.</p>
Full article ">Figure 12
<p>Fast solution velocity profile (solid blue), <math display="inline"><semantics> <mrow> <mi>v</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math> vs. <span class="html-italic">u</span>. Six different <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law velocity profiles are also plotted. It is clearly seen that the <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law approximation is a good one for <math display="inline"><semantics> <mrow> <mn>0.7</mn> <mo>≲</mo> <mi>β</mi> <mo>≲</mo> <mn>1.2</mn> </mrow> </semantics></math>. See text for details.</p>
Full article ">Figure 13
<p><math display="inline"><semantics> <mi>δ</mi> </semantics></math>-slow solution velocity profile (solid blue), <math display="inline"><semantics> <mrow> <mi>v</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math> vs. <span class="html-italic">u</span>. Six different <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law velocity profiles are also plotted. For values around <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>=</mo> <mn>0.7</mn> </mrow> </semantics></math>, the profiles can be considered similar, but it can be clearly concluded that for <math display="inline"><semantics> <mrow> <mi>β</mi> <mo>&gt;</mo> <mn>0.7</mn> </mrow> </semantics></math>, the <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law profile cannot fit the m-CAK hydrodynamical <math display="inline"><semantics> <mi>δ</mi> </semantics></math>-slow solution.</p>
Full article ">Figure 14
<p><math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>-slow solution velocity profile (solid blue), <math display="inline"><semantics> <mrow> <mi>v</mi> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </semantics></math> vs. <span class="html-italic">u</span>. Six different <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law velocity profiles are also plotted. It can be clearly concluded that no <math display="inline"><semantics> <mi>β</mi> </semantics></math>-law profile can fit the m-CAK hydrodynamical <math display="inline"><semantics> <mi mathvariant="sans-serif">Ω</mi> </semantics></math>-slow solution.</p>
Full article ">Figure 15
<p>Velocity profiles as a function of the inverse radial coordinate <math display="inline"><semantics> <mrow> <mi>u</mi> <mo>=</mo> <mo>−</mo> <msub> <mi>R</mi> <mo>*</mo> </msub> <mo>/</mo> <mi>r</mi> <mo>=</mo> <mo>−</mo> <mn>1</mn> <mo>/</mo> <mover accent="true"> <mi>r</mi> <mo>^</mo> </mover> </mrow> </semantics></math> for four models. The hydrodynamic results from <span class="html-small-caps">Hydwind</span> are shown in solid blue lines and the analytical solutions are shown by dashed lines. The stellar and line force parameters for the models are given in Araya et al. [<a href="#B60-galaxies-11-00068" class="html-bibr">60</a>].</p>
Full article ">Figure 16
<p>Velocity profiles of <math display="inline"><semantics> <mi>ϵ</mi> </semantics></math> Ori as a function of <math display="inline"><semantics> <mrow> <mo form="prefix">log</mo> <mo>(</mo> <mi>r</mi> <mo>/</mo> <msub> <mi>R</mi> <mo>*</mo> </msub> <mo>−</mo> <mn>1</mn> <mo>)</mo> </mrow> </semantics></math> in a region near to the stellar surface. The solid blue line shows the numerical hydrodynamic result and the analytical solution is shown by a dashed line. The dot symbol indicates the position of the sonic (or critical) point. The difference between both curves is around one thermal speed.</p>
Full article ">Figure 17
<p>Real and complex regions where the line acceleration expression given by Villata [<a href="#B58-galaxies-11-00068" class="html-bibr">58</a>] can be found. These regions are delimited by the values of the line force parameters <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>δ</mi> </semantics></math>.</p>
Full article ">Figure 18
<p>Hydrodynamic models in the <math display="inline"><semantics> <msub> <mi>T</mi> <mi>eff</mi> </msub> </semantics></math>-<math display="inline"><semantics> <mrow> <mo form="prefix">log</mo> <mspace width="0.166667em"/> <mi>g</mi> </mrow> </semantics></math> plane. Blue dots represent the converged solutions. Grey solid lines are the evolutionary tracks for stars of <math display="inline"><semantics> <mrow> <mn>7</mn> <mo> </mo> <msub> <mi>M</mi> <mo>⊙</mo> </msub> </mrow> </semantics></math> to <math display="inline"><semantics> <mrow> <mn>60</mn> <mo> </mo> <msub> <mi>M</mi> <mo>⊙</mo> </msub> </mrow> </semantics></math> without rotation [<a href="#B72-galaxies-11-00068" class="html-bibr">72</a>], and black lines represent the zero-age main sequence (ZAMS) and the terminal age main sequence (TAMS).</p>
Full article ">Figure 19
<p>Velocity profiles as a function of the inverse radial coordinate <math display="inline"><semantics> <mrow> <mi>u</mi> <mo>=</mo> <mo>−</mo> <msub> <mi>R</mi> <mo>*</mo> </msub> <mo>/</mo> <mi>r</mi> <mo>=</mo> <mo>−</mo> <mn>1</mn> <mo>/</mo> <mover accent="true"> <mi>r</mi> <mo>^</mo> </mover> </mrow> </semantics></math> for the model R24 from Curé et al. [<a href="#B27-galaxies-11-00068" class="html-bibr">27</a>]. The hydrodynamic result from <span class="html-small-caps">Hydwind</span> is shown in the solid blue line and the analytical solution is a dashed black line.</p>
Full article ">
20 pages, 930 KiB  
Review
Wide-Angle-Tail (WAT) Radio Sources
by Christopher P. O’Dea and Stefi A. Baum
Galaxies 2023, 11(3), 67; https://doi.org/10.3390/galaxies11030067 - 12 May 2023
Cited by 14 | Viewed by 2750
Abstract
We review the properties of Wide-Angle-Tail (WAT) radio sources. The WAT radio sources are powerful, bent radio sources typically associated with the dominant galaxy in a cluster or group. For the purpose of this review, we define the radio morphology properties of WATs [...] Read more.
We review the properties of Wide-Angle-Tail (WAT) radio sources. The WAT radio sources are powerful, bent radio sources typically associated with the dominant galaxy in a cluster or group. For the purpose of this review, we define the radio morphology properties of WATs as (1) a sudden jet-tail transition, (2) overall bending of the tails to one side, and (3) non-parallel tails. The mechanism for the rapid jet-tail transition is uncertain but it seems to occur near the transition from the host ISM to ICM. The jet-tail transition may make the jets easier to bend. The narrow range in radio luminosity can be understood if there is a minimum luminosity required to allow the jets to propagate undisturbed for tens of kpc and a maximum luminosity required to allow the jet disruption mechanism to act. WATs are typically hosted by the brightest cluster galaxies in clusters which are currently merging. Thus, WATs can be used as tracers of merging clusters. The merging produces large-scale bulk motions in the ICM which can provide sufficient ram pressure to bend the jets. We suggest that although the Lorentz force may not bend the jets in WATs, it may be relevant in other sources, e.g., protostellar jets. Full article
(This article belongs to the Special Issue New Perspectives on Radio Galaxy Dynamics)
Show Figures

Figure 1

Figure 1
<p>The WAT 3C465 in A2634. JVLA 1.5 GHz image with a resolution of 1.37 × 1.13 arcsec taken from [<a href="#B43-galaxies-11-00067" class="html-bibr">43</a>].</p>
Full article ">Figure 2
<p>The ratio of jet length (longer/shorter) to average jet length in WATs. Data are taken from O’Donoghue et al. [<a href="#B5-galaxies-11-00067" class="html-bibr">5</a>].</p>
Full article ">Figure 3
<p>The WAT 1233+168 in the merging cluster A1569. Contours of the VLA image are shown on the Chandra X-ray image. The arrow shows the bisection of the angle between the two tails of the WAT, showing that the elongated X-ray emission extends down the middle of the radio tails. The Figure, including all annotations, is from Tiwari and Singh [<a href="#B115-galaxies-11-00067" class="html-bibr">115</a>].</p>
Full article ">
8 pages, 904 KiB  
Communication
Advanced Life Peaked Billions of Years Ago According to Black Holes
by David Garofalo
Galaxies 2023, 11(3), 66; https://doi.org/10.3390/galaxies11030066 - 11 May 2023
Cited by 1 | Viewed by 2183
Abstract
The link between black holes and star formation allows for us to draw a connection between black holes and the places and times when extraterrestrial intelligences (ETIs) had a greater chance of emerging. Within the context of the gap paradigm for black holes, [...] Read more.
The link between black holes and star formation allows for us to draw a connection between black holes and the places and times when extraterrestrial intelligences (ETIs) had a greater chance of emerging. Within the context of the gap paradigm for black holes, we show that denser cluster environments that led to gas-rich mergers and copious star formation were places less compatible on average with the emergence of ETIs compared to isolated elliptical galaxies by almost two orders of magnitude. The probability for ETIs peaked in these isolated environments around 6 billion years ago and cosmic downsizing shifted the likelihood of ETIs emerging to galaxies with weak black hole feedback, such as in spiral galaxies, at late times. Full article
Show Figures

Figure 1

Figure 1
<p>Isolated environments: Following a gas-rich merger, the accretion disk (in blue) settles in a counterrotating configuration around a rapidly spinning black hole. This configuration produces a narrow jet whose effect is to drive star formation rates to higher values (<b>left panel</b>). This enhancement lasts at most 8 million years after which the tilted disk fails to produce a jet, so star formation is unaffected (<b>middle panel</b>). The star formation rate dies down as the cold gas reservoir comes to an end and a dead quasar is formed after over a billion years (<b>right panel</b>).</p>
Full article ">Figure 2
<p>Richest environments: Following a gas-rich merger, the accretion disk (in blue) settles in a counterrotating configuration around a rapidly spinning black hole. This configuration produces a narrow jet whose effect is to drive star formation rates to higher values (<b>left panel</b>). This enhancement lasts for a few hundred million years after which the disk becomes tilted, a new jet direction is generated, and jet feedback suppresses star formation and heats the interstellar medium (ISM) with X-rays (<b>middle panel</b>). After billions of years, the gas supply dies down and the ISM of the dead quasar is suffused in a hot X-ray halo (<b>right panel</b>).</p>
Full article ">
17 pages, 1987 KiB  
Article
Faint Galaxy Number Counts in the Durham and SDSS Catalogues
by John H. Marr
Galaxies 2023, 11(3), 65; https://doi.org/10.3390/galaxies11030065 - 7 May 2023
Cited by 1 | Viewed by 1840
Abstract
Galaxy number counts in the K-, H-, I-, R-, B- and U-bands from the Durham Extragalactic Astronomy and Cosmology catalogue could be well-fitted over their whole range using luminosity function (LF) parameters derived from the SDSS at [...] Read more.
Galaxy number counts in the K-, H-, I-, R-, B- and U-bands from the Durham Extragalactic Astronomy and Cosmology catalogue could be well-fitted over their whole range using luminosity function (LF) parameters derived from the SDSS at the bright region and required only modest luminosity evolution with the steepening of the LF slope (α), except for a sudden steep increase in the B-band and a less steep increase in the U-band at faint magnitudes that required a starburst evolutionary model to account for the excess faint number counts. A cosmological model treating Hubble expansion as an Einstein curvature required less correction at faint magnitudes than a standard ΛCDM model, without requiring dark matter or dark energy. Data from DR17 of the SDSS in the g, i, r, u and z bands over two areas of the sky centred on the North Galactic Cap (NGC) and above the South Galactic Cap (SGC), with areas of 5954 and 859 sq. deg., respectively, and a combined count of 622,121 galaxies, were used to construct bright galaxy number counts and galaxy redshift/density plots within the limits of redshift 0.4 and mag 20. Their comparative densities confirmed an extensive void in the Southern sky with a deficit of 26% out to a redshift z ≤ 0.15. Although not included in the number count data set because of its incompleteness at fainter magnitudes, extending the SDSS redshift-number count survey to fainter and more distant galaxies with redshift ≤ 1.20 showed a secondary peak in the number counts with many QSOs, bright X-ray and radio sources, and evolving irregular galaxies with rapid star formation rates. This sub-population at redshifts of 0.45–0.65 may account for the excess counts observed in the B-band. Recent observations from the HST and James Webb Space Telescope (JWST) have also begun to reveal a high density of massive galaxies at high redshifts (z>7) with high UV and X-ray emissions, and future observations by the JWST may reveal the assembly of galaxies in the early universe going back to the first light in the universe. Full article
Show Figures

Figure 1

Figure 1
<p>A graphic showing the observed volume/square degree with redshift.</p>
Full article ">Figure 2
<p>Building a theoretical number count curve (red line) with composite Schechter curves at increasing redshifts (<math display="inline"><semantics> <mrow> <msup> <mi>M</mi> <mo>*</mo> </msup> <mo>=</mo> <mo>−</mo> <mn>23.0</mn> <mo>;</mo> <mi>α</mi> <mo>=</mo> <mo>−</mo> <mn>1.5</mn> </mrow> </semantics></math>). The transition from a Euclidean slope of 0.6 to the <math display="inline"><semantics> <mi>α</mi> </semantics></math>-dependent slope is clearly shown (dashed lines).</p>
Full article ">Figure 3
<p>Schechter curves for 5 bands at uniform redshift (<math display="inline"><semantics> <mrow> <mi>z</mi> <mo>=</mo> <mn>0.11</mn> </mrow> </semantics></math>), after Blanton et al. [<a href="#B16-galaxies-11-00065" class="html-bibr">16</a>].</p>
Full article ">Figure 4
<p>SDSS galaxy density in redshift bins of 0.0125. NGC (red bars) and SGC (blue bars). <math display="inline"><semantics> <mrow> <msub> <mi>B</mi> <mi>J</mi> </msub> <mo>≤</mo> <mn>20</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Galaxy number counts for the five SDSS observational bands (crosses). Overlain are the bright <span class="html-italic">B</span>-band counts (circles ± bin ranges) and the SDSS-derived <span class="html-italic">B</span> curve (solid red line).</p>
Full article ">Figure 6
<p><span class="html-italic">K</span>-magnitude plots with no evolution (dashed line, <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mo>−</mo> <mn>1.03</mn> </mrow> </semantics></math>) and pure luminosity evolution (PLE) from <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mo>−</mo> <mn>1.03</mn> </mrow> </semantics></math> to <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mo>−</mo> <mn>1.70</mn> </mrow> </semantics></math> (solid line) in the selected GR model, both with <math display="inline"><semantics> <mrow> <msup> <mi>M</mi> <mo>*</mo> </msup> <mo>=</mo> <mo>−</mo> <mn>23</mn> <mo>+</mo> <mn>5</mn> <msub> <mo>log</mo> <mn>10</mn> </msub> </mrow> </semantics></math>h. The error bars reflect the maximum range of each bin; the absence of error bars implies the bin has a single member. The GR model fully overlays a <math display="inline"><semantics> <mi mathvariant="normal">Λ</mi> </semantics></math>CDM model with <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">Ω</mi> <mi>m</mi> </msub> <mo>=</mo> <mn>0.3</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">Ω</mi> <mi mathvariant="normal">Λ</mi> </msub> <mo>=</mo> <mn>0.7</mn> </mrow> </semantics></math>. The Euclidean 0.6 slope is also shown (dashed black line). Adapted from composite means of 42 <span class="html-italic">K</span>-band surveys from Durham number count archive [<a href="#B1-galaxies-11-00065" class="html-bibr">1</a>].</p>
Full article ">Figure 7
<p>All six colour bands: GR model with PLE for <span class="html-italic">K</span>-, <span class="html-italic">H</span>-, <span class="html-italic">I</span>- and <span class="html-italic">R</span>-bands and starburst GR model for <span class="html-italic">B</span>- and <span class="html-italic">U</span>-bands. Composite means of all surveys (error bars omitted for clarity). (Data from Durham number count archive [<a href="#B1-galaxies-11-00065" class="html-bibr">1</a>].)</p>
Full article ">Figure 8
<p><span class="html-italic">B</span>-band magnitude counts binned ±0.25 m (composite means of 35 <math display="inline"><semantics> <msub> <mi>B</mi> <mi>J</mi> </msub> </semantics></math>-band surveys) overlain with the GR starburst model (blue line) with starburst evolution from <math display="inline"><semantics> <mrow> <mn>0.3</mn> <mo>&lt;</mo> <mi>z</mi> <mo>&lt;</mo> <mn>1.2</mn> </mrow> </semantics></math>, and the <math display="inline"><semantics> <mi mathvariant="normal">Λ</mi> </semantics></math>CDM model (red line). Bars indicate the range of observations in each bin. Adapted from Durham galaxy number count archive [<a href="#B1-galaxies-11-00065" class="html-bibr">1</a>].</p>
Full article ">Figure 9
<p>Change in slopes of <span class="html-italic">B</span>-band magnitude counts compared to <span class="html-italic">K</span>-mag counts.</p>
Full article ">Figure 10
<p>Log–log plot of the observed emission wavelenths vs. redshift for the six optical bands (solid lines) with reference lines for each band (dashed lines). Diamonds show the cut-off redshifts used to generate the number count plots for each band.</p>
Full article ">
14 pages, 37354 KiB  
Article
Large-Scale Ejecta of Z CMa—Proper Motion Study and New Features Discovered
by Tiina Liimets, Michaela Kraus, Lydia Cidale, Sergey Karpov and Anthony Marston
Galaxies 2023, 11(3), 64; https://doi.org/10.3390/galaxies11030064 - 4 May 2023
Cited by 1 | Viewed by 1425
Abstract
Z Canis Majoris is a fascinating early-type binary with a Herbig Be primary and a FU Orionis-type secondary. Both of the stars exhibit sub-arcsecond jet-like ejecta. In addition, the primary is associated with the extended jet as well as with the large-scale outflow. [...] Read more.
Z Canis Majoris is a fascinating early-type binary with a Herbig Be primary and a FU Orionis-type secondary. Both of the stars exhibit sub-arcsecond jet-like ejecta. In addition, the primary is associated with the extended jet as well as with the large-scale outflow. In this study, we investigate further the nature of the large-scale outflow, which has not been studied since its discovery almost three and a half decades ago. We present proper motion measurements of individual features of the large-scale outflow and determine their kinematical ages. Furthermore, with our newly acquired deep images, we have discovered additional faint arc-shaped features that can be associated with the central binary. Full article
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)
Show Figures

Figure 1

Figure 1
<p>(<b>Left</b>): Schematic view of the various extended nebular features surrounding Z CMa. The stars in the field of view are drawn as black filled circles, large-scale outflow features are marked with grayish areas or circles with black contures. The dashed line is drawn at the 60° position angle. Individual numbers refer to radial velocities (Poetzel et al. [<a href="#B5-galaxies-11-00064" class="html-bibr">5</a>]). Field of view is <math display="inline"><semantics> <mrow> <msup> <mn>9</mn> <mo>′</mo> </msup> <mo>.</mo> <mn>5</mn> <mo>×</mo> <msup> <mn>5</mn> <mo>′</mo> </msup> <mo>.</mo> <mn>6</mn> </mrow> </semantics></math>. The base of the figure is depicted from Poetzel et al. [<a href="#B5-galaxies-11-00064" class="html-bibr">5</a>]. Reproduced with permission © ESO. (<b>Right</b>): Schematic view of the sub-arcsecond features around Z CMa presented in true proportions. The lengths of the micro-jets are taken from Whelan et al. [<a href="#B6-galaxies-11-00064" class="html-bibr">6</a>] and those for the streamer are taken from Dong et al. [<a href="#B2-galaxies-11-00064" class="html-bibr">2</a>]. The dashed line represents the position angle of the large-scale outflow. Field of view is <math display="inline"><semantics> <mrow> <msup> <mn>1</mn> <mrow> <mo>″</mo> </mrow> </msup> <mo>.</mo> <mn>6</mn> <mo>×</mo> <msup> <mn>2</mn> <mrow> <mo>″</mo> </mrow> </msup> <mo>.</mo> <mn>7</mn> </mrow> </semantics></math>. A white square representing the same size is drawn on the left panel at the position of the central binary inside the comma nebula. On both panels, north is up and east is to the left.</p>
Full article ">Figure 2
<p>(<b>Left</b>): [S <span class="html-small-caps">ii</span>] image of Z CMa acquired with GMOS attached to Gemini-South. The white square shape feature, indicated with a white arrow, is an artifact from vignetting of the guiding probe. (<b>Right</b>): Insets of the resolved features in the [S <span class="html-small-caps">ii</span>] image. The FOVs of the smallest insets are <math display="inline"><semantics> <mrow> <msup> <mn>20</mn> <mrow> <mo>″</mo> </mrow> </msup> <mo>×</mo> <msup> <mn>15</mn> <mrow> <mo>″</mo> </mrow> </msup> </mrow> </semantics></math> each. On all images, north is up, east is to the left, and the intensity is in log scale to improve the contrast. See text for more details.</p>
Full article ">Figure 3
<p>(<b>Left</b>): H<math display="inline"><semantics> <mi>α</mi> </semantics></math> image of Z CMa taken in 2002. GMOS FOV is shown for comparison. (<b>Right</b>): Insets of the features resolvable in the H<math display="inline"><semantics> <mi>α</mi> </semantics></math> image and which are outside GMOS FOV. On all images, north is up, east is to the left, and the intensity is in log scale to improve the contrast. See text for more details.</p>
Full article ">Figure 4
<p>Pan-STARRS RGB image presenting the faint extended features around Z CMa. Red channel corresponds to <span class="html-italic">z</span> filter, green channel corresponds to <span class="html-italic">i</span>, and blue channel corresponds to <span class="html-italic">g</span> filter. The intensity is in log scale to improve the contrast. North is up and east is to the left. FOV <math display="inline"><semantics> <mrow> <msup> <mn>9</mn> <mo>′</mo> </msup> <mo>×</mo> <msup> <mn>9</mn> <mo>′</mo> </msup> </mrow> </semantics></math>. See text for more details.</p>
Full article ">
16 pages, 433 KiB  
Article
Language Models for Multimessenger Astronomy
by Vladimir Sotnikov and Anastasiia Chaikova
Galaxies 2023, 11(3), 63; https://doi.org/10.3390/galaxies11030063 - 1 May 2023
Cited by 4 | Viewed by 3493
Abstract
With the increasing reliance of astronomy on multi-instrument and multi-messenger observations for detecting transient phenomena, communication among astronomers has become more critical. Apart from automatic prompt follow-up observations, short reports, e.g., GCN circulars and ATels, provide essential human-written interpretations and discussions of observations. [...] Read more.
With the increasing reliance of astronomy on multi-instrument and multi-messenger observations for detecting transient phenomena, communication among astronomers has become more critical. Apart from automatic prompt follow-up observations, short reports, e.g., GCN circulars and ATels, provide essential human-written interpretations and discussions of observations. These reports lack a defined format, unlike machine-readable messages, making it challenging to associate phenomena with specific objects or coordinates in the sky. This paper examines the use of large language models (LLMs)—machine learning models with billions of trainable parameters or more that are trained on text—such as InstructGPT-3 and open-source Flan-T5-XXL for extracting information from astronomical reports. The study investigates the zero-shot and few-shot learning capabilities of LLMs and demonstrates various techniques to improve the accuracy of predictions. The study shows the importance of careful prompt engineering while working with LLMs, as demonstrated through edge case examples. The study’s findings have significant implications for the development of data-driven applications for astrophysical text analysis. Full article
(This article belongs to the Special Issue The New Era of Real-Time Multi-Messenger Astronomy)
Show Figures

Figure 1

Figure 1
<p>An example message from ATel. Named entities are marked with color: <span style="color: #BAD7DE">■</span> object name; <span style="color: #F4FFAA">■</span> the type of the object or physical phenomena; <span style="color: #C6DDC4">■</span> event type. Taken from Ref. [<a href="#B1-galaxies-11-00063" class="html-bibr">1</a>].</p>
Full article ">Figure 2
<p>An example of description prompt for the type of the object or physical phenomena.</p>
Full article ">Figure 3
<p>An example of a “few−shot” prompt for an event type. T. and A. denote input text and expected answer correspondingly. Named entities are marked with color: <span style="color: #BAD7DE">■</span> object name; <span style="color: #F4FFAA">■</span> the type of the object or physical phenomena; <span style="color: #C6DDC4">■</span> event type.</p>
Full article ">Figure 4
<p>An example of explanation inserted in “few−shot” prompt for the event type. A. denotes the expected answer.</p>
Full article ">Figure 5
<p>An example of prompt obtained by APE method for the type of the object or physical phenomena.</p>
Full article ">Figure 6
<p>An example of a prompt for coming up with motivation for an answer.</p>
Full article ">Figure 7
<p>An example of output with three entities.</p>
Full article ">Figure 8
<p>Steps of the embedding ranking pipeline (from up to bottom). Step I. Calculate the top three answers and their embeddings for a given text message and a prompt. Step II. Predict an embedding of a correct answer for a given text message. Step III. Calculate the L2 distance between answer embedding and predicted embedding. The answer with the lowest L2 is assumed to be correct.</p>
Full article ">Figure 9
<p>Comparison of FFN validation accuracy during training depending on a loss function and a hidden dimension size. Shown on the example of event type entity and fine-tuned <tt>pearsonkyle/gpt2-exomachina</tt> embeddings. The MSELoss is represented by the solid line, ContrastiveLoss by the dashed line, and TripletMarginLoss by the line with a dot.</p>
Full article ">
20 pages, 2221 KiB  
Review
Opacities and Atomic Diffusion
by Georges Alecian and Morgan Deal
Galaxies 2023, 11(3), 62; https://doi.org/10.3390/galaxies11030062 - 25 Apr 2023
Cited by 1 | Viewed by 1819
Abstract
Opacity is a fundamental quantity for stellar modeling, and it plays an essential role throughout the life of stars. After gravity drives the collapse of interstellar matter into a protostar, the opacity determines how this matter is structured around the stellar core. The [...] Read more.
Opacity is a fundamental quantity for stellar modeling, and it plays an essential role throughout the life of stars. After gravity drives the collapse of interstellar matter into a protostar, the opacity determines how this matter is structured around the stellar core. The opacity explains how the radiation field interacts with the matter and how a major part of the energy flows through the star. It results from all the microscopic interactions of photons with atoms. Part of the momentum exchange between photons and atoms gives rise to radiative accelerations (specific to each type of atom), which are strongly involved in a second-order process: atomic diffusion. Although this process is a slow one, it can have a significant impact on stellar structure and chemical composition measurements. In this review, we discuss the way opacities are presently computed and used in numerical codes. Atomic diffusion is described, and the current status of the consideration of this process is presented. Full article
(This article belongs to the Special Issue The Structure and Evolution of Stars)
Show Figures

Figure 1

Figure 1
<p>An example of monochromatic opacity. The case of calcium from OP data for a temperature of <math display="inline"><semantics> <msup> <mn>10</mn> <mn>5</mn> </msup> </semantics></math> K and density of <math display="inline"><semantics> <mrow> <mn>7.8</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>7</mn> </mrow> </msup> </mrow> </semantics></math> g·cm<math display="inline"><semantics> <msup> <mrow/> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </semantics></math>.</p>
Full article ">Figure 2
<p>Profiles of the logarithm of the Rosseland mean opacity according to the logarithm of the temperature for <math display="inline"><semantics> <mrow> <mn>1.0</mn> </mrow> </semantics></math> (blue) and <math display="inline"><semantics> <mrow> <mn>1.7</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> (orange) models at <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.55</mn> </mrow> </semantics></math>. Higher temperatures correspond to deeper layers inside the models. The dashed green curve represents the opacity profile for the same model but using the Rosseland mean opacity not computed from the monochromatic opacities with the exact mixture. The vertical dotted line shows the bottom of the mixed region, where the iron accumulation is induced by radiative acceleration, i.e., where iron is one of the main contributors to the opacity.</p>
Full article ">Figure 3
<p>Ratio of the SCO-RCG (<math display="inline"><semantics> <msub> <mi>κ</mi> <mrow> <mi>R</mi> <mo>,</mo> <mi>S</mi> <mi>R</mi> </mrow> </msub> </semantics></math>) and OP (<math display="inline"><semantics> <msub> <mi>κ</mi> <mrow> <mi>R</mi> <mo>,</mo> <mi>O</mi> <mi>P</mi> </mrow> </msub> </semantics></math>) Rosseland means vs. log T for the stellar mixture (thick line) and for Ni alone (thin line) for a <math display="inline"><semantics> <mrow> <mn>9.5</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> main-sequence model without atomic diffusion. The dotted line represents a ratio of one. Courtesy of Alain Hui-Bon-Hoa, adapted from [<a href="#B11-galaxies-11-00062" class="html-bibr">11</a>].</p>
Full article ">Figure 4
<p>Profiles of the logarithm of the mass fraction according to the radius of several elements for a <math display="inline"><semantics> <mrow> <mn>1.0</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> model corresponding to ages of 3.48 Gyr. The model is computed with the Cesam2k20 stellar evolution code [<a href="#B52-galaxies-11-00062" class="html-bibr">52</a>,<a href="#B57-galaxies-11-00062" class="html-bibr">57</a>,<a href="#B58-galaxies-11-00062" class="html-bibr">58</a>]. The track is not calibrated on the Sun. The grey area represents the convective region. The vertical dash line represents the boundary between convective and radiative zones.</p>
Full article ">Figure 5
<p>The same legend as <a href="#galaxies-11-00062-f004" class="html-fig">Figure 4</a> for <math display="inline"><semantics> <mrow> <mn>1.7</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> models at <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.45</mn> </mrow> </semantics></math> (<b>top panel</b>) and <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math> (<b>bottom panel</b>). The ages are 0.98 Gyr and 1.76 Gyr, respectively. The grey and green areas represent the convective and turbulent mixing regions, respectively. In the bottom panel, the turbulent mixing region is located inside the convective envelope.</p>
Full article ">Figure 6
<p>Same legend as <a href="#galaxies-11-00062-f005" class="html-fig">Figure 5</a> for <math display="inline"><semantics> <mrow> <mn>5.0</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> models. The ages are 57 Myr and 96 Myr, respectively.</p>
Full article ">Figure 7
<p>Chromium cloud in an atmosphere <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mi>eff</mi> </msub> <mo>=</mo> <mn>12</mn> <mo>,</mo> <mn>000</mn> </mrow> </semantics></math> K, and <math display="inline"><semantics> <mrow> <mo form="prefix">log</mo> <mi>g</mi> <mo>=</mo> <mn>4.0</mn> </mrow> </semantics></math>, assuming an anisotropic wind. The top panel shows the positions (point 1 and point 2) for probing the atmosphere (Hammer projection of the full stellar surface). The right panel presents the chromium abundance stratifications at the last time step, the blue line corresponding to point 1, the pink line to point 2. The horizontal solid black line represents the initial solar abundance. The vertical dashed lines show the depth limits of the slab shown in the top panel.</p>
Full article ">Figure 8
<p>Rosseland mean opacity profiles of <math display="inline"><semantics> <mrow> <mn>1.4</mn> </mrow> </semantics></math> M<math display="inline"><semantics> <msub> <mrow/> <mo>⊙</mo> </msub> </semantics></math> models with <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.6</mn> </mrow> </semantics></math> (solid curves), <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.4</mn> </mrow> </semantics></math> (dashed curves), and <math display="inline"><semantics> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mn>0.2</mn> </mrow> </semantics></math> (dotted curves) from [<a href="#B52-galaxies-11-00062" class="html-bibr">52</a>]. The blue and red curves represent, respectively, the models without and with radiative accelerations taken into account in the diffusion velocity. The solid dashed and dotted vertical lines represent the positions of the bottom of the surface convection zone for the model without radiative accelerations for the same value of <math display="inline"><semantics> <msub> <mi>X</mi> <mi>c</mi> </msub> </semantics></math> as the opacity profiles (for clarity, they are not represented for the model with grad). Adapted from <a href="#galaxies-11-00062-f005" class="html-fig">Figure 5</a> of [<a href="#B52-galaxies-11-00062" class="html-bibr">52</a>].</p>
Full article ">
38 pages, 9925 KiB  
Article
Key Science Goals for the Next-Generation Event Horizon Telescope
by Michael D. Johnson, Kazunori Akiyama, Lindy Blackburn, Katherine L. Bouman, Avery E. Broderick, Vitor Cardoso, Rob P. Fender, Christian M. Fromm, Peter Galison, José L. Gómez, Daryl Haggard, Matthew L. Lister, Andrei P. Lobanov, Sera Markoff, Ramesh Narayan, Priyamvada Natarajan, Tiffany Nichols, Dominic W. Pesce, Ziri Younsi, Andrew Chael, Koushik Chatterjee, Ryan Chaves, Juliusz Doboszewski, Richard Dodson, Sheperd S. Doeleman, Jamee Elder, Garret Fitzpatrick, Kari Haworth, Janice Houston, Sara Issaoun, Yuri Y. Kovalev, Aviad Levis, Rocco Lico, Alexandru Marcoci, Niels C. M. Martens, Neil M. Nagar, Aaron Oppenheimer, Daniel C. M. Palumbo, Angelo Ricarte, María  J. Rioja, Freek Roelofs, Ann C. Thresher, Paul Tiede, Jonathan Weintroub and Maciek Wielgusadd Show full author list remove Hide full author list
Galaxies 2023, 11(3), 61; https://doi.org/10.3390/galaxies11030061 - 24 Apr 2023
Cited by 55 | Viewed by 6653
Abstract
The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply [...] Read more.
The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will uniquely enable a wealth of transformative new discoveries related to black hole science, extending from event-horizon-scale studies of strong gravity to studies of explosive transients to the cosmological growth and influence of supermassive black holes. Here, we present the key science goals for the ngEHT and their associated instrument requirements, both of which have been formulated through a multi-year international effort involving hundreds of scientists worldwide. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution of EHT and ngEHT sites around the globe. Sites that have joined EHT campaigns are shown in white (see [<a href="#B2-galaxies-11-00061" class="html-bibr">2</a>]), additional ngEHT Phase-1 sites are shown in cyan, and ngEHT Phase-2 sites are shown in green. Three of the EHT sites have joined since its initial observing campaign in 2017: the 12 m Greenland Telescope (GLT; [<a href="#B49-galaxies-11-00061" class="html-bibr">49</a>]), the 12 m Kitt Peak Telescope (KP), the Northern Extended Millimeter Array (NOEMA) composed of twelve 15 m dishes. Several other existing or upcoming sites that plan to join EHT/ngEHT observations are shown in yellow: the 37 m Haystack Telescope (HAY; [<a href="#B41-galaxies-11-00061" class="html-bibr">41</a>]), the 21 m Yonsei Radio Observatory of the Korea VLBI Network (KVN-YS; [<a href="#B42-galaxies-11-00061" class="html-bibr">42</a>]), the 15 m Africa Millimetre Telescope (AMT; [<a href="#B43-galaxies-11-00061" class="html-bibr">43</a>]), and the 12 m Large Latin American Millimeter Array (LLA; [<a href="#B44-galaxies-11-00061" class="html-bibr">44</a>]). For additional details on the planned ngEHT specifications, see ngEHT Collaboration [<a href="#B50-galaxies-11-00061" class="html-bibr">50</a>].</p>
Full article ">Figure 2
<p>Comparison of image angular resolutions and timescales accessible to the EHT and ngEHT and the associated scientific opportunities. For M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> and Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math>, the ranges of angular resolution and timescale needed to study the three primary domains–fundamental physics, accretion, and jet launching–are indicated with the tilted shaded regions. These shaded regions are centered on the resolution-timescale for each source determined by the speed of light (<math display="inline"><semantics> <mrow> <mi>c</mi> <mi>t</mi> <mo>=</mo> <mi>D</mi> <mi>θ</mi> </mrow> </semantics></math>). Snapshot images require an array to form images on these timescales or shorter; average images require an array to form images over significantly longer timescales; movies require that an array can form images of the full range of timescales from snapshots to averages. The primary difference in M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> and Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> is the factor of <math display="inline"><semantics> <mrow> <mo>∼</mo> </mrow> </semantics></math>1500 difference in the SMBH mass, which sets the system timescale. In contrast, the relevant angular scales in these systems are determined by the mass-to-distance ratio, which only differs by <math display="inline"><semantics> <mrow> <mo>∼</mo> </mrow> </semantics></math>20% for these two SMBHs. The approximate resolution-timescale pair to study each of the ngEHT Key Science Goals is indicated with the inset labeled boxes. Goals associated with Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> or M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> are colored in blue or purple, respectively.</p>
Full article ">Figure 3
<p>Range of observing frequency and angular resolution for selected current and upcoming facilities, from radio to the infrared. The ngEHT can achieve an imaging angular resolution that is significantly finer than any other planned facility or experiment. The ngEHT also envisions simultaneous multi-band observations, extending from 86 to 345 GHz, which will significantly expand the frequency coverage of currently published EHT data (black filled region). Figure adapted from Selina et al. [<a href="#B51-galaxies-11-00061" class="html-bibr">51</a>].</p>
Full article ">Figure 4
<p>BH images display a series of distinctive relativistic features such as the BH apparent “shadow” (e.g., [<a href="#B27-galaxies-11-00061" class="html-bibr">27</a>]), “inner shadow” (e.g., [<a href="#B59-galaxies-11-00061" class="html-bibr">59</a>]), and “photon ring” (e.g., [<a href="#B60-galaxies-11-00061" class="html-bibr">60</a>]).</p>
Full article ">Figure 5
<p>EHT representative average image of Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> using data from 7 April 2017 [<a href="#B11-galaxies-11-00061" class="html-bibr">11</a>]. The white circle in the lower-right shows a <math display="inline"><semantics> <mrow> <mn>20</mn> </mrow> </semantics></math>μas beam that gives the approximate EHT resolution. The overlaid annuli show the predicted ranges of the Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> critical curve using measurements of resolved stellar orbits using the VLTI (blue; [<a href="#B72-galaxies-11-00061" class="html-bibr">72</a>]) and Keck (red; [<a href="#B73-galaxies-11-00061" class="html-bibr">73</a>]); the ranges are dominated by the potential variation in size with spin, <math display="inline"><semantics> <msub> <mi>d</mi> <mi>sh</mi> </msub> </semantics></math> = (9.6–10.4) <math display="inline"><semantics> <msub> <mi>θ</mi> <mi mathvariant="normal">g</mi> </msub> </semantics></math> [<a href="#B25-galaxies-11-00061" class="html-bibr">25</a>,<a href="#B74-galaxies-11-00061" class="html-bibr">74</a>]. The green annulus shows the estimated range (<math display="inline"><semantics> <mrow> <mo>±</mo> <mn>1</mn> <mi>σ</mi> </mrow> </semantics></math>) of the critical curve using EHT measurements, which is consistent with these predictions [<a href="#B14-galaxies-11-00061" class="html-bibr">14</a>]. However, because of the limited baseline coverage of the EHT, key image features such as the azimuthal brightness around the ring and the depth and shape of the central brightness depression are only weakly constrained with current observations.</p>
Full article ">Figure 6
<p>Accessing signatures of the event horizon with the ngEHT. Each panel shows an image on a logarithmic scale, with an inset shown with a linear scale. The left panel shows a time-averaged simulated image of M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math>, which shows a prominent photon ring and inner shadow. The right panel shows a reconstructed ngEHT image using the Bayesian VLBI analysis package <tt>Comrade.jl</tt> [<a href="#B75-galaxies-11-00061" class="html-bibr">75</a>] applied to simulated ngEHT phase-1 observations. The ngEHT provides both the angular resolution and dynamic range required to identify the deep brightness depression produced by the inner shadow in this simulated image.</p>
Full article ">Figure 7
<p>Summary of spin signatures in polarized images of time-averaged GRMHD simulations. In each panel, color indicates brightness and ticks show linear polarization direction. Rows show time-averaged primary (<b>top</b>) and secondary (<b>bottom</b>) images from MAD GRMHD simulations of M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math>; columns show varying BH spin, ranging from a rapidly spinning BH with a retrograde accretion flow (<b>left</b>) to a non-spinning BH (<b>center</b>) to a rapidly spinning BH with a prograde accretion flow (<b>right</b>). The angular radius of the black hole, <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>/</mo> <mi>D</mi> </mrow> </semantics></math>, is identical in each panel. The polarization pattern becomes more radial at higher spin, as frame dragging enforces toroidal magnetic fields near the horizon. In retrograde flows, the spirals pattern reverses handedeness over radius, indicating the transition from the prograde rotation within the ergosphere to the retrograde flow at larger radii. The handedness flips across sub-images, leading to depolarization in the photon ring of the full image (see [<a href="#B83-galaxies-11-00061" class="html-bibr">83</a>,<a href="#B84-galaxies-11-00061" class="html-bibr">84</a>]). By studying the polarized structure and its radial evolution, the ngEHT can estimate the spin of M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> and Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> and quantify the effects of frame dragging. Adapted from Palumbo [<a href="#B85-galaxies-11-00061" class="html-bibr">85</a>].</p>
Full article ">Figure 8
<p>Conceptual illustration of the science cases explored within the “Black holes and their cosmic context” SWG: BH growth, binary BHs and gravitational waves, and MWL studies of BHs and jets. Credits from left to right: Perimeter Institute, NASA’s Goddard Space Flight Center/Jeremy Schnittman and Brian P. Powell, J. C. Algaba for the EHT Collaboration [<a href="#B15-galaxies-11-00061" class="html-bibr">15</a>]. Composition: Thalia Traianou, IAA-CSIC.</p>
Full article ">Figure 9
<p>SMBH population studies with the ngEHT. (<b>left</b>) Black contours show the estimated cumulative number density of SMBHs as a function of shadow diameter and 230 GHz flux density. Colored contours indicate threshold values at which the ngEHT Phase-1 could plausibly measure the SMBH mass (red), spin (green), and shadow (blue) in a superresolution regime. Reproduced from Pesce et al. [<a href="#B129-galaxies-11-00061" class="html-bibr">129</a>]. (<b>right</b>) Estimated 230 GHz compact flux density and BH shadow diameter for a subset of bright VLBI-detected SMBHs in the ETHER database. Colored lines again indicate the approximate measurement thresholds for the ngEHT Phase-1 array to measure the BH mass, spin, and shadow as shown on the left. Adapted from Ramakrishnan et al. [<a href="#B131-galaxies-11-00061" class="html-bibr">131</a>].</p>
Full article ">Figure 10
<p>Studying accretion and jet dynamics with the ngEHT. (<b>left</b>) A frame from a simulated movie of M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> [<a href="#B175-galaxies-11-00061" class="html-bibr">175</a>]. (<b>right</b>) Azimuthal (<b>top</b>) and radial (<b>bottom</b>) brightness variations in a reconstructed movie of M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> using ngEHT Phase-1 coverage. The top panel shows how azimuthal variations around the black dashed circle track orbital dynamics near the BH, evident here as diagonal striations with sub-Keplerian angular velocity. The bottom panel shows how radial variations along the white dashed line will reveal the SMBH-jet connection and measure acceleration within the jet-launching region. Initial ngEHT monitoring campaigns (light blue vertical bands) will span 3 months per year with a dense (sub-week) observing cadence; for comparison, current EHT campaigns (dark vertical bands) only span <math display="inline"><semantics> <mrow> <mo>∼</mo> <mn>2</mn> <mspace width="0.166667em"/> <mi>weeks</mi> </mrow> </semantics></math> per year, which is insufficient to measure the dynamics of the accretion disk or jet.</p>
Full article ">Figure 11
<p>Example ngEHT reconstructions for Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> (top two rows) and M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> (bottom two rows), using submissions for the second ngEHT Analysis Challenge [<a href="#B34-galaxies-11-00061" class="html-bibr">34</a>]. For each source, upper panels show ground truth movie frames, and lower panels show example reconstructions. The M87<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> ground truth movie is a GRMHD simulation generated with H-AMR [<a href="#B245-galaxies-11-00061" class="html-bibr">245</a>] and ray-traced with <tt>ipole</tt> [<a href="#B246-galaxies-11-00061" class="html-bibr">246</a>]; the reconstructed movie was produced using <tt>resolve</tt> [<a href="#B225-galaxies-11-00061" class="html-bibr">225</a>]. The Sgr A<math display="inline"><semantics> <msup> <mrow/> <mo>*</mo> </msup> </semantics></math> simulation is a semi-analytic accretion flow with a shearing hot spot [<a href="#B97-galaxies-11-00061" class="html-bibr">97</a>,<a href="#B247-galaxies-11-00061" class="html-bibr">247</a>]; the reconstructed movie was produced using <tt>StarWarps</tt> [<a href="#B248-galaxies-11-00061" class="html-bibr">248</a>]. Panels are reproduced from Roelofs et al. [<a href="#B34-galaxies-11-00061" class="html-bibr">34</a>].</p>
Full article ">Figure 12
<p>Representative subset of the ngEHT Science Traceability Matrix (STM). Daggers (†) indicate threshold science goals. The STM is used to guide the array design and to inform decisions about the multi-phase deployment.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop