Kinematic properties of the TW Hya association

Vadim V. Bobylev

arxiv: 2605.00687 · v1 · submitted 2026-05-01 · 🌌 astro-ph.GA

Kinematic properties of the TW Hya association

Vadim V. Bobylev This is my paper

Pith reviewed 2026-05-09 18:56 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords TW Hya associationkinematic analysisstellar associationsOgorodnikov-Milne modelproper motionsvelocity fielddynamical agegalactic kinematics

0 comments

The pith

The TW Hya association expands with coefficient K_xyz of 103 km s^{-1} kpc^{-1}, implying a dynamical age of 9.7 million years.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper conducts a kinematic analysis of the young stellar association TW Hya by estimating components of the displacement matrix in the Ogorodnikov-Milne linear model, both graphically and through direct solution of the kinematic equations. It confirms the association's volume expansion and derives a dynamical age from that rate while also reporting first estimates of rigid-body rotation around two galactic axes from the graphical approach. A sympathetic reader would care because these properties constrain how recently the group formed and how its stars are moving relative to the galaxy as a whole. The work shows that different analysis methods can give conflicting signals for the rotation components.

Core claim

The association exhibits volume expansion with K_xyz = 103 ± 9 km s^{-1} kpc^{-1}, which corresponds to a dynamical age t = 9.7 ± 0.8 Myr. Graphical analysis yields proper rigid-body rotation parameters ω around galactic axes x and y in the range 50-70 km s^{-1} kpc^{-1} with errors of 14-19 km s^{-1} kpc^{-1}, but solving the kinematic equations using proper motions alone finds all three rotation components consistent with zero: (ω_x, ω_y, ω_z) = (4, 7, 11) ± (5, 5, 5) km s^{-1} kpc^{-1}.

What carries the argument

The Ogorodnikov-Milne linear model of the stellar velocity field, used to extract the expansion coefficient K_xyz and the rigid-body rotation parameters ω from observed proper motions and radial velocities.

If this is right

The TW Hya association has a dynamical age of 9.7 million years derived from its measured expansion rate.
Volume expansion at 103 km s^{-1} kpc^{-1} is established as a dominant kinematic feature of the group.
Rigid-body rotation parameters around galactic x and y axes reach 50-70 km s^{-1} kpc^{-1} when estimated graphically.
The linear velocity-field model applies to the association's spatial scale.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The confirmed expansion implies the association is dispersing and will spread further over the next few million years.
Method-dependent rotation results indicate that higher-precision data could clarify whether any net rotation is present.
These kinematic parameters could help match the association to possible parent molecular clouds or other nearby young groups.
The age estimate provides an independent check on isochronal ages derived from stellar evolution models for the same stars.

Load-bearing premise

The chosen sample of stars consists of genuine association members free from significant contamination or selection bias, and the linear Ogorodnikov-Milne approximation remains valid across the spatial extent of the association.

What would settle it

A new catalog of radial velocities and proper motions for additional candidate members that yields an expansion coefficient statistically inconsistent with 103 km s^{-1} kpc^{-1} or an age estimate outside the 9.7 ± 0.8 Myr range.

Figures

Figures reproduced from arXiv: 2605.00687 by Vadim V. Bobylev.

**Figure 1.** Figure 1: Dependences of velocities U, V, W on coordinates x, y, z. We will also find the time interval that has passed from the beginning of the expansion of the star system to the presentmoment: t = 977.5/Kxyz = 9.1 ± 1.1 Myr. Based on the found values (6) using relations (5), we obtain the following values of the three angular velocities of rotation: ωx = M − 32 = 61 ± 18 km s−1 kpc−1 , ωy = M − 13 = 68 ± 17 km s… view at source ↗

read the original abstract

