Astrobites RSS

J0416

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Magnifying the Early Episodes of Star Formation: Super-star clusters at Cosmological Distances
Authors: E. Vanzella et al.
First Author’s Institution: INAF–Osservatorio Astronomico di Bologna
Status: Submitted to ApJL, open access

Have another look at the cover image, which depicts the Hubble Frontier Fields of the galaxy cluster MACS J0416. It always amazes me to see the manifestation of gravitational lensing in deep Hubble images — light from very-far-away galaxies being magnified and stretched into arcs by the strong gravity of the quite-far-away galaxy clusters. The gravity of the galaxy clusters acts as a “natural telescope” that focuses light to reveal background galaxies, which otherwise are too faint to be seen.

Astrophysicists have been puzzling over the mystery of reionization. How did reionization occur and what sources caused it? To try to answer these questions we need to know the origins and the properties of the early, far-away galaxies that were responsible. Recently, it was found that the huge number of faint galaxies may provide enough photons to reionize the universe. The technique of gravitational lensing comes in very handy because it allows far-away and faint objects to be observed!

Directly observing galaxies during reionization (with redshift z > 6) is hard. They are extremely faint and the strong characteristic spectral lines lie outside the limits of our detectors (no worries, JWST will come to rescue!) One way astrophysicists get around this problem is to study objects with slightly lower redshifts at z ~ 3, which are the younger analogs of the sources that reionized the universe. Today’s paper follows this approach.

Typically gravitational lensing reveals far-away galaxies. Today’s story is extraordinary: the authors managed to unravel two star clusters at redshift z = 3.2 by the lensing technique, and derived important hints about the ionization history of the universe.

Hubble insets

Figure 1. HST color image of the galaxy cluster MACS J0416 (middle section of the cover image.) The insets show the six images of the object ID14 (marked a to f). The annotated numbers are the magnitudes of each component. [Adapted from Vanzella et al. 2017]

OK, let’s get into the beautiful observations. Figure 1 shows the the middle section of the cover image, again centered at the galaxy cluster J0416. The insets (Image 1, 2, and 3) are the multiple images of the object ID14 generated by the gravity of J0416. Image 1 is further magnified into four additional images — ID14a, b, c, and d — by the elliptical galaxy pair E1 and E2 (ID14 is therefore called a doubly lensed system). Each ID14 image has two components, marked “1” and “2”.

MUSE spectra

Figure 2. Spectra of ID14 taken by the MUSE instrument of the Very Large Telescope. The colors of the spectra denote contributions from different images (black is sum of ID14a, b, and c; red is ID14b and c; blue is ID14a only). The main features are the strong metal lines and the weak Lyman-alpha line. [Adapted from Vanzella et al. 2017]

Spectra of ID14 are shown in Figure 2, with different colors representing contributions from different images. The magenta spectrum is taken from a Lyman-alpha emitting region ~2.1 kpc away from ID14 at the same redshift (magenta ellipse, Figure 3). There are two main points to take away: first, there are multiple strong high ionization lines (from highly ionized atoms He+, C2+, C3+, O2+) characteristic for energetic ionized environments; second, there is a weak Lyman-alpha emitting region not too far from ID14.

Figure 3. Image showing the Lyman-alpha emitting cloud (magenta ellipse) near ID14a, b, and c (black arc). The separation is estimated to be ~2 kpc. [Adapted from Vanzella et al. 2017]

By analyzing the images in Figure 1, the authors find that the source ID14 comprises two compact systems with sizes of ~30 pc each, separated by ~300 pc. Also, the line ratios measured from the spectra (Figure 2) are consistent with a stellar population. Further modeling of additional spectra gives a mass estimate of 106–107 solar masses. These suggest that ID14 consists of two ancient, compact, young, and massive star clusters — commonly referred to as super star clusters.

It is intriguing to see star clusters so far away. What’s more? ID14 also hints at the structure of ionizing radiation in the early universe! With the observed Hβ spectral line (not shown, see Figure 3 of the original paper), the Lyman-alpha line is predicted to be >150 times brighter than currently observed. Such deficiency can be explained by (1) dust absorption and (2) Lyman-alpha photons being scattered out of the observer’s line of sight by irregular distributions of gas. The authors proposed a plausible picture where ionizing radiation escapes the star clusters and hits the neutral cloud nearby, where we see the Lyman-alpha emission in fluorescence (Figure 3). This finding suggests that direction-dependent visibility of ionizing radiation observed on galactic scales could also prevail on the scale of star clusters.

We astrophysicists are cosmic detectives. By combining our advanced telescopes with the natural gravitational lenses, we are able to grasp information that would otherwise be out of reach. Today’s story highlights that reionization is really an incredibly complex problem, one that connects tiny star clusters to the scale of the cosmos. More and better observations will further constrain the properties of the ionizing sources and help us uncover the process of reionization!

About the author, Benny Tsang:

I am a graduate student at the University of Texas at Austin working with Prof. Milos Milosavljevic. Using Texas-sized supercomputers and computer simulations, I focus on understanding the effects of radiation from stars when massive star clusters are being assembled. When I am not staring at computer screens, you will find me running around Austin, exploring this beautiful city.

volcano

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: A Volcanic Hydrogen Habitable Zone
Authors: Ramses Ramirez and Lisa Kaltenegger
First Author’s Institution: Cornell University
Status: Published in ApJL, open access

The search for life beyond the solar system has long focused on the habitable zone (HZ). This is the region around a star where a planet with the right properties could maintain liquid water on its surface for a substantial period of time. The classical inner edge of the HZ is set using the runaway greenhouse effect, in which a positive feedback loop causes oceans to evaporate, creating an oven-like world similar to Venus. The classical outer edge of the HZ is set using the maximum greenhouse effect from carbon dioxide, which is the distance at which adding carbon dioxide to a planet’s atmosphere starts cooling the planet (due to scattering the light or condensation). There have been many other calculations of the HZ edges using different assumptions, such as a nearly desert planet and planets with different masses. In this paper, the authors try to use volcanoes to expand the edges of the HZ. They calculate the HZ edges for atmospheres with significant amounts of hydrogen gas produced by volcanoes, another powerful greenhouse gas.

Hydrogen-Induced Greenhouse Warming

