Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consequence”

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,171 records · Page 65

CONTRACTING AND ERUPTING COMPONENTS OF SIGMOIDAL ACTIVE REGIONS

It has recently been noted that solar eruptions can be associated with the contraction of coronal loops that are not involved in magnetic reconnection processes. In this paper, we investigate five coronal eruptions originating from four sigmoidal active regions, using high-cadence, high-resolution narrowband EUV images obtained by the Solar Dynamic Observatory (SDO). The magnitudes of the flares associated with the eruptions range from GOES class B to class X. Owing to the high-sensitivity and broad temperature coverage of the Atmospheric Imaging Assembly (AIA) on board SDO, we are able to identify both the contracting and erupting components of the eruptions: the former is observed in cold AIA channels as the contracting coronal loops overlying the elbows of the sigmoid, and the latter is preferentially observed in warm/hot AIA channels as an expanding bubble originating from the center of the sigmoid. The initiation of eruption always precedes the contraction, and in the energetically mild events (B- and C-flares), it also precedes the increase in GOES soft X-ray fluxes. In the more energetic events, the eruption is simultaneous with the impulsive phase of the nonthermal hard X-ray emission. These observations confirm that loop contraction is an integrated process in eruptions with partially opened arcades. The consequence of contraction is a new equilibrium with reduced magnetic energy, as the contracting loops never regain their original positions. The contracting process is a direct consequence of flare energy release, as evidenced by the strong correlation of the maximal contracting speed, and strong anti-correlation of the time delay of contraction relative to expansion, with the peak soft X-ray flux. This is also implied by the relationship between contraction and expansion, i.e., their timing and speed.

79 ASTRONOMY AND ASTROPHYSICS↗

THE X-FACTOR IN GALAXIES. II. THE MOLECULAR-HYDROGEN-STAR-FORMATION RELATION

There is ample observational evidence that the star formation rate (SFR) surface density, {Sigma}{sub SFR}, is closely correlated with the surface density of molecular hydrogen, {Sigma}{sub H{sub 2}}. This empirical relation holds both for galaxy-wide averages and for individual {approx}>kpc sized patches of the interstellar medium, but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed {Sigma}{sub H{sub 2}}-{Sigma}{sub SFR} relation using a set of cosmological, galaxy formation simulations with a peak resolution of {approx}100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on {approx}100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H{sub 2}. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and {Sigma}{sub H{sub 2}}-{Sigma}{sub SFR} relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low ({approx}< 20 M {sub Sun} pc{sup -2}) H{sub 2} surface densities on sub-kpc scales. Our models also predict that the inferred {Sigma}{sub H{sub 2}}-{Sigma}{sub SFR} relation steepens at high H{sub 2} surface densities as a result of the surface density dependence of the CO/H{sub 2} conversion factor. Finally, we show that on sub-kpc scales most of the scatter of the relation is a consequence of discreteness effects of the star formation process. In contrast, variations of the CO/H{sub 2} conversion factor are responsible for most of the scatter measured on super-kpc scales.

79 ASTRONOMY AND ASTROPHYSICS↗

ON ABSORPTION BY CIRCUMSTELLAR DUST, WITH THE PROGENITOR OF SN 2012aw AS A CASE STUDY

We use the progenitor of SN 2012aw to illustrate the consequences of modeling circumstellar dust using Galactic (interstellar) extinction laws that (1) ignore dust emission in the near-IR and beyond, (2) average over dust compositions, and (3) mischaracterize the optical/UV absorption by assuming that scattered photons are lost to the observer. The primary consequences for the progenitor of SN 2012aw are that both the luminosity and the absorption are significantly overestimated. In particular, the stellar luminosity is most likely in the range 10{sup 4.8} < L {sub *}/L {sub Sun} < 10{sup 5.0} and the star was not extremely massive for a Type IIP progenitor, with M {sub *} < 15 M {sub Sun }. Given the properties of the circumstellar dust and the early X-ray/radio detections of SN 2012aw, the star was probably obscured by an ongoing wind with M-dot {approx}10{sup -5.5} to 10{sup -5.0} M {sub Sun} yr{sup -1} at the time of the explosion, roughly consistent with the expected mass-loss rates for a star of its temperature (T{sub *} {approx_equal} 3600{sup +300} {sub -200} K) and luminosity. In the spirit of Galactic extinction laws, we supply simple interpolation formulae for circumstellar extinction by dusty graphitic and silicate shells as a function of wavelength ({lambda} {>=} 0.3 {mu}m) and total (absorption plus scattering) V-band optical depth ({tau}{sub V} {<=} 20). These do not include the contributions of dust emission, but provide a simple, physical alternative to incorrectly using interstellar extinction laws.