A kinematic analysis of the young stellar association TWHya has been performed. The components of the displacement matrix in the Ogorodnikov-Milne linear model have been estimated both graphically and by solving the basic kinematic equations. The association's volume expansion with a coefficient of $K_{xyz}=103\pm9$ km s$^{-1}$ kpc$^{-1}$ was confirmed, which yields a dynamical age estimate of $t=9.7 \pm0.8$ Myr. Using the graphical method, estimates of the association's proper rigid-body rotation parameters $\omega$ around the galactic axes x and y have been obtained for the first time, with velocity values in the range of 50-70 km s$^{-1}$ kpc$^{-1}$ and errors in their determination of 14-19 km s$^{-1}$ kpc$^{-1}$. However, these values are not confirmed by another method. For example, when solving kinematic equations only using proper motions, all three components of rigid body rotation do not differ significantly from zero, $(\omega_x,\omega_y,\omega_z)=(4,7,11)\pm(5,5,5)$ km s$^{-1}$ kpc$^{-1}$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

TW Hya paper confirms expansion with both methods but shows clear inconsistency on rotation parameters that undercuts the new claims.

read the letter

Bobylev's paper on the TW Hya association confirms the group's volume expansion using the Ogorodnikov-Milne model. Both the graphical approach and direct solution of the kinematic equations recover K_xyz near 103 km s^{-1} kpc^{-1}, which gives a dynamical age of 9.7 Myr. That part lines up with earlier results on this nearby young group and shows the two techniques agree on the expansion term. The paper also presents what it calls first estimates of rigid-body rotation around the galactic x and y axes from the graphical method, in the 50-70 km s^{-1} kpc^{-1} range. Those numbers appear in the abstract as new. The analysis uses proper motions and the linear displacement matrix, which is standard for this kind of local kinematics work. The sample is drawn from known members of the association, and the methods are laid out plainly. The soft spot is the rotation result itself. The paper states outright that solving the kinematic equations with proper motions alone returns omega components of only 4, 7, and 11 km s^{-1} kpc^{-1} with uncertainties of 5, fully consistent with zero. This mismatch between the two techniques is not minor; it directly questions how stably the linear model captures the velocity field. The dynamical age is simply the inverse of the fitted expansion coefficient with no external anchor, which is a common limitation but reduces its weight. Sample biases or higher-order terms could be at play, though the text does not test them. This is narrow work aimed at specialists tracking nearby associations. It adds some numbers for comparison but does not resolve the rotation issue or broaden to galactic dynamics. The honest presentation of both methods is a plus. It deserves peer review so referees can examine the data handling and model assumptions.

Referee Report

2 major / 1 minor

Summary. The paper performs a kinematic analysis of the TW Hya young stellar association using the Ogorodnikov-Milne linear model. Components of the displacement matrix are estimated both graphically and by solving the basic kinematic equations. The volume expansion coefficient is reported as K_xyz = 103 ± 9 km s^{-1} kpc^{-1} (confirmed by both methods), yielding a dynamical age t = 9.7 ± 0.8 Myr. Graphical estimates of rigid-body rotation parameters give ω_x and ω_y in the range 50-70 km s^{-1} kpc^{-1} (with errors 14-19 km s^{-1} kpc^{-1}), presented as first-time results, but these are not confirmed by the kinematic-equation method using proper motions, which yields (ω_x, ω_y, ω_z) = (4, 7, 11) ± (5, 5, 5) km s^{-1} kpc^{-1}, consistent with zero.

Significance. Confirmation of the expansion term would supply a dynamical age for the association that can be compared against other indicators such as isochrones or lithium depletion. The work adds to the set of nearby young groups with measured velocity gradients. However, the large method-to-method discrepancy in the rotation components indicates that the linear approximation may not be uniformly applicable, limiting the reliability of all derived parameters including K_xyz.

major comments (2)

[Abstract] Abstract: the graphical and kinematic-equation methods return mutually inconsistent values for the rigid-body rotation components (50-70 vs. consistent with zero). Because both methods are applied to the same linear displacement matrix, this discrepancy raises the possibility that the reported K_xyz expansion coefficient is also sensitive to the choice of technique, sample weighting, or unmodeled higher-order velocity terms.
[Abstract] Abstract: the dynamical age t = 9.7 ± 0.8 Myr is obtained directly as the reciprocal of the fitted K_xyz without reference to any independent age anchor or external validation. This makes the age estimate circular with respect to the kinematic fit itself.

