Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error correction.”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81

Patient motion effects on the quantification of regional myocardial blood flow with dynamic PET imaging

Purpose: Patient motion is a common problem during dynamic positron emission tomography (PET) scans for quantification of myocardial blood flow (MBF). The purpose of this study was to quantify the prevalence of body motion in a clinical setting and evaluate with realistic phantoms the effects of motion on blood flow quantification, including CT attenuation correction (CTAC) artifacts that result from PET–CT misalignment. Methods: A cohort of 236 sequential patients was analyzed for patient motion under resting and peak stress conditions by two independent observers. The presence of motion, affected time-frames, and direction of motion was recorded; discrepancy between observers was resolved by consensus review. Based on these results, patient body motion effects on MBF quantification were characterized using the digital NURBS-based cardiac-torso phantom, with characteristic time activity curves (TACs) assigned to the heart wall (myocardium) and blood regions. Simulated projection data were corrected for attenuation and reconstructed using filtered back-projection. All simulations were performed without noise added, and a single CT image was used for attenuation correction and aligned to the early- or late-frame PET images. Results: In the patient cohort, mild motion of 0.5 ± 0.1 cm occurred in 24% and moderate motion of 1.0 ± 0.3 cm occurred in 38% of patients. Motion in the superior/inferior direction accounted for 45% of all detected motion, with 30% in the superior direction. Anterior/posterior motion was predominant (29%) in the posterior direction. Left/right motion occurred in 24% of cases, with similar proportions in the left and right directions. Computer simulation studies indicated that errors in MBF can approach 500% for scans with severe patient motion (up to 2 cm). The largest errors occurred when the heart wall was shifted left toward the adjacent lung region, resulting in a severe undercorrection for attenuation of the heart wall. Simulations also indicated that the magnitude of MBF errors resulting from motion in the superior/inferior and anterior/posterior directions was similar (up to 250%). Body motion effects were more detrimental for higher resolution PET imaging (2 vs 10 mm full-width at half-maximum), and for motion occurring during the mid-to-late time-frames. Motion correction of the reconstructed dynamic image series resulted in significant reduction in MBF errors, but did not account for the residual PET–CTAC misalignment artifacts. MBF bias was reduced further using global partial-volume correction, and using dynamic alignment of the PET projection data to the CT scan for accurate attenuation correction during image reconstruction. Conclusions: Patient body motion can produce MBF estimation errors up to 500%. To reduce these errors, new motion correction algorithms must be effective in identifying motion in the left/right direction, and in the mid-to-late time-frames, since these conditions produce the largest errors in MBF, particularly for high resolution PET imaging. Ideally, motion correction should be done before or during image reconstruction to eliminate PET-CTAC misalignment artifacts.

60 APPLIED LIFE SCIENCES↗

First measurements of error fields on W7-X using flux surface mapping

Error fields have been detected and quantified using the flux surface mapping diagnostic system on Wendelstein 7-X (W7-X). A low-field '${\rlap{-}\ \iota} =1/2$ ' magnetic configuration (${\rlap{-}\ \iota} =\iota /2\pi $ ), sensitive to error fields, was developed in order to detect their presence using the flux surface mapping diagnostic. In this configuration, a vacuum flux surface with rotational transform of n/m = 1/2 is created at the mid-radius of the vacuum flux surfaces. If no error fields are present a vanishingly small n/m = 5/10 island chain should be present. Modeling indicates that if an n = 1 perturbing field is applied by the trim coils, a large n/m = 1/2 island chain will be opened. This island chain is used to create a perturbation large enough to be imaged by the diagnostic. Phase and amplitude scans of the applied field allow the measurement of a small $\sim 0.04$ m intrinsic island chain with a ${{130}^{\circ}}$ phase relative to the first module of the W7-X experiment. Lastly, these error fields are determined to be small and easily correctable by the trim coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Smoothing tool design and performance during subaperture glass polishing