79 ASTRONOMY AND ASTROPHYSICS↗

A THEORETICAL COLOR-VELOCITY CORRELATION FOR SUPERNOVAE ASSOCIATED WITH GAMMA-RAY BURSTS

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These allow us to make self-consistent predictions for the orientation dependence of both color evolution and spectral features. We compare models with different degrees of asphericity and metallicity and find significant observable consequences of both. In aspherical models, we find spectral and light curve features that vary systematically with observer orientation. In particular, we find that the early-phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on observer orientation, including the velocity of the Si II 6355 A line. Consequently, our models predict a correlation between line velocity and color that could assist the identification of supernovae associated with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree of asphericity, the metallicity, and the epoch of observation but we find that it is always present and acts in the same direction.

79 ASTRONOMY AND ASTROPHYSICS↗

THE STABILITY OF WEAKLY COLLISIONAL PLASMAS WITH THERMAL AND COMPOSITION GRADIENTS

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport phenomena, and long-term evolution of weakly collisional, magnetized plasmas which are stratified in temperature. The insights gained via these studies have led to a significant improvement of our understanding of the processes that determine the physical evolution and observational properties of the intracluster medium (ICM) permeating galaxy clusters. These studies have been carried out under the assumption that the ICM is a homogeneous medium. This, however, might not be a good approximation if heavy elements are able to sediment in the inner region of the galaxy cluster. Motivated by the need to obtain a more complete picture of the dynamical properties of the ICM, we analyze the stability of a weakly collisional, magnetized plane-parallel atmosphere which is stratified in both temperature and composition. This allows us to discuss for the first time the dynamics of weakly collisional environments where heat conduction, momentum transport, and ion-diffusion are anisotropic with respect to the direction of the magnetic field. We show that depending on the relative signs and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous media. We also find that there are new modes which are driven by heat conduction and particle diffusion. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could be alleviated if the mean molecular weight gradient is steep enough, i.e., ({nabla}{mu})/{mu} > ({nabla}T)/T. This study constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of heavy elements, and its consequences for the long-term evolution and observational signatures of the ICM in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

RED SUPERGIANT STARS AS COSMIC ABUNDANCE PROBES. II. NLTE EFFECTS IN J-BAND SILICON LINES

Medium-resolution J-band spectroscopy of individual red supergiant stars is a promising tool to investigate the chemical composition of the young stellar population in star-forming galaxies. As a continuation of recent work on iron and titanium, detailed non-LTE (NLTE) calculations are presented to investigate the influence of NLTE on the formation of silicon lines in the J-band spectra of red supergiants. Substantial effects are found resulting in significantly stronger absorption lines of neutral silicon in NLTE. As a consequence, silicon abundances determined in NLTE are significantly smaller than in local thermodynamic equilibrium (LTE) with the NLTE abundance corrections varying smoothly between -0.4 dex and -0.1 dex for effective temperatures between 3400 K and 4400 K. The effects are largest at low metallicity. The physical reasons behind the NLTE effects and the consequences for extragalactic J-band abundance studies are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

TYPE Ia SINGLE DEGENERATE SURVIVORS MUST BE OVERLUMINOUS

