---
title: 'ROTFIT: Stellar Parameter Estimation Tool'
url: https://www.emergentmind.com/topics/rotfit
type: topic
---

# ROTFIT: Stellar Parameter Estimation Tool

ROTFIT is a spectrum-matching code for deriving stellar atmospheric parameters and kinematic diagnostics from observed spectra. Across its documented implementations, it is used to determine effective temperature $T_{\rm eff}$, surface gravity $\log g$, metallicity $[\mathrm{Fe}/\mathrm{H}]$, radial velocity (RV or $V_r$), and projected rotational velocity $v\sin i$ by comparing continuum-normalized target spectra with libraries of real-star or synthetic templates that have been degraded to the instrumental resolution, rotationally broadened, shifted in velocity, and ranked by a goodness-of-fit statistic. The method has been applied to high-resolution echelle spectroscopy, X-Shooter spectra of young stellar objects, low-resolution LAMOST-Kepler spectra, and medium-resolution LAMOST MRS spectra of field stars and open clusters, often in conjunction with spectral subtraction for activity and lithium diagnostics [1606.09149] [1011.1550] [1703.01251] [2205.04757].

## 1. Origins, scope, and scientific role

ROTFIT is described as a spectrum-matching code originally developed by Frasca et al. (2003, 2006) and subsequently adapted to multiple data sets and resolutions [1606.09149]. In the Kepler asteroseismic context, it was used on high-resolution spectroscopic observations to derive $T_{\rm eff}$, $\log g$, $[\mathrm{Fe}/\mathrm{H}]$, and $v\sin i$ for F–K targets, alongside RVs obtained by cross-correlation [1011.1550]. In surveys of young stars and young stellar objects, it was configured with synthetic BT-Settl libraries to infer $T_{\rm eff}$, $\log g$, veiling, RV, and $v\sin i$ from X-Shooter spectra [1703.01251] [1705.10075]. In LAMOST-based work, dedicated implementations were developed for low-resolution LAMOST-Kepler spectra and later for medium-resolution dual-arm LAMOST MRS spectra [1606.09149] [2205.04757].

The scientific role of ROTFIT is broader than parameter estimation alone. In the LAMOST-Kepler survey it supported the construction of a catalog containing atmospheric parameters, RV, and an estimate of $v\sin i$, and for cool stars it also enabled measurement of H$\alpha$ and Ca II IRT chromospheric fluxes [1606.09149]. In the medium-resolution LAMOST analysis of Kepler targets, it was used to determine atmospheric parameters, activity level, lithium atmospheric content, and binarity, including the identification of double-lined binaries, stars with variable RVs, lithium-rich giants, and emission-line objects [2205.04757]. In the Pleiades and NGC 1647 studies, the same framework underpinned cluster RV distributions, metallicity estimates, lithium-based ages, and the detection of SB1 and SB2 systems [2504.09031] [2604.11765].

A recurrent feature of ROTFIT is that it operates empirically when real-star template libraries are used and physically through model atmospheres when synthetic grids are used. This duality is explicit in the literature: some implementations rely on ELODIE, Indo–US, or FRESCO reference spectra, while others adopt BT-Settl synthetic spectra [1011.1550] [1606.09149] [1705.10075]. This suggests that the code is best understood as a flexible fitting framework rather than a single fixed pipeline.

## 2. Core fitting formalism

The central operation of ROTFIT is $\chi^2$ minimization between an observed spectrum and a set of prepared templates. In the LAMOST MRS Kepler implementation, the goodness-of-fit statistic is defined as
$$
\chi^2 \;=\;\sum_{\lambda}
\frac{\bigl[F_{\rm obs}(\lambda)-F_{\rm temp}(\lambda;\,T_{\rm eff},\log g,[\mathrm{Fe}/\mathrm{H}],v\sin i)\bigr]^2}
{\sigma^2(\lambda)}.
$$
The target spectrum is continuum-normalized, shifted to the rest frame using an RV estimate from cross-correlation, and then compared with templates spanning parameter space; the best solutions are combined through weighted means, with weight $w=1/\chi^2$ for the ten best template–rotation combinations in the LAMOST MRS Kepler analysis [2205.04757].