minor comments (1)

[Abstract] Abstract: the graphical-method errors are stated as a range (14-19 km s^{-1} kpc^{-1}) rather than component-specific uncertainties; individual values for ω_x and ω_y should be tabulated.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive comments on our manuscript. We provide point-by-point responses to the major comments below, and we plan to incorporate revisions to address the concerns raised.

read point-by-point responses

Referee: [Abstract] Abstract: the graphical and kinematic-equation methods return mutually inconsistent values for the rigid-body rotation components (50-70 vs. consistent with zero). Because both methods are applied to the same linear displacement matrix, this discrepancy raises the possibility that the reported K_xyz expansion coefficient is also sensitive to the choice of technique, sample weighting, or unmodeled higher-order velocity terms.

Authors: The manuscript already notes that the rigid-body rotation parameters from the graphical method are not confirmed by the kinematic-equation method using proper motions. While this discrepancy exists for the rotation components, the expansion coefficient K_xyz is consistent between both methods at 103 ± 9 km s^{-1} kpc^{-1}. This suggests that the expansion term is reliably determined despite the differences in rotation estimates. Nevertheless, we agree that this highlights potential limitations of the linear model. In the revised manuscript, we will expand the discussion to include a more thorough analysis of the method discrepancies and the applicability of the Ogorodnikov-Milne model to this association. revision: partial
Referee: [Abstract] Abstract: the dynamical age t = 9.7 ± 0.8 Myr is obtained directly as the reciprocal of the fitted K_xyz without reference to any independent age anchor or external validation. This makes the age estimate circular with respect to the kinematic fit itself.

Authors: We agree that the dynamical age is derived solely from the kinematic expansion coefficient as t = 1/K_xyz. This is the standard definition for a dynamical age in expanding groups and is presented as such. To avoid any implication of independence, we will revise the abstract and relevant sections to explicitly state that this is a model-dependent dynamical age estimate and provide comparisons with other age indicators from the literature, such as isochrone fitting or lithium depletion boundaries, where available. revision: yes

Circularity Check

1 steps flagged

Dynamical age is direct reciprocal of fitted K_xyz with no independent anchor

specific steps

fitted input called prediction [Abstract]
"The association's volume expansion with a coefficient of K_{xyz}=103±9 km s^{-1} kpc^{-1} was confirmed, which yields a dynamical age estimate of t=9.7 ±0.8 Myr."

K_xyz is obtained by fitting the linear velocity field model to the data. The age t is then stated as the direct reciprocal (t ≈ 1/K_xyz in consistent units), so the numerical value and uncertainty are algebraically determined by the fit itself rather than by any independent chronological anchor or external validation.

full rationale

The paper's central result on volume expansion is obtained by fitting the Ogorodnikov-Milne displacement matrix to the stellar velocities. The reported dynamical age is then obtained solely by algebraic inversion of that fitted coefficient. This is a standard transformation but qualifies as a fitted input presented as a derived estimate. No other load-bearing claims (rotation parameters, model confirmation) reduce to self-definition or self-citation chains; the paper explicitly notes the method discrepancy for rotation and does not force consistency. The derivation chain therefore contains one partial circularity but remains largely self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central results depend on fitting a linear kinematic model to stellar velocity data, with the age derived directly from the expansion parameter.

free parameters (2)

K_xyz = 103 ±9 km s^{-1} kpc^{-1}
Volume expansion coefficient fitted from the kinematic data using the Ogorodnikov-Milne model.
omega_x, omega_y, omega_z = 50-70 (graphical); 4,7,11 (equations)
Rigid-body rotation components estimated from graphical method and kinematic equations.

axioms (1)

domain assumption The Ogorodnikov-Milne linear model is appropriate for describing the velocity field of the TW Hya association.
Invoked to estimate displacement matrix components, expansion, and rotation parameters.

pith-pipeline@v0.9.0 · 5507 in / 1368 out tokens · 92946 ms · 2026-05-09T18:56:41.861286+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

