2023-09 September - 12 December Science Readings

Yeah, well for several months, I've done damned-all. Try and get back into harness. I've been thinning out the Arχiv paper listings, even if I couldn't bear to get to work on them. So let's see what remains in the pile.

2023-12-20

Over a week without doing anything. Try to get back into the habit. The HD110067 paper stuck it's nose out at me because I recognised the HD number. Which isn't normal, even for me. It's the system which was recently reported as having 6 planets at a "perfect spacing". Which is a rather problemtic description. Anyway, follow the link into the main record.

New Year

Get a few more done.

Contents

Articles studied this September to December - some of which might go to Slashdot.
Forming Asteroid Moons
MOND Count
Star Formation in the Cartwheel Galaxy
December HTML things
The "Cosmic Train Wreck" Galaxy, A520
The abundance discrepancy in ionized nebulae: which are the correct abundances?
6-planet systemis a multi-star system.
Dating of a Latin Astrolabe.
Early Galaxies initial Initial Mass Function
End of document

2023-09 September to 12 December Science Readings.

Revisiting Dimorphos formation: A pyramidal regime perspective and application to Dinkinesh’s satellite

https://arxiv.org/pdf/2311.07271.pdf

Last year NASA hit an asteroid moon with a spacecraft to test the efficiency of momentum transfer into asteroids. Being an (astronomically) small body orbiting in an (astronomically) weak gravitational field, the (astronomically) small amount of energy the NASA can afford to pump into such a system could produce detectable effects. The "mainstream media" didn't generally understand the point of the exercise, which was not to test if there would be an effect - everyone involved was sure there would be an effect - but to test how efficiently momentum was transferred to the body. That was a very open question, because astronomers were (and still are) pretty sure that many PHAs (Potentially Hazardous Asteroids) are "rubble piles" of many smaller bodies (from dust, up to mountains) held together by their own gravity, and were expected to be very weak. So - if hit by (say) the force of a megatonne nuke, it was very unclear if 90% of the energy would go into accelerating 99% of the asteroid's mass, or if that 90% would go into accelerating 1% of the asteroid's mass, and only 10% MT go into accelerating the main part of the asteroid away for an impact situation. That's quite an important question if we were ever faced with solving the problem in the coming year.