The RV step is commonly handled through cross-correlation. In one formulation,
$$
{\rm CCF}(v) \;=\; A\,\exp\Bigl[-\tfrac12\bigl(v - v_0\bigr)^2/\sigma_{\rm CCF}^2\Bigr]
\;\longrightarrow\; v_0\equiv{\rm RV},
$$
where a Gaussian is fit to the CCF peak and the centroid is taken as the RV [2205.04757]. Closely related CCF formulations appear in the high-resolution Kepler asteroseismic work, the low-resolution LAMOST-Kepler pipeline, the Lupus X-Shooter analysis, and the RasTyc survey [1011.1550] [1606.09149] [1703.01251] [1801.00671].

Rotational broadening is incorporated by convolving each template with a Gray rotational profile. For the LAMOST MRS Kepler work, the profile is written schematically as
$$
G(\Delta\lambda) \propto \left[1 - (\Delta\lambda/\Delta\lambda_D)^2\right]^{1/2}
\quad {\rm for}\ |\Delta\lambda|<\Delta\lambda_D,
$$
with $\Delta\lambda_D = (\lambda_0\,v\sin i)/c$ [2205.04757]. In X-Shooter applications, the Gray kernel is used with linear limb darkening and $v\sin i$ scanned over an explicit trial grid [1703.01251] [2006.07921]. In high-resolution echelle work, trial broadenings are applied order by order, and the final $v\sin i$ is the weighted mean over orders or the value minimizing the summed $\chi^2$ [1011.1550] [1801.00671].

Weighted combination schemes are a further defining element. In the Molenda-Żakowicz et al. Kepler asteroseismic implementation, per-order parameters are averaged using weights proportional to $f_i\cdot\chi_{ij}^{-2}\cdot\sigma_i^{-2}$, where $f_i$ is an information-content factor proportional to total line absorption in the order; the order-wise results are then combined into a global stellar value [1011.1550]. In the RasTyc analysis, order weights are similarly tied to total line absorption and best-fit $\chi^2$ [1801.00671]. In several LAMOST MRS studies, blue- and red-arm solutions are combined as weighted means with $w=1/\sigma^2$ [2504.09031] [2604.11765].

## 3. Template libraries and observational regimes

ROTFIT has been deployed with both empirical and synthetic template libraries, and the choice of library is strongly tied to spectral type, wavelength coverage, and instrumental configuration.

For high-resolution spectroscopy of Kepler asteroseismic targets, two libraries were used in parallel: the ELODIE archive, degraded from $R\approx 42\,000$ to $R=21\,000$, and a FRESCO library of 122 slowly rotating stars calibrated on the same instrumental setup [1011.1550]. In the RasTyc survey, 270 high-S/N, slowly rotating ELODIE spectra with parameters from the PASTEL catalog were employed and convolved to match SARG, FRESCO, or AURELIE resolution [1801.00671]. The low-resolution LAMOST-Kepler pipeline instead used a wide and homogeneous collection of real-star spectra from the Indo–US library, degraded to $R\approx 1\,800$ and re-sampled onto the LAMOST wavelength scale [1606.09149].

For X-Shooter studies of pre-main sequence stars and young stellar objects, ROTFIT was configured with BT-Settl synthetic spectra. In the Lupus work, the grid spans solar metallicity, $T_{\rm eff}=2000$ to $6000\,\mathrm{K}$ in steps of $100\,\mathrm{K}$, and $\log g=0.5$ to $5.5$ dex in steps of $0.5$ dex, with selected UVB and VIS windows analyzed after masking accretion and chromospheric lines [1703.01251]. In the Class III X-Shooter template study, BT-Settl models were again adopted, degraded to match each arm, and analyzed segment by segment with veiling fixed to zero [1705.10075]. For ISO-ChaI 52, the VIS arm alone was used, with BT-Settl templates at solar metallicity spanning $T_{\rm eff}\simeq 2400$–$7000\,\mathrm{K}$ and $\log g=3.5$–$5.5$ dex [2006.07921].

The medium-resolution LAMOST MRS implementations are characterized by dual-arm fitting. In the Kepler-field analysis, the code was modified to accept blue-arm spectra covering 495–535 nm and red-arm spectra covering 630–680 nm at $R\simeq 7\,500$ [2205.04757]. The template library consists of 388 high-S/N, slowly rotating real-star spectra from ELODIE covering $T_{\rm eff}$, $\log g$, and $[\mathrm{Fe}/\mathrm{H}]$ space, convolved to the MRS resolution and resampled onto the MRS wavelength grid [2205.04757]. The Pleiades work used the same 388 ELODIE spectra and a similar dual-arm strategy [2504.09031]. In NGC 1647, cool stars with $T_{\rm eff}\leq 7000\,\mathrm{K}$ were fitted with the same 388 inactive ELODIE templates, whereas hot stars with $T_{\rm eff}\gtrsim 7000\,\mathrm{K}$ were fitted to a grid of synthetic BT-Settl models with $T_{\rm eff}=6000$–$12000\,\mathrm{K}$, $\log g=3.0$–$5.0$ dex, and $[\mathrm{Fe}/\mathrm{H}]$ fixed at 0 [2604.11765].