An atmosphere with significant greenhouse warming due to hydrogen is difficult to maintain, because hydrogen gas escapes atmospheres quickly. However, early in the Earth and Martian geological histories, volcanic outgassing of hydrogen may have exceeded atmospheric escape of hydrogen. This means that both planets might have had significant amounts of hydrogen in their atmosphere. Different conditions in the mantle could make this hydrogen outgassing occur over a much longer timescale and therefore give the planet a longer hydrogen-induced greenhouse effect. Using this volcanic outgassing as their source of hydrogen, the authors used a 1D atmospheric climate model to compute the edges of the HZ for an atmosphere composed of nitrogen, water vapor, carbon dioxide, and hydrogen for stars with temperatures between 2,600 K and 10,000 K. A variety of atmospheric concentrations were tested up to 50% hydrogen (30% hydrogen is the highest concentration they could reasonably acquire by assuming different geologies, but 50% was included as an extreme outlier). Because these models depended on so many variables, many assumptions were necessary in the model too, such as plate tectonics, the carbon-silicate cycle, an oxygen-reduced (i.e., oxygen-poor) mantle, and a constant albedo.

New Habitable Zone Results

Adding hydrogen into planetary atmospheres moved both the inner and outer edges of the HZ outward. The outer edge moved farther than the inner edge, which widened the HZ. The incident stellar flux (the amount of energy hitting the planet per second per square meter) needed to maintain liquid water on the surface at the outer edge of the HZ decreased by 25%, 44%, and 52% when the atmosphere was 5%, 30%, and 50% hydrogen, respectively. This moved the classical HZ edge from 1.67 AU to 1.94 AU, 2.23 AU, and 2.4 AU, respectively. The HZ expanded much more for the hotter stars than the cooler stars. The inner edge of the HZ, on the other hand, shifted only a tiny bit: 0.1% outward for 1% hydrogen and 4% outward for 50% hydrogen.

Figure 1: The outer edge of the habitable zone. The x-axis is the amount of energy received from the star per second per area, and the y-axis is the temperature of the star. The stellar temperature is important because it changes the distribution of energy hitting the planet (e.g., a higher proportion of the incident energy is in the infrared as the stellar temperature decreases). The dashed line is the classical outer edge of the HZ. The solid line is the empirical outer edge using evidence suggesting that early Mars had liquid oceans. The red lines are the outer edges of the HZ for atmospheres with different concentrations of hydrogen. For reference, the Sun’s effective temperature is 5,780 K (where Mars is).

Conclusions

It is expected that terrestrial planets are born with oxygen-reduced mantles. A planet with a reduced mantle is more likely to have an extended period of hydrogen outgassing and therefore a longer hydrogen-induced greenhouse effect. Over time, however, the mantles become oxidized. Some research has suggested that smaller planets’ mantles (like Mars’s) stay reduced, while larger planets’ mantles oxidize quickly. This suggests that hydrogen outgassing might only be relevant for smaller planets. On the other hand, more massive planets can hold onto hydrogen more easily due to their gravity and higher likelihood to have a strong magnetic field. The relationship between planetary mass and the effectiveness of hydrogen outgassing on habitability remains unclear.

In our own solar system, Earth’s mantle may have become oxidized only about 100 million years after formation. Mars’s mantle, though, may have stayed reduced for a billion years. Two meteorites (called ALH84001 and NWA Black Beauty) from 4 billion years ago support this idea. Therefore, hydrogen outgassing from volcanoes could have contributed to a warm, wet, early Mars.

About the author, Joseph Schmitt:

I’m a 5th year graduate student at Yale University. My main research is on the discovery, characterization, and statistics of exoplanets. I’m also one of the science leads on the citizen science project Planet Hunters, a website where the general public can join the search for exoplanets.

Titan

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Compositional Similarities and Distinctions Between Titan’s Evaporitic Terrains
Authors: S.M. MacKenzie and Jason W. Barnes
First Author’s Institution: University of Idaho
Status: Published in ApJ, open access

Titan, Saturn’s largest moon, is the only solar system object other than the Earth to have a thick atmosphere and standing surface liquid. When the Cassini spacecraft began observing Titan, it even discovered lakes and seas dotting the northern hemisphere. Don’t fire up your rocket just yet, though — because Titan is so cold, the lakes and seas are filled with liquid methane and ethane rather than water.

Cassini

Artist’s illustration of the Cassini mission at Saturn. [NASA]

Titan’s thick, methane-rich atmosphere makes it difficult to observe the surface at visible wavelengths. Luckily, there are several windows in the near-infrared through which light can pass and reveal the surface. Seven of these windows overlap with the wavelength range covered by Cassini’s Visual and Infrared Mapping Spectrometer (VIMS). By looking at how the brightness of the surface changes with wavelength, we can learn about the composition of the surface material. The cover photo above depicts a three-color map of Titan’s surface made with VIMS. The pinkish regions show where the surface reflects strongly at 5 microns.

The 5-micron-bright regions are found circling lakes in the northern hemisphere, in dry lake beds in both hemispheres, and in the desert-like equatorial regions. The bright rings around the lakes are believed to be evaporites—solid material left behind after the liquid in which it was dissolved evaporates. This explains the presence of the bright material surrounding the lakes and the dry lake beds, but what about the desert? Linking 5-micron-bright regions in what is today a desert to the bright rings around the lakes could provide evidence that the equatorial regions of Titan were once covered with liquid.

In this paper, the authors searched for a compositional link between the bright regions in the desert and the evaporites around the lakes and seas. They used an absorption feature at 4.92 microns in order to investigate whether or not the 5-micron-bright material in each of these regions is the same. The 4.92-micron absorption feature has been observed previously in the desert region, but no one has been able to definitively say what compound causes it. Because of this, finding the same feature in the desert and around the lakes can indicate that the regions are geologically similar, but can’t yet tell us about the chemical makeup of the material.

VIMS data

A non-projected version of the VIMS map of Titan shown above. [JPL/NASA/Univ. of Arizona/CNRS/LPGNantes]

The authors used Principal Component Analysis (PCA) to isolate the weak 4.92-micron absorption feature. PCA is a mathematical method that separates the individual components that make up an observed signal. In this case, the main contributors to the signal (i.e. the “principal components”) could be changes in the surface reflectivity, instrumental noise, or compositional variations. Once the components have been separated, the unwanted contributors can be removed. As a result, PCA can be used to isolate a signal that is much smaller than the background noise. (PCA is also used in the direct detection of exoplanets and is described in more detail here.) After applying PCA, the authors observed the 4.92-micron absorption feature in both the desert and around the lakes, strengthening the hypothesis that the desert once had liquid. However, they also found that not all of the lake regions had the absorption feature, and some of the regions that did have it didn’t have it in every observation. They suggested that material with a crystalline structure that reflects light more strongly at some angles or transient effects like methane rain could cause the absorption feature to appear intermittently.