During subaperture tool grinding and polishing, overlaps of the tool influence function can result in undesirable mid-spatial frequency (MSF) errors in the form of surface ripples, which are often corrected using a smoothing polishing step. Here, in this study, flat multi-layer smoothing polishing tools are designed and tested to simultaneously (1) reduce or remove MSF errors, (2) minimize surface figure degradation, and (3) maximize the material removal rate. A time-dependent convergence model in which spatial material removal varies with a workpiece-tool height mismatch, combined with a finite element mechanical analysis to determine the interface contact pressure distribution, was developed to evaluate various smoothing tool designs as a function of tool material properties, thicknesses, pad textures, and displacements. An improvement in smoothing tool performance is achieved when the gap pressure constant, $\overline{h}$ (which describes the inverse rate at which the pressure drops with a workpiece-tool height mismatch), is minimized for smaller spatial scale length surface features (namely, MSF errors) and maximized for large spatial scale length features (i.e., surface figure). Five specific smoothing tool designs were experimentally evaluated. A two-layer smoothing tool using a thin, grooved IC1000 polyurethane pad (with a high elastic modulus, E pad =360 MPa ), thicker blue foam (with an intermediate modulus, E foam =5.3 MPa ) underlayer, and an optimized displacement (d t =1 mm ) provided the best overall performance (namely, high MSF error convergence, minimal surface figure degradation, and high material removal rate).

47 OTHER INSTRUMENTATION↗

LAF-Net: A Deep Residual and Cross-Attention Framework for Day-Ahead Load Forecasting: Preprint

Accurate day-ahead load forecasting is essential for reliable power system operations and market efficiency. System operators such as the Midcontinent Independent System Operator (MISO) rely on forecasts from multiple vendors, yet combining them effectively remains a persistent challenge due to vendor-specific biases. This paper presents a novel LSTM-Attention Fusion Network with Error Representation (LAF-Net) that enhances day-ahead hourly load forecasting through deep residual learning and multi-modal cross-attention. The proposed model builds a historical error memory from past vendor performance and dynamically queries it with future hour context to generate adaptive, hour-specific trust weights for each vendor. A bounded residual correction further refines forecasts by mitigating systematic and temporally localized errors. Tested on real MISO LBA data with multi-vendor forecasts, LAF-Net consistently outperforms the best vendor baseline across all 38 LBAs, achieving more than a 40% reduction in system-level mean absolute error (MAE) during peak load hours relative to the best vendor baseline.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SU-E-T-299: Small Fields Profiles Correction Through Detectors Spatial Response Functions and Field Size Dependence Analysis

Purpose: To calculate the spatial response function of various radiation detectors, to evaluate the dependence on the field size and to analyze the small fields profiles corrections by deconvolution techniques. Methods: Crossline profiles were measured on a Novalis Tx 6MV beam with a HDMLC. The configuration setup was SSD=100cm and depth=5cm. Five fields were studied (200×200mm2,100×100mm2, 20×20mm2, 10×10mm2and 5×5mm2) and measured were made with passive detectors (EBT3 radiochromic films and TLD700 thermoluminescent detectors), ionization chambers (PTW30013, PTW31003, CC04 and PTW31016) and diodes (PTW60012 and IBA SFD). The results of passive detectors were adopted as the actual beam profile. To calculate the detectors kernels, modeled by Gaussian functions, an iterative process based on a least squares criterion was used. The deconvolutions of the measured profiles were calculated with the Richardson-Lucy method. Results: The profiles of the passive detectors corresponded with a difference in the penumbra less than 0.1mm. Both diodes resolve the profiles with an overestimation of the penumbra smaller than 0.2mm. For the other detectors, response functions were calculated and resulted in Gaussian functions with a standard deviation approximate to the radius of the detector in study (with a variation less than 3%). The corrected profiles resolve the penumbra with less than 1% error. Major discrepancies were observed for cases in extreme conditions (PTW31003 and 5×5mm2 field size). Conclusion: This work concludes that the response function of a radiation detector is independent on the field size, even for small radiation beams. The profiles correction, using deconvolution techniques and response functions of standard deviation equal to the radius of the detector, gives penumbra values with less than 1% difference to the real profile. The implementation of this technique allows estimating the real profile, freeing from the effects of the detector used for the acquisition.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Towards Automated Reasoning Chains for Verification of LLM-Generated Scientific Code