A consistent pre-processing step across these regimes is resolution matching by Gaussian convolution. One LAMOST MRS formulation gives
$$
W = \sqrt{1/R_{\rm MRS}^2 - 1/R_{\rm ELODIE}^2},
$$
for degrading ELODIE templates to $R_{\rm MRS}\simeq 7\,500$ before resampling [2205.04757]. In the Pleiades analysis the same idea is written in velocity form, $W=c\cdot\sqrt{1/R_{\rm MRS}^2-1/R_{\rm ELODIE}^2}$ [2504.09031]. In NGC 1647, the true LAMOST MRS resolving power was measured from Th-Ar lamps as $R_{\rm blue}=8200\pm 200$ and $R_{\rm red}=8400\pm 300$, and all templates were convolved to these values before fitting [2604.11765].

## 4. Derived quantities beyond atmospheric parameters

Although the core products of ROTFIT are atmospheric parameters, RV, and $v\sin i$, several studies use the same fitted template as a photospheric reference for subtraction analyses.

In the low-resolution LAMOST-Kepler survey, for stars with $T_{\rm eff}<6000\,\mathrm{K}$ the nearest non-active template, shifted to the stellar RV and broadened, is subtracted from the target to derive residual H$\alpha$ emission and Ca II IRT fluxes [1606.09149]. The residual equivalent width is measured over defined integration windows and converted to surface flux through
$$
F_{\rm line}=F_{\rm cont}(T_{\rm eff},\log g)\cdot EW^{\rm res},
$$
with the activity index
$$
R'_{\rm H\alpha}=\frac{F_{\rm H\alpha}}{\sigma_{\rm SB}T_{\rm eff}^4}.
$$
This application led to the identification of 442 chromospherically active stars and one likely accreting object [1606.09149].

In the LAMOST MRS Kepler-field work, H$\alpha$ and Li I 6708 equivalent widths were calculated for cool stars with $T_{\rm eff}<6500\,\mathrm{K}$ [2205.04757]. The analysis yielded 327 active stars on the basis of H$\alpha$ flux and detections of the Li I 6708 line for 1657 stars, both giants and stars on the main sequence [2205.04757]. Among the giants, 195 Li-rich stars were found, 161 of them reported there for the first time; among main-sequence stars, a discrete age classification was performed based on atmospheric lithium abundance and upper envelopes of a few open clusters [2205.04757].

The same template-subtraction logic is used in cluster studies. In the Pleiades analysis, late-type stars with $T_{\rm eff}<6500\,\mathrm{K}$ were characterized through H$\alpha$ and Li I-6708 net equivalent widths by subtraction of inactive photospheric templates, enabling studies of activity, lithium depletion, and flares [2504.09031]. In NGC 1647, H$\alpha$ and Li I-6708 net equivalent widths were measured for solar-type and cooler stars, supporting inference of a lithium-based cluster age of $203\pm 27\,\mathrm{Myr}$ [2604.11765].

In young-star spectroscopy, the photospheric template also anchors accretion and chromospheric diagnostics. In the Lupus X-Shooter study, the best-fit BT-Settl spectrum, rotationally broadened and veiled, was subtracted from the observed spectrum to recover net emission in H$\alpha$, H$\beta$, Ca II K and IRT, and Na I D [1703.01251]. Residual line equivalent widths were converted to surface fluxes using model continuum fluxes,
$$
F_{\rm line}=EW_{\rm line}\times F_C^{\rm model}(T_{\rm eff},\log g),
$$
and these diagnostics were combined with $\log g$, RV, and Li I 6708 equivalent width for membership assessment [1703.01251].

## 5. Calibration, uncertainty estimation, and validation

Uncertainty estimation in ROTFIT depends on implementation, but several patterns recur: internal errors from the scatter among orders or segments, formal errors from the curvature of the $\chi^2$ minimum, and external checks against independent surveys or blue–red arm comparisons.