[1]

J. J. Armstrong, J. C. Tan, N. J.Wright, et al., Monthly Not ices Royal Astron. Soc. 543 (3), 2349 (2025)

work page 2025
[2]

Blaauw, Annual Rev

A. Blaauw, Annual Rev. Astron. Astrophys. 2, 213 (1964)

work page 1964
[3]

V. V. Bobylev and A. T. Bajkova, Astronomy Letters 49 (7), 4 10 (2023)

work page 2023
[4]

V. V. Bobylev and A. T. Bajkova, Astrophysical Bulletin 79 (3), 473 (2024)

work page 2024
[5]

de la Reza, E

R. de la Reza, E. Jilinski, and V. G. Ortega, Astron. J. 131 ( 5), 2609 (2006)

work page 2006
[6]

de la Reza, C

R. de la Reza, C. A. O. Torres, G. Quast, et al., Astrophys. J . 343, L61 (1989). 13

work page 1989
[7]

J. K. Donaldson, A. J. Weinberger, J. Gagn´ e, et al., Astro phys. J. 833 (1), article id. 95 (2016)

work page 2016
[8]

Ducourant, R

C. Ducourant, R. Teixeira, P. A. B. Galli, et al., Astron. a nd Astrophys. 563, id. A121 (2014)

work page 2014
[9]

Gagn´ e, O

J. Gagn´ e, O. Roy-Loubier, J. K. Faherty, et al., Astrophy s. J. 860 (1), article id. 43 (2018)

work page 2018
[10]

Gregorio-Hetem, J

J. Gregorio-Hetem, J. R. D. Lepine, G. R. Quast, et al., As tron. J. 103, 549 (1992)

work page 1992
[11]

K. L. Luhman, Astron. J. 165 (6), id. 269 (2023)

work page 2023
[12]

V. V. Makarov and C. Fabricius, Astron. and Astrophys. 36 8, 866 (2001)

work page 2001
[13]

E. E. Mamajek, Astrophys. J. 634 (2), 1385 (2005)

work page 2005
[14]

A. M. Mel’nik and A. K. Dambis, Monthly Notices Royal Astr on. Soc. 472 (4), 3887 (2017)

work page 2017
[15]

A. M. Mel’nik and A. K. Dambis, Astronomy Reports 62 (12), 998 (2018)

work page 2018
[16]

Miret-Roig, J

N. Miret-Roig, J. Alves, S. Ratzenb¨ ock, et al., Astron. and Astrophys. 694, id. A60 (2025)

work page 2025
[17]

Nagananda, L

N. Nagananda, L. Vican, B. Zuckerman, et al., Open Journa l of Astrophysics 7, id. 80 (2024)

work page 2024
[18]

K. F. Ogorodnikov, Dynamics of Stellar Systems (Pergamo n, Oxford, 1965)

work page 1965
[19]

Olivares, N

J. Olivares, N. Miret-Roig, P. A. B. Galli, and H. Bouy, As tron. and Astrophys. 699, id. A122 (2025)

work page 2025
[20]

Posch, J

L. Posch, J. Alves, N. Miret-Roig, et al., Astron. and Ast rophys. 693, id. A175 (2025)

work page 2025
[21]

Preibisch and H

T. Preibisch and H. Zinnecker,Astron. J. 117 (5), 2381 (1 999)

work page
[22]

Sch¨ onrich, J

R. Sch¨ onrich, J. Binney, and W. Dehnen, Monthly Notices Royal Astron. Soc. 403 (4), 1829 (2010)

work page 2010
[23]

Vallenari et al

A. Vallenari et al. (Gaia Collab.), Astron. and Astrophy s. 674, id. A1 (2023)

work page 2023
[24]

N. J. Wright, New Astronomy Reviews 90, article id. 10154 9 (2020)

work page 2020
[25]

Zuckerman and E

B. Zuckerman and E. E. Becklin, Astrophys. J. 406, L25 (19 93). 14

work page

[1] [1]

