Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inference”

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,045 records · Page 58

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗

Three-dimensional simulations of turbulent convective mixing in ONe and CO classical nova explosions

Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures of ~(1-4) × 10 8 K. During these events, about 10 -3 -10 -7 M ⊙ , enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the high concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. We performed three-dimensional simulations of mixing at the core-envelope interface during nova outbursts with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We also show that the presence of an ONe-rich substrate, as in “neon novae”, yields higher metallicity enhancements in the ejecta than CO-rich substrates (i.e., non-neon novae). Finally, a number of requirements and constraints for such 3D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum efficiency investigations of type-II InAs/GaSb midwave infrared superlattice photodetectors

We present in this paper a comparison between different type-II InAs/GaSb superlattice (T2SL) photodiodes and focal plane array (FPA) in the mid-wavelength infrared domain to understand which phenomenon drives the performances of the T2SL structure in terms of quantum efficiency (QE). Our measurements on test photodiodes suggest low minority carrier diffusion length in the “InAs-rich” design, which penalizes carriers' collection in this structure for low bias voltage and front side illumination. This analysis is completed by a comparison of the experimental data with a fully analytic model, which allows to infer a hole diffusion length shorter than 100 nm. In addition, measurements on a FPA with backside illumination are finally presented. Results show an average QE in the 3–4.7 μm window equal to 42% for U{sub bias} = −0.1 V, 77 K operating temperature and no anti-reflection coating. These measurements, completed by modulation transfer function and noise measurements, reveal that the InAs-rich design, despite a low hole diffusion length, is promising for high performance infrared imaging applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantum learning of classical stochastic processes: The completely positive realization problem