What causes some lake regions to have the absorption feature while others don’t? The authors suggested that the material that causes the 4.92-micron absorption feature could be just one of several solids that are left behind as the lakes evaporate away. Whether or not a lake rim has the absorption feature could be a function of how far the evaporation has progressed. As evaporation proceeds, materials that are more soluble precipitate out in sequence. We could see a 5-micron-bright evaporite ring without the absorption feature if the lake hasn’t evaporated enough for the material causing the absorption to precipitate out. The authors even have a suggestion for why this might happen to some lakes in the northern hemisphere but not others—lakes closest to the north pole might experience more methane rainfall than more southern lagoons, periodically halting the evaporation before the absorbing material can crystallize.

Although the authors posit many explanations for the mysterious behavior of the 4.92-micron absorption feature, they can’t yet settle on one cause. It’s not surprising that Titan, an inhospitable but strangely familiar world with complex geology and weather systems, presents a challenge to astronomers. In the future, by modeling how Titan’s climate changes over time, we can hope to learn more about what causes the distribution of evaporites on Titan’s surface.

About the author, Kerrin Hensley:

I am a second year graduate student at Boston University, where I study the upper atmospheres and ionospheres of Venus and Mars. I’m especially interested in how the ionospheres of these planets change as the Sun proceeds through its solar activity cycle and what this can tell us about the ionospheres of planets around other stars. Outside of grad school, you can find me rock climbing, drawing, or exploring Boston.

hot Jupiter atmosphere

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Atmospheric circulation of hot Jupiters: dayside–nightside temperature differences
Authors: Thaddeus D. Komacek and Adam P. Showman
First Author’s Institution: University of Arizona
Status: Published in ApJ, open access

There is an old sci-fi movie “The Chronicles of Riddick” that’s set on a bizarre planet. One scene I still remember depicted the dawn on the planet. As the sun was rising (yes, somehow not tidally-locked…), the frozen surface from the night suddenly became so boiling hot that it burned anything into ash within seconds. Fortunately, on Earth, the day–night temperature difference is much milder — it’s usually less than 30°C, because the atmosphere mitigates the temperature and the Earth rotates quickly (like making evenly grilled chicken by turning it).

Hot Jupiters, on the other hand, are extreme worlds. They orbit very close to their host stars (< 0.1 AU) and are locked by the tidal force into synchronous rotation, with the same side always facing their stars. This makes for interesting atmospheric dynamics. In today’s astrobite, we take a look at these exotic worlds. The authors examined what essentially controls the day–night temperature differences and compare their theory to numerical simulations (so-called general circulation models or GCM).

The Hotter It Is, the Faster It Cools

Figure 1. Normalized day–night brightness temperature difference vs. equilibrium temperature from observation of transiting hot Jupiters.

Figure 1 summarizes the observed day–night temperature difference vs. the equilibrium temperature of several most-studied hot Jupiters. The value shown on the y-axis is the normalized day–night temperature difference. The equilibrium temperature shown on the x-axis essentially tells us how close the planets are to their stars. There is clearly a trend, and the reasoning can be understood via the Stefan–Boltzmann law: the hotter (more strongly irradiated) the planet, the faster it cools (radiates). The winds cannot carry hot gas to the planet’s nightside if the heat is radiated into space too quickly. This would lead to a larger day–night temperature difference. On the other hand, if the atmospheric circulation acts efficiently enough, the winds can smear out the temperature variations by redistributing heat across the planet. Another way to look at it is that hotter planets have shorter radiative timescales. Astrophysicists like to think in terms of “timescales” while dealing with different competing processes. The process with a shorter timescale occurs faster, and it therefore dominates over the others. The authors developed a simple analytical theory that shows the dependence of day–night temperature difference, as a function of several factors: day–night equilibrium temperature difference, the timescale of radiation, and the timescale of drag. These factors, together with planetary parameters, are the key input of GCM simulation and allow a direct comparison.

Scale Analysis

Scale analysis is a powerful tool for understanding which mechanisms play important roles. The authors use scale analysis to simplify the problem and construct the theory. For example, by looking at the continuity equation, we can readily understand that  U/L ~ W/H; i.e., the ratio of the horizontal wind speed (U) to the vertical wind speed (W) is about equal to that of the planet radius (L) to the pressure scale height (H). For a typical hot Jupiter, this ratio is ~100. In other words, the observable atmosphere is thin and the horizontal flow is much faster than the vertical motion.

Pen & Paper vs. Computer

Applying scale analysis to the set of primitive equations describing the atmosphere, the authors worked out analytical expressions for various quantities, like wind speed and day–night temperature difference. Before comparing with the GCM results, let’s briefly follow their analysis based on the governing law of the gas movement — the conservation of momentum. The gas motion has four contributing factors: gas flowing in or out (advection), planet’s rotation (Coriolis force), drag force (friction), and how the pressure changes spatially (pressure gradient). The pressure gradient is the driving force of the circulation, which has to be balanced by the others. Depending on whether drag, Coriolis force, or advection dominates, the atmosphere may be in different regimes and behave differently. When rotation is important compared to advection, the pressure gradient is balanced by rotation (weak-drag) or by drag force (strong-drag), but when rotation is negligible compared to advection (slow rotation or near the equator), the pressure gradient is balanced by advection (weak-drag) or by drag force (strong-drag).

Figure 2. Overall day–night temperature difference (colors) as a function of radiative timescale and drag timescale, showing a comparison of the GCM results (left) and theory (right). The black line in the right panels indicates the transition from drag dominance to rotation dominance.

The theoretical prediction is compared to GCM results in Figure 2, showing the day–night difference for various values of drag timescales and radiative timescales for the rotation-important regime. Firstly, the day–night difference decreases as radiative timescale increases, in accordance with the trend observed. Secondly, the day–night difference has no dependence on the drag timescale until the drag timescale is shorter than ~105 s (the black line in the righthand panels). This indicates the transition as discussed above: the weak-drag regime lies above the black line and the strong-drag regime under it. Lastly, a term that corresponds to the timescale of wave propagation appears in the analytical solutions in all regimes. These waves are similar to the waves in the tropics on Earth, which are strongly connected to the climate. The relative value of the wave timescale essentially controls the day–night temperature difference. It implies the wave propagation is crucial to mediate the day–night temperature difference (read more about this here).