In the single-degenerate (SD) channel of a Type Ia supernovae (SNe Ia) explosion, a main-sequence (MS) donor star survives the explosion but it is stripped of mass and shock heated. An essentially unavoidable consequence of mass loss during the explosion is that the companion must have an overextended envelope after the explosion. While this has been noted previously, it has not been strongly emphasized as an inevitable consequence. We calculate the future evolution of the companion by injecting 2-6 Multiplication-Sign 10{sup 47} erg into the stellar evolution model of a 1 M{sub Sun} donor star based on the post-explosion progenitors seen in simulations. We find that, due to the Kelvin-Helmholtz collapse of the envelope, the companion must become significantly more luminous (10-10{sup 3} L{sub Sun }) for a long period of time (10{sup 3}-10{sup 4} yr). The lack of such a luminous ''leftover'' star in the LMC supernova remnant SNR 0609-67.5 provides another piece of evidence against the SD scenario. We also show that none of the stars proposed as the survivors of the Tycho supernova, including Tycho G, could plausibly be the donor star. Additionally, luminous donors closer than {approx}10 Mpc should be observable with the Hubble Space Telescope starting {approx}2 yr post-peak. Such systems include SN 1937C, SN 1972E, SN 1986G, and SN 2011fe. Thus, the SD channel is already ruled out for at least two nearby SNe Ia and can easily be tested for a number of additional ones. We also discuss similar implications for the companions of core-collapse SNe.

79 ASTRONOMY AND ASTROPHYSICS↗

THERMAL AND DYNAMICAL PROPERTIES OF GAS ACCRETING ONTO A SUPERMASSIVE BLACK HOLE IN AN ACTIVE GALACTIC NUCLEUS

We study stability of gas accretion in active galactic nuclei (AGNs). Our grid-based simulations cover a radial range from 0.1 to 200 pc, which may enable linking the galactic/cosmological simulations with small-scale black hole (BH) accretion models within a few hundreds of Schwarzschild radii. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive BH of 10{sup 8} M{sub Sun }. Our theoretical estimates and simulations show that for the X-ray luminosity, L{sub X} {approx} 0.008 L{sub Edd}, the gas is thermally and convectively unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then nonlinearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For L{sub X} = 0.015 L{sub Edd} or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For L{sub X} < 0.015 L{sub Edd}, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. who used smoothed particle hydrodynamic (SPH) simulations to tackle the same problem. Here, however, because we use a grid-based code to solve equations in one dimension and two dimensions, we are able to follow the gas dynamics at much higher spacial resolution and for longer time compared with the three-dimensional SPH simulations. One of the new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outward within the hot outflow. Therefore, these simulations may serve as an attractive model for the so-called narrow-line region in AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

SOLAR ENERGETIC PARTICLE MODULATIONS ASSOCIATED WITH COHERENT MAGNETIC STRUCTURES

In situ observations of solar energetic particles (SEPs) often show rapid variations of their intensity profile, affecting all energies simultaneously, without time dispersion. A previously proposed interpretation suggests that these modulations are directly related to the presence of magnetic structures with a different magnetic topology. However, no compelling evidence of local changes in magnetic field or in plasma parameters during SEP modulations has been reported. In this paper, we performed a detailed analysis of SEP events and we found several signatures in the local magnetic field and/or plasma parameters associated with SEP modulations. The study of magnetic helicity allowed us to identify magnetic boundaries, associated with variations of plasma parameters, which are thought to represent the borders between adjacent magnetic flux tubes. It is found that SEP dispersionless modulations are generally associated with such magnetic boundaries. Consequently, we support the idea that SEP modulations are observed when the spacecraft passes through magnetic flux tubes, filled or devoid of SEPs, which are alternatively connected and not connected with the flare site. In other cases, we found SEP dropouts associated with large-scale magnetic holes. A possible generation mechanism suggests that these holes are formed in the high solar corona as a consequence of magnetic reconnection. This reconnection process modifies the magnetic field topology, and therefore, these holes can be magnetically isolated from the surrounding plasma and could also explain their association with SEP dropouts.

79 ASTRONOMY AND ASTROPHYSICS↗

THE IMPACT OF METALLICITY ON THE RATE OF TYPE Ia SUPERNOVAE