Among several tasks in Machine Learning, a specially important one is the problem of inferring the latent variables of a system and their causal relations with the observed behavior. A paradigmatic instance of this is the task of inferring the hidden Markov model underlying a given stochastic process. This is known as the positive realization problem (PRP), [L. Benvenuti and L. Farina, IEEE Trans. Autom. Control 49(5), 651–664 (2004)] and constitutes a central problem in machine learning. The PRP and its solutions have far-reaching consequences in many areas of systems and control theory, and is nowadays an important piece in the broad field of positive systems theory. We consider the scenario where the latent variables are quantum (i.e., quantum states of a finite-dimensional system) and the system dynamics is constrained only by physical transformations on the quantum system. The observable dynamics is then described by a quantum instrument, and the task is to determine which quantum instrument — if any — yields the process at hand by iterative application. We take as a starting point the theory of quasi-realizations, whence a description of the dynamics of the process is given in terms of linear maps on state vectors and probabilities are given by linear functionals on the state vectors. This description, despite its remarkable resemblance with the hidden Markov model, or the iterated quantum instrument, is however devoid of any stochastic or quantum mechanical interpretation, as said maps fail to satisfy any positivity conditions. The completely positive realization problem then consists in determining whether an equivalent quantum mechanical description of the same process exists. We generalize some key results of stochastic realization theory, and show that the problem has deep connections with operator systems theory, giving possible insight to the lifting problem in quotient operator systems. Our results have potential applications in quantum machine learning, device-independent characterization and reverse-engineering of stochastic processes and quantum processors, and more generally, of dynamical processes with quantum memory [M. Guţă, Phys. Rev. A 83(6), 062324 (2011); M. Guţă and N. Yamamoto, e-print http://arxiv.org/abs/1303.3771 (2013)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of preheat and mix on the fuel adiabat of an imploding capsule

We demonstrate the effect of preheat, hydrodynamic mix and vorticity on the adiabat of the deuterium-tritium (DT) fuel in fusion capsule experiments. We show that the adiabat of the DT fuel increases resulting from hydrodynamic mixing due to the phenomenon of entropy of mixture. An upper limit of mix, M clean =M DT ≥ 0:98 is found necessary to keep the DT fuel on a low adiabat. We demonstrate in this study that the use of a high adiabat for the DT fuel in theoretical analysis and with the aid of 1D code simulations could explain some aspects of 3D effects and mix in capsule implosion. Furthermore, we can infer from our physics model and the observed neutron images the adiabat of the DT fuel in the capsule and the amount of mix produced on the hot spot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using multiple neutron time of flight detectors to determine the hot spot velocity

An important diagnostic value of a shot at the National Ignition Facility is the resultant center-of-mass motion of the imploding capsule. This residual velocity reduces the efficiency of converting laser energy into plasma temperature. A new analysis method extracts the effective hot spot motion by using information from multiple neutron time-of-flight (nToF) lines-of-sight (LoSs). This technique fits a near Gaussian spectrum to the nToF scope traces and overcomes reliance on models to relate the plasma temperature to the mean energy of the emitted neutrons. This method requires having at least four nToF LoSs. The results of this analysis will be compared to an approach where each LoS is analyzed separately and a model is employed to infer the mean energy of the emitted neutrons

47 OTHER INSTRUMENTATION↗

Simple and sophisticated models of Taylor’s cylinder impact

We present a study based on the simplest possible models for Taylor’s cylinder impact problem, in addition to examina-tion of using convolutional neural networks (CNNs) to map cylinder profiles to strength calibrations. We find that the approximate treatments of Taylor and Hawkyard compare well with hydrodynamic simulations using an equivalent assumption of constant flow stress. The CNN models prove to be well suited to successfully infer parameterizations of the Preston-Tonks-Wallace model of plastic deformation based on the deformed profile of an impacted cylinder.

36 MATERIALS SCIENCE↗

SPECTRAL ENERGY DISTRIBUTIONS OF YOUNG STARS IN IC 348: THE ROLE OF DISKS IN ANGULAR MOMENTUM EVOLUTION OF YOUNG, LOW-MASS STARS

Theoretical work suggests that a young star's angular momentum content and rotation rate may be strongly influenced by magnetic interactions with its circumstellar disk. A generic prediction of these 'disk-locking' theories is that a disk-locked star will be forced to co-rotate with the Keplerian angular velocity of the inner edge of the disk; that is, the disk's inner-truncation radius should equal its co-rotation radius. These theories have also been interpreted to suggest a gross correlation between young stars' rotation periods and the structural properties of their circumstellar disks, such that slowly rotating stars possess close-in disks that enforce the star's slow rotation, whereas rapidly rotating stars possess anemic or evacuated inner disks that are unable to brake the stars and instead the stars spin up as they contract. To test these expectations, we model the spectral energy distributions (SEDs) of 33 young stars in IC 348 with known rotation periods and infrared excesses indicating the presence of circumstellar disks. For each star, we match the observed SED, typically sampling 0.6-8.0 {mu}m, to a grid of 200,000 pre-computed star+disk radiative transfer models, from which we infer the disk's inner-truncation radius. We then compare this truncation radius to the disk's co-rotation radius, calculated from the star's measured rotation period. We do not find obvious differences in the disk truncation radii of slow rotators versus rapid rotators. This holds true both at the level of whether close-in disk material is present at all, and in analyzing the precise location of the inner disk edge relative to the co-rotation radius among the subset of stars with close-in disk material. One interpretation is that disk locking is unimportant for the IC 348 stars in our sample. Alternatively, if disk locking does operate, then it must operate on both the slow and rapid rotators, potentially producing both spin-up and spin-down torques, and the transition from the disk-locked state to the disk-released state must occur more rapidly than the stellar contraction timescale.

79 ASTRONOMY AND ASTROPHYSICS↗

EXTENDED SUBMILLIMETER EMISSION OF THE GALACTIC CENTER AND NEAR-INFRARED/SUBMILLIMETER VARIABILITY OF ITS SUPERMASSIVE BLACK HOLE

The innermost tens of parsecs of our Galaxy are characterized by the presence of molecular cloud complexes surrounding Sgr A*, the radiative counterpart of the supermassive black hole ({approx}4 x 10{sup 6} M{sub sun}) at the Galactic center. We seek to distinguish the different physical mechanisms that dominate the molecular clouds at the Galactic center, with special emphasis on the circumnuclear disk (CND). We also want to study the energy flow and model the variable emission of Sgr A*. Our study is based on NIR and submillimeter (sub-mm) observations. Using sub-mm maps, we describe the complex morphology of the molecular clouds and the circumnuclear disk, along with their masses (of order 10{sup 5}-10{sup 6} M{sub sun}), and derive also the temperature and spectral index maps of the regions under study. We conclude that the average temperature of the dust is 14 {+-} 4 K. The spectral index map shows that the 20 and 50 km s{sup -1} clouds are dominated by dust emission. Comparatively, in the CND and its surroundings the spectral indices decrease toward Sgr A* and range between about 1 and -0.6. These values are mostly explained with a combination of dust, synchrotron, and free-free emission in different ratios. The presence of non-thermal emission also accounts for the apparent low temperatures derived in these areas, indicating their unreliability. The Sgr A* light curves show significant flux density excursions in both the NIR and sub-mm domains. We have defined a classification system to account for the NIR variability of Sgr A*. Also, we have modeled on the NIR/sub-mm events. From our modeling results we can infer a sub-mm emission delay with respect to the NIR; we argue that the delay is due to the adiabatic expansion of the synchrotron source components.

79 ASTRONOMY AND ASTROPHYSICS↗

FIELD TOPOLOGY ANALYSIS OF A LONG-LASTING CORONAL SIGMOID

We present the first field topology analysis based on nonlinear force-free field (NLFFF) models of a long-lasting coronal sigmoid observed in 2007 February with the X-Ray Telescope on Hinode. The NLFFF models are built with the flux rope insertion method and give the three-dimensional coronal magnetic field as constrained by observed coronal loop structures and photospheric magnetograms. Based on these models, we have computed horizontal maps of the current and the squashing factor Q for 25 different heights in the corona for all six days of the evolution of the region. We use the squashing factor to quantify the degree of change of the field line linkage and to identify prominent quasi-separatrix layers (QSLs). We discuss the major properties of these QSL maps and devise a way to pick out important QSLs since our calculation cannot reach high values of Q. The complexity in the QSL maps reflects the high degree of fragmentation of the photospheric field. We find main QSLs and current concentrations that outline the flux rope cavity and that become characteristically S-shaped during the evolution of the sigmoid. We note that, although intermittent bald patches exist along the length of the sigmoid during its whole evolution, the flux rope remains stable for several days. However, shortly after the topology of the field exhibits hyperbolic flux tubes (HFT) on February 7 and February 12 the sigmoid loses equilibrium and produces two B-class flares and associated coronal mass ejections (CMEs). The location of the most elevated part of the HFT in our model coincides with the inferred locations of the two flares. Therefore, we suggest that the presence of an HFT in a coronal magnetic configuration may be an indication that the system is ready to erupt. We offer a scenario in which magnetic reconnection at the HFT drives the system toward the marginally stable state. Once this state is reached, loss of equilibrium occurs via the torus instability, producing a CME.

79 ASTRONOMY AND ASTROPHYSICS↗

STELLAR POPULATIONS OF LYMAN BREAK GALAXIES AT z {approx_equal} 1-3 IN THE HST/WFC3 EARLY RELEASE SCIENCE OBSERVATIONS

We analyze the spectral energy distributions (SEDs) of Lyman break galaxies (LBGs) at z {approx_equal} 1-3 selected using the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) UVIS channel filters. These HST/WFC3 observations cover about 50 arcmin{sup 2} in the GOODS-South field as a part of the WFC3 Early Release Science program. These LBGs at z {approx_equal} 1-3 are selected using dropout selection criteria similar to high-redshift LBGs. The deep multi-band photometry in this field is used to identify best-fit SED models, from which we infer the following results: (1) the photometric redshift estimate of these dropout-selected LBGs is accurate to within few percent; (2) the UV spectral slope {beta} is redder than at high redshift (z > 3), where LBGs are less dusty; (3) on average, LBGs at z {approx_equal} 1-3 are massive, dustier, and more highly star forming, compared to LBGs at higher redshifts with similar luminosities (0.1L* {approx}< L {approx}< 2.5L*), though their median values are similar within 1{sigma} uncertainties. This could imply that identical dropout selection technique, at all redshifts, finds physically similar galaxies; and (4) the stellar masses of these LBGs are directly proportional to their UV luminosities with a logarithmic slope of {approx}0.46, and star formation rates are proportional to their stellar masses with a logarithmic slope of {approx}0.90. These relations hold true-within luminosities probed in this study-for LBGs from z {approx_equal} 1.5 to 5. The star-forming galaxies selected using other color-based techniques show similar correlations at z {approx_equal} 2, but to avoid any selection biases, and for direct comparison with LBGs at z > 3, a true Lyman break selection at z {approx_equal} 2 is essential. The future HST UV surveys, both wider and deeper, covering a large luminosity range are important to better understand LBG properties and their evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DEMOGRAPHICS OF THE GALAXIES HOSTING SHORT-DURATION GAMMA-RAY BURSTS