In addition, the circulation and temperature structures in the GCM are displayed in Figure 3. Again, the radiative timescale predominantly determines the day–night temperature difference, until the drag becomes strong enough and the equatorial jet vanishes, exhibiting a symmetric circulation pattern.   

Figure 3: How temperature (colors) and wind (vectors) in GCM simulations vary with drag timescales (the y direction) and radiative timescales (the x direction).

This simple scaling theory helps us to understand the leading role of the dynamics and tells us a useful story about the physics behind it. It is always nice to see successful pen-and-paper solutions that help decipher the intricate results of computer simulations!

About the author, Shang-Min Tsai:

I am a 3rd year PhD student at the University of Bern and part of the Exoplanets & Exoclimes Group led by Prof. Kevin Heng. We are developing various open-source tools to study exoplanets. I work on modelling of atmospheric chemistry and dynamics. When I am not coding or debugging, I enjoy basketball and playing board games.

M67

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Galactic Doppelganger: The chemical similarity among field stars and among stars with a common birth origin
Authors: M. Ness, H-W. Rix, David W. Hogg et al.
First Author’s Institution: Max Planck Institute for Astronomy, Germany
Status: Submitted to ApJ, open access

Stars are not born alone. Stellar formation occurs inside giant molecular clouds, which collapse when the internal gas pressure cannot overcome gravity. This is known as Jeans instability. The increase in density caused by the collapse leads to a fragmentation of the cloud, forming thus not only one, but many stars from the same original material. Consequently stars that are born together have the same initial composition. In addiction, these sibling stars are subjected to each others’ gravitational pull, so the formation process usually results in a stellar cluster rather than isolated stars.

Some stars can be ejected from the cluster in many-body interactions. Moreover, most clusters don’t survive long, since they are disrupted due to gravitational interaction with the surroundings. This implies that siblings are not always found together. They still preserve, though, one common birthmark: their similar chemical composition. Thus one can in principle identify population members by comparing stellar abundances. This is known as chemical tagging, and it is widely used to characterize the components of the Milky Way and constrain its assembly history (see this bite and this bite).

However, there are a few considerations to take into account. First, because of differences in the evolutionary process, some spread in abundances is expected among sibling stars. Second, there’s always the chance that two stars that are not born together have the same current abundances by pure chance — so-called doppelgangers. These two difficulties were usually ignored in chemical tagging, since no estimate of their magnitude existed. Today’s paper aims to put an end to this limitation, using precise abundance determinations to estimate intra-cluster dispersion and doppelganger rates in the Milky Way’s thin disk.

Data and abundance estimates

To obtain a large sample of stars and achieve high-precision abundance measurements, the authors used data from APOGEE, a project that is part of the Sloan Digital Sky Survey. APOGEE takes spectra of mainly giant stars and allows abundance estimates for tens of chemical elements. APOGEE’s pipeline estimates each star’s physical parameters, effective temperature (Teff) and surface gravity (log(g)), and also element abundances. The authors’ first step was to correct the abundances for systematic uncertainties (such as those arising from dependence on stellar physical parameters) with a training sample containing thousands of field stars.

These corrections were applied to test data, consisting of spectra of about a hundred stars in seven well-studied open clusters. Fig. 1 shows the estimated abundances for Messier 67 (shown at the top of the page) as an example. As expected, there isn’t much spread in abundances for the stars that are part of the cluster. Using the same method, abundances were estimated for other 90 open clusters. The obtained uncertainties were typically 20% – 50% lower than obtained in other determinations.

Figure 1: Abundance estimate for 20 elements in M67 stars, coloured by effective temperature. The grey points are the training data. The values at the top of each subplot are mean and standard deviation of the estimates. All elements are measures with respect to Fe, except [Fe/H].

Figure 1: Abundance estimate for 20 elements in M67 stars, coloured by effective temperature. The grey points are the training data. The values at the top of each subplot are mean and standard deviation of the estimates. All elements are measures with respect to Fe, except [Fe/H]. [Ness et al. 2017]

Siblings or doppelgangers?

The results point out that the typical dispersion of abundances among sibling stars is comparable to the measurement uncertainties. This implies that clusters are unsurprisingly nearly homogeneous in their abundances. So there is no need to worry about intrinsic abundance dispersion in chemical tagging analysis, at least in the case of large-scale surveys.

Knowing that this dispersion is negligible for the APOGEE’s uncertainties, the authors next estimated the similarity between stars in the same cluster, and between stars in the field. Do stars within a cluster always look alike, in terms of abundance? Do stars in the field always look different? The answer to both these questions is no. The authors showed that by estimating the distribution of abundance differences, through computation of the chi-square for the measured abundances for each pair of stars, both intra-cluster and in the field. This was first done considering stars with similar Teff and log(g), and then also pairs with similar [Fe/H], a measure of metallicity.

As can be seen in Fig. 2, the distributions for intra-cluster pairs and field pairs are very different. The former has a median very close to the number of degrees of freedom, as is expected for a chi-square distribution when a good fit is obtained, confirming that most siblings are a good fit to one another, and thus very similar. There exists, however, some spread, indicating that siblings can be significantly different. The distribution for the field pairs, on the other hand, has a much larger spread and much larger values. That reflects the fact that field stars are usually not similar. However, about 1% of field pairs are as similar as siblings. The estimated probability that two stars chosen at random were born together is, according to the authors, much smaller than that (about 0.003%). Thus a significant fraction of these pairs is not in fact siblings, but only doppelgangers: very similar, but unrelated.

Figure 2: Top plot shows the chi-square distribution of abundance differences for pairs with similar Teff and log(g). The black histogram represents intra-cluster pairs, and the red-dashed shows field pairs. The distributions are very different, but not disjoint: some field pairs are as similar as siblings. Bottom panel shows analogous estimates, but the [Fe/H] is set to the similar value of the clusters M67 and NGC6819. Their intra-cluster distributions are shown in the blue dash-dot histogram for comparison. Even for similar metallicity, most field stars are distinguishable, but there's still a considerable number of doppelgangers.

Figure 2: Top panel shows the chi-square distribution of abundance differences for pairs with similar Teff and log(g). The black histogram represents intra-cluster pairs, and the red-dashed shows field pairs. The distributions are very different, but not disjoint: some field pairs are as similar as siblings. In the bottom panel pairs have also similar [Fe/H], close to the value of the clusters M67 and NGC6819. Their intra-cluster distributions are shown in the blue dash-dot histogram for comparison. Even for similar metallicity, most field stars are distinguishable, but there’s still a considerable number of doppelgangers. [Ness et al. 2017]