Anyway, that was the "why" of the mission, and a reasonably useful reesult was obtained (IIRC, about 1/3 of the momentum applied stayed in the "rubble pile, resulting in the moon's orbit being changed by about 30s - which was easily detectable. But this paper in't about that - it's about how such asteroid moons can form.

Available scenarios include : forming together (in which case, why didn't the last parts of the accumulating debris go onto the larger body, with it's stronger gravity field), capture (which Darwin (George) showed is actually quite difficult, unless you've got a "third body" to take away some of the energy and angular momentum from the encounter), impact on the main asteroid (as for the Earth's Moon, most likely ; but why then are asteroid moons so common?). Or, as this paper proposes :

[An asteroid can] deposit material into orbit via landslides.

Now, there's an idea. It doesn't sound right. If the energy comes from the movement of material under gravity, how can that material get into orbit? The trick here is that this was happening while the primary was having it's spinn increased by the Yarkovsky–O'Keefe–Radzievskii–Paddack (YORP) effect, and that effectively reduces the weight of materials near it's surface (far from the centre of spin). And that is enough to allow some of the materials disturbed by such movements.

Which sounds very weird. But that is how weak the gravitational fields are on the "small bodies" of the solar system.

The paper has more work about the time sequence of events they infer leading to the formation of this moon, and in this shape. We're probably going to see a number more examples of this as the "Lucy" mission works through it's list of asteroid fly-bys.

Back to List.

MOND Count

Formation and Growth of the First Supermassive Black Holes in MOG

A while ago, I did a count of papers posted to Arχiv over the last few years of papers on the subject of "MOND" (MOdified Newtonian Dynamics") because someone had posted (again) the false claim that scientists don't look at alternatives to the main theories of science. They do get attention - just not a lot. Previously I'd searched for "Mordehai Milgrom", "MOND" and "Non-Newtonian Gravity". "MOG" is another search term to watch out for. My previous searches have looked at yyyy-12-31 as a cut-off date, so I'll carry on with that.

Hmmm, moderately interesting - at current publishing rates, Mordecai Milgrom looks likely to retire (or stop publishing - same thing, really) about 2030~31.

Back to List.

Formation of the Cartwheel Galaxy

Star formation history of the post-collisional Cartwheel galaxy using Astrosat/UVIT FUV images

The Cartwheel Galaxy is one of the prettier - and weirder - sights on the sky. To the telescopic eye, one can see an outer ring of bright, bluish stars, and an inner ring of slightly redder stars around the margin where approximately the edge of the nucleus would be, and a set of slightly curved "spokes" faintly link the two main structures.

view of the Cartwheel galaxy and it's two, nearer, companions 'G1' and 'G2', about one diameter from the centre of the CW galaxy.
It's a pretty thing, shouting a really loud "How?" across the universe across to us.

The general understanding of this galaxy's structure is that another galaxy passed, nearly centrally, through the "Cartwheel" galaxy. In 1996 radio telescopy detected a discontinuous band of neutral hydrogen emission stretching from the Cartwheel towards a nearby galaxy, "G3", which is interpreted as the impacting or "bullet" galaxy. Critically, this emission structure extends beyond the other nearby galxies, "G1" and "G2". A wider view of the system than the "pretty one" shows the relationship of the galaxies: (this is an image from the DSS - Deep Sky Survey - but as a negative image) :

annotated wider view of the system, negative, with 'G3' galaxy about 5 times the Cartwheel galaxy's diameter out from it's centre

This paper reports studies on UV emissions and brightness to estimate the star-forming rates of the Cartwheel's brighter regions, and infer the time since their formation. So we now know that the galaxies were in collision at about the start of Earth's Carboniferous period, 300 million years ago. With some modest error bars.

Back to List.

December HTML things

link

What got me prompted to try to get back to work was seeing a reference to an unusual bit of HTML CSS - the "rotate" property. With this, you can ... go on, guess ... rotate text (well, strictly, a division of a page) about a point.

So, naturally, I wondered, will it work in Blogger's sandbox? Let's try.

I'll start by using my handy-dandy bit of table template :

Heading Column
a Α

And now, I wrap it up in the "div-Wol-Rotate" division I've defined up in the page's style header.

div-Wol-Rotate { height: 100px; width: 200px; background-color: lightgrey; rotate: 1 1 0 60deg;

Not working. 3 axes? no Angle? No (and the bgnd-colour isn't working either) Ah, whoops, missing trailing ";" Size? Nope. Still doesn't seem to be working. Quick dive over to the "try it" sandbox at W3schools.com ... and the code above works there, but not here. Deduction : Blogger doesn't like rotation. Oh well

I picked up this question from reading Elf Sternberg's blog, about recreating the 1994 alt.sex FAQ website in modern coding. Now, I don't know how often Mr Sternberg gets people linking to the programming bits of his website, but just so he knows, from the above - people don't only read his website one-handed.

I can't claim to have fully understood the "rotate" property - the height, width and definition of the rotation axis interact in some way I need to look at to shear the image before the rotation is applied. I need to work on that. but it doesn't seem to work here, so I'll shelve it for now.

Fiddled with the CSS for HR elements though.

A bit more mucking around with tables in IMF article.

Back to List.

The Cosmic Train Wreck Galaxy, Abell-520

How Zwicky already ruled out modified gravity theories without dark matter,

This 2017 paper just jumped out at me while doing the counts of "MOND" and "MOG" and "dark matter" papers while checking out the question above.

It has interesting bits. Who could resist a phrase like "The Cosmic Train Wreck Galaxy"? Not me - that's for sure. It sounds like something "fun" happened, in the sense of "may you live in interesting times", and "best observed using a telescope". Even better that it's in a context of overheated extravagant astronomers, not staid, restrained stage actors who couldn't emote if you pulled the stick out of their back, leaving splinters. Talk dirty to me, like a central 3-4 × 1013M mass clump with mass-to-light ratio 800 × M/LR☉. i don't know how dirty, but that's definitely dirty talk! (There's probably a typo - "LR☉" for "L" ; I'm sure I'll find a resolution somewhere.)

OK - how dirty are we? 800 Solar masses for each Solar luminosity? Let's see - Sol is heavier than about 90% of stars, but the decrease in luminosity with decreasing star mass is steep. A half-M star has a luminosity of 0.063 × L (about 1/16th) while at the other end of the mass scale, the luminosity ramps far more steeply than the star mass. A 20 × M star shies with luminosity 53600 × L. So, calculating the mean brightness of matter on a gelactic scale is .. hairy. It's not clear to my limited maths if it'll be a figure above or below 1 × M/L. Probably not wildly off from 1 × . So 800 × is a galaxy with remarkably little light for it's mass. Something seems to have wrecked the galaxy's mechanism for turning mass into luminosity. No wonder people noticed them (once they'd done their accounting).

The "Train-Wreck Galaxy" has 4 references linked ( 13 M. Jee, A. Mahdavi, Hoekstra et al. The Astrophysical Journal 747(2), 96 (2012).
14 D. Clowe, M. Markevitch, Bradaˇc et al. The Astrophysical Journal 758(2), 128 (2012).
15 M. J. Jee, H. Hoekstra, A. Mahdavi, and A. Babul The Astrophysical Journal 783(2), 78 (2014).
16 Q. Wang, M. Markevitch, and S. Giacintucci arXiv preprint arXiv:1603.05232 (2016).
) but the rather more accessible Wiki article suggests that later observations have picked up considerably more luminosity, reducing the seeming unusualness of the galaxy cluster by quite a lot. It's still odd, but less odd.

The other list of cosmological oddities is worth looking at too. As they phrase it, so Lambda;CDM is rightfully a good effective theory. But is it a fundamental theory? Many puzzling observations make this nonevident. :

  1. The CDM particle, the WIMP, if it exists, keeps on hiding itself 6 more years since its “moment of truth” ;
  2. One observes 19 quasars with spins aligned with their hosts large-scale structures on a scale of almost 1 Gpc (ref.10) and dozens of radiojets from AGNs aligned on a scale of 30 Mpc (ref.11) ;
  3. A ring (actually, a spiral) of 9 gamma ray bursts extends over nearly 2 Gpc (ref.12) ;
  4. the "train wreck galaxy" that caught my attention above (which may not be so weird) ;
  5. in the cluster A3827 the offset between baryonic and dark mass (ref.17) is an order of magnitude ‘too large”(ref.18) ;
  6. Puzzles in galaxies include: the brightness fluctuations in the Twin Quasar allow a DM interpretation in terms of a large halo of rogue planets in the lensing galaxy (ref.19) ;
  7. the observed satellites of the Galaxy lie in a plane, not in random ΛCDM directions (ref.20) ;
  8. the predicted transition for the most massive galaxies to transform from their initial halo assembly at redshifts z = 8 − 4 to the later baryonic evolution known from star-forming galaxies and quasars is not observed (ref.21) ;
  9. the galaxy power spectrum deduced from SDSS-III observations fits well to the stretched exponential exp[−(k/kb)1/2] from turbulence (ref.22) ;
  10. Various further arguments can be found in our investigations (refs.23,24) ;

That's an interesting list of "problems" with the ΛCDM model, while acknowledging that it is still a good model. No joy there for Archimedes Plutonium and the Electric Universe crowd of the deluded. Frankly, while they're interesting, they're not interesting enough. Just worth me noting down, to counterblast people claiming "censorship" and "suppression" of unorthodox thinking.

Why did Zwicky think you needed dark matter, not MOND/MOG, in the 1940s?

On a third reading, I still don't get the argument the authors are making. They seem to think that, for the cluster they study (Abell 1689 ; the mention of Abell 520 is just in the "problems with cosmology" section of the introduction ; it doesn't help that this is a fairly early draft, with several arguments re-written in the text.) their calculation under several generic "modified gravity" models, still need to add some non-visible mass to the system to match the observed lensing [strong ("SL" with multiple arcs for identifiable background galaxies) and weak ("WL" with shape distortion of individual background galaxies)]. Their preferred (why?, not clear) DM candidate is neutrinos. Particularly "heavy" neutrinos with a mass of ≦ 2.0 eV. Most experimental estimates for conventional neutrinos gives them a mass of ≈ 0.1 eV - a substantial difference, but some particle physics ideas produce "sterile" neutrinos in this mass range.

I think they need to make their point on Zwicky clearer - the title is attractive, attention-getting. But I can't follow their reasoning.

Back to List.

Element Abundances

The abundance discrepancy in ionized nebulae: which are the correct abundances?

A long-standing question I've had is : what is the average composition of a planet (asteroid, "small body"), versus this particular surface, or an ore body. There are a lot of slightly different answers in the literature, and not a lot of clues pointing to the most-accepted (and most detailed) answers. Slightly different, but related, answers often given for this question include :

  • "What Is the Composition Of ("WICO") the Sun?" (not necessarily the same as the average planets)
  • "WICO the Sun, if you lose the H and He?" (ditto)
  • "WICO a CI meteorite? Is it the same as their parent bodies before impact heating?"
  • "WICO the Earth (including the core)?"
  • "WICO the Earth's (or any other planet's) surface?"

They're all quite similar questions, with not quite the same answers. And the differences may (or may not) be significant when you're trying to consider the changes implied by formation processes. Obviously, I'm looking for a full, 92 element composition. How it's expressed is an issue to be wary of - a short while ago I had one of those "drowning goldfish" moments when someone was talking about "atmospheric carbon di-oxide being absorbed by the lime (calcium oxide) present in basalts". It took a time to realise that he was reading the traditional geochemist's way of reporting compositions as "oxides" (e.g. 40% SiO2, 10% CaO, 10% Al2O3 …), and taking that to mean that a basalt contains grain or molecules of actual CaO, with all the chemical properties of a lump of quicklime. Which ... well, I can see why this relic from gravimetric "wet chemical analysis" of the 18th century is confusing. But it obviously can lead to profound misunderstandings. Gotta watch that banana skin.

I'm slightly surprised that my confused correspondent didn't notice that the trace elements are normally reported as "ppm of [Element]", rather than in the "oxides" presentation. I suspect I'm going to have to address this in more detail one day. It may explain a lot of confusion about geochemistry.

So, interesting paper title, given that context. Does it deliver?

OK - first thing is to see the paper as a contribution to a programme called "Planetary Nebulae: a Universal Toolbox in the Era of Precision Astrophysics" - which isn't terribly optimistic for my "geochemistry" interests. Still useful enough though - what a star spits out this gigayear is likely to affect the composition of next Gyr's molecular cloud, hence stars and attendant planets.

From TFA : "However, the heavy-element abundances derived from collisional excited lines (CELs) and recombination lines (RLs) do not align." - which is what they mean by an "abundance discrepancy". Going back to the "oxides" presentation discussion above, the technology of measurement influences the precise data collected, and how it's expressed. Specifically, since they're measuring line strengths in spectra, I'd be surprised if they presented their abundances as "oxides".

Actually, it's quite a big effect. From the paper's introduction :

However, there is a problem, and it is a big one: we are not certain about the correct abundance of heavy elements like O, C, N, Ne, collectively known as “metals”. Since the pioneering works of Bowen & Wyse (1939) and Wyse (1942), it has been known that the abundances of these elements determined from the bright collisionally excited lines (CELs) (e.g., [O III] 5007, 4959) are systematically lower than those inferred from the faint recombination lines (RLs) (e.g., O II 4649, 4650). The ratio between both estimates is known as the ”Abundance Discrepancy Factor” (ADF), and it has been found to be around a factor of 2-4 in H II regions (García-Rojas & Esteban 2007) but can reach values of more than 500 in some PNe (Wesson et al. 2003).

A discrepancy of 2~4×, but up to 500× - yep, that's a problem. But, it's an astronomer's problem, not really a geochemist's (or planetary scientist's) problem, so I think I'll just file that as a detail to watch out for.

Annnnd … the paper is discussing the various ideas for explaining those discrepancies, but not presenting an example composition for any of their target objects. Ho hum. That just leaves me with the materials I've already collected in my "compositions" collection ( http://ArXiv.org/pdf/2105.01661v1.pdf which gives Solar photosphere and CI chrondite compositions for 83 elements ; 54 terrestrial rocks from "Science Snippets" on my HDD ; … ; ; ; ) and anything I add to it. Back to the grindstone. I still don't have estimates for the average composition of the (local) universe. Granted that the outer arms of the Milky way are second (or third) generation stars, enriched by at least one (possibly two or even three) previous generations of stars burning, exploding (for the larger ones) and contaminating the (primordial) interstellar medium ("ISM", H, He, trace Li) with "metals" which change the structure and nuclear processes of the stars. Probably I'd be looking for a local Milky Way composition, a Solar system one, and maybe a composition for an estimate first generation ISM.

Back to List.

6-planet system is a multi-star system.

HD 110067 is a wide hierarchical triple system

A couple of weeks ago there was a moderate amount of news fuss over the announced discovery of a planetary system containing 6 planets whose orbital periods were at small-integer ratios (2:3 ; 3:4, etc). News coverage was all over the place. Which was all very well and good, but the source paper was a deal more nuanced than the write-ups in the popular press (well colour me sideways and brush me with a Blackadder-shaped paint pot!) Most was made of the small-integer orbital period ratios - which to me (and most people) talks of the system having gone through a phase of planetary migration that carved these ratios - probably when there was some degree of "viscosity" in the system to damp interactions to some degree. (all that means is - considerable remaining gas, or many "small bodies", being present to diffuse energy and angular momentum less violently than by collisions or close approaches). The more interesting thing, to me, is how a relatively sparse set of observations were translated into "six planets". Although the reviewers at Nature must have challenged this, I'm sure they were convinced by the data in a way I can't (and am not interested) to challenge. I accept that they've got at least some of the planets adequately located.

This original paper was published in Nature, (here, titled, "A resonant sextuplet of sub-Neptunes transiting the bright star HD 110067", which contains substantial "Supplementary Information"), but the paper was also released on Arχiv without the access restrictions of Nature (but also without the "SI"). You should be able to judge the paper on it's own merits, but the SI is probably in large part a response to the reviewer's criticisms (or their anticipated criticisms - same difference really). I'll refer to this as the "Nature paper", with the titular paper of the "triple star" system being referred to as the "Arχiv paper".

The Nature paper reports several long-period observations of the target star (HD110067 - "HD"= the Henry Draper Catalogue and a reference number) using the TESS (Transiting Exoplanet Survey Satellite) observatory in two sessions, two years apart. The first observing session recorded 4 transits (3 of which had a common spacing - indicating one planet ; one unconstrained transit). The second session recorded 9 more transits, which confirmed the first candidate (HD110067 b, period 9.113678 days) and a second (HD110067 c, period 13.673694 days)candidate observed in both sessions. Leaving 8 unexplained transit events - which was sufficient justification for observing periods on CHEOPS (CHaracterising ExOPlanets Satellite) of fractions of a day around possible repetitions of one or more of those potential objects. In considering potential observing times, the negative evidence - that neither of the TESS sessions had recorded transits except at the noted times - which constrained the possible orbital periods associated with those transits. This series of observations confirmed HD110067 d (period of 20.519 days). Those periods (bc ; cd) were both close to a "Mean Motion Resonance" (MMR) ratio of 3:2, justifying further data collection.

This is where the "chain" of logic starts to get a bit thin. There was adequate justification for targetting more observation time, and some of that justification is in the "SI" of the Nature paper. Another planet was discovered by analysis of one of the transits, for which there was an absence of a second transit in the other TESS data set. This was identified as HD110067 e (period 30.7931 days, in a 3/2 MMR with planet d) - but only on the basis of two transits. That left two unexplained transits in the TESS data. With the assumption that these represent more planets, which continue the chain of low-integer period ratios. Which is how they proposed two additonal planets, which should have been seen in the original TESS data. But the periods they should have been there were also when the Moon or Earth were also near the satellite's line of sight (TESS orbits between a little above geocentric orbit, and a little below Lunar orbit), so originally these records weren't analysed because of glare. When this study gave adequate reason for study, ... there were the expected transits, in the (rather noisy) data. So that gave confirmation of HD110067 f (41.05854 days) and HD110067 g (54.76992 days). Which is sufficient for the reviewers and for science - but it's a lot thinner than the stereotype of multiple transits per planet.

That is how future planet identifications are going to go. If we're looking for evidence of "Earth-like" planets, you're not going to get the stereotypical multiple tranists at a rate faster than 1 observation per planetary "year". So inevitably, we're going to see increasing numbers of single-transit detections, and longer chains of inference from raw observation to deduced planet (and their properties).

Some work was also done using ground-based observations, which revealed a fair amount of magnetig and radial-velocity noise in the system, but some support for the planetary system too.

The inferred planet properties are nothing wildly unusual. They're close to the star (necessary, to find multiple transits of the first few members of the chain), at average radii (of 0.0793, 0.1039, 0.1362, 0.1785, 0.2163, 0.2621 AU) all within Mercury's orbit compared to the Sun. The deduced surface temperatures (from the star's temperature of 5266 K, and the orbital distances) are high, but not extreme (800, 699, 602, 533, 489, 440 K); a short distance out from the detected planets would be the "habitable zone" of the star. The star is estimated to be about 8 billion years old - about half way through it's estimated 17 billion year lifespan, from it's 0.798 solar-mass mass. The planetary masses (5.69, <6.3, 8.52, <3.9, 5.04, <8.4) are in a range not found in the Solar system, but actually very common in exoplanet studies.

The triple-star aspect

The second, "Arχiv" paper is also an event likely to recur. While the first group had been doing their transit-based study, other data was present in already-taken databases. The Nature paper took the primary as being a single star, it had been classified as a "wide binary" since at least 2001 (references in the Arχiv paper) using positional records from 1893 to 2015 (that must be a later revision of the "Washington Double Star Catalogue" than the 2001 reference). To my moderate surprise, the listed companion is also in the HD Catalogue, and in the original numbering system (number 110106), which implies it is quite bright, quite close, and could be spectroscopically investigated at the turn of the last century. That clearly shows that the two stars are indeed, a wide binary, separated by 415 seconds of arc, which at the star's distance (32.25pc, 105.2ly, from Gaia data) translates to some 13394 AU - which is very comparable to the spearation of Proxima Centauri and Alpha Centauri A+B. To further the similarity of the two systems, both Alpha Centauri and HD110106 are double stars - though Alpha Centauri can be resolved by telescope, while HD110106 hasn't (yet).

Since these secondary stars are well out of the line of sight to the primary, with it's resonant planetary system, they don't interfere with the detection of the primary's planets. The presence of a wide binary companion (s) do raise questions about the planet-forming process around the primary, so it's a safe bet that more investigations of this system are planned. You heard it here first.

Back to List. 2023-12-20

Dating of a Latin Astrolabe

https://arxiv.org/pdf/2311.17966.pdf quote mark

This is just a piece of fun. Ha-ha, but serious fun.

An astrolabe is an astronomical observation instrument consisting of a plate (the "mater"), a targeting sight bar (the "alidade"), a "plate" engraved with the coordinates of the celestial sphere, a movable metal net ("rete") which indicates the position of various stars, and a suspension point and vertical bar ("rule", which keeps it aligned with local vertical). It can be used for surveying, for measuring the altitude of heavenly bodies (and hence, local time), and of course identifying heavenly bodies (particularly since they normally had a celestial chart engraved on one face). Mediæval users suggested that there were a thousand uses, many of which were listed in a manual on the device written by Geoffrey Chaucer.

One of these days I'm going to find one in a second hand shop (I've got a sextant somewhere a working toy) and learn how to use it.

By pointing the sight ("alidade") at an identified star, and aligning the net ("rete") to match the star's indicator to it's altitude, one can work out the local time from the celestial grids enmgraved on the plate ("mater"). Other ways of solving two knowns for an unknown can also be done with the instrument - including (for the many Islamic users) locating the direction of Mecca in order to pray towards it (the "quiblah". [Totally tangential, I wonder if that was the inspiration behind Rowling's "Potterverse" equivalent of the "Fortean times", the "Quibbler".]

That's the mechanics. The astronomy, by which the astrolabe's manufacture (or design, or calibration) date was determined is that the precession of the equinoxes (rotation of the Earth's spin axis compared to the celestial sphere) moved the position of the stars relative to the celestial grids, allowing the determination of the manufacture (or design, or calibration) to approximately 1550 (CE ; all dates are CE, not AH).

The reason I specify "manufacture (or design, or calibration)" is that once soneone has made one of these, then the additional work to make a second is relatively small, separating "design" and "manufacture" events. The way of indicating star positions (extensions of the "rete" as long pointers, with the star being at the tip of the pointer) is potentially amenable to being fine-tuned with a pair of pliers, so a 1550 design might still be modifiable for use in 1650 (or a 1450 design for use in 1550). In fact, if the "rete" were made by casting in bronze, the raw casting would very likely need some cleaning up, and adjustment because castings - particularly reticulate ("net-like" ; same Latin root) castings often distort as the cool from molten bronze to room temperature. It's very unlikely that the raw casting would have been flat to the plane of the "mater" without some tweaking.

The pictures in the paper aren't high enough resolution to see if the "pointers" have been tweaked. But I bet they have. (when I posted the full-resolution image, I looked more closely - the "rete" is made of several pieces, which look to be cast and have then been soldered together. Yep - that almost certainly needed "calibration" after soldering together. Whether that was higher precision than this naked eye instrument would have needed ... maybe not. But it's at the limit.

Just a fun bit of sideline stuff. For me. Isn't this a thing of beauty?

Figure 6: The front of the astrolabe, with a plate, the rete, the alidade (diagonal bar), the rule (vertical bar).Caption : Figure 6: The front of the astrolabe, with a plate, the rete, the alidade (diagonal bar), the rule (vertical bar).
Back to List.2023-12-21

Initial IMF

Nebular dominated galaxies in the early Universe with top-heavy stellar initial mass functions

That's a convoluted name for a relatively simple idea, with an awful lot of consequences.

As JWST has rapidly improved our state of knowledge of galaxies in the early universe, there has been an increasing tension between the inferred properties of those galaxies. What we can observe are luminosities at different wavelengths. But to convert those observations into a model of stellar populations, ages, brightnesses and interactions, we need a way to convert those luminosites into a number of stars, at a selection of brightnesses. This is not simple, because the brightness of a star does not relate simply to it's mass. The relationship is far from linear.

Which, frankly, has been known since the start of astrophotometry - the (reasonably) accurate measurement of star brightness - in the 1870s. When that was combined with positional measurement (revealing parallax, for relatively close stars, and from that the distance to the stars), and analysis of eclipsing binaries (which gives disc sizes, relative brightness in a pair, and the absolute masses of the stars in a pair), a relationship between a star's mass and it's brightness (when on the main sequence - where stars spend most of their lifetimes) emerges. But it's not a simple relationship.

In my handy-dandy pile of astronomical data and calculations, I put together a handful of ranges of mass to cover the (main sequence) brightness of stars in the Milky Way.

Stellar Mass - Brightness relation (referenced to Sol)
Star Mass (M) Luminosity (L) Main Sequence lifetime (T, Gyr)Comment
0.1 0.001 1000 Quite approximate
0.2 0.006 333 Luminosity ≃ 0.23 × (mass ratio)2.3
0.3 0.014 214
0.5 0.063 79 Luminosity ≃ (mass ratio)4.0
0.7 0.24 29
0.9 0.7 13
1.0 1.0 10.0 (solar estimate)
1.3 2.9 4.48
1.8 9.4 1.86
2.0 16.0 1.25 Luminosity ≃ 1.5 × (mass ratio)3.5
2.5 37.1 0.6739
5.0 419.3 0.1192 Gyr
7.5 1733 43 Myr
10 4743 21.1 Myr
15 19600 7.5 Myr
20 56670 3.73 Myr, Luminosity ≃ 3200 × (mass ratio)
30 84000 3.57 Myr
40 152000 2.63 Myr
50 240000 2.08 Myr
75 550000 1.37 Myr
100 980000 1.02 Myr
I have two references for this : arxiv 1710.11134.pdf and "The Astronomical Journal, 154:115 2017, Table 1" ; it's likely that the exact numbers will change with time, but not by much.

Already the lower rows of that table are decidedly uncertain. Massive stars (more than a few solar masses) typically lose more than a few percent of their mass during their evolution (through "solar wind", strong coronal mass ejections and other processes) meaning thir luminosity and lifetime will vary considerably from the initial expectation, which is based on nuclear physics. But you can clearly see that the total luminosity varies a lot more than the stellar mass. 100 solar masses of material will shine (approximately) as brightly as 100 Sol if in individual stars, but nearly 10,000 times brighter if collected into one body.

A significant uncertainty in the nuclear physics part of this model is how much non-primordial (H, He, a touch of Li) matter there is in the mix. Since nuclei like C and N take part in catalytic series of nuclear reactions at lower temperatures than without them, even small amounts of such nuclei can considerably alter the rate of energy release form the cores of medium-size stars, and hence the mass- luminosity relationship.

When all the data you have about a distant object is it's total luminosity and distance (from it's redshift), you don't have enough data to directly infer the mass of the body (galaxy) from that data. You have to, somehow, incorporate some relationship between luminosity and mass - such as the table given above. The conventional way of doing this - in the absence of any better information - is to assume that the remote object has a similar relationship to that present in the Milky Way. It's an assumption, but it's an explicit assumption, and very open to challenge. This paper challenges that assumption for high red-shift galaxies.

(It is worth noting that, because the evolution of bright stars releases non-primordial nuclei (heavier than Li) into the interstellar medium of a galaxy, it is expected that the brightness-mass relationship will change with time - and that always has been expected. Exactly how the mass-luminosity relationship changes with time remains a matter of debate though.)

This paper looks at the spectrum of the light from three of these high red-shift galaxies (GS-NDG-9422, at red-shift z = 5.943 [reference 14]; the "Lynx arc" [from "Massive Star Formation in a Gravitationally Lensed H II Galaxy at z = 3.357. reference 15 from the paper] ; and A2744-NDG-ZD4 at red-shift z = 7.88, reference 16) and sees a distinct break in the spectrum that they infer as meaning there is a lot of gas heated to tens of thousands of Kelvin, producing significant "free-bound" electron interactions. Those temperatures are a lot hotter than the hydrogen ionisation event which gave the Cosmic Microwave Background (CMB) radiation at a temperature of approximately 3000 K, allowing the inference of considerable numbers of stars of high brightness that "pump up" the galaxy's gas temperature to these levels. Which means that the "Initial Mass Function" of the galaxies is richer in heavy stars than the Milky way of today. And in turn, that means that the inferred mass of the galaxies - which had been stressing models of the evolution of the early universe would be considerably lower.

Having a "top heavy" initial mass function in the early universe will require a lot of re-calibration of models of the evolution of galaxies. But at least there is now some data to suggest in which direction, and by how much, to adjust the models.

I see that Dr Becky Smethurst has done a video-blog on this paper on her "Night Sky News" channel.

(The heaviest species detected in the spectra is argon (lines at λ4711 and λ4740 &Angstrom;, which would suggest some "silicon burning" nuclear reactions.)

End of Document
Back to List.

And that's New Year, so on to the next page.

2023-08 August Science Readings

August Science Readings

Well I'm only 5 days into the month as I start, this time.

Articles studied this August - some of which might go to Slashdot.
ToDo List Black Body radiation/a>
Latter-day Gamma-ray Coordinate Network
New evidence of plant food processing in Italy before 40ka
End of document

2023-08 August Science Readings.

And I've got to the end of the month with very little done. Other stuff in life. Several things got into the pipeline, but I haven't got time for.

To Do List - Black Body radiation & Mean Free Path

link

Once again a question of the form "If this [whatever] gets so-much hotter, how will it's colour (black body spectrum) change?"

OK, it's obviously a question for whatever equation produces the "hump" (log-log) or "spike" (linear-linear) graph in all the physics (astrophysics) text books. So I know qualitatively what will happen. But I want to calculate how much it'll change. Which calls for a spreadsheet.

    Specification
  1. Put in a base temperature, get a colour (peak frequency). Put in a second temperature, get a second colour.
  2. Tabulate temperatures and colours. One datum is that 2.83 K gives a signal in the microwaves - CMB. 21cm Hydrogen? Orange is arounf 1600k
  3. subsidiary : wavelengths for colours - it's got to be tabiulated somewhere. If only well-known emission/ absorbtion lines like sodium-D = yellow. Obviously this is going to be rather arbitrary.
  4. I shouoldn't neeed to plot the BB radiation curves, but I'd like to. Two temperatures, can OO(Calc) do "fill between"? Can Google Calc? Obviously ties into item 2.

Just listening to the radio, about getting stars started, and the concept of Mean Free Path reared it's ugly head again. Need to look at that too.


Latter-day Gamma-ray Coordinate Network

link

For a number of years (since ... when I was on CI$, so pre-2000) prompt reporting from space-based Gamma Ray detectors has used a mailing list to distribute alerts of spike in GR detections, and by inference, the occurrence of a gamma-ray burst somewhere on the sky. That system has been deprecated as larger numbers of "high energy events" are being monitored, from gravity-wave detectors (3 systems), neutrino detectors (3 operating, several in construction) gamma- and x-ray space telescopes amd other systems. That's annoying, because the simplicity of checking my email has been replaced with needing to register on a NASA website, download and install Python, compile and install several programmes (I'm not sure how many), and then get really informative responses :

 topic=gcn.classic.text.AMON_ICECUBE_COINC, offset=None
b'Subscribed topic not available: gcn.classic.text.AMON_ICECUBE_COINC: Broker: Unknown topic or partition'
topic=gcn.classic.text.FERMI_GBM_TRANS, offset=None
b'Subscribed topic not available: gcn.classic.text.FERMI_GBM_TRANS: Broker: Unknown topic or partition'

Which is as useful as something not very useful.

As so often, the documentation seem to know that all users will know everything about what and how a "streaming protocol" is, and how to use one. Which ... well they call it "Kafka", and the name is well-chosen. I know how K. felt.

OK, now I'm getting some "content" - I left the terminal with the python code running while doing other stuff :

topic=gcn.classic.text.FERMI_GBM_ALERT, offset=917
b'TITLE:           GCN/FERMI NOTICE\n
  NOTICE_DATE:     Mon 07 Aug 23 14:38:00 UT\n
  NOTICE_TYPE:     Fermi-GBM Alert\n
  RECORD_NUM:      1\n
  TRIGGER_NUM:     713111879\n
  GRB_DATE:        20163 TJD;   219 DOY;   23/08/07\n
  GRB_TIME:        52674.82 SOD {14:37:54.82} UT\n
  TRIGGER_SIGNIF:  6.7 [sigma]\n
  TRIGGER_DUR:     0.064 [sec]\n
  E_RANGE:         2-2 [chan]   23-47 [keV]\n
  ALGORITHM:       26\n
  DETECTORS:       0,0,1, 0,0,1, 0,0,0, 0,0,0, 0,0,\n
  LC_URL:          http://heasarc.gsfc.nasa.gov/FTP/fermi/data/gbm/triggers/2023/bn230807610/quicklook/glg_lc_medres34_bn230807610.gif\n
  COMMENTS:        Fermi-GBM Trigger Alert.  \n
  COMMENTS:        This trigger occurred at longitude,latitude = 236.15,-12.07 [deg].  \n
  COMMENTS:        The LC_URL file will not be created until ~15 min after the trigger.  \n'

Which isn't much help, but I'd also received an email with the same content. The email is more USABLE.


New evidence of plant food processing in Italy before 40ka

https://www.sciencedirect.com/science/article/abs/pii/S0277379123002093?dgcid=coauthor Quaternary Science Reviews Volume 312, 15 July 2023, 108161

(Prepared a while ago, offline. Stalled.)

If anyone actually thought about the implications behind the hype about the "Palæolithic Diet", none of this would come as a surprise. But since thinking about things is antithetic to the interests of the "influencers" behind the "Palæolithic Diet", then it's unlikely to get much traction from them.

Abstract : Evidence of plant food processing is a significant indicator of the human ability to exploit environmental resources. The recovery of starch grains associated with use-wear on Palaeolithic grinding tools offers proof of a specific technology for making flour among Pleistocene hunter-gatherers. Let’s get this clear – this is HUNTER GATHERERs making flour – and by implication, breads, porrages, gruels, etc. Just because they’re called “HUNTER GATHERERs” doesn’t mean that they’re living on mammoth steaks and bronto-burgers exclusively. As modern studies of modern HGs suggests, upwards of 50% of their calories come from the GATHER part of the lifestyle, and very often, it's gathering by the women-folk. Collecting tubers with a waen on the tit is probably less population-risky than mammoth-hunting while giving suck. Less appealing ot the "I wanna eat a mammoth-burger" crowd of Andrew Tate wannabes.

Continiing FTFAbstract : The recovery of starch grains on a Mousterian grindstone at Bombrini suggests that the last Neanderthals not only consumed and processed plants but also made flour 43 - 41,000 years ago. Starch grains attributable to Triticeae on Protoaurignacian grindstones at both sites testify that Sapiens were processing wild cereals at least 41,500 - 36,500 years ago when they expanded into Eurasia, long before the dawn of agriculture. Does that need expansion? The sites are in Italy, Bombrini cave overlooks the Mediterranean, near the Monaco border ; Castelciveta is in Campania, well inland, and is sealed by the Campanian ignimbrite from the Phlegraean Fields supervolcano in the outskirts of Naples (erupted 39,220 ~ 39,705 BCE). One of the caves (I didn't note which) has plant-preparation tools at two significantly different levels (ages), giving three sites at two locales. There’s no particular reason to believe that proto-Aurignacian in Bombrini cave is close to the same date as at Castelcivita cave ; it may be earlier, overlapping, or later. As always, other actual implements may not have been identified, as always, though the archaeologists have used a fairly broad set of criteria, andone of the tools was identified as such during excavation, allowing immediate "sterile" (dig sites aren't steril ; nor is soil) collection.

The dating at Castelcivetta has a latest-possible date in the sealing Campanian Ignimbrite deposit. The Bombrini specimens had “areas - including an apex - covered by carbonate incrustations formed during their permanence in the cave.” (OK ; clearly not edited by a native-English speaker. Bear in mind. I recently met "permanencia" in my Spanish as somewhat equivalent to "period of residence". Trivial point.) And they got U-series dates from that carbonate, giving a latest possible date there. Why the Campanian Ignimbrite is considered to mark the end of human occupation of Campania for a considerable period is left as an exercise for the reader.

From my PoV, it is interesting to see the morphology of the detected starch grains on the tools. That sort of material wasn't covered in my mineralogy microscopy. I note in particular the pseudo-isotropic bisectrix ficures of the starch grains under XPL ... which implies that their spherical shape makes them quite strong converging lenese in a generally plane (not convergent) polarisation field. I'll have to try to show that to the Microscopy Club if they ever meet again.

Fig. 4. Starch grains and phytoliths from the grindstones of Riparo Bombrini and Grotta di Castelcivita, at bright-field microscope and at polarizing light microscope… note specifically the extinction croses in the microcrystalline starch grains.

The body of the paper had a few worthwhile highlights too :

  • Half of each sample of B-A2 and B-M1 was also subjected to heavy liquid separation using zinc chloride, according to Mariotti Lippi et al. (2015). Why? What were they expecting to find? We’d use ZnCl2 to make density columns between about 1.8 and 2.5 SG, so that covers a lot of territory. Obviously looking ... ah, if the ZnCl2 is fairly dense, say 1.8 SG, it would float off starch (organic)grains at ~1 SG, but drop out minerals like calcite (2.7) quartz (2.6), and phyllosilicates (clay-ish, 2.0~2.6 SG) all while keeping the starches in a low-osmosis potential fluid.
  • In Africa a cobble used for grinding plant materials is mentioned from the Early to Middle Stone Age site of Sai Island, Sudan (Van Peer et al., 2003 “The Early to Middle Stone Age transition and the emergence of modern human behaviour at site 8-B-11, Sai Island, Sudan.” J. Hum. Evol. 45 (2), 187 – 193).” Which gets a Spock-Fascinating.GIF from me, not least because it's very old and a long way from the "Fertile Crescent" associated with the origin of agriculture.
  • New evidence of processed plant food is illustrated by pulse remains from the Late Middle Palaeolithic to Upper Palaeolithic at Shanidar cave (Iraq) and Franchthi cave (Greece) (Kabukcu et al., 2022 [Cooking in caves: palaeolithic carbonised plant food remains from Franchthi and Shanidar. Antiquity 2023 Vol. 97 (391): 12–28 https://doi.org/10.15184/aqy.2022.143 ]” I know the name of Shanidar – isn’t that the "crippled Neanderthal" cave? These two sites pretty well "bracket" the Fertile Crescent", but again, much eariler than the conventional "origin of agriculture".
  • From Kabukcu, above, Almost all sites from these regions dating to the Middle and Upper Palaeolithic and the Epipalaeolithic/ Mesolithic periods, for example, provide evidence for the use of wild almonds, which contain high levels of cyanogenic metabolites that can produce hydrogen cyanide. [...] Several other plants also feature prominently in the regional archaeobotanical record, including tannin-rich wild pistachios (terebinth), wild pulses (some containing neuro-toxic compounds) and astringent wild mustards. Most of these plants require several preparation steps to leach out unpalatable and/or toxic compounds prior to consumption. The long-term and widespread use of almonds, terebinths and pulses therefore suggests that Palaeolithic foragers developed processing technologies and associated food preparation practices that enabled their routine safe consumption. (Me : That's going to stay right OFF the Palæolithic Diet menus.) … new evidence concerning the long-term histories of Palaeolithic plant food use and associated food preparation practices from two multi-period sites: Franchthi Cave (Greece) and Shanidar Cave (Iraqi Kurdistan). We focus on the analysis of amorphous, charred plant aggregates retrieved from flotation samples from the two sites; … Franchthi Cave is located in the Argolid peninsula of southern mainland Greece. It was excavated between 1969 and 1976 by T.W. Jacobsen of Indiana University and M.H. Jameson of Pennsylvania University, […] Occupation at the site spans the Upper and Final Palaeolithic, Mesolithic and Neolithic (c. 38,000–6,000 cal BP) […] Shanidar Cave, “located on the western flanks of the Zagros Mountains of Iraqi Kurdistan, was originally excavated between 1951 and 1960 […] Since 2015, a team led by Graeme Barker has conducted systematic excavations at the site (Reynolds et al. 2015), during which the fragments analysed in this study were collected. Five charred plant aggregates were recovered from Upper Palaeolithic (Baradostian) and one further fragment from the Middle Palaeolithic (Mousterian) deposits. Various levels there 43 – 30 kyr BP, 54.4 – 46.05 kyr BP, 75 – 70 kyr BP […] based on their broad stratigraphic association with the well-known Neanderthal flower burial and the recently discovered Shanidar Z articulated skeletal remains, dated to c. 73 kyr BP (Pomeroy et al. 2017, 2020).
    All of the charred food remains were further examined under a Meiji MT6500 darkfield/ brightfield incident light microscope (magnification ×50–500) and subsequently mounted on SEM aluminium stubs and gold sputter coated (to a thickness of 20nμ) to allow for more detailed observation nµ ?? nm, surely? Beyond the Eastern Mediterranean and South-west Asia, archaeobotanical studies at sites such as Niah Cave (Sarawak, Borneo) have revealed evidence for the processing of the highly toxic Dioscorea (yam) and Pangium edule nuts from as early as 50 kyr ago, underscoring the complexity and deep ancestry of such food preparation practices (Barker et al. 2007; Barton et al. 2016). “Niah” rings bells for me. Not hobbits. But … just a few bones, though Palaeolithic.

All, uh, grist to the "Palæolithic Diet" menu's non-existant non-meat part. Not that it ever had any connection ot archaeology.


And that's all I've got time for this (last) month.
Back to List.

2023-07 July Science Readings


2023-07 July Science Readings

Well, I didn't do damned-all last month. Other things.

There seems to be a lot of work on stellar dimmings going on at the moment.

Articles studied this July- some of which might go to Slashdot.
A first great dimming.Betelgeuse Misbehaviour
Beware of Hubris : World Is Not Enough - WINE
Discovery of Gaia17bpp, a Giant Star with The Deepest and Longest Known [Stellar] Dimming Event
Astrophysics observatories as seismographs
Great Dimming But RW Cephei, not Betelgeuse
Potential Supernova Candidates
Exponential distance relation (aka Titius-Bode "Law") in extra solar planetary systems
End of document

Betelgeuse Misbehaviour

The evolutionary stage of Betelgeuse inferred from its pulsation periods

(Further discussion from Molńar, Meridith and Leung appended to original article.)

About a year ago, I posted on the unusual behaviour (misbehaviour, even) of Betelgeuse over the last few years, with the "Great Dimming" superimposed on it's semi-regular brightness variations. My "laughing price" thesis was that by extreme extrapolation of diameter estimates for the star, one could "predict" [sic, my tyops] simplistically, the expected "zero diameter [tyop, "] date range for Betelgeuse is 2030 to about 2040. That was from a 2009 paper, well before the "Great Dimming". Totally unreasonable from any statistical point of view, and I knew it as I calculated it.

So, another paper has come out, also pointing in the same sort of date range, but predicting a supernova not a disappearence up the stellar jacksy. Meh, I can live with that. They've probably approached the question more seriously than I did.

These authors are using a slightly more realistic stellar model than mine, and also new interpretations of old data which suggest (Neuhäuser R., Torres G., et al, 2022, MNRAS, 516, 693) that in the pre-instrumental period, Betelgeuse was yellower (and so, hotter) then than it is today. They put this - the colour change ; the semi-periodic pulsations - together to come up with a revised model for Betelgeuses fundamental parameters thar are a bit different to generally accepted models, but not wildly so. And which, these authors suggest, means that Betelgeuse is teetering on the point of experiencing a core collapse supernova in the imminent future - within a few tens of years.

Well, that'd be fun, for almost everybody. An end most devoutly to be hoped for. And on an "I told you so, young whipper-snapper Grasshopper!" timescale.

Now for the inevitable back-pedalling :

  • Modelling like this is pretty uncertain. Sure, we've now got a mass estimate that is 19 solar masses (M hereafter) compared to previous mass estimates of 11~14 M. That has a lot of consequence for the stars expected main sequence (MS) lifetime (16~9 Myr for the lower mass range, 4Myr for the higher mass range). The post-MS lifetimes also become much shorter as the star mass goes up, and the variability and the dust clouds and mass outflows visible in and around Betelgeuse rather suggest it's already in it's post-MS phase.
  • It's almost unrelated to the "Great Dimming", which is most likely (IMHO) a previous dust emission, thinning and dispersing. That could easily happen on multiple century time scales, and could also affect the stars colour (as discussed by Neuhäuser et al, 2022, marked above). Which would weaken at least one major plank in this paper's estimate of Betelgeuise's mass. Since the invention of the telescope, we may never have seen Betelgeuse/s un-dimmed state.
  • Dr Becky had doubts too, expressed on her "Night Sky News" vlog (timestamps 06:50 ~ 14:30) - but they're pretty similar to mine. I'd add to her "general principles" concerns that we also don't know the metallicity of the core of Betelgeuse terribly well, and thats … somewhere between likely to have an effect, and unlikely to have no effect. We think that the heat transfer from the core outwards is by radiation, not by convection, so the surface (photosphere) composition is decoupled from the core composition. How decoupled ... aye, there's the rub. (Smaller stars than the Sun, "M dwarfs", may be "fully convective", where convection stirs the core material into the main body of the star, but this takes time, and in the latter stages of evolution, things happen in the core on a matter of years, while mixing may be on much longer timescales.
  • She was probably wrong about one thing though - the last supernova in the Milky Way wasn't "Kepler's Supernova" of 1604, but a previously unremarked star in Sagittarius, which went supernova between about 1890 and 1908 producing the radio feature G1.9+0.3. That was invisible from Earth behind dust clouds, but the supernova remenant (SNR) is obvious (and visibly growing) in radio wavelengths. Kepler's was the last naked-eye supernova in the Milky Way. It's very likely that humans will see Betelgeuse going supernova when it goes, because it's so close (450~650 ly, 150~200 pc) that it'll be naked-eye even during daylight hours (unless it happens when in conjunction with the Sun - approximately May to August, depending somewhat on the observer's latitude). Whether it happens before humans lose telescope technology is a more open question.
    Of course, all this attention on Betelgeuse deflects attention from other stars which have a good chance of being the next Milky Way supernova. There's Eta Carinae, trembling away like a Hollywood child actor's lower lip; a list of 31 candidates within 800 ly of Earth (OK, Deneb/ α Cygni is 802 ly, plus or minus error bars), and a lot of other runners and riders. (After writing this, a paper concerning recent strong dimming in the "hypergiant star RW Cephei" has added another, distant, name to the "runners and riders" list. See discussion below.
  • It's almost inevitable in discussions like this that we have to re-stress that all the times mentioned in this topic are in Earth's reference frame, not Betelgeuse's reference frame. Yes, of course the star could have gone SN already, and the light (and neutrinos) from the event are on their way to us. But we won't know that until they get here. making the point moot. Astronomers understand this well enough that since before SN 1987A people have been plotting to observe "light echoes" from SNs (and ordinary novæ) bouncing off dust clouds behind or to one side of the (super)nova, including looking for the arival (at Earth) of light from flares and outbursts of the progenitor star before the "main event". (Eta Carinae has been a notable success in such searches, IIRC. I didn't make notes.)
    (Including this caveat won't mean that nobody re-flogs the bones of this particular flayed horse. Suffer, dead Dobbin, suffer!)

Routinely, when discussing "the big one", along with the dead horse flogging over reference frames (really guys, we get it. Why is it always guys who flog this flensed Equus?) there are questions along the lines of "are we all going to die?"

The answer is "almost certainly not". When we get a "big one", it'll be the 3rd or 4th in a millennium. One going off about every 250 years, that we see. Seriously, outside astronomers, almost nobody notices. Horses chowing grass might riase head from sward for a few seconds, then go back to grazing and worrying about internet trolls and temporal reference frames. Late night taxi drivers might ask "what's that light?". At the right time of the year, daytime people might repeat the question. And almost everyone will forget about it in a few days or weeks. Alart form astronomy nerds, who will be boring people to absolute death about it. Boredom by an astronomy nerd (or bludgoning of the same nerds) is likely to be the biggest death toll consequent on the next "big one".

Specifically, for Betelgeuse, there is one potential "everybody is dead" scenario which is very unlikely. Supernovae remain potential sources of Gamma Ray Bursts ("GRB", possibly also sourced by "compact body" (neutron stars and black holes) mergers, now visible to our gravity-wave telescopes), and one of those going off 600-odd ly away could be a potential bad hair day for a lot of people. Depending on how long it lasts, Homo sapiens mighte be depending on the crews of steel-hulled ships, oil rigs & submarines to bounce back from being hit by such. Which might be game over for the species, and certainly a rough time. But probably not either. Living organisms in the open might die. But all seeds below the soil's surface, nope. Ditto some organisms living in cave entrances facing in the right direction should get enough protection. It's very unlikely to be an "end of life" situation - a mass extinction scenario more likely (hint : 6 of the 7 biggest mass extinctions have obvious non-GRB causes). Even if it happens.

But for that - the Earth to be hit by the beamed energy from a GRB - to happen, we would need to be very closely aligned with the progenitor's rotation axis. Various imaging campaigns from the mid-1990s to the present see noticeable variation in Betelgeuse's brightness from epoch (date of observation) to epoch, which suggest that we're not on the stars rotation axis (the longer-established overall brightness variations argue similarly). A 1998 paper ("Spatially Resolved Hubble Space Telescope Spectra of the Chromosphere of alpha-Orionis". The Astronomical Journal. 116 (5): 2501–2512. Bibcode:1998AJ....116.2501U. doi:10.1086/300596. S2CID 117596395.) suggested that we're about 20° off the star's axis. Which is probably enough to protect us from the "relativistic beaming" of a full-on GRB, and reduce the radiation from a pulsar to about 1 second in 5 - which the atmosphere will reduce further. I just don't see Betelgeuse being much of a (realistic) GRB threat.

Besides, Betelgeuse is probably too small ; if you really want to worry about a GRB, worry about Eta Carinae - but that is around 10 times further away, which factor alone would reduce the radiation dose about 100-fold. And we're 30~45° from the axis of that system. So that's a big fat nothing-burger too.

That's enough Betel-juicing for now.

That was last night

This morning's IArXiv (summary of today's ArXiv postings, weighted by my responses to previous IArXiv listings; though I'm actually somewhat behind on drinking from this particular hosepipe) brings this comment paper : "Comment on the feasibility of carbon burning in Betelgeuse: a response to “The evolutionary stage of Betelgeuse inferred from its pulsation periods”, arXiv:2306.00287", which is various credentialled astronomers responding from Hungary and America. I've dropped the link to Dr Becky, since she's posted on this tool. What are they saying?

These people assert (to be evaluated) that Saio et al's interpretation of Betelgeuse's brightness variation as solely a consequence of pulsations in the diameter of the star is incorrect. "However, the angular diameter measurements of the star are in conflict with the stellar radius required by their models". So, are these objections justified?

I recall thinking, as I read Saio et al's estimate of the radius of Betelgeuse - derived from their pulsation model - as being "higher than normally suggested". But an inflated radius for red giant and red supergiant stars is, as the name suggests, rather normal, so I wasn't wildly upset about it.

From The Friendly Paper (FTFP), on the subject of the star's diameter, "However, measurements do not necessarily detect the visible-light photosphere of the star directly. They are affected by factors like limb darkening, spots, molecular layers and circumstellar dust." Which, yes, that's a large part of why I wasn't terribly upset by the difference between Saio et al's radius and previous estimates. These authors present a considerable compilation of "Visual", "V-band" (an astronomical wavelength-range filter) and "SMEI" ("Solar Mass Ejection Imager" ? Why would a solar telescope be pointed at Betelgeuse? From the first light photos, it can image the Milky Way, but … ) measurements that show Betelgeuse's 380~430 day variation very well, and less distinctly, the 2200 day period. These authors interpret the "Great Dimming" as a sporadic, out of phase event (such as release of a dusty cloud from the star's surface). Which might be an important point, if Saio et al's had relied on it. They mention it, say it was about a minimum of the "~400 day vairation", and then leave it aside as just another event in that cycle.

So basically the dispute here is whether the 380~430 signal is the star's radial pulsation, or if the 2200 day pulsation fits that rôle. The "Great Dimming" doesn't really come into it - it's a random event on top of the informative variation cycles.

The new authors then "work outwards" from the photosphere into interactions of the photosphere and the star's atmosphere. Which they say, prefectly correctly, are likely to have a significant effect on the star's external appearence - as the whole "Great Dimming", and observation of other stars would agree. However ...

A major point of Saio et al's work was to introduce into the discussion "Neuhäuser et al. (2022) …, marked above] interpretation of pre-telescopic records which reveal that this star exhibited significantly higher temperatures [up to] two millennia ago." (as I quoted in the first version of this article). And Molńar, Meridith and Leung don't look at that, they look only at essentially atmospheric effects, while Saio et al's work primarily addresses the internal heat engine that drives the photospheric effects which Molńar, Meridith and Leung are looking at.

Molńar, Meridith and Leung may have valid points about the "atmospheric" effects, but they don't really invalidate Saio et al's discussion of the internal power-source of the star. We're just going to have to wait to find out.

Now, does Dr Becky have owt to say?

See also the discussion of RW Cephei below. I can't get away from the topic.


The World Is Not Enough - WINE - spacecraft ; hubris for those lacking in Classical scholarship.

https://arxiv.org/pdf/2306.03776.pdf

Boringly, the paper title is "Thermal Extraction of Volatiles from Lunar and Asteroid Regolith in Axisymmetric Crank-Nicholson Modeling", but I just love the idea of "the World Is Not Enough (WINE) spacecraft concept". Clearly someone has either not heard of the concept of hubris, and how various pantheons of gods have punished humans showing this characteristic. The designers will have to be tunr widdershins at the right moment, and perform their functional sacrifices with appropriate chanting burning of incense and prayer to the Gods of spaceflight.

What are they actually doing? Well, with Philip T. Metzger (on Twitter, @DrPhilTill on Twitter, while it's still alive) involved, you can guess its' something to do with the regolith on space bodies. In this specific case, they'ere looking at the (theoretical) efficiency of extracting volatiles (in general, water in particular) from regoliths on the Moon, asteroids and (potentially) Mars.

Previous work has shown (unsurprisingly) that in a vacuum, as you heat a sample of regolith the volatiles move away from the heat source, but the details depend on the vacuum, the temperatures, temperature gradients and many other details. A particular construction of mining machine could collect voatiles efficiently in one nevironment, and very inefficiently in another. So, there is a need to do multiple tests in various environments. Which is where the WINE (hubristic name expanded above) spacecraft comes in. In their words, WINE will be a small spacecraft, approximately 27U in CubeSat dimensions (3 by 3 by 3 cubes), with legs for walking short distances and a steam propulsion system for hopping multiple kilometers. […] WINE will drive a corer into the regolith to extract water and perform science and prospecting measurements on the regolith. The water will be extracted thermally by heating the material in the corer. Vapor will travel into a collection chamber where it is frozen onto a cold finger. After multiple coring operations have collected enough water, the tank will be heated to high pressure and vented through a nozzle to produce hopping thrust Which is a fairly ballsy set of optimisms.

They then follow with a lot of mathematical modelling, which really just serves to show that we don't have a large amount of experimental data on physical and chemical properties of regoliths, apart from a small amount of Apollo data density/ strength obtained by using core penetrometers (tube, hammer, measuring stick, notepaper) to collect soil samples. Some properties (thermal conductivity, density) were measured after several days transport and re-warming in the returning spacecraft, plus a sea landing. Yeah, I'm not going to put too much faith on those values without some "ground truthing" measurements. (These authors don't mention it, but if the "mole" of Mars Insight had worked, it would have provided at least one log of some of these parameters in Martian soil. possibly including sub-surface ice.)

The table of products form the LCROSS impact mission is interesting.

Volatiles in LCROSS Ejecta (from https://arxiv.org/pdf/2306.03776.pdf)
Compound Symbol Concentration (wt%)
Water H20 5.50
Hydrogen sulfide H2S 1.73
Sulfur dioxide Si [sic, typo] SO2 0.61
Ammonia NH3 0.32
Carbon dioxide CO2 0.29
Ethylene C2H4 0.27
Methanol CH3OH 0.15
Methane CH4 0.03
Hydroxyl OH 0.0017
Carbon monoxide CO 0.000003
Calcium Ca 0.0000008
Hydrogen gas H2 0.0000007
Mercury Hg 0.0000006
Magnesium Mg 0.0000002

I hadn't particularly paid attention to that previously. But having 1.73% of the ejecta mass as hydrogen sulphide is ... challenging. If you were to pump this mix into a chemical plant to make atmosphere, rocket fuel, whatever, you're going to want to either dump that stuff, or move the sulphur to a differentl, less hazardous, oxidation state. The simple process of oxidising it to sulphurous acid (H2SO2 or 3, which is relatively innocuous) would help, but would account for something like 1/3 of the potential oxygen output from electrolysing the water. Hmmm. That's a problem. Maybe just bring the sulphur up from oxidation state -1 to state 0 would be better. "Flowers of sulphur" is a damned sight less horrible a chemical to have near living quarters then H2S. Someone is going to have to look at that.

There's a lot more work in this paper, modelling temperature changes in lunar soils through a day-night cycle, considering temperature profiles in the corer of the "WINE" spacecraft proposal. Lots of theory. Really, some of this needs ground-truth, and I'm sure this paper will be fairly heavily drawn-upon in arguments over which next missions to launch to the Moon (or various asteroids, or Mars) to evaluate their soils physical and chemical properties. There are a lot of free parameters in these models which need pinning down to something more closely constrained by reality.

Not a very "fun" paper, but probably quite important for actually building humanity into a space-dwelling civilisation, if not necessarily a Mars dwelling one. Who would want to live at the bottom of a hole? (Known Space reference.


Stars with really long dimming events

Discovery of Gaia17bpp, a Giant Star with the Deepest and Longest Known Dimming Event

Do you remember where you were when Tabetha Boyajian said "Where's The FLux?" about star KIC 8462852 (KIC - Kepler Input Catalog)? That was a lot of a surprise (and still remains somewhat puzzling), but with the rapidly increasing completeness, cadence and sensitivity of astronomical survey programmes the number of - and extremes of magnitude and duration of - dimmings of stars recorded have beceom a commoner thing.

This is a new record-breaker. Until the next record-breaker. This star has spent about 6.5 years (2000-some days) in it's dimmed state - the start is somewhat uncertain between surveys which serendipitously covered this field in 2012 and 2013 - before which the star exhibited no obvious brightness changes. However with up to 11 or 12 year gaps in observation in the 1950s and 60s, other dimmings are plausible.

The interpretation - as for a number of other dimming events - is that the star's light has been dimmed by passing through an approximately neutrally-coloured disk of material orbiting a companion (or possibly independent, if this dimming is a one-off event). The modelled disc is at a small angle to our line of sight - if the disc had been inclined closer to face-on, the depth of dimming would have been less.

Similar explanations have been put forward for other "dimming" stars such as Epsilon Aurigae, though it's dimming recurs on a period of some 27 years, with each dimming lasting around 2 years.

We're going to see more of these. Lots more of these. Following the same arguments as justified the Kepler Observatory, around 1% of multiple star systems where one component has a significant absoebing disc will be detectable from Earth, given sufficient numbers of observations. Similar biases to close-in dust disk carrying companions will exist.

This model doesn't explain "Tabby's Star" KIC 8462852 - or at best, not very well. There are hints of periodicity in that one's dimmings, but there is a lot more stochastic noise.


Astrophysics observatories as Seismographs

Status of the GINGER (Gyroscopes IN GEneral Relativity) project.

After the successes of gravity-wave telescopes (LIGO, to a lesser extent KANGA and Virgo ; and now the various *PTA "NANAOgrav discoveries.), people are naturally looking to turn it up to 11. s. GINGER (Gyroscopes IN GEneral Relativity) is a proposal to build a number of ring laser gyroscopes (a technology that has been used in surveying since almost the invention of the lasers) to "measure general relativity effects and Lorentz Violation in the gravity sector". Which I'm sure would be terribly exciting for physicists and astronomers everywhere. But the other stipulation, that the devices be "rigidly connected to the Earth", also means that these would both be subject to non-trivial levels of seismic noise, and would also serve as a remarkably broadband seismograph. The seismic noise they've got a technique for, as deployed by LIGO Hanford and LIGO Livingstone, by having two (or more) more-or-less identical machines at considerable separation, and having some degree of appropriate coincidence detection and rejection to differentiate between GW signals and seismic signals.

That is going to be more challenging with two sensors at one location.

schematic of two rectangular frames inclined with respect to each other ; human-for-scale indicates the frames to be about two person-heights on edge, roughly square

(The proximal part of "GINGER", time wise, is about building two "technology demonstrators" in the existing deep underground lab at Grand Sasso, which is currently used for low-radiation-background experiments, neutrino detection, dark matter telescopes, that sort of thing.) They'd have to interface with more conventional seismic networks to reject the "understood" noise sources, which would leave the not-understood signals in the noise. Or for that matter, in the not-previously-"understood" signal. But I'm sure that's do-able, before building a second site. Or fully interfacing with one of the various other high-sensitivity ring gyro systems under development.

That's viewing the seismic noise as noise. But of course, it's also going to be an exquisitely sensitive seismograph. If they're looking at sensitivities in the order of 1 part in 109 to 1011 that would imply a bandwidth down to around millisecond or microsecond signals (since they're looking towards astronomical sources from an Earht that rotates in approximately 105 s. That's a pretty good step up for seismography. And the physicists pay for it! What's not to like?

The project already has, they report (references in original paper), collapborations with high-precision observatories in Germany (G9 at Wettzell and ROMY in Bavaria), New Zealand (U.o.Canturbury) and China (HUST, Wuhan)). That Wuhan connection is going to bring the anti-vaxx freaks out of the woodwork, though I'd be moderately amused to see how they get from Length-of-Day and Earthquake measurements to Bill Gates' 5G mind control chips in the COVID vaccines. After all, it can't be coincidence. Those pre-existing collaborations are proof positive of nefarious and long-standing plans.

I see, checking those co-project's references, that several of the paper's titles include reference to siesmic and seismological results. Evidently they too have seen this ... interdisciplinary link. So they're probably touching the seismology community for some contributions of "folding" already.

As a footnote, this is a contender for the DOOFAAS project. If you like it, mail that project's instigator. It is not known if ACRONYM software was used to construct this ... "name". Or if it is natural talent on someone's behalf.


Another Great Dimming - not Betelgeuse

The Great Dimming of the hypergiant star RW Cephei

It seems as if this month's recurring theme is the dimming of Beteleuse, but many other intrinsically bright stars also have irregular dimmings, as well as more regular periodicities in brightness. RW Cephei is another such. Although it is a lot brighter intrinsically than Betelgeuse, it's not so well known because it's further away. We can't actually tell how far away to within a factor of 2 - estimates on Wiki give it at 3,416 or 6,666 pc (11,100 - 21700 ly) - an eighth to a quarter of the way around the Galaxy. It's not clear why the distance is uncertain, but the presence of an extended envelope of emissions around the star probably doesn't help. That's 10 to 20 times the distance to Betelgeuse, so it's not a bright star in our sky.

Cepheus is a circumpolar constellation from the UK, but it's not terribly well known because it's not got many bright stars. As seen on the sky, it is to the west (anticlockwise, widdershins, centred on Polaris) from the bright "W" asterism of Cassiopeia. RW Cephei is in the Cassiopeia corner of the constellation, though actually closer to the border with the dim constellation of Lacerta. It's well below naked-eye visibility during the current dimming, but at it's brightest should just be visible under excellent conditions.

Despite being a lot further away from us than Betelgeuse (which has been subject to interferometric disc-size measurement since the 1920s, and more recent actual disc imaging), RW Cephei is such a large star thet it too has a disc capable of being imaged. This paper reports on the imaging results. THe observations were carried out with the CHARA telescope array, comprised of six 1m telescopes with delay lines so that act as an interferometer with a resolving power of 200 µas. These results couples with the distance estimates given above indicate a stellar radius of 900 to 1760  R, making it one of the largest stars measured, to date.

The report also includes AAVSO data for RW Cephei's visual brightness (with Betelgeuse for comparison) as figure 1 :

Caption :

If you want, you can get updated data from here. At the report's writing, the dimming event had lasted about 1100 days (compared to about 200 days for Betelgeuse's recent dimming and 2000+ days for Gaia17bpp described above). At points in the deepest part of the dimming, the magnitude was being reported as 8 (visual), compared to the normal range of brightness of 6.0 to 7.6 . The most recent (late July 2023) AAVSO data has the brightness back to about M 7.0, well within the normal range of variation.

But that's not the exciting thing in this report. What this has, which most reports like this don't have, is that with the telescopes set up for interferometry, the observers could map variations in the brightness of the star's surface and reconstruct the image of the star's surface. Which has been done before (with Betelgeuse, for example). But this is on a star an eighth to a quarter of the way across the galaxy, and a considerably more luminous star too. RW Cephei is about 300,000 times the Sun's luminosity (300,000 L) while Betelgeuse is less than half that (90000 - 105000 L.

The first level of reconstruction produces these images (fig 5 upper) :

But those images combine both changes in colour and intensity, which include the effects of limb darkening. If you account for that effect, you get a better impression of what is happening to the star itself. (fig 5 lower, and caption)

The colour and/ or shape variations of these images probably reflect variations in the temperature of the star's surface, which for a large star is thought to be dominated by convection cells bringing heat from the interior to the surface. But the largest stars have a competing process of mass loss with both material jetting off the surface into space (which implies cooling, with spectroscopic consequences, which can be seen) and also the release of dust to form obscuring clouds between us and the staR - a process thought to be happening with Betelgeuse.

Section 3 of the paper reports near-infrared spectroscopy taken during the dimming. Through most of the spectrum the intensity is lower than in archived (normal brightness) IR spectra, but the change is greater at shorter wavelengths. With some well established physics this can be turned into mean surface temperatures


And now the end-of-month tidy-up. Only 3 overhangs.


Red Supergiant Candidates for Multimessenger Monitoring of the Next Galactic Supernova

https://arxiv.org/pdf/2307.08785.pdf

Science wise, this is just a catalogue of these candidates. But it’s also a public bet. 677 public bets. The justification is to shorten the list of candidates for pointing “light buckets” at, in the event of an early detection of a likely supernova by neutrino telescopy. If (“if”) one or several of the operating neutrino telescopes (French ANTARES, Antarctican ICECUBE, Russian BDUNT, Japan’s SuperKamiokaNDE …) detect a burst of neutrinos incoming, and can get a directional fix (some are more directional than others) then we may have an approximate detection.

The neutrino burst from a star's core collapse (or potentially, a gravitational wave signal) should escape from the core of a star in seconds (about 5 seconds for the Sun) to the surface where we can see it, while the explosion shock wave can take hours to days to emerge (“shock breakout”). Thus, a prompt observing campaign could capture brightness, spectroscopic or even compositional data on the star as close as possible before the actual supernova. Hence, a catalogue of “usual suspects” could help with prioritising candidates. It’s a numbers game, but with a reasonable chance of prompt payoff. The 1987A supernova in the Large Magellanic Cloud was accompanied by a burst of 20 (or maybe 25 - 1 of the 4 detector instruments involved is discordant) neutrinos and antineutrinos 2 to 3 hours before the first optical detection of the star's brightening. That's a very useful amount of warning.

Probably major observatories are already incorporating this into their “Target Of Opportunity” (TOO) decision process. To a significant degree, time and date are going to affect each observatory’s listing, day by day and hour by hour, but also which instruments are deployed, the slewing time from [current target] to [TOO] to [next target] … it’s not something the night shift operators would want to have to do “on the fly” without prior planning.

The criteria used for inclusion in the catalogue are luminosity (absolute, obviously, which requires a reasonably good distance from the Gaia stellar catalogue), temperature (requiring at least multiple filter measurements, if not full-blown spectroscopy, so not all potential targets have this data available, to this date ; that'll change), and chemical contamination (also requiring spectroscopy).

Obviously there’s a lot more detail in the selection process. But the core idea is there.

The catalogue is online on GitHub , but is also already incorporated into VizieR, a standard astronomical database collection. (I need to improve my VizieR-fu! Ditto the “Jupyter notebook” mentioned on GitHub, which I’ve heard of but never needed to use.)

I’d already got a catalogue of 112 RSGs and B binaries in my astronomy workbook, so I’ve added this lot. Plus a list of 31 potential progenitors from https://arxiv.org/pdf/2004.02045.pdf - which list is of suspect stars within a kiloparsec (3200 ly, approx) whereas this more spectroscopic list has nothing closer than about 250 pc (it gives parallaxes in mas, not pc). The lists use different sets of names, so are hard to compare directly, but at least three stars (HD 17958 (HR 861 in Cassiopeia, Gaia DR3 ID 467907038749283000, DR2 ID 467907038747132000, 2MASS ID J02562466+6419563), HD 80108 (HR 3692 in Vela, DR3 ID 5423960064637140000, DR2 ID 5423960064637140000, 2MASS ID J09162303-4415564) and HD 205349 (HR 8248 in Cygnus, DR3 ID 1971358279140130000, DR2 ID 1971358279140130000, 2MASS ID J21331788+4551142 ) appear on both lists.

Ye gods, and by gods I mean Cthulhu and the FSM , I'd forgotten how bad star names were. It was bad enough with AB Doradus and ZZ UMajor a decade ago, but these Gaia ones are a factor of several worse. If only they were comprehensive catalogues. But that's not really feasible this side of 102,023 CE. Sorry, 1,102,023 CE, because we'd need to geet to the other end of the galaxy. and get the signal back. That's assuming an average human expansion rate of 0.1C, and a speciation rate of zero - both of which are rather implausible.

If you can find any more "duplicates", you're welcome to. And you deserve a chocolate biscuit.


Exponential distance relation (aka Titius-Bode "Law") in extra solar planetary systems

https://arxiv.org/abs/2307.06070

Titus and Bode’s “Law” was an exercise in applied numerology which became popular once the geometry of the Solar system became calculable, using Newton’s dynamics, plus measurements of the Astronomical Unit (AU) by observation of parallaxes, and specifically transits of Venus across the Sun’s disc. Per Wiki, the first statement was in 1715, with variations and various authors (including Titus and Bode, and others) until 1772. When Uranus was discovered in 1781, which fitted reasonably closely to the “Law”, it became more popular. And when Ceres (now classified as a dwarf planet) was discovered in 1801, further numerology ensued. And again, with the discovery of Neptune, in “not-quite” the right place, more flogging of the never-alive horse happened. Pluto, also not being in quite the right place, further encouraged the numerologists. With 3 (or 4) parameters (for different ways of expressing the idea) and 7, 8, or 9 data points (for planets, dwarf or not), anyone could join in the fun of trying to make a better-fitting model. However, no “universal” solution of free parameters has been found – every system has different parameters.

If the Solar system had formed by natural processes (potentially a question in the 1700s, not so these-centuries) then it’s appealing to expect that there should be some relationship between the orbits of the planets. Unfortunately, “appealing to shaved apes on a mud-ball” isn’t a particularly compelling argument to planets in other stellar systems or scientists anywhere. It would be harsh to call it “not even wrong”, but it’s not very well founded in physical reality for a number of reasons :

  • as generally presented, the Titus-Bode “Law” looks for relation in the semi-major axes of bodies orbits ; but basic Newtonian gravitation (Einsteinian too, not that it matters) shows that the biggest forces between planets would occur when the outer planet is at perihelion (nearest point in the orbit to the Sun) and the inner planet at aphelion (furthest point from the Sun), when they’re in the appropriate phase relationship. And on a significant time scale (I’m a geologist – the mega-year is a convenient unit) phase relationships do change.
  • orbits evolve with time. They evolve significantly. Today’s measurements of orbital parameters are good for predicting what they will be next year, and not bad for predicting a million years hence. But a billion years … nope, you’ve got no real choice but to calculate it numerically because those 10th decimal places really do add up over time. If you do a lot of calculation of models there is, for example, about a 1% chance that Earth will be hit by Mercury before the Sun goes red giant. That’s a consequence of Mercury – the second least circular planetary orbit, if you include Pluto as a planet – occasionally approaching Venus relatively closely, which increases Mercury’s eccentricity leading to closer approaches … and 1% of the time, Mercury gets ejected from the Solar system, but finds Earth in the right place at the wrong time. (Don’t worry, the Earth’s oceans boil and the Sun goes red giant with pretty much 100% probability. Even if another star hits the Sun (less than a billion-to one chance, even when the Andromeda galaxy hits the Milky Way), that’s only going to bring the red giant date forward.)
  • there seems to be quite a lot of randomness in the development of planets. If our best model of planetary formation is right, little things (dust grains, sand grains, dirt balls, asteroids, big asteroids) meet to form bigger things until they run out of bigger things crossing their orbits. At all stages similarly-sized bodies are colliding and merging– which is a sensitively random process. When people tried modelling the “Moon-forming impact” in the late 1980s (when supercomputer power became affordable), they discovered that it was exquisitely sensitive to the impact factor (distance between proto-Earth centre and projected path of the impactor), and to the rotation rate of both bodies. I’ve got a compendium of evidences for “late giant impacts” in the Solar system, and the only planet without good reasons for believing there to have been a “late giant impact” is Saturn – the one with the recent debris rings and the “Death Star Moon” Mimas
    Wikimedia image of satellite Mimas with crater Herschel, almost a quarter of the satellite's diameter, resembling the Death Star of Star Wars
    ). Yeah, not a lot of evidence for giant impacts and hierarchical growth around Saturn.
  • Would you really expect a Titus-Bode like "Law" in the satellites of a planet within a stellar system? That’s quite unintuitive to me – surely the gravitational influences of the other planets would be have an effect, particularly on the outermost satellites where the statistically greatest influence on a “Law” would happen, per satellite. You could view the satellites of a planet in a multi-planet system as being somewhat similar to the planets orbiting one star in a (widely-separated) multiple-star system.

So, yeah, it’s maybe not so obvious that there should be a Titus-Bode-like "Law" for planetary systems. Certainly not to me.

FTFA, the authors find that an exponential (“Titus-Bode-like”) model fits the 32 reasonably-sized (i.e. 5 or more planets) planetary systems known to date. They use “R²” and “Median” comparison statistics – (I know “R²”- Pearson’s correlation coefficient. Their “Median” test test though is the https://en.m.wikipedia.org/wiki/Average_absolute_deviation which is a moderately well-known description of clustering in a data set around it’s central point. Here they’re using the mean as a measure of central tendency: ) Mean Average Error ( MAE ) = 1 N i | y i y i ^ | Mean Average Error (MAE) = 1 over N sum from{i} abs{ {y_i} - {hat{y}_i} } (y-hat being the sample mean).

Historically, the Titus-Bode “Law” expression started from Earth’s orbit, and went both in- and out- wards. These authors (like many others) simplify matters by working from the smallest known orbit and just working outwards. Which makes for a simpler expression :

r ( n ) = ae 2 λ n n = n n = 1,2,3 r(n) = ae^{2 %lambda n} phantom { n = n } n = 1,2,3 dotslow

The authors are careful to exclude multi-star systems from their consideration. That suggests (to me) that the effects of the lighter star being closer to the inner planets some times - but not others - would be an issue.

By generating naïve artificial planetary systems and analysing them similarly to the actual star systems they get average regression values of 0.905 while for the Kepler poly-planet set it’s 0.966 – appreciably better. It’s unsettling that they find the (generalised) Titus-Bode "Law" to be a good descriptor of planetary systems, when I’m so disparaging of it. But they agree with me that the lack of a physical basis is problematical.

An interesting sideline is that by introducing a physically realistic constraint on placing of their “pseudo planets” in their pseudo-systems (separating planets by several Hill radii https://en.m.wikipedia.org/wiki/Hill_sphere ), their pseudo-systems approach the statistics for Titus-Bode models. That’s … suggestive. But far from convincing.

The authors spend some pages discussing the use of the Titus-Bode "Law" (originally, and in more recent formulations) to predict the positions (well, semi-major axes ; equivalently periods) of as-yet undiscovered planets. As always, the example of the discovery of Ceres is touted, and indeed, the search that yielded Ceres does actually fit the original Titus-Bode "Law". The surprising thing is that the search didn’t find Neptune, with over three times the angular diameter (though something like 1/10th of the overall brightness, so maybe not that surprising). However, Ceres itself, and the whole asteroid belt adds up to really, really little. People get this wrong. They think the asteroid belt is much more significant than it really is. By mass it is about 0.004 that of Mars (proportionally less than half the size of Earth's Moon), and 0.0004 (count the zeros! I did.) of the mass of Earth. The lack of mass in the asteroid belt has long been a problem for people trying to work out a reasonable scheme for building the Solar system from a circumstellar disc (such as we see around protostars and young stars today). Considering that, really we should do calculations about the Titus-Bode “Law” excluding both Ceres and Pluto, since they’re so small and almost massless that they can’t greatly influence anything of significant size at significant range.

While I’m worrying about Ceres’ mass, the next section of the paper comes along … suggesting the Titus-Bode “Law” may in some way be related to the age of a planetary system, as if it had been “running in” like a coarsely-machined engine. Problematically for this thesis, what you need for shuffling the components of a planetary system around is, bluntly, mass. But by including Ceres in their considerations as a “planet”, they’re also saying that, in their opinion, mass is unimportant. Fortunately, their correlation is no correlation. “Correlation is not causality” may be a mantra of statistics lecturers (even if the students tend to forget it), but “No correlation means no causality” is a considerably more firmly based claim. The authors dug this pit for themselves, then chose to not climb into it. Well done.

But then the next section they choose to throw themselves back into this pit. After a not-unreasonable discussion of the mutual influence of neighbouring planets on each other’s periods and spacing (see above comments about mass), they go on to say It is in fact well-known that “small” objects, like comets, small asteroids, debris, etc. do not follow, singularly, any particular pattern in the size of their orbits, and can be found at any distance (allowed by classical mechanics) from the central body. and then they go on to consider Ceres to be significant in the dynamics of the Solar system. Dodging a bullet is a skill, but getting hit by it after dodging it is a really difficult skill.

There may be something to the “Harmonic Resonances (HR) method”. But not by using this argument. Not if they’re considering Ceres to be throwing Mars around. Their subsequent point is that the "HR method" generally uses distinct “small number” ratios between neighbouring planetary neighbours, so the system effectively has approximately two free parameters per planet – which would practically guarantee a fairly good fit. That's well-founded. So, why pay much attention to the "method"?

Apparently the journal Icarus refuses to accept papers purporting to improve the Titus-Bode “Law”. I can understand why. Archaeology journals probably don’t accept many papers on von Daniken blurb or, latterly, Hancockian ink-waste. There may be something behind the (“a”, generalised) Titus-Bode “Law”, but it still has no theoretical basis, and it is (IMO) very problematic that real planetary systems (at least, the Solar system, and several young extra-Solar systems which have large excesses of dust) seem to have had significant but essentially stochastic (random) events. You might end up with a Titus-Bode-like "Law" from some filtering or damping of an initially random system, but that is a connection that hasn’t yet, to my knowledge, been made. And this paper doesn’t do it either. It’s probably worthwhile trying it, but this sort of numerology isn’t likely to do it. Numerical modelling of evolving systems in silicio is more likely (IMHO) to get somewhere in the right direction.

One thing that worries me is how all this depends on counting planets (gravitationally significant bodies) out from the star. If you’ve got a “hot Jupiter” in a 3-day orbit so close to the star that it's tidally locked, there probably isn’t room for a stable inner un-detected planet. But in the Solar system the innermost planet is also the smallest, and therefore the hardest to detect (by most methods). So surely, when looking at “adding planets” to an observationally based model, the case of an undiscovered small innermost planet is one that should be considered before messing around adding planets further out.

And no, Ceres still is not a gravitationally significant planet. Pluto is over twice the diameter and 13-odd times the mass, and I’m not screaming for it’s re-inclusion as a planet. (For the record, I’d have gone for a gravitational rounding + orbiting the star(s) definition rather than the IAU's “clearing the orbit” method. But I’m a geologist not an astronomer. And I’m definitely not going to get into that fight.)


And finally I'm caught up!

End of Document
Back to List.

Popular Posts