Variable magnetic field and adaptive mixing-length: reproducing Li abundances and constraining rotational evolution of solar-type stars in clusters

This study employs rotating stellar models with dynamically varying magnetic field strength and mixing-length parameters to successfully reproduce observed lithium abundances and general rotational trends in solar-type stars, though it currently overestimates the Sun's present-day rotation rate and magnetic field, highlighting the need for additional angular momentum loss mechanisms.

The Gist

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Picture: Why Do Stars "Forget" Their Lithium?

Imagine you have a cup of coffee with a little bit of sugar (Lithium) at the bottom. If you stir the coffee vigorously, the sugar dissolves and mixes everywhere. If you leave it alone, the sugar stays at the bottom.

In stars, Lithium is like that sugar. It gets destroyed (burned up) if it gets mixed down deep into the hot core of the star. Astronomers have been puzzled for decades by a mystery: Why does our Sun have so little Lithium left, while many other young stars of the same age have a lot more?

This paper is like a team of astrophysicists trying to fix the "recipe" for how stars spin and mix to finally solve this mystery.

The Old Recipe vs. The New Recipe

For a long time, scientists modeled stars using a "static" recipe. They assumed two things never changed during a star's life:

1. The Magnetic Field: They assumed the star's magnetic field was a fixed, constant strength (like a lightbulb that never changes brightness).
2. The Mixing Efficiency: They used a fixed number to describe how well the star's outer layers churned and mixed (like a blender set to a single, unchangeable speed).

The Problem: Stars aren't static. They spin, they shrink, they grow, and they change their magnetic fields as they age. Using a fixed recipe was like trying to bake a cake by assuming the oven temperature never changes, even though the cake is rising and cooling down.

The New Approach:
The authors of this paper decided to make their model adaptive. They introduced two "smart" variables that change dynamically as the star evolves:

• Variable Magnetic Field: Instead of a fixed lightbulb, they modeled a magnetic field that acts like a dimmer switch. It gets very strong when the star spins fast (when it's young) and gets weaker as the star slows down (when it's old).
• Adaptive Mixing: Instead of a fixed blender speed, they let the "mixing length" (how much the star churns) change based on the star's temperature and size, just like a real blender might adjust its speed depending on how thick the mixture is.

How They Did It: The "Star Simulator"

The team used a super-computer program called MESA (which is like a video game engine for stars, but much more serious). They ran simulations of 1-solar-mass stars (stars exactly like our Sun) from their birth (as a baby star) all the way to their current age (4.57 billion years).

They compared their new "adaptive" models against real data from:

• The Sun: We know exactly how much Lithium it has and how fast it spins.
• 64 Open Clusters: These are groups of stars born at the same time, acting like a "time machine." By looking at young clusters, we see stars at the Sun's "childhood." By looking at old clusters, we see stars at the Sun's "old age."

The Results: A Mixed Bag of Success

The Good News (The Lithium Puzzle Solved!):
The new model was a huge success at predicting Lithium.

• The model predicted that a Sun-like star should have a Lithium abundance of 1.12.
• The actual measured value for our Sun is 1.1.
• The Analogy: It's like trying to hit a bullseye on a dartboard. The old models were throwing darts that landed on the outer rings. The new model threw the dart and hit the bullseye almost perfectly. By letting the magnetic field and mixing change over time, they finally figured out exactly how much Lithium gets burned up.

The Bad News (The Spin Problem):
While they nailed the Lithium, they missed the mark on rotation speed.

• The model predicted the Sun should be spinning at 4.72 km/s.
• The real Sun is spinning at 2.0 km/s.
• The Analogy: Imagine you built a car that gets perfect gas mileage (Lithium), but the speedometer says it's going 100 mph when it's actually only going 40 mph. The car is spinning too fast in the simulation.

The Magnetic Mystery:
The model also predicted the Sun should have a magnetic field of about 37 Gauss, but the real Sun averages only 1 Gauss.

• Why? The model assumes the magnetic field is generated purely by the star's spin. It seems to be missing a "brake" or a "leak" that slows the spin down faster in reality. The authors suspect they are missing a mechanism called "disk locking" (where a baby star is held back by the disk of gas and dust it was born from) or some other internal magnetic coupling that slows the spin down more efficiently than their model accounts for.

The Takeaway

Think of this paper as a major upgrade to the "Star Operating System."

1. We fixed the Lithium bug: By making the magnetic field and mixing "smart" and changeable, the computer finally understands why the Sun has so little Lithium.
2. We found a new glitch: The computer thinks the Sun is spinning too fast and has too much magnetic power. This tells us there is still a missing piece of the puzzle—likely something that happens when the star is a baby (like being held back by a cosmic leash).

In short: The authors successfully updated the rules of the game to explain the "sugar" (Lithium) in the coffee, but they realized they still need to figure out exactly how the "stirring" (rotation) slows down so perfectly over billions of years. This sets the stage for the next generation of star models to include even more physics, like the interaction between baby stars and their birth disks.

Technical Summary

Here is a detailed technical summary of the paper "Variable magnetic field and adaptive mixing-length: reproducing Li abundances and constraining rotational evolution of solar-type stars in clusters."

1. Problem Statement

The paper addresses two persistent anomalies in the modeling of solar-type stars:

1. Lithium (Li) Depletion: Standard stellar evolution models struggle to simultaneously reproduce the observed surface lithium abundance of the Sun ($A(\text{Li}) \approx 1.1$ dex) and the rotational evolution of stars in open clusters. Specifically, models often fail to predict the correct amount of Li depletion observed in young clusters and the Sun.
2. Fixed Physical Parameters: Traditional models rely on fixed values for critical free parameters, specifically the magnetic field strength ($B$) and the mixing-length parameter ($\alpha_{\text{MLT}}$) used in Mixing-Length Theory (MLT) to describe convection. The authors argue that assuming these parameters remain constant throughout a star's evolution is an unrealistic simplification that fails to capture the dynamic interplay between rotation, magnetic activity, and internal mixing.