The metallicity of a star strongly affects both its evolution and the properties of the stellar remnant that results from its demise. It is generally accepted that stars with initial masses below {approx}8 M{sub Sun} leave behind white dwarfs and that some sub-population of these lead to Type Ia supernovae (SNe Ia). However, it is often tacitly assumed that metallicity has no effect on the rate of SNe Ia. We propose that a consequence of the effects of metallicity is to significantly increase the SN Ia rate in lower-metallicity galaxies, in contrast to previous expectations. This is because lower-metallicity stars leave behind higher-mass white dwarfs, which should be easier to bring to explosion. We first model SN Ia rates in relation to galaxy masses and ages alone, finding that the elevation in the rate of SNe Ia in lower-mass galaxies measured by Lick Observatory SN Search is readily explained. However, we then see that models incorporating this effect of metallicity agree just as well. Using the same parameters to estimate the cosmic SN Ia rate, we again find good agreement with data up to z Almost-Equal-To 2. We suggest that this degeneracy warrants more detailed examination of host galaxy metallicities. We discuss additional implications, including for hosts of high-z SNe Ia, the SN Ia delay time distribution, super-Chandrasekhar SNe, and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

CARBON DEFLAGRATION IN TYPE Ia SUPERNOVA. I. CENTRALLY IGNITED MODELS

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, but is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of {approx}100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed-within the range expected from turbulent combustion-or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g., {sup 56}Ni of 0.48-0.56 M{sub Sun }) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

79 ASTRONOMY AND ASTROPHYSICS↗

INTERNAL GRAVITY WAVES IN MASSIVE STARS: ANGULAR MOMENTUM TRANSPORT

We present numerical simulations of internal gravity waves (IGW) in a star with a convective core and extended radiative envelope. We report on amplitudes, spectra, dissipation, and consequent angular momentum transport by such waves. We find that these waves are generated efficiently and transport angular momentum on short timescales over large distances. We show that, as in Earth's atmosphere, IGW drive equatorial flows which change magnitude and direction on short timescales. These results have profound consequences for the observational inferences of massive stars, as well as their long term angular momentum evolution. We suggest IGW angular momentum transport may explain many observational mysteries, such as: the misalignment of hot Jupiters around hot stars, the Be class of stars, Ni enrichment anomalies in massive stars, and the non-synchronous orbits of interacting binaries.

79 ASTRONOMY AND ASTROPHYSICS↗

INTO THE LAIR: GRAVITATIONAL-WAVE SIGNATURES OF DARK MATTER

The nature and properties of dark matter (DM) are both outstanding issues in physics. Besides clustering in halos, the universal character of gravity implies that self-gravitating compact DM configurations-predicted by various models-might be spread throughout the universe. Their astrophysical signature can be used to probe fundamental particle physics, or to test alternative descriptions of compact objects in active galactic nuclei. Here, we discuss the most promising dissection tool of such configurations: the inspiral of a compact stellar-size object and consequent gravitational-wave (GW) emission. The inward motion of this ''test probe'' encodes unique information about the nature of the supermassive configuration. When the probe travels through some compact region we show, within a Newtonian approximation, that the quasi-adiabatic inspiral is mainly driven by DM accretion and by dynamical friction, rather than by radiation reaction. When accretion dominates, the frequency and amplitude of the GW signal produced during the latest stages of the inspiral are nearly constant. In the exterior region we study a model in which the inspiral is driven by GW and scalar-wave emission, described at a fully relativistic level. Resonances in the energy flux appear whenever the orbital frequency matches the effective mass of the DM particle, corresponding to the excitation of the central object's quasinormal frequencies. Unexpectedly, these resonances can lead to large dephasing with respect to standard inspiral templates, to such an extent as to prevent detection with matched filtering techniques. We discuss some observational consequences of these effects for GW detection.

79 ASTRONOMY AND ASTROPHYSICS↗