We present observations of the afterglows and host galaxies of three short-duration gamma-ray bursts (GRBs): 100625A, 101219A, and 110112A. We find that GRB 100625A occurred in a z = 0.452 early-type galaxy with a stellar mass of Almost-Equal-To 4.6 Multiplication-Sign 10{sup 9} M{sub Sun} and a stellar population age of Almost-Equal-To 0.7 Gyr, and GRB 101219A originated in a star-forming galaxy at z = 0.718 with a stellar mass of Almost-Equal-To 1.4 Multiplication-Sign 10{sup 9} M{sub Sun }, a star formation rate of Almost-Equal-To 16 M{sub Sun} yr{sup -1}, and a stellar population age of Almost-Equal-To 50 Myr. We also report the discovery of the optical afterglow of GRB 110112A, which lacks a coincident host galaxy to i {approx}> 26 mag, and we cannot conclusively identify any field galaxy as a possible host. From afterglow modeling, the bursts have inferred circumburst densities of Almost-Equal-To 10{sup -4}-1 cm{sup -3} and isotropic-equivalent gamma-ray and kinetic energies of Almost-Equal-To 10{sup 50}-10{sup 51} erg. These three events highlight the diversity of galactic environments that host short GRBs. To quantify this diversity, we use the sample of 36 Swift short GRBs with robust associations to an environment ({approx}1/2 of 68 short bursts detected by Swift to 2012 May) and classify bursts originating from four types of environments: late-type ( Almost-Equal-To 50%), early-type ( Almost-Equal-To 15%), inconclusive ( Almost-Equal-To 20%), and ''host-less'' (lacking a coincident host galaxy to limits of {approx}> 26 mag; Almost-Equal-To 15%). To find likely ranges for the true late- and early-type fractions, we assign each of the host-less bursts to either the late- or early-type category using probabilistic arguments and consider the scenario that all hosts in the inconclusive category are early-type galaxies to set an upper bound on the early-type fraction. We calculate most likely ranges for the late- and early-type fractions of Almost-Equal-To 60%-80% and Almost-Equal-To 20%-40%, respectively. We find no clear trend between gamma-ray duration and host type. We also find no change to the fractions when excluding events recently claimed as possible contaminants from the long GRB/collapsar population. Our reported demographics are consistent with a short GRB rate driven by both stellar mass and star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Microlensing of quasar ultraviolet iron emission