Implications for chemical tagging

Chemical tagging alone shouldn’t be used to characterize a population; certainly not in the Milky Way disk, even with the high data quality achieved by the authors. The relatively high doppelganger rate found by the authors, combined with the fact that siblings are not always chemically alike, makes chemical tagging uncertain. Other parameters, such as velocity information, should be combined with the abundance analysis. Keep that in mind!

About the author, Ingrid Pelisoli:

I am a second year PhD student at Universidade Federal do Rio Grande do Sul, in Brazil. I study white dwarf stars and (try to) use what we learn about them to understand more about the structure and evolution of our Galaxy. When I am not sciencing, I like to binge-watch sci-fi and fantasy series, eat pizza, and drink beer.

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Modeling the Historical Flux of Planetary Impactors
Authors: David Nesvorny, Fernando Roig, William Bottke
First Author’s Institution: Southwest Research Institute
Status: Accepted to AJ, open access

Asteroids have a bad reputation. They may have wiped out the dinosaurs, and they have threatened the survival of humanity in many terrible movies.

Until today’s featured paper, asteroids were also blamed for the Late Heavy Bombardment (LHB, for short) — a period of time in the early solar system when the Moon, the Earth, and the other rocky planets were hit with an unusually high number of impactors from a variety of material in space. It lasted from the Sun’s formation 4.6 billion years (Gyr) ago until the impact that created the Orientale crater on the Moon 800 million years later (3.8 Gyr ago). Planetary scientists working on the Apollo missions in the 1960s first hypothesized the LHB when astronauts brought back impact melt rocks from various craters that unexpectedly all dated back to this time. The idea was later extended to include the rest of the inner solar system when planetary scientists found similar cratering histories on each of the rocky planets. These LHB-era impacts are believed to have been caused by some combination of asteroids (from the asteroid belt), comets (from the Kuiper Belt and beyond), and leftover material from the formation of the inner rocky planets. However, it is still not clear which of these three sources was the main culprit.

To address one part of this issue, Nesvorny et al. — the authors of today’s paper — ask: Were there enough impacting asteroids to account for all of the craters on the Moon from the Late Heavy Bombardment?

How Did Asteroids End Up Impacting the Moon?

The leading idea for the source of asteroid impactors during the LHB is that they did not come from the main asteroid belt that exists today. Instead, they originated in what used to be the inner part of the belt — called the “E-belt” — that spanned from 1.7 to 2.1 AU, but became extinct (“E” is for Extinct!) after interacting with the planets in our solar system before they settled into their current orbits.

During this time period, the planets were not located where you think of them today. In the Nice model, Jupiter started out at 5.35 AU (compared to 5.2 AU today) and Saturn at 8.40 AU (compared to 9.5 AU today). Then, as Jupiter migrated inward and Saturn moved outward — both by about 0.04 AU — they moved from almost being in resonance to exactly in resonance. At this point, Saturn began to orbit around the Sun two times for every one time that Jupiter completed an orbit. This perfect alignment created complete chaos! Jupiter and Saturn quickly migrated towards their current locations, and in the process, they flung Uranus and Neptune much further away from the Sun, while also ejecting asteroids out of the asteroid belt. Many of these asteroids fled for the inner solar system, where they could then impact the Moon, the Earth, and the other inner planets.

Nice model

Figure 1. Nice Model: Evolution of the early solar system. When Jupiter (J) and Saturn (S) reached a 2:1 resonance (dotted line), they ejected asteroids out of the asteroid belt and pushed Uranus and Neptune away. There is a decent amount of evidence that something like this occurred (albeit, there are competing models), although it is unknown whether planet migration during the Nice model is directly responsible for the LHB. [Adapted from Rivkin et al. 2010]

Counting Asteroid Impacts in the Nice Model

Nesvorny et al. test whether asteroids could be responsible for the LHB by simulating the orbital evolution of a full belt of asteroids that stretches from 1.7 to 3.3 AU with a wide range of low eccentricities (0.0 to 0.4) and low inclinations (0 to 20 degrees). Concurrently, they also integrate the orbits of all of the planets (except Mercury to save computation time) according to how they evolve in the Nice model. They run this simulation from not too long before Jupiter and Saturn fall into resonance to the present day.

At the end of the simulation, they calibrate the fraction of asteroids that impact the Moon into an actual number of impacts to expect in the real solar system by comparing the fraction of surviving asteroids in the entire belt (in the simulation) to the number of asteroids in the real asteroid belt. Better yet, the authors count only the real asteroids above a certain size (e.g. 10 km, and 130 km) to figure out whether there are enough large impactors in the simulation to account for all of the large craters on the actual Moon.

As Figure 2 shows, Nesvorny et al. find that there are not enough large asteroids to explain the number of known impact craters on the Moon. Thus, asteroids from the asteroid belt could not have been responsible for the majority of the impacts during the LHB.

Number of impacts

Figure 2. Total number of impacts on the Moon above a given size (left: 50 km, 130km; right: 10 km, 20 km) for three simulations. Left: The Imbrium crater was created by a 130-km impactor, but not even one asteroid that size hit the Moon throughout its history in the simulations. Right: There are 200 craters with a 150-km diameter or more on the Moon, each of which was created by a 10-km impactor. However, fewer than 20 asteroids that size hit the Moon. [Nesvorny et al. 2017]

Other Ideas

Recent work has also explored whether comets could explain the LHB-era impacts on the Moon, but found that the rocks returned from the Apollo missions do not have the right chemical make-up — specifically, oxygen isotope ratios — to match up with comets. This leaves leftover material from when the Earth, Mars, and other inner planets formed as the most likely remaining explanation for the impacts during the Late Heavy Bombardment.

The lead author, David Nesvorny, mentions that he has submitted another paper (with Alessandro Morbidelli and several other planetary scientists) that validates the rocky planet leftovers as the main cause of the Late Heavy Bombardment. Hopefully when that paper is published, I can write another Astrobite that more conclusively finishes this important part of the story.

About the author, Michael Hammer:

I am a 2nd-year graduate student at the University of Arizona, where I am working with Kaitlin Kratter on studying planetary dynamics and planet-disk interactions through numerical simulations. I am from Queens, NYC.

Tibet Air Shower Array

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Northern sky Galactic Cosmic Ray anisotropy between 10-1000 TeV with the Tibet Air Shower Array
Authors: The Tibet AS Collaboration
First Author’s Institution: Yale University
Status: Accepted to ApJ, open access