THE BEHAVIOR OF TRANSVERSE WAVES IN NONUNIFORM SOLAR FLUX TUBES. I. COMPARISON OF IDEAL AND RESISTIVE RESULTS

Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius of the tube, i.e., the so-called thin boundary approximation. Here, we develop a general analytic method to compute the dispersion relation and the eigenfunctions of ideal MHD waves in pressureless flux tubes with transversely nonuniform layers of arbitrary thickness. Results for kink waves are produced and compared with fully numerical resistive MHD eigenvalue computations in the limit of small resistivity. We find that the frequency and resonant damping rate are the same in both ideal and resistive cases. The actual results for thick nonuniform layers deviate from the behavior predicted in the thin boundary approximation and strongly depend on the shape of the nonuniform layer. The eigenfunctions in ideal MHD are very different from those in resistive MHD. The ideal eigenfunctions display a global character regardless of the thickness of the nonuniform layer, while the resistive eigenfunctions are localized around the resonance and are indistinguishable from those of ordinary resistive Alfvén modes. Consequently, the spatial distribution of wave energy in the ideal and resistive cases is dramatically different. This poses a fundamental theoretical problem with clear observational consequences.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic and gravitational disk-star interactions: an interdependence of PMS stellar rotation rates and spin-orbit misalignments

The presence of giant gaseous planets that reside in close proximity to their host stars, i.e., hot Jupiters, may be a consequence of large-scale radial migration through the protoplanetary nebulae. Within the framework of this picture, significant orbital obliquities characteristic of a substantial fraction of such planets can be attributed to external torques that perturb the natal disks out of alignment with the spin axes of their host stars. Therefore, the acquisition of orbital obliquity likely exhibits sensitive dependence on the physics of disk-star interactions. Here, we analyze the primordial excitation of spin-orbit misalignment of Sun-like stars in light of disk-star angular momentum transfer. We begin by calculating the stellar pre-main-sequence rotational evolution, accounting for spin-up due to gravitational contraction and accretion as well as spin-down due to magnetic star-disk coupling. We devote particular attention to angular momentum transfer by accretion, and show that while generally subdominant to gravitational contraction, this process is largely controlled by the morphology of the stellar magnetic field (that is, specific angular momentum accreted by stars with octupole-dominated surface fields is smaller than that accreted by dipole-dominated stars by an order of magnitude). Subsequently, we examine the secular spin-axis dynamics of disk-bearing stars, accounting for the time-evolution of stellar and disk properties, and demonstrate that misalignments are preferentially excited in systems where stellar rotation is not overwhelmingly rapid. Moreover, we show that the excitation of spin-orbit misalignment occurs impulsively through an encounter with a resonance between the stellar precession frequency and the disk-torquing frequency. Cumulatively, the model developed herein opens up a previously unexplored avenue toward understanding star-disk evolution and its consequences in a unified manner.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on cosmology from the cosmic microwave background power spectrum of the 2500 deg{sup 2} SPT-SZ survey