We measure the differential microlensing of the UV Fe II and Fe III emission line blends between 14 quasar image pairs in 13 gravitational lenses. We find that the UV iron emission is strongly microlensed in four cases with amplitudes comparable to that of the continuum. Statistically modeling the magnifications, we infer a typical size of r{sub s}∼4√(M/M{sub ⊙}) light-days for the Fe line-emitting regions, which is comparable to the size of the region generating the UV continuum (∼3-7 light-days). This may indicate that a significant part of the UV Fe II and Fe III emission originates in the quasar accretion disk.

79 ASTRONOMY AND ASTROPHYSICS↗

New near-infrared period–luminosity–metallicity relations for RR Lyrae stars and the outlook for GAIA

We present results of the analysis of 70 RR Lyrae stars located in the bar of the Large Magellanic Cloud (LMC). Combining the spectroscopically determined metallicity of these stars from the literature with precise periods from the OGLE III catalog and multi-epoch K{sub s} photometry from the VISTA survey of the Magellanic Clouds system, we derive a new near-infrared period–luminosity–metallicity (PL{sub K{sub s}}Z) relation for RR Lyrae variables. In order to fit the relation we use a fitting method developed specifically for this study. The zero-point of the relation is estimated two different ways: by assuming the value of the distance to the LMC and by using Hubble Space Telescope parallaxes of five RR Lyrae stars in the Milky Way (MW). The difference in distance moduli derived by applying these two approaches is ∼0.2 mag. To investigate this point further we derive the PL{sub K{sub s}}Z relation based on 23 MW RR Lyrae stars that had been analyzed in Baade–Wesselink studies. We compared the derived PL{sub K{sub s}}Z relations for RR Lyrae stars in the MW and LMC. Slopes and zero-points are different, but still consistent within the errors. The shallow slope of the metallicity term is confirmed by both LMC and MW variables. The astrometric space mission Gaia is expected to provide a huge contribution to the determination of the RR Lyrae PL{sub K{sub s}}Z relation; however, calculating an absolute magnitude from the trigonometric parallax of each star and fitting a PL{sub K{sub s}}Z relation directly to period and absolute magnitude leads to biased results. We present a tool to achieve an unbiased solution by modeling the data and inferring the slope and zero-point of the relation via statistical methods.