With the rise of Large Language Model (LLM) generated code, including in domains like scientific computing, ensuring not only syntactical, but also mathematical correctness, has become a critical task. Traditional formal methods approaches often struggle with the ambiguity of floating-point code, and full symbolic execution is extremely costly and limited. We propose a chain-of-reasoning approach that iteratively lifts basic semantics from code into the SPIRAL system and then establishes numerical equivalency to the desired mathematical operation. Here, we leverage the ample mathematical knowledge already formalized in SPIRAL to enable the system to recognize not just different implementations of the same algorithm but fully separate approaches to solving the given problem. The chain establishes tight error bounds on the output of given code with respect to the true continuous solution it approximates, quantifying all sources of error. We demonstrate this approach by establishing the correctness of a pseudospectral solver for a simple 1-dimensional Poisson problem.

Oschatz, Quentin [Carnegie Mellon University,Pitts↗

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Image-Guided Radiotherapy for Breast Cancer Patients: Surgical Clips as Surrogate for Breast Excision Cavity

Purpose: To determine the use of surgical clips as a surrogate for localization of the excision cavity and to quantify the stability of the clips' positions during the course of external beam radiotherapy for breast cancer patients, using cone beam computed tomography (CBCT) scans. Methods and Materials: Twenty-one breast cancer patients with surgical clips placed in the breast excision cavity were treated in a supine position with 28 daily fractions. CBCT scans were regularly acquired for a setup correction protocol. Retrospectively, the CBCT scans were registered to the planning CT scans, using gray-value registration of the excision cavity region and chamfer matching of the clips. Subsequently, residual setup errors (systematic [{Sigma}] and random [{sigma}]) of the excision cavity were estimated relative to the clips' registration. Finally, the stability of the clips' positions were quantified as the movement of each separate clip according to the center of gravity of the excision cavity. Results: When clips were used for online setup corrections, the residual errors of the excision cavity were {Sigma}{sub left-right} = 1.2, {sigma}{sub left-right} = 1.0; {Sigma}{sub cranial-caudal} = 1.3, {sigma}{sub cranial-caudal} = 1.2; and {Sigma}{sub anterior-posterior} = 0.7, {sigma}{sub anterior-posterior} = 0.9 mm. Furthermore, the average distance (over all patients) between the clips and centers of gravity of the excision cavities was 18.8 mm (on the planning CT) and was reduced to 17.4 mm (measured on the last CBCT scan). Conclusion: Clips move in the direction of the center of gravity of the excision cavity, on average, 1.4 mm. The clips are good surrogates for locating the excision cavity and providing small residual errors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Software-controlled, highly automated intrafraction prostate motion correction with intrafraction stereographic targeting: System description and clinical results

Purpose: A new system for software-controlled, highly automated correction of intrafraction prostate motion,'' intrafraction stereographic targeting'' (iSGT), is described and evaluated. Methods: At our institute, daily prostate positioning is routinely performed at the start of treatment beam using stereographic targeting (SGT). iSGT was implemented by extension of the SGT software to facilitate fast and accurate intrafraction motion corrections with minimal user interaction. iSGT entails megavoltage (MV) image acquisitions with the first segment of selected IMRT beams, automatic registration of implanted markers, followed by remote couch repositioning to correct for intrafraction motion above a predefined threshold, prior to delivery of the remaining segments. For a group of 120 patients, iSGT with corrections for two nearly lateral beams was evaluated in terms of workload and impact on effective intrafraction displacements in the sagittal plane. Results: SDs of systematic ({Sigma}) and random ({sigma}) displacements relative to the planning CT measured directly after initial SGT setup correction were <0.5 and <0.8 mm, respectively. Without iSGT corrections, effective {Sigma} and {sigma} for the 11-min treatments would increase to {Sigma}{sub eff} < 1.1 mm and {sigma}{sub eff} < 1.2 mm. With the iSGT procedure with an action level of 4 mm, effective positioning errors were reduced to {Sigma}{sub eff} < 0.8 mm and {sigma}{sub eff} < 1.0 mm, with 23.1% of all fractions requiring a correction. Computer simulations demonstrated that with an action level of 2 mm, the errors would have been reduced to {Sigma}{sub eff} < 0.6 mm and {sigma}{sub eff} < 0.7 mm, requiring corrections in 82.4% of the fractions. Because iSGT is highly automated, the extra time added by iSGT is <30 s if a correction is required. Conclusions: Without increasing imaging dose, iSGT successfully reduces intrafraction prostate motion with minimal workload and increase in fraction time. An action level of 2 mm is recommended.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Investigating Tropical Versus Extratropical Influences on the Southern Hemisphere Tropical Edge in the Unified Model

Abstract Since the late 1970s, observations have shown a widening of the tropical Hadley cell (HC) circulation. State‐of‐the‐art climate models reproduce the general trend along with a projected continuous expansion. Discrepancies in expansion rates of observation‐ and model‐based studies have been attributed to differences in applied methods, natural variability and model shortcomings. Furthermore, the driving influence of tropical or extratropical processes on these changes is not well understood. All of this highlights the dynamical mechanisms and the region of origin controlling the tropical width are still insufficiently understood. Here we examine the influence of systematic model biases of the atmosphere‐only Unified Model (UM) onto the simulation of the Southern Hemisphere (SH) tropical edge. We utilize nudged experiments with prescribed sea surface temperatures, where potential temperature and horizontal winds are relaxed back to ERA‐Interim reanalysis for a 20‐year period in selected regions. Correcting model biases in the tropics and extratropics separately allows us to dissect the dominant remote impacts of present model errors onto the SH tropical edge simulation. The experiments are applied to established tropical width metrics ranging from near‐surface to upper‐level metrics capturing the poleward flank of the HC. We find both regions work remotely to reduce errors in the UM fields and location of the tropical edge. Surprisingly, correcting the extratropical biases, south of 45°S, more consistently improves the tropical width across the metrics and seasons than nudging the tropics (10°N–10°S). These findings demonstrate the substantial role of extratropical influences in locating the SH tropical edge.

54 ENVIRONMENTAL SCIENCES↗

Validation of a pretreatment delivery quality assurance method for the CyberKnife Synchrony system

Purpose: To evaluate the geometric and dosimetric accuracies of the CyberKnife Synchrony respiratory tracking system (RTS) and to validate a method for pretreatment patient-specific delivery quality assurance (DQA). Methods: An EasyCube phantom was mounted on the ExacTrac gating phantom, which can move along the superior–inferior (SI) axis of a patient to simulate a moving target. The authors compared dynamic and static measurements. For each case, a Gafchromic EBT3 film was positioned between two slabs of the EasyCube, while a PinPoint ionization chamber was placed in the appropriate space. There were three steps to their evaluation: (1) the field size, the penumbra, and the symmetry of six secondary collimators were measured along the two main orthogonal axes. Dynamic measurements with deliberately simulated errors were also taken. (2) The delivered dose distributions (from step 1) were compared with the planned ones, using the gamma analysis method. The local gamma passing rates were evaluated using three acceptance criteria: 3% local dose difference (LDD)/3 mm, 2%LDD/2 mm, and 3%LDD/1 mm. (3) The DQA plans for six clinical patients were irradiated in different dynamic conditions, to give a total of 19 cases. The measured and planned dose distributions were evaluated with the same gamma-index criteria used in step 2 and the measured chamber doses were compared with the planned mean doses in the sensitive volume of the chamber. Results: (1) A very slight enlargement of the field size and of the penumbra was observed in the SI direction (on average <1 mm), in line with the overall average CyberKnife system error for tracking treatments. (2) Comparison between the planned and the correctly delivered dose distributions confirmed the dosimetric accuracy of the RTS for simple plans. The multicriteria gamma analysis was able to detect the simulated errors, proving the robustness of their method of analysis. (3) All of the DQA clinical plans passed the tests, both in static and dynamic conditions. No statistically significant differences were found between static and dynamic cases, confirming the high degree of accuracy of the Synchrony RTS. Conclusions: The presented methods and measurements verified the mechanical and dosimetric accuracy of the Synchrony RTS. Their method confirms the fact that the RTS, if used properly, is able to treat a moving target with great precision. By combining PinPoint ion chamber, EBT3 films, and gamma evaluation of dose distributions, their DQA method robustly validated the effectiveness of CyberKnife and Synchrony system.

60 APPLIED LIFE SCIENCES↗

MEASURING HIGH-PRECISION ASTROMETRY WITH THE INFRARED ARRAY CAMERA ON THE SPITZER SPACE TELESCOPE

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope currently offers the greatest potential for high-precision astrometry of faint mid-IR sources across arcminute-scale fields, which would be especially valuable for measuring parallaxes of cold brown dwarfs in the solar neighborhood and proper motions of obscured members of nearby star-forming regions. To more fully realize IRAC's astrometric capabilities, we have sought to minimize the largest sources of uncertainty in astrometry with its 3.6 and 4.5 μm bands. By comparing different routines that estimate stellar positions, we have found that Point Response Function (PRF) fitting with the Spitzer Science Center's Astronomical Point Source Extractor produces both the smallest systematic errors from varying intra-pixel sensitivity and the greatest precision in measurements of positions. In addition, self-calibration has been used to derive new 7th and 8th order distortion corrections for the 3.6 and 4.5 μm arrays of IRAC, respectively. These corrections are suitable for data throughout the mission of Spitzer when a time-dependent scale factor is applied to the corrections. To illustrate the astrometric accuracy that can be achieved by combining PRF fitting with our new distortion corrections, we have applied them to archival data for a nearby star-forming region, arriving at total astrometric errors of ∼20 and 70 mas at signal to noise ratios of 100 and 10, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

An AFM Dynamic Contact Model with Finite Thickness Correction to Study Micro-Rheology of Biological Cells

One of the major sources of error in Atomic Force Microscopy based Force Spectroscopic studies is bottom substrate effect, which comes due to assumption of infinite sample thickness. In the present work, anew experimentally validated mathematical correction for finite thickness of thesample is incorporated into a dynamic indentation model. This model is an extension of Sneddon’s static indentation contact model developed for conical tip indenters by incorporating the tip geometric effects, correction to hydrodynamic drag, viscoelastic nature of the sample and its finite thickness. The new model is applied to extract the viscoelastic properties of primary mesenchymal stem cells (MSC’s) extracted from human bone marrow and HCT-116 aggressive epithelial colorectal cancer cells. The micro-rheological studies indicated the presence of profound influence of infinite thickness assumption depending upon indentation depth and sample thickness at the point of indentation. The present model nullified this error and provided a better estimate of the rheological properties. The tests were conducted over a wide frequency bandwidth (0.5 to 126 Hz) and frequency dependent material properties were analyzed using power law structural damping model. It was found that MSC’s are stiffer and less viscous in nature compared to cancer cells.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Band-gap corrected density functional theory calculations for InAs/GaSb type II superlattices

We performed pseudopotential based density functional theory (DFT) calculations for GaSb/InAs type II superlattices (T2SLs), with bandgap errors from the local density approximation mitigated by applying an empirical method to correct the bulk bandgaps. Specifically, this work (1) compared the calculated bandgaps with experimental data and non-self-consistent atomistic methods; (2) calculated the T2SL band structures with varying structural parameters; (3) investigated the interfacial effects associated with the no-common-atom heterostructure; and (4) studied the strain effect due to lattice mismatch between the two components. This work demonstrates the feasibility of applying the DFT method to more exotic heterostructures and defect problems related to this material system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A REVISED EFFECTIVE TEMPERATURE SCALE FOR THE KEPLER INPUT CATALOG

We present a catalog of revised effective temperatures for stars observed in long-cadence mode in the Kepler Input Catalog (KIC). We use Sloan Digital Sky Survey (SDSS) griz filters tied to the fundamental temperature scale. Polynomials for griz color-temperature relations are presented, along with correction terms for surface gravity effects, metallicity, and statistical corrections for binary companions or blending. We compare our temperature scale to the published infrared flux method (IRFM) scale for V{sub T}JK{sub s} in both open clusters and the Kepler fields. We find good agreement overall, with some deviations between (J - K{sub s} )-based temperatures from the IRFM and both SDSS filter and other diagnostic IRFM color-temperature relationships above 6000 K. For field dwarfs, we find a mean shift toward hotter temperatures relative to the KIC, of order 215 K, in the regime where the IRFM scale is well defined (4000 K to 6500 K). This change is of comparable magnitude in both color systems and in spectroscopy for stars with T{sub eff} below 6000 K. Systematic differences between temperature estimators appear for hotter stars, and we define corrections to put the SDSS temperatures on the IRFM scale for them. When the theoretical dependence on gravity is accounted for, we find a similar temperature scale offset between the fundamental and KIC scales for giants. We demonstrate that statistical corrections to color-based temperatures from binaries are significant. Typical errors, mostly from uncertainties in extinction, are of order 100 K. Implications for other applications of the KIC are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic characterization of phase detection algorithms in phase-shifting interferometry

Phase-shifting interferometry (PSI) is the preferred non-contact method for profiling sub-nanometer surfaces. Based on monochromatic light interference, the method computes the surface profile from a set of interferograms collected at separate stepping positions. Errors in the estimated profile are introduced when these positions are not located correctly. In order to cope with this problem, various algorithms that minimize the effects of certain types of stepping errors (linear, sinusoidal, etc.) have been developed. Despite the relatively large number of algorithms suggested in the literature, there is no unified way of characterizing their performance when additional unaccounted random errors are present. Here, we suggest a procedure for quantifying the expected behavior of each algorithm in the presence of independent and identically distributed (i.i.d.) random stepping errors, which can occur in addition to the systematic errors for which the algorithm has been designed. As a result, the usefulness of this method derives from the fact that it can guide the selection of the best algorithm for specific measurement situations.

42 ENGINEERING↗

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TOWARD A NEW KIND OF ASTEROSEISMIC GRID FITTING

Recent developments in instrumentation (e.g., in particular the Kepler and CoRoT satellites) provide a new opportunity to improve the models of stellar pulsations. Surface layers, rotation, and magnetic fields imprint erratic frequency shifts, trends, and other non-random behavior in the frequency spectra. As our observational uncertainties become smaller, these are increasingly important and difficult to deal with using standard fitting techniques. To improve the models, new ways to compare their predictions with observations need to be conceived. In this paper, we present a completely probabilistic (Bayesian) approach to asteroseismic model fitting. It allows for varying degrees of prior mode identification, corrections for the discrete nature of the grid, and most importantly implements a treatment of systematic errors, such as the 'surface effects'. It removes the need to apply semi-empirical corrections to the observations prior to fitting them to the models and results in a consistent set of probabilities with which the model physics can be probed and compared. As an example, we show a detailed asteroseismic analysis of the Sun. We find a most probable solar age, including a 35 {+-} 5 million year pre-main-sequence phase, of 4.591 billion years, and initial element mass fractions of X{sub 0} = 0.72, Y{sub 0} = 0.264, Z{sub 0} = 0.016, consistent with recent asteroseismic and non-asteroseismic studies.

79 ASTRONOMY AND ASTROPHYSICS↗