One of the big unsolved mysteries in particle astrophysics is the origin and acceleration mechanisms of cosmic rays, or charged particles that are constantly bombarding the Earth. It is impossible to tell which astronomical source is the origin of any one particular cosmic ray. Because cosmic rays are charged particles, they are deflected on their way to Earth by our galactic magnetic field. Due to this effect, one might expect cosmic rays to arrive isotropically (in equal numbers in every direction), with slight anisotropy because of diffusion effects.  However, that is not what is observed. Instead, there is a large-scale anisotropy with an energy-dependent amplitude. This causes problems for traditional diffusion models. Studying this anisotropy is important in learning more about cosmic rays.

Today’s bite uses roughly five years of data from the Tibet Air Shower experiment, located at 4300 meters above sea level in Tibet, to provide an update in the study of the cosmic-ray anisotropy. The area of the sky they looked at is slightly larger than in their previous papers, and when combined with data from the IceCube experiment in Antarctica, allows us to have a complete picture of the entire sky in the energy range of a few hundred TeV.

New Results

new results

Left: The anisotropy, as seen in 5 different energy bins. From top to bottom, the median energy is 15, 50, 100, 300 and 1000 TeV. Right: The 1D projection of the plot on the left. The blue curve is the first harmonic fit to the data. Note how the phase changes with energy. [Amenomori et al. 2017]

Maps were binned according to energy and analyzed. A summary of results can be found below.

  • In the map corresponding to a median energy of ~300 TeV, there are two regions that can be seen by eye: one area where there is an excess of cosmic rays, and another area where there is a deficit. Only the excess is statistically significant after taking trials into account. (For a description of what “trials” are in a statistical context, check out this page about the “look-elsewhere effect“). This result is consistent with what IceCube sees in a similar energy range. The authors looked to see if this anisotropy is due to the Compton-Getting effect (i.e. if it was caused by the motion of our solar system around the center of our Galaxy). The Compton-Getting effect predicts that we would see more cosmic rays coming from the direction that the Earth is moving toward, with a deficit coming from the other direction. They concluded that this particular excess is not related to the Compton-Getting effect.
  • The anisotropy appears to be energy dependent. The maps in the 15–50 TeV range have completely different features than the maps at a few hundred TeV/1 PeV. (See the figure above for an illustration)
  • These effects are not related to seasonal variations in the performance of the detector (atmospheric effects can affect experiments such as the Tibet Air Shower array, but it was shown to be negligible in this case).

So what are the implications of this? Well, for one, studies like these give us some hints as to the origins of cosmic rays. For example, at a few hundred TeV, the most significant excess is coming from the direction of the Galactic center. However, the observed energy dependence shows that there is still a lot to learn about how cosmic rays diffuse. Changes in the anisotropy with energy may imply that the cosmic ray propagation parameters are evolving.

About the author, Kelly Malone:

I am a fourth year physics graduate student at Penn State University studying gamma-ray astrophysics. I previously received bachelor’s degrees in physics and astronomy from UMass Amherst in 2013.

elliptical galaxy

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we repost astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Morphology and the Color-Mass Diagram as Clues to Galaxy Evolution at z~1
Authors: Meredith C. Powell, C. Megan Urry, Carolin N. Cardamone, et al.
First Author’s Institution: Yale University
Status: Published in ApJ, open access

Introduction

In astronomy, the relationship between the color and the brightness of objects has been used as a basic classification tool for decades. For stars, color is a measure of temperature; the bluer a star, the hotter it is. Measuring the relationship between observed colors and brightnesses of stars helps astronomers figure out the evolutionary sequences of stars.

It turns out that you can do almost the exact same thing for galaxies. Except this time, the color of a galaxy is a measure of how quickly it produces stars, called its star formation rate. This is because hot (and therefore blue) stars are extremely bright, despite being relatively rare, so they can make a whole galaxy appear blue. And such hot stars don’t live for long, so if a galaxy has these stars it must have formed them very recently. On the other hand, the brightness of a galaxy is roughly correlated with its mass, since a more massive galaxy will generally have more stars and therefore be brighter than a less-massive galaxy.

The color-mass diagram

In the same way we use color-magnitude diagrams to determine how stars evolve, we can use a plot of galaxy color (star formation rate) against mass (brightness) to figure out how galaxies evolve.

Figure 1. A cartoon of a low-redshift color-mass diagram. On the y-axis, “number of stars forming” is the star formation rate, which is correlated with color; on the x-axis, “number of existing stars” is the mass, which is correlated with brightness. [CANDELS Collaboration]

Figure 1 (which is a cartoon that I’ve shamelessly borrowed from this previous Astrobite) shows a few basic trends that emerge from this color-mass plot. There’s a “main sequence” of blue star-forming galaxies, where most galaxies lie. It’s thought that galaxies move from this main sequence to the “red and dead cloud,” made up of red galaxies that aren’t actively forming many stars. Not many galaxies are located in the region between these two areas, sometimes called the “green valley,” suggesting that galaxies don’t stay there for a very long time — so the processes that “quench” (turn off) star formation in blue galaxies must be relatively fast.

These basic trends appear on the color-mass diagram in the nearby universe at z~0 (z is what astronomers call “redshift,” which is both a measure of distance and of the universe’s age). Today’s paper aims to understand the trends on the color-mass diagram at higher redshifts (z~1), when the universe was younger — less than half its current age!

Today’s paper

More specifically, the authors want to answer a few questions:

  1. Do overall trends in the color-mass diagram look different at high-redshift (z~1) than at low-redshift (z~0)?
  2. Do elliptical (what the authors call “bulge-shaped”) and spiral (“disk-shaped”) galaxies follow the same trends on the color-mass diagram?

If the answer to question (1) is “yes,” then we know that the processes that shut off star formation have changed from redshift z~1 to now. And because most bulge-shaped galaxies are thought to be the results of galaxy mergers (although this is up for some debate), answering question (2) would help us figure out how mergers affect galaxy evolution.

To address these questions, the authors first used Hubble Space Telescope images (specifically, the GOODS and CANDELS surveys) to measure the colors of about 5,000 galaxies. Next, they matched the galaxy spectral energy distributions to theoretical models in order to infer their masses. Finally, they fit the galaxy images to different profile shapes to distinguish between “bulge” and “disk” galaxies.

Armed with this data, let’s see if the questions above were answered.

screen-shot-2017-01-22-at-11-38-52-pm