79 ASTRONOMY AND ASTROPHYSICS↗

CONSTRUCTING A FLEXIBLE LIKELIHOOD FUNCTION FOR SPECTROSCOPIC INFERENCE

We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically addresses the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line “outliers.” By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high-resolution V-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate-resolution K-band spectrum of Gliese 51, an M5 field dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating multilevel Markov Chain Monte Carlo using machine learning models

Here, this work presents an efficient approach for accelerating multilevel Markov Chain Monte Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning models. While conventional techniques for large-scale Bayesian inference often substitute computationally expensive high-fidelity models with machine learning models, thereby introducing approximation errors, our approach offers a computationally efficient alternative by augmenting high-fidelity models with low-fidelity ones within a hierarchical framework. The multilevel approach utilizes the low-fidelity machine learning model (MLM) for inexpensive evaluation of proposed samples thereby improving the acceptance of samples by the high-fidelity model. The hierarchy in our multilevel algorithm is derived from geometric multigrid hierarchy. We utilize an MLM to accelerate the coarse level sampling. Training machine learning model for the coarsest level significantly reduces the computational cost associated with generating training data and training the model. We present an MCMC algorithm to accelerate the coarsest level sampling using MLM and account for the approximation error introduced. We provide theoretical proofs of detailed balance and demonstrate that our multilevel approach constitutes a consistent MCMC algorithm. Additionally, we derive the expression for cost reduction due to machine learning model to facilitate cost analysis of the hierarchical sampling algorithm. Our technique is demonstrated on a standard benchmark inference problem in groundwater flow, where we estimate the probability density of a quantity of interest using a four-level MCMC algorithm. Our proposed algorithm accelerates multilevel sampling by a factor of two while achieving similar accuracy compared to sampling using the standard multilevel algorithm.

97 MATHEMATICS AND COMPUTING↗

Modulations of the cosmic muon signal in ten years of Borexino data

We have measured the flux of cosmic muons in the Laboratori Nazionali del Gran Sasso at 3800\,m\,w.e. to be (3.432±0.003)∙10 -4 m -2 s -1 based on ten years of Borexino data acquired between May 2007 and May 2017. A seasonal modulation with a period of (366.3±0.6)d and a relative amplitude of (1.36±0.04)% is observed. The phase is measured to be (181.7±0.4)d, corresponding to a maximum at the 1st of July. Using data inferred from global atmospheric models, we show the muon flux to be positively correlated with the atmospheric temperature and measure the effective temperature coefficient α T =0.90±0.02. The origin of cosmic muons from pion and kaon decays in the atmosphere allows to interpret the effective temperature coefficient as an indirect measurement of the atmospheric kaon-to-pion production ratio rK/π=0.11 +0.11 -0.07 for primary energies above 18TeV. We find evidence for a long-term modulation of the muon flux with a period of ~3000d and a maximum in June 2012 that is not present in the atmospheric temperature data. A possible correlation between this modulation and the solar activity is investigated. The cosmogenic neutron production rate is found to show a seasonal modulation in phase with the cosmic muon flux but with an increased amplitude of (2.6±0.4)%.

79 ASTRONOMY AND ASTROPHYSICS↗

X-RAY DETERMINATION OF THE VARIABLE RATE OF MASS ACCRETION ONTO TW HYDRAE

Diagnostics of electron temperature (T{sub e} ), electron density (n{sub e} ), and hydrogen column density (N{sub H}) from the Chandra High Energy Transmission Grating spectrum of He-like Ne IX in TW Hydrae (TW Hya), in conjunction with a classical accretion model, allow us to infer the accretion rate onto the star directly from measurements of the accreting material. The new method introduces the use of the absorption of Ne IX lines as a measure of the column density of the intervening, accreting material. On average, the derived mass accretion rate for TW Hya is 1.5 Multiplication-Sign 10{sup -9} M{sub Sun} yr{sup -1}, for a stellar magnetic field strength of 600 G and a filling factor of 3.5%. Three individual Chandra exposures show statistically significant differences in the Ne IX line ratios, indicating changes in N{sub H}, T{sub e} , and n{sub e} by factors of 0.28, 1.6, and 1.3, respectively. In exposures separated by 2.7 days, the observations reported here suggest a five-fold reduction in the accretion rate. This powerful new technique promises to substantially improve our understanding of the accretion process in young stars.

79 ASTRONOMY AND ASTROPHYSICS↗