We explore extensions to the ΛCDM cosmology using measurements of the cosmic microwave background (CMB) from the recent SPT-SZ survey, along with data from WMAP7 and measurements of H {sub 0} and baryon acoustic oscillation (BAO). We check for consistency within ΛCDM between these data sets, and find some tension. The CMB alone gives weak support to physics beyond ΛCDM, due to a slight trend relative to ΛCDM of decreasing power toward smaller angular scales. While it may be due to statistical fluctuation, this trend could also be explained by several extensions. We consider running of the primordial spectral index (dn{sub s} /dln k), as well as two extensions that modify the damping tail power (the primordial helium abundance Y{sub p} and the effective number of neutrino species N {sub eff}) and one that modifies the large-scale power due to the integrated Sachs-Wolfe effect (the sum of neutrino masses ∑m {sub ν}). These extensions have similar observational consequences and are partially degenerate when considered simultaneously. Of the six one-parameter extensions considered, we find CMB to have the largest preference for dn{sub s} /dln k with –0.046 < dn{sub s} /dln k < –0.003 at 95% confidence, which strengthens to a 2.7σ indication of dn{sub s} /dln k < 0 from CMB+BAO+H {sub 0}. Detectable dn{sub s} /dln k ≠ 0 is difficult to explain in the context of single-field, slow-roll inflation models. We find N {sub eff} = 3.62 ± 0.48 for the CMB, which tightens to N {sub eff} = 3.71 ± 0.35 from CMB+BAO+H {sub 0}. Larger values of N {sub eff} relieve the mild tension between CMB, BAO, and H {sub 0}. When the Sunyaev-Zel'dovich selected galaxy cluster abundances (SPT{sub CL}) data are also included, we obtain N {sub eff} = 3.29 ± 0.31. Allowing for ∑m {sub ν} gives a 3.0σ detection of ∑m {sub ν} > 0 from CMB+BAO+H {sub 0} +SPT{sub CL}. The median value is (0.32 ± 0.11) eV, a factor of six above the lower bound set by neutrino oscillation observations. All data sets except H {sub 0} show some preference for massive neutrinos; data combinations including H {sub 0} favor nonzero masses only if BAO data are also included. We also constrain the two-parameter extensions N {sub eff} + ∑m {sub ν} and N {sub eff} + Y{sub p} to explore constraints on additional light species and big bang nucleosynthesis, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Time evolution of a viscous protoplanetary disk with a free geometry: Toward a more self-consistent picture

Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as r {sup –1} or aspect ratios of about 0.03-0.23. Accretion rates are also recurrently found around 10{sup –8}-10{sup –6} M {sub ☉} yr{sup –1} for disks that have already evolved. Several models have been developed in order to recover these properties. However, most of them usually simplify the disk geometry if not its mid-plane temperature. This has major consequences for modeling the disk evolution over millions of years and consequently planet migration. In the present paper, we develop a viscous evolution hydrodynamical numerical code that simultaneously determines the disk photosphere geometry and the mid-plane temperature. We then compare our results of long-term simulations with similar simulations of disks with a constrained geometry along the Chiang and Goldreich prescription (d lnH/d lnr = 9/7). We find that the constrained geometry models provide a good approximation of the disk surface density evolution. However, they differ significantly regarding the temperature-time evolution. In addition, we find that shadowed regions naturally appear at the transition between viscously dominated and radiation-dominated regions that falls in the region of planetary formation. We show that χ (photosphere height to pressure scale height ratio) cannot be considered a constant, which is consistent with the findings of Watanabe and Lin. Comparisons with observations show that all disks naturally evolve toward a shallow surface density disk (Σ∝r {sup –1}). The mass flux across the disk typically stabilizes in about 1 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

Supervised detection of anomalous light curves in massive astronomical catalogs

The development of synoptic sky surveys has led to a massive amount of data for which resources needed for analysis are beyond human capabilities. In order to process this information and to extract all possible knowledge, machine learning techniques become necessary. Here we present a new methodology to automatically discover unknown variable objects in large astronomical catalogs. With the aim of taking full advantage of all information we have about known objects, our method is based on a supervised algorithm. In particular, we train a random forest classifier using known variability classes of objects and obtain votes for each of the objects in the training set. We then model this voting distribution with a Bayesian network and obtain the joint voting distribution among the training objects. Consequently, an unknown object is considered as an outlier insofar it has a low joint probability. By leaving out one of the classes on the training set, we perform a validity test and show that when the random forest classifier attempts to classify unknown light curves (the class left out), it votes with an unusual distribution among the classes. This rare voting is detected by the Bayesian network and expressed as a low joint probability. Our method is suitable for exploring massive data sets given that the training process is performed offline. We tested our algorithm on 20 million light curves from the MACHO catalog and generated a list of anomalous candidates. After analysis, we divided the candidates into two main classes of outliers: artifacts and intrinsic outliers. Artifacts were principally due to air mass variation, seasonal variation, bad calibration, or instrumental errors and were consequently removed from our outlier list and added to the training set. After retraining, we selected about 4000 objects, which we passed to a post-analysis stage by performing a cross-match with all publicly available catalogs. Within these candidates we identified certain known but rare objects such as eclipsing Cepheids, blue variables, cataclysmic variables, and X-ray sources. For some outliers there was no additional information. Among them we identified three unknown variability types and a few individual outliers that will be followed up in order to perform a deeper analysis.

79 ASTRONOMY AND ASTROPHYSICS↗