2. Methodology

The authors developed a new grid of rotating stellar evolution models using MESA (Modules for Experiments in Stellar Astrophysics). The core innovation lies in making two key parameters time-dependent rather than fixed:

• Adaptive Mixing-Length ($\alpha_{\text{MLT}}$): Instead of a fixed value (e.g., 1.82), $\alpha_{\text{MLT}}$ is calculated dynamically based on the star's effective temperature ($T_{\text{eff}}$) and surface gravity ($\log g$). This follows a calibration derived from 3D hydrodynamic atmosphere simulations (Sonoi et al. 2019). This approach links convective efficiency directly to the star's thermal structure, which evolves with rotation and age.
• Variable Magnetic Field ($B$): The magnetic field strength is not a preset constant. It is calculated using a dynamo prescription (Gallet & Bouvier 2013) where $B$ depends on the Rossby number ($R_o$), which is a function of the rotation period and convective turnover time ($\tau_{\text{conv}}$). This allows $B$ to evolve naturally as the star spins down and its internal structure changes.
• Magnetic Braking (MB): The models implement a magnetic braking routine that removes angular momentum ($J$) based on the variable magnetic field and mass loss rates. The torque is calculated using the Alfvén radius, which scales with the magnetic field strength.
• Rotational Instabilities: The models include various rotation-induced instabilities (Solberg-Hoiland, Eddington-Sweet, Goldreich-Schubert-Fricke, Secular Shear, and Dynamical Shear) to handle angular momentum transport, while explicitly excluding the Tayler-Spruit dynamo to isolate the effects of the new prescriptions.
• Data Validation: The models were tested against observational data from 64 Open Clusters (OCs) obtained via the Gaia-ESO Survey (GES) and Gaia data releases (GDR3). The sample focused on solar analogues ($1 M_{\odot}$) across a wide age range (from pre-main sequence to the current solar age).

3. Key Contributions

• Dynamic Parameterization: The study successfully implements a framework where both $\alpha_{\text{MLT}}$ and $B$ are free, time-dependent variables derived from stellar state parameters, reducing the number of arbitrary free parameters in stellar modeling.
• Reproduction of Solar Li Abundance: The models successfully reproduce the Sun's current lithium abundance ($A(\text{Li}) \approx 1.12$ dex) within observational uncertainties ($1.1 \pm 0.1$ dex), a significant improvement over models with fixed parameters.
• Physical Consistency: The study demonstrates that the variability of $\alpha_{\text{MLT}}$ during the Pre-Main Sequence (PMS) phase is crucial. The adaptive $\alpha_{\text{MLT}}$ is lower during early phases (leading to less Li destruction) and converges to the solar-calibrated value ($\approx 1.76-1.78$) at the present age.
• Comparison with Previous Work: The paper provides a rigorous comparison with previous studies (Caballero et al. 2020; Gossage et al. 2021), highlighting how variable parameters alter the predicted Li depletion and rotational spin-down trends compared to fixed-parameter models.

4. Results

• Lithium Abundance: The models with variable $B$ and $\alpha_{\text{MLT}}$ successfully match the observed Li abundance of the Sun and solar-like stars in clusters. The adaptive mixing length reduces Li destruction during the PMS, allowing the star to retain enough Li to match solar observations despite high initial rotation rates.
• Rotational Evolution Discrepancy: While the models reproduce Li abundances, they over-predict the current equatorial rotation speed of the Sun.
  • Model Prediction: $v \approx 4.72$ km/s.
  • Observed Solar Value: $v \approx 2.0$ km/s.
  • Model Prediction: Mean surface magnetic field $B \approx 36.9$ G.
  • Observed Solar Value: $B \approx 1$ G.
• Magnetic Field Evolution: The models predict magnetic field strengths reaching $\sim 1$ kG during the PMS (consistent with T Tauri stars) and declining to $\sim 40$ G at the solar age. This is significantly higher than the observed solar mean field, suggesting the magnetic braking efficiency or saturation physics in the current formalism is too strong or lacks a saturation mechanism.
• Cluster Trends: The models qualitatively reproduce the rotational spin-down trends across different cluster ages but fail to quantitatively match the tight rotation-Li relation in all clusters without additional physics.

5. Significance and Limitations

• Significance: The paper establishes that variable mixing efficiency is a critical, previously underutilized factor in explaining lithium depletion. It proves that a "one-size-fits-all" $\alpha_{\text{MLT}}$ is insufficient for solar-type stars. The approach reduces the reliance on ad-hoc tuning, moving toward a more physically grounded stellar evolution theory.
• Limitations:
  • Missing Physics: The discrepancy in the Sun's current rotation rate and magnetic field strength suggests the models are missing key angular momentum loss mechanisms, specifically disk locking (during the PMS) and internal magnetic coupling (between the radiative core and convective envelope).
  • Magnetic Saturation: The current dynamo prescription may overestimate magnetic field strengths in the saturation regime for slow rotators.
• Future Work: The authors plan to integrate disk locking and more sophisticated internal angular momentum transport mechanisms to resolve the rotation rate discrepancy and refine the magnetic field saturation physics.

In conclusion, this work represents a significant step forward in stellar modeling by demonstrating that time-dependent magnetic fields and adaptive convection parameters are essential for reproducing the chemical and rotational history of solar-type stars, even if further refinements are needed to fully match the Sun's current state.