Figure 2. The color-mass diagram at z~1. Note that the y-axis is reversed from Figure 1, so that bluer (higher SFR) galaxies are lower and redder galaxies are higher. This means that the “main sequence” is marked by the contours in the bottom left, while the “red and dead cloud” is marked by the contours in the upper right. Purple points are active galactic nuclei (AGNs). [Powell et al. 2017]

  1. Do overall color-mass trends change with redshift?
    Figure 2 shows the color-mass diagram at z~1. Like Figure 1 (which is a cartoon of the z~0 color-mass diagram), it shows that the “main sequence” and “red and dead cloud” are mostly distinct. The main difference is that there are fewer “red and dead” galaxies than at z~0. This is probably because there hasn’t been as much time since the peak of star formation in the universe (which was at redshifts z~1–3) for galaxies to move away from the main sequence.
  2. Do bulges and disks follow the same color-mass trends?
    Disk galaxies follow pretty much the same trend as the gray points in Figure 2, although the “main sequence” and “red and dead cloud” are not as distinct. In contrast, bulge galaxies (see Figure 3) are mostly clustered in the “red and dead cloud,” with very few galaxies between the cloud and the “main sequence.” This means that disks probably evolve gradually, but bulge galaxies (which are the results of mergers) undergo a much more abrupt quenching process. (Note that this agrees with the results reported in this Astrobite for low redshifts!)

Figure 3. Color-mass diagrams for bulge galaxies. The galaxies are colored by specific star formation rate, which is just star formation rate divided by mass. [Powell et al. 2017]

So we’ve answered all our questions! But the authors note that there are still some unclear results. For instance, they also compare the color-mass trends of active galactic nuclei (AGNs; the purple points in Figure 2) with the trends of inactive galaxies. This is an important question in galaxy evolution, since it’s not clear whether feedback from AGNs helps to trigger or quench star formation. The authors find that AGNs might have some effect on quenching, particularly in mergers, but more data is needed to be certain.

The color-mass diagram has been around for a while — but as today’s paper shows, it’s still an incredibly useful diagnostic tool that can teach us a thing or two about galaxy evolution.

About the author, Mia de los Reyes:

I’m a postgrad student at the University of Cambridge, where I work with Rob Kennicutt to study star formation in galaxies (mostly I make lots of graphs). I did my undergrad at NC State University, which is where I started making graphs. Now, whenever I’m not making graphs, I can usually be found rock climbing or eating baked goods.

NGC 2392

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we will be reposting astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Masses of the Planetary Nebula Central Stars in the Galactic Globular Cluster System from HST Imaging and Spectroscopy
Authors: George H. Jacoby, Orsola De Marco, James Davies, I. Lotarevich, Howard E. Bond, J. Patrick Harrington, Thierry Lanz
First Author’s Institution: Lowell Observatory
Status: Accepted to ApJ, open access

Have you ever looked at something and wondered, “How did that get there?!” Has that something ever been a planetary nebula? Astronomers are scratching their heads over four planetary nebulae that have turned up in the unlikeliest of locations: globular clusters.

We know of thousands of planetary nebulae in the Milky Way and can even study planetary nebulae in other galaxies. So the first question is: why shouldn’t we see them in globular clusters? To understand this, we need to know a little bit more about how planetary nebulae are formed.

Planetary nebulae (PNe, singular PN) form during a fleeting phase in the life cycle of low- and intermediate-mass stars (a good range to remember is 0.8 to 8 times the mass of the Sun; our Sun is likely to become a PN in ~5 billion years). PNe consist of a central white dwarf wreathed in shells of hot diffuse gas ejected during the asymptotic giant branch phase late in the star’s evolution. The gas is ionized by ultraviolet photons from the central star, giving rise to some of the most colorful and beautiful objects in the Universe (see Figure 1). The term “planetary nebula” is an unfortunate misnomer; an overzealous William Herschel, fresh from his discovery of Uranus in 1781, thought the fuzzy greenish object he spied through a telescope bore a striking resemblance to his recent planetary find.

m57_nasa_hst

Figure 1. M57, otherwise known as the Ring Nebula, is an elliptical planetary nebula that can be observed even with a small telescope. [NASA/HST]

In order to create a visible PN, the central star must become hot enough to produce ultraviolet photons to ionize the nebula before it drifts away from the central star and dissipates. Stars with masses less than ~0.8 solar masses are thought not to form visible PNe because by the time the central star becomes hot enough to ionize the nebula, the nebular material has already dispersed. Given the age of the globular clusters investigated in the work (11.6 – 13.2 billion years) and stellar evolution models, stars in these clusters with masses of ~0.8 solar masses should be departing the main sequence. In other words, the turn-off mass is approximately equal to the lower limit for PN formation. If the stars departing the main sequence in globular clusters aren’t massive enough to form PNe, how did the PNe we observe come to be?

The leading theories for how we might be able to see PNe in globular clusters involve interactions between two stars. First, the observed PNe might be the descendants of blue stragglers — main sequence stars with masses higher than the cluster turn-off mass which are thought to be the result of two stars merging. The resultant PN would appear to be the evolutionary product of a single massive star, but the PN central star would be more massive than a typical white dwarf in the cluster. Alternatively, the PNe could result from post-common-envelope binaries. (You can learn more about post-common-envelope binaries and how they relate to PNe here.) The common envelope accelerates the star’s transition between asymptotic giant branch star and white dwarf. As a result, the PN central star can have a mass equal to or less than the typical cluster white dwarf mass.

In this paper, the authors analyzed Hubble Space Telescope observations of the four known globular cluster PNe in order to determine the most likely formation scenario for these objects. Optical images of the four target objects can be seen in Figure 2.

pne

Figure 2. The four PNe investigated in this paper, as seen by the Hubble Space Telescope. From left to right: Ps 1, IRAS 18333, JaFu 1, JaFu 2. [Jacoby et al. 2017]

The authors used stellar evolutionary tracks — models of how the temperature and luminosity evolve after a star leaves the main sequence — to determine the masses of the central stars. Combining the derived central star masses with secondary information such as the morphology of the individual nebulae, they conclude:

  1. Two PNe (first and third in Figure 2) most likely resulted from a merger or mass transfer. However, the masses cannot yet be determined precisely enough to distinguish between the formation scenarios described above.
  2. One PN (far right in Figure 2) shows only weak evidence for a binary interaction. If PNe in globular clusters arise from single stars, it would require a re-evaluation of established evolutionary timescales.
  3. The last object (second from the left in Figure 2) is so bizarre that the authors questioned its membership in the PN class altogether!

With only half of the known four globular cluster PNe requiring some form of binary interaction, we can’t yet invoke binaries as the cause of globular cluster PNe. Despite this, it’s still important to understand the role that binary and multiple star systems play in PN formation because so many stars in the Galaxy fall into the PN progenitor mass range. Do binary systems help shape many PN? Or can single stars put on a spectacular end-of-the-line display without assistance from a companion? Earthlings, five billion years from now, will be waiting to find out.

About the author, Kerrin Hensley:

I am a second year graduate student at Boston University, where I study the upper atmospheres and ionospheres of Venus and Mars. I’m especially interested in how the ionospheres of these planets change as the Sun proceeds through its solar activity cycle and what this can tell us about the ionospheres of planets around other stars. Outside of grad school, you can find me rock climbing, drawing, or exploring Boston.

pulsar bow shock

Editor’s note: Astrobites is a graduate-student-run organization that digests astrophysical literature for undergraduate students. As part of the partnership between the AAS and astrobites, we will be reposting astrobites content here at AAS Nova once a week. We hope you enjoy this post from astrobites; the original can be viewed at astrobites.org!

Title: Hubble Space Telescope detection of the millisecond pulsar J2124–3358 and its far-ultraviolet bow shock nebula
Authors: B. Rangelov, G. G. Pavlov, O. Kargaltsev, A. Reisenegger, S. Guillot, M. van Kerkwijk & C. Reyes
First Author’s Institution: George Washington University
Status: Accepted to ApJ, open access

Pulsars — the rapidly rotating, highly magnetized neutron stars that beam radiation from their magnetic axes — are as mysterious as they are exotic. They’re most often observed at radio frequencies using single-dish telescopes, and they’re sometimes glimpsed in X-ray and gamma-ray bands. Far rarer are pulsar observations at “in-between” frequencies, such as ultraviolet (UV), optical, and infrared (IR) (collectively, UVOIR); in fact, only about a dozen pulsars have been detected this way. However, their study in this frequency range has proved enlightening, as we will see in today’s post.

A pulsar too hot to handle

While one would expect a neutron star to cool with age if an internal heating mechanism does not operate throughout its lifetime, observations of the millisecond pulsar J0437–4715 (an interesting object in its own right) yielded surprising results. In a 2016 study, far-UV observations revealed the 7-billion-year-old pulsar to have a surface temperature of about 2 × 105 K — about 35 times the temperature of the Sun’s photosphere. This finding inspired Rangelov et al. to observe another millisecond pulsar, J2124–3358 (a 3.8-billion-year-old pulsar with a spin period of 4.93 ms), in the far-UV and optical bands using the Hubble Space Telescope (HST).

Because so few pulsars have been studied in these frequency ranges, their spectral energy distributions (SEDs) in this regime are poorly understood. Generally speaking, the spectra of normal, rotation-powered pulsars reveal a nonthermal (not dependent on temperature) component in optical and X-rays caused by electrons and positrons in the pulsar magnetosphere. In the far-UV, some pulsars show a thermal (blackbody) component in their spectra, thought to come from the surface of the cooling object. Analysis of the team’s HST images revealed an SED that is best modeled by a combined nonthermal and thermal spectral fit, with nonthermal emission dominating at optical wavelengths and thermal emission appearing in the far-UV (see Figure 1). If their interpretation is correct, this implies a surface temperature for J2124–3358 that is between 0.5 × 105 and 2.1 × 105 K, which is very much in line with the temperature of J0437–4715. If this proves to be the case, these two measurements will strongly suggest the presence of a heating mechanism in millisecond pulsars. However, various fits using only nonthermal components in the far-UV are still valid, so it is impossible to make an absolute determination of the correct fit.

There are quite a few heating mechanisms that could be invoked to explain these objects’ high temperatures, ranging from the release of stored strain energy from the pulsar’s crust to dark matter annihilation in the pulsar’s interior. More spectral coverage of J2124–3358 is necessary to both check the validity of the nonthermal and thermal combined fit and to get closer to determining more specifically the heating mechanism in play.

Figure 1: thermal and nonthermal combined fit to HST far-UV/optical data for J2124

Figure 1: Thermal (red dashed) and nonthermal (blue dashed) combined spectral fit to HST far-UV/optical data for J2124–3358. The black line signifies the sum of both components. Because there is uncertainty about the nature of the nonthermal component, two possible spectral slopes are shown. [Rangelov et al. 2016]

A (bow) shocking find in the far-UV

Images of J2124–3358 also show the presence of a bow shock, which is an arc-shaped shock that occurs when an object is moving faster than the interstellar medium (ISM) sound speed. J2124–3358 was known before this study to be accompanied by such a shock in H-alpha (Hydrogen transition from n=3 to n=2) filters, for which plenty of neutral hydrogen is required. As a result of the HST observations, J2124–3358 was found to have an (albeit fainter) far-UV shock coincident with the H-alpha shock (see Figure 2). This is only the second such object (after J0437–4715) to show a far-UV bow shock. It is absolutely possible that many pulsars cause bow shocks that don’t emit in H-alpha, but do in other wavelength regimes. Studying these more carefully will yield information about the nature of the ISM.

In order to learn more about the heating mechanisms operating in these objects as well as the bow shocks that sometimes accompany them, many more pulsars will need to be studied using various optical, UV, and IR filters. Studies in the far-UV are only possible with Hubble, so it will be a long time before a sufficient number of objects will be studied at these frequencies in order to make solid conclusions about the nature of such interesting phenomena.

Figure 1: New observations from this study using the HST at three different wavelengths are shown in the top (left and right) and bottom left images. The shock is clearly visible in the far-UV using the F125LP filter. The bottom right image shows a previous H-alpha observation of the same pulsar. Figure 1 in the paper.Figure 2: New observations of J2124–3358 from this study using the HST at three different wavelengths are shown in the top (left and right) and bottom left images. The shock is clearly visible in the far-UV using the F125LP filter. The bottom right image shows a previous H-alpha observation of the same pulsar. [Rangelov et al. 2016]

About the author, Thankful Cromartie:

I am a graduate student at the University of Virginia and completed my B.S. in Physics at UNC-Chapel Hill. As a member of the NANOGrav collaboration, my research focuses on millisecond pulsars and how we can use them as precise tools for detecting nanohertz-frequency gravitational waves. Additionally, I use the world’s largest radio telescopes to search for new millisecond pulsars. Outside of research, I enjoy video games, exploring the mountains, traveling to music festivals, and yoga.

1 … 41 42 43 44 45 46