J. J. Armstrong, J. C. Tan, N. J.Wright, et al., Monthly Not ices Royal Astron. Soc. 543 (3), 2349 (2025)

work page 2025

[2] [2]

Blaauw, Annual Rev

A. Blaauw, Annual Rev. Astron. Astrophys. 2, 213 (1964)

work page 1964

[3] [3]

V. V. Bobylev and A. T. Bajkova, Astronomy Letters 49 (7), 4 10 (2023)

work page 2023

[4] [4]

V. V. Bobylev and A. T. Bajkova, Astrophysical Bulletin 79 (3), 473 (2024)

work page 2024

[5] [5]

de la Reza, E

R. de la Reza, E. Jilinski, and V. G. Ortega, Astron. J. 131 ( 5), 2609 (2006)

work page 2006

[6] [6]

de la Reza, C

R. de la Reza, C. A. O. Torres, G. Quast, et al., Astrophys. J . 343, L61 (1989). 13

work page 1989

[7] [7]

J. K. Donaldson, A. J. Weinberger, J. Gagn´ e, et al., Astro phys. J. 833 (1), article id. 95 (2016)

work page 2016

[8] [8]

Ducourant, R

C. Ducourant, R. Teixeira, P. A. B. Galli, et al., Astron. a nd Astrophys. 563, id. A121 (2014)

work page 2014

[9] [9]

Gagn´ e, O

J. Gagn´ e, O. Roy-Loubier, J. K. Faherty, et al., Astrophy s. J. 860 (1), article id. 43 (2018)

work page 2018

[10] [10]

Gregorio-Hetem, J

J. Gregorio-Hetem, J. R. D. Lepine, G. R. Quast, et al., As tron. J. 103, 549 (1992)

work page 1992

[11] [11]

K. L. Luhman, Astron. J. 165 (6), id. 269 (2023)

work page 2023

[12] [12]

V. V. Makarov and C. Fabricius, Astron. and Astrophys. 36 8, 866 (2001)

work page 2001

[13] [13]

E. E. Mamajek, Astrophys. J. 634 (2), 1385 (2005)

work page 2005

[14] [14]

A. M. Mel’nik and A. K. Dambis, Monthly Notices Royal Astr on. Soc. 472 (4), 3887 (2017)

work page 2017

[15] [15]

A. M. Mel’nik and A. K. Dambis, Astronomy Reports 62 (12), 998 (2018)

work page 2018

[16] [16]

Miret-Roig, J

N. Miret-Roig, J. Alves, S. Ratzenb¨ ock, et al., Astron. and Astrophys. 694, id. A60 (2025)

work page 2025

[17] [17]

Nagananda, L

N. Nagananda, L. Vican, B. Zuckerman, et al., Open Journa l of Astrophysics 7, id. 80 (2024)

work page 2024

[18] [18]

K. F. Ogorodnikov, Dynamics of Stellar Systems (Pergamo n, Oxford, 1965)

work page 1965

[19] [19]

Olivares, N

J. Olivares, N. Miret-Roig, P. A. B. Galli, and H. Bouy, As tron. and Astrophys. 699, id. A122 (2025)

work page 2025

[20] [20]

Posch, J

L. Posch, J. Alves, N. Miret-Roig, et al., Astron. and Ast rophys. 693, id. A175 (2025)

work page 2025

[21] [21]

Preibisch and H

T. Preibisch and H. Zinnecker,Astron. J. 117 (5), 2381 (1 999)

work page

[22] [22]

Sch¨ onrich, J

R. Sch¨ onrich, J. Binney, and W. Dehnen, Monthly Notices Royal Astron. Soc. 403 (4), 1829 (2010)

work page 2010

[23] [23]

Vallenari et al

A. Vallenari et al. (Gaia Collab.), Astron. and Astrophy s. 674, id. A1 (2023)

work page 2023

[24] [24]

N. J. Wright, New Astronomy Reviews 90, article id. 10154 9 (2020)

work page 2020

[25] [25]

Zuckerman and E

B. Zuckerman and E. E. Becklin, Astrophys. J. 406, L25 (19 93). 14

work page