For the Kepler asteroseismic analysis at $R\approx 21\,000$, validation against ELODIE-based and FRESCO-based solutions showed agreement to within $\sigma(T_{\rm eff})\approx 100\,\mathrm{K}$, $\sigma(\log g)\approx 0.1$ dex, and $\sigma([\mathrm{Fe}/\mathrm{H}])\approx 0.1$ dex, while independent synthetic-spectrum fitting for a subset agreed with ROTFIT to better than the adopted grid spacing of $\pm 125\,\mathrm{K}$, $\pm 0.25$ dex, and $\pm 2\,\mathrm{km\,s^{-1}}$ [1011.1550]. In the X-Shooter Class III analysis, typical internal uncertainties are quoted as $\sigma(T_{\rm eff})\sim 70$–$200\,\mathrm{K}$, $\sigma(\log g)\sim 0.1$–$0.3$ dex, $\sigma(v\sin i)\sim 1$–$9\,\mathrm{km\,s^{-1}}$, and $\sigma(\mathrm{RV})\sim 2$–$3\,\mathrm{km\,s^{-1}}$ [1705.10075]. In the Lupus study, typical internal errors are $\sigma_{T_{\rm eff}}\sim 70\,\mathrm{K}$, $\sigma_{\log g}\sim 0.21$ dex, $\sigma_{\rm RV}\sim 2.3\,\mathrm{km\,s^{-1}}$, and $\sigma_{v\sin i}\sim 4\,\mathrm{km\,s^{-1}}$ [1703.01251].

The low-resolution LAMOST-Kepler pipeline reported internal repeat-based precisions of $\sigma(T_{\rm eff})\approx 66\,\mathrm{K}$, $\sigma(\log g)\approx 0.046$ dex, and $\sigma([\mathrm{Fe}/\mathrm{H}])\approx 0.055$ dex, but external accuracies against high-resolution literature values were about $14\,\mathrm{km\,s^{-1}}$ in RV, $3.5\%$ in $T_{\rm eff}$, $0.3$ dex in $\log g$, and $0.2$ dex in $[\mathrm{Fe}/\mathrm{H}]$ [1606.09149]. This distinction between internal consistency and external accuracy is central to interpreting ROTFIT outputs in low-resolution surveys.

The medium-resolution LAMOST MRS analyses introduced explicit blue–red arm cross-validation. In the Kepler-field study, blue–red differences of $65\pm 175\,\mathrm{K}$ in $T_{\rm eff}$, $-0.08\pm 0.26$ in $\log g$, $0.00\pm 0.12$ in $[\mathrm{Fe}/\mathrm{H}]$, $-3.8\pm 3.4\,\mathrm{km\,s^{-1}}$ in $v\sin i$, and $-0.1\pm 2.3\,\mathrm{km\,s^{-1}}$ in RV imply per-arm uncertainties of about $120\,\mathrm{K}$, $0.18$ dex, $0.09$ dex, $2.4\,\mathrm{km\,s^{-1}}$, and $1.6\,\mathrm{km\,s^{-1}}$, respectively [2205.04757]. Comparison with APOGEE DR16 for about 2500 stars gave rms dispersions of about $160\,\mathrm{K}$ in $T_{\rm eff}$, $0.30$ dex in $\log g$, and $\sigma\sim 6.1\,\mathrm{km\,s^{-1}}$ in RV including genuine RV variables, or $\sigma\sim 3.4\,\mathrm{km\,s^{-1}}$ if variables were removed, corresponding to an accuracy of about $1.3\,\mathrm{km\,s^{-1}}$ on a clean sample [2205.04757].

The Pleiades study used the rms of blue–red differences divided by $\sqrt{2}$ as an external empirical error estimate and found $\sigma_{\rm RV}\approx 2.2\,\mathrm{km\,s^{-1}}$, $\sigma_{T_{\rm eff}}\approx 143\,\mathrm{K}$ for $T_{\rm eff}>4500\,\mathrm{K}$, $\sigma_{\log g}\approx 0.15$ dex, $\sigma_{[\mathrm{Fe}/\mathrm{H}]}\approx 0.09$ dex, and $\sigma_v\approx 5$–$6\,\mathrm{km\,s^{-1}}$ for $v\sin i\geq 8\,\mathrm{km\,s^{-1}}$ [2504.09031]. In NGC 1647, blue–red arm comparisons yielded rms $\Delta V_r\approx 6.9\,\mathrm{km\,s^{-1}}$ overall, or about $3.6\,\mathrm{km\,s^{-1}}$ for $v\sin i<70\,\mathrm{km\,s^{-1}}$, corresponding to a single-arm precision of about $2.5\,\mathrm{km\,s^{-1}}$; for cool stars, the rms scatters were about $250\,\mathrm{K}$ in $T_{\rm eff}$, $0.15$ dex in $\log g$, $0.18$ dex in $[\mathrm{Fe}/\mathrm{H}]$, and $23\,\mathrm{km\,s^{-1}}$ in $v\sin i$, reduced to about $15\,\mathrm{km\,s^{-1}}$ for $v\sin i<100\,\mathrm{km\,s^{-1}}$ [2604.11765].

RV zero-point calibration is especially important in LAMOST MRS work. In the Kepler-field analysis, plate-by-plate zero-point offsets of up to about $6\,\mathrm{km\,s^{-1}}$ before May 2018 were corrected by comparison to Gaia RVs of non-variable stars on the same plate [2205.04757]. The Pleiades work states that early MRS runs before May 2018 were calibrated with Sc lamps and later ones with Th–Ar, causing $\lesssim 6\,\mathrm{km\,s^{-1}}$ RV offsets; corrections were applied plate by plate and arm by arm using contemporaneous RV standards or Gaia DR3 stars with $|\Delta \mathrm{RV}|<0.5\,\mathrm{km\,s^{-1}}$ [2504.09031]. NGC 1647 likewise required correction for nightly wavelength-calibration zero-point shifts using Th-Ar lamp runs or DR12 recipes [2604.11765].

## 6. Resolution limits, biases, and methodological constraints

A major practical limitation of ROTFIT is the resolving-power floor on $v\sin i$. In the medium-resolution LAMOST MRS Kepler analysis, Monte Carlo tests with BT-Settl spectra at S/N = 50–100 showed that $v\sin i<8\,\mathrm{km\,s^{-1}}$ cannot be reliably resolved, and values below that threshold are reported as “$<8\,\mathrm{km\,s^{-1}}$” [2205.04757]. The Pleiades and NGC 1647 studies adopt the same lower limit for LAMOST MRS [2504.09031] [2604.11765]. In X-Shooter work, the minimum detectable $v\sin i$ is about $8\,\mathrm{km\,s^{-1}}$ for the VIS slit at $R\simeq 8400$ and about $6\,\mathrm{km\,s^{-1}}$ at $R\simeq 17\,000$, as verified by Monte Carlo tests [1703.01251]. At low LAMOST resolution, projected rotation is only robustly measurable for $v\sin i\gtrsim 120\,\mathrm{km\,s^{-1}}$ [1606.09149].

Template-grid structure also imposes systematic effects. The LAMOST MRS Kepler study notes that template-grid non-uniformity leads to clustering of returned parameters at nodes, especially at low $[\mathrm{Fe}/\mathrm{H}]$ and $\log g$, and that no interpolation or regularization between templates produces discrete jumps in output [2205.04757]. The Pleiades paper describes analogous “boxy” behavior in parameter space because ROTFIT picks the single best non-interpolated template, producing clusters of solutions at discrete grid points [2504.09031]. In the NGC 1647 analysis, the hot-star template grid is restricted to solar metallicity and coarse $\log g$ steps, and early-type results carry larger uncertainties [2604.11765].

Metallicity calibration is another recurrent issue. In the low-resolution LAMOST-Kepler survey, comparison with APOKASC and SAGA revealed systematic biases, leading to the relation
$$
[\mathrm{Fe}/\mathrm{H}]_{\rm corr}=2.16\,[\mathrm{Fe}/\mathrm{H}]_{\rm LAMOST}+0.17,
\qquad [\mathrm{Fe}/\mathrm{H}]>-1.0,
$$
while the raw ROTFIT values were retained in the catalog with advice that users apply the correction as needed [1606.09149]. In the medium-resolution LAMOST MRS Kepler analysis, comparison with APOGEE showed overestimation at $[\mathrm{Fe}/\mathrm{H}]<0$ and underestimation at $[\mathrm{Fe}/\mathrm{H}]>0.2$, leading to
$$
[\mathrm{Fe}/\mathrm{H}]_{\rm corr}
=
2.57\cdot[\mathrm{Fe}/\mathrm{H}]_{\rm ROTFIT}
-0.01,
\qquad {\rm valid\ for}\ [\mathrm{Fe}/\mathrm{H}]_{\rm ROTFIT}>-1.0.
$$
The same study states that better template coverage at $[\mathrm{Fe}/\mathrm{H}]<-0.5$ and $>+0.3$ dex would help [2205.04757].

Other limitations are regime-specific. In the X-Shooter Class III analysis, BT-Settl synthetic spectra are said to show small mismatches in certain molecular bands, notably TiO and VO, at $\log g\approx 3.5$–4.0; this introduces systematic uncertainties of about $100\,\mathrm{K}$ in $T_{\rm eff}$ and about $0.2$ dex in $\log g$ for the coolest objects beyond the quoted statistical errors [1705.10075]. In the Lupus implementation, veiling below $r\lesssim 0.2$ is not significant, and under-luminous objects with edge-on disks require special treatment when placed on the HR diagram [1703.01251]. In NGC 1647, double-lined spectroscopic binaries are not handled by ROTFIT in its single-template mode and require dedicated SB2 analysis after identification through double peaks in the CCF [2604.11765].

## 7. Representative applications and evolution of the framework

ROTFIT has been used in diverse observational programs, and its scientific outputs illustrate the evolution of the code from stellar parameter estimation to broader survey analysis.

In Kepler-related spectroscopy, the high-resolution asteroseismic study derived RV, $v\sin i$, $T_{\rm eff}$, $\log g$, and metallicity for 44 targets, discovered three double-lined spectroscopic binaries and two suspected single-lined systems, and emphasized the need for verification and refinement of atmospheric parameters in the Kepler Input Catalog [1011.1550]. The later low-resolution LAMOST-Kepler application produced atmospheric parameters and RVs for 61,753 spectra of 51,385 stars, identified RV variables, ultrafast rotators, and emission-line objects, and found 442 chromospherically active stars [1606.09149]. The medium-resolution LAMOST Kepler-field analysis then extended the same line of work to 16,300 spectra, deriving RV and atmospheric parameters for 14,300 spectra of 7443 stars and reporting 327 active stars, 1657 stars with Li I 6708 detections, and a fraction of Li-rich giants of about 4%, higher than expected [2205.04757].

In young-star studies, ROTFIT was integrated with X-Shooter analyses of pre-main sequence and accreting systems. The Lupus survey used it to determine $T_{\rm eff}$, $\log g$, veiling, RV, and $v\sin i$ for 102 YSO candidates, rejecting 13 as non-members primarily through discrepant RV or very low $\log g$ [1703.01251]. The Class III template work used ROTFIT-derived stellar parameters as part of a library of reduced, flux-calibrated, and telluric-corrected spectra for disk-less pre-main sequence stars [1705.10075]. In the ISO-ChaI 52 study, ROTFIT on the X-Shooter VIS spectrum delivered $T_{\rm eff}=3240\pm 70\,\mathrm{K}$, $\log g=4.1\pm 0.3$ dex, $\mathrm{RV}=+4.3\pm 1.0\,\mathrm{km\,s^{-1}}$, $v\sin i=9.9\pm 0.6\,\mathrm{km\,s^{-1}}$, and veiling $r<0.2$ in the red VIS, anchoring the subsequent SED and disk-warp analysis [2006.07921].

In open-cluster spectroscopy with LAMOST MRS, the framework became closely tied to rotation, activity, lithium depletion, and binarity. The Pleiades analysis applied ROTFIT to 1581 spectra of 283 stars, finding an RV distribution that peaks at $5.0\,\mathrm{km\,s^{-1}}$ with a dispersion of $1.4\,\mathrm{km\,s^{-1}}$, an average metallicity of $[\mathrm{Fe}/\mathrm{H}]=-0.03\pm 0.06$, 39 possible SB1 and ten SB2 systems, and a lithium-isochrone age of $118\pm 6\,\mathrm{Myr}$ [2504.09031]. The NGC 1647 work used ROTFIT for 158 cluster members, found an RV distribution peaking at $-5.3\,\mathrm{km\,s^{-1}}$ with a dispersion of $1.6\,\mathrm{km\,s^{-1}}$, an average metallicity of $[\mathrm{Fe}/\mathrm{H}]=-0.08\pm 0.08$ dex, four double-lined spectroscopic systems, and a lithium-based age of $203\pm 27\,\mathrm{Myr}$ [2604.11765].

Taken together, these implementations show a stable methodological core—cross-correlation for RV, rotational broadening, $\chi^2$-based template matching, and weighted parameter aggregation—combined with instrument-specific calibration, library selection, and post-fit spectral subtraction. A plausible implication is that ROTFIT’s enduring utility derives less from a single fixed parameterization than from its adaptability to different resolutions, spectral types, and survey objectives, provided that template coverage, instrumental resolution, and zero-point calibration are explicitly controlled.

Source: https://www.emergentmind.com/topics/rotfit