The Kepler mission observed many thousands of red giants. The long time series, some as long as the mission itself, have allowed us to study red giants with unprecedented detail. Given that red giants are intrinsically luminous, and hence can be observed from very large distances, knowing the properties of red giants, in particular ages, is of immense value for studies of the formation and evolution of the Galaxy, an endeavor known as "Galactic archaeology". In this article we review what we have learned about red giants using asteroseismic data. We start with the properties of the power spectrum and move on to internal structure and dynamics of these stars; we also touch upon unsolved issues in red-giant asteroseismology and the prospects of making further progress in understanding these stars.

The formation of circumstellar disks is investigated using three-dimensional resistive magnetohydrodynamic simulations, in which the initial prestellar cloud has a misaligned rotation axis with respect to the magnetic field. We examine the effects of (i) the initial angle difference between the global magnetic field and the cloud rotation axis ($\theta_0$) and (ii) the ratio of the thermal to gravitational energy ($\alpha_0$). We study $16$ models in total and calculate the cloud evolution until $\sim \! 5000$ yr after protostar formation. Our simulation results indicate that an initial non-zero $\theta_0$ ($> 0$) promotes the disk formation but tends to suppress the outflow driving, for models that are moderately gravitationally unstable, $\alpha_0 \lesssim 1$. In these models, a large-sized rotationally-supported disk forms and a weak outflow appears, in contrast to a smaller disk and strong outflow in the aligned case ($\theta_0 = 0$). Furthermore, we find that when the initial cloud is highly unstable with small $\alpha_0$, the initial angle difference $\theta_0$ does not significantly affect the disk formation and outflow driving.

In the oscillation spectra of giant stars, nonradial modes may be seen to undergo avoided crossings, which produce a characteristic "mode bumping" of the otherwise uniform asymptotic p- and g-mode patterns in their respective echelle diagrams. Avoided crossings evolve very quickly relative to typical observational errors, and are therefore extremely useful in determining precise ages of stars, particularly in subgiants. This phenomenon is caused by coupling between modes in the p- and g-mode cavities that are near resonance with each other. Most theoretical analyses of the coupling between these mode cavities rely on the JWKB approach, which is strictly speaking inapplicable for the low-order g-modes observed in subgiants, or the low-order p-modes seen in very evolved red giants. We present both a nonasymptotic prescription for isolating the two mode cavities, as well as a perturbative (and also nonasymptotic) description of the coupling between them, which we show to hold good for the low-order g- and p-modes in these physical situations. Finally, we discuss how these results may be applied to modelling subgiant stars and determining their global properties from oscillation frequencies. We also make our code for all of these computations publicly available.

EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.

Transport of angular momentum in stellar interiors is currently not well understood. Asteroseismology can provide us with estimates of internal rotation of stars and thereby advances our understanding of angular momentum transport. We can measure core-rotation rates in red-giant stars and we can place upper bounds on surface-rotation rates using measurements of dipole ($l=1$) modes. Here, we aim to determine the theoretical sensitivity of modes of different spherical degree towards the surface rotation. Additionally, we aim to identify modes that can potentially add sensitivity at intermediate radii. We used asteroseismic rotational inversions to probe the internal stellar rotation profiles in red-giant models from the base of the red-giant branch up to the luminosity bump. We used the inversion method of multiplicative optimally localised averages (MOLA) to assess how well internal and surface rotation rates can be recovered from different mode sets and different synthetic rotation profiles. We confirm that dipole mixed modes are sufficient to set constraints on the average core-rotation rates in red giants. However, surface-rotation rates estimated with only dipole mixed modes are contaminated by the core rotation. We show that the sensitivity to the surface rotation decreases from the base of the red-giant branch until it reaches a minimum at 0.6-0.8$L_\text{bump}$ due to a glitch in the buoyancy frequency. Thereafter a narrow range of increased surface sensitivity just below the bump luminosity exists. Quadrupole and octopole modes have more sensitivity in the outer parts of the star. If observed, quadrupole and octopole modes enable us to distinguish between differential and solid body rotation in the convection zone. To obtain accurate estimates of rotation rates at intermediate radii, acoustic oscillation modes with a spherical degree of $l\approx10$ are needed.

This paper presents the results of a set of radiative hydrodynamic (RHD) simulations of convection in the near-surface regions of a rapidly rotating star. The simulations use microphysics consistent with stellar models, and include the effects of realistic convection and radiative transfer. We find that the overall effect of rotation is to reduce the strength of turbulence. The combination of rotation and radiative cooling creates a zonal velocity profile in which the motion of fluid parcels near the surface is independent of rotation. Their motion is controlled by the strong up and down flows generated by radiative cooling. The fluid parcels in the deeper layers, on the other hand, are controlled by rotation.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

A determination of the mass function (MF) of stellar clusters can be quite dependent on the range of measured masses, the fitting technique, and the analytic function that is being fit to the data. Here, we use HST/WFPC2 data of NGC 1711, a stellar cluster in the Large Magellanic Cloud, as a test case to explore a range of possible determinations of the MF from a single dataset. We employ the analytic modified lognormal power-law (MLP) distribution, a hybrid function that has a peaked lognormal-like body and a power-law tail at intermediate and high masses. A fit with the MLP has the advantage that the resulting best-fit function can be either a hybrid function, a pure lognormal, or a pure power law, in different limits of the function. The completeness limit for the observations means that the data contains masses above $\sim 0.90\,M_{\odot}$. In this case, the MLP fits yield essentially a pure power-law MF. We demonstrate that the nonlinear regression/least-squares approach is not justified since the underlying assumptions are not satisfied. By using maximum likelihood estimation, which is independent of binning, we find a best-fit functional form $dN/d\ln m \propto m^{-\alpha}$, where $\alpha = 1.72 \pm 0.05$ or $1.75 \pm 0.05$ for two different theoretical isochrone models, respectively. Furthermore, we explore the possibility of systematic errors in the determination of the power-law index due to the depth of the observations. When we combine the observational data with artificially generated data from the lognormal Chabrier IMF for masses below $0.90\, M_{\odot}$, the best fit MLP is a hybrid function but with a steeper asymptotic slope i.e., $\alpha = 2.04 \pm 0.07$. This illustrates the systematic uncertainties in commonly used MF parameters that can depend on the range of data that is fitted.

This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfv\'enic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17{\deg}. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10{\mu}G to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5.

The Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission aiming to search for exoplanets that transit bright stars. The high-quality photometric data of TESS are excellent for the asteroseismic study of solar-like stars. In this work, we present an asteroseismic analysis of the red-giant star HD~222076 hosting a long-period (2.4 yr) giant planet discovered through radial velocities. Solar-like oscillations of HD~222076 are detected around $203 \, \mu$Hz by TESS for the first time. Asteroseismic modeling, using global asteroseismic parameters as input, yields a determination of the stellar mass ($M_\star = 1.12 \pm 0.12\, M_\odot$), radius ($R_\star = 4.34 \pm 0.21\,R_\odot$), and age ($7.4 \pm 2.7\,$Gyr), with precisions greatly improved from previous studies. The period spacing of the dipolar mixed modes extracted from the observed power spectrum reveals that the star is on the red-giant branch burning hydrogen in a shell surrounding the core. We find that the planet will not escape the tidal pull of the star and be engulfed into it within about $800\,$Myr, before the tip of the red-giant branch is reached.

We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framework to demonstrate that the interplay between a mass-dependent and a time-dependent step-like dropout rate from accretion leads to IMFs that exhibit multiple power laws for an exponential mass growth. While the mass-dependent accretion and its dropout is intrinsic to each star, the time-dependent dropout might be tied to a specific history such as the rapid consumption of nebular material by nearby stars or the sweeping away of some material by shock waves. The time-dependent dropout folded into the mass-dependent process of star formation is shown to have a significant influence on the IMFs.

The accretion phase of star formation is investigated in magnetically-dominated clouds that have an initial subcritical mass-to-flux ratio. We employ nonideal magnetohydrodynamic simulations that include ambipolar diffusion and ohmic dissipation. During the early prestellar phase the mass-to-flux ratio rises toward the critical value for collapse, and during this time the angular momentum of the cloud core is reduced significantly by magnetic braking. Once a protostar is formed in the core, the accretion phase is characterized by the presence of a small amount of angular momentum but a large amount of magnetic flux in the near-protostellar environment. The low angular momentum leads to a very small (or even nonexistent) disk and weak outflow, while the large magnetic flux can lead to an interchange instability that rapidly removes flux from the central region. The effective magnetic braking in the early collapse phase can even lead to a counter-rotating disk and outflow, in which the rotation direction of the disk and outflow is opposite to that of the infalling envelope. The solutions with a counter-rotating disk, tiny disk, or nonexistent disk (direct collapse) are unique outcomes that are realized in collapse from magnetically-dominated clouds with an initial subcritical mass-to-flux ratio.

Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot' measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraints from pressure modes to infer the CBM properties of Kepler and CoRoT main-sequence and subgiant oscillators, as well as in two radial velocity targets (Procyon A and $\alpha$ Cen A). Collectively these targets allow us to identify how measurement precision, stellar spectral type, and overshoot implementation impact the asteroseismic solution. With these new measures we find that the `effective overshoot' for most stars is in line with physical expectations and calibrations from binaries and clusters. However, two F-stars in the CoRoT field (HD 49933 and HD 181906) still necessitate high overshoot in the models. Due to short mode lifetimes, mode identification can be difficult in these stars. We demonstrate that an incongruence between the radial and non-radial modes drives the asteroseismic solution to extreme structures with highly-efficient CBM as an inevitable outcome. Understanding the cause of seemingly anomalous physics for such stars is vital for inferring accurate stellar parameters from TESS data with comparable timeseries length.

The search for twins of the Sun and Earth relies on accurate characterization of stellar and exoplanetary parameters: i.e., ages, masses, and radii. In the modern era of asteroseismology, parameters of solar-like stars are derived by fitting theoretical models to observational data, which include measurements of their oscillation frequencies, metallicity [Fe/H], and effective temperature Teff. Combining this information with transit data furthermore yields the corresponding parameters for their exoplanets. While [Fe/H] and Teff are commonly stated to a precision of ~0.1 dex and ~100 K, the impact of errors in their measurement has not been studied in practice within the context of the parameters derived from them. Here we use the Stellar Parameters in an Instant (SPI) pipeline to estimate the parameters of nearly 100 stars observed by Kepler and Gaia, many of which are confirmed planet hosts. We adjust the reported spectroscopic measurements of these stars by introducing faux systematic errors and artificially increasing the reported uncertainties, and quantify the differences in the resulting parameters. We find that a systematic error of 0.1 dex in [Fe/H] translates to differences of only 4%, 2%, and 1% on average in the resulting stellar ages, masses, and radii, which are well within their uncertainties (~11%, 3.5%, 1.4%) as derived by SPI. We also find that increasing the uncertainty of [Fe/H] measurements by 0.1 dex increases the uncertainties by only 0.01 Gyr, 0.02 M_sun, and 0.01 R_sun, which are again well below their reported uncertainties (0.5 Gyr, 0.04 M_sun, 0.02 R_sun). The results for Teff at 100 K are similar. Stellar parameters from SPI are unchanged within uncertainties by errors of up to 0.14 dex or 175 K, and are even more robust to errors in Teff than the seismic scaling relations. Consequently, the parameters for their exoplanets are robust as well.

Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate its potential for investigating the Galaxy by carrying out the first asteroseismic ensemble study of red giant stars observed by TESS. We use a sample of 25 stars for which we measure their global asteroseimic observables and estimate their fundamental stellar properties, such as radius, mass, and age. Significant improvements are seen in the uncertainties of our estimates when combining seismic observables from TESS with astrometric measurements from the Gaia mission compared to when the seismology and astrometry are applied separately. Specifically, when combined we show that stellar radii can be determined to a precision of a few percent, masses to 5-10% and ages to the 20% level. This is comparable to the precision typically obtained using end-of-mission Kepler data

The frequency, $\nu_{\rm max}$, at which the envelope of pulsation power peaks for solar-like oscillators is an important quantity in asteroseismology. We measure $\nu_{\rm max}$ for the Sun using 25 years of Sun-as-a-Star Doppler velocity observations with the Birmingham Solar-Oscillations Network (BiSON), by fitting a simple model to binned power spectra of the data. We also apply the fit to Sun-as-a-Star Doppler velocity data from GONG and GOLF, and photometry data from VIRGO/SPM on the ESA/NASA SOHO spacecraft. We discover a weak but nevertheless significant positive correlation of the solar $\nu_{\rm max}$ with solar activity. The uncovered shift between low and high activity, of $\simeq 25\,\rm \mu Hz$, translates to an uncertainty of 0.8 per cent in radius and 2.4 per cent in mass, based on direct use of asteroseismic scaling relations calibrated to the Sun. The mean $\nu_{\rm max}$ in the different datasets is also clearly offset in frequency. Our results flag the need for caution when using $\nu_{\rm max}$ in asteroseismology.

We report the discovery of a warm sub-Saturn, TOI-257b (HD 19916b), based on data from NASA's Transiting Exoplanet Survey Satellite (TESS). The transit signal was detected by TESS and confirmed to be of planetary origin based on radial-velocity observations with the Minerva-Australis telescope array. An analysis of the TESS photometry, the Minerva-Australis, FEROS, and HARPS radial velocities, and the asteroseismic data of the stellar oscillations reveals that TOI-257b has a mass of $M_P=0.134^{+0.023}_{-0.022}$$\rm{M_J}$ ($42.6^{+7.3}_{-7.0}$$\rm{M_{\oplus}}$), a radius of $R_P=0.626^{+0.013}_{-0.012}$$\rm{R_J}$ ($7.02^{+0.15}_{-0.13}$$\rm{R_{\oplus}}$), and an orbit with eccentricity $0.242^{+0.040}_{-0.065}$ and period $18.38827\pm0.00072$$\rm{days}$. TOI-257b orbits a bright ($\mathrm{V}=7.570$mag) somewhat evolved late F-type star with $M_*=1.390\pm0.046$$\rm{M_{\odot}}$, $R_*=1.888\pm0.033$$\rm{R_{\odot}}$, $T_{\rm eff}=6075\pm90$$\rm{K}$, and $v\sin{i}=11.3\pm0.5$km/s. Additionally, we statistically validate a second non-transiting sub-Saturn mass planet on a $\sim71$day orbit using the radial velocity data. This system joins the ranks of a small number of exoplanet host stars that have been characterized with asteroseismology. Warm sub-Saturns are rare in the known sample of exoplanets, and thus the discovery of TOI-257b is important in the context of future work studying the formation and migration history of similar planetary systems.

The evolution of low-mass stars into red giants is still poorly understood. During this evolution the core of the star contracts and, simultaneously, the envelope expands -- a process known as the `mirror'. Additionally, there is a short phase where the trend for increasing luminosity is reversed. This is known as the red-giant-branch bump. We explore the underlying physical reasons for these two phenomena by considering the specific entropy distribution in the star and its temporal changes. We find that between the luminosity maximum and luminosity minimum of the bump there is no mirror present and the star is fully contracting. The contraction is halted and the star regains its mirror when the hydrogen-burning shell reaches the mean molecular weight discontinuity. This marks the luminosity minimum of the bump.

We revisit the analysis of the bright multiplanet system K2-93, discovered with data taken by the K2 mission. This system contains five identified planets ranging in size from sub-Neptune to Jupiter size. The K2 data available at the discovery of the system only showed single transits for the three outer planets, which allowed weak constraints to be put on their periods. As these planets are interesting candidates for future atmospheric studies, a better characterization of the host star and tighter constraints on their orbital periods are essential. Using new data from the K2 mission taken after the discovery of the system, we perform an asteroseismic characterization of the host star. We are able to place strong constraints on the stellar parameters and obtain a value for the stellar mass of $1.22^{+0.03}_{-0.02}\, \rm M_{\odot}$, a stellar radius of $1.30\pm 0.01\, \rm R_{\odot}$, and an age of $2.07^{+0.36}_{-0.27}$ Gyr. Put together with the additional transits identified for two of the three outer planets, we constrain the orbital periods of the outer planets and provide updated estimates for the stellar reflex velocities induced by the planets.

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

An increasingly popular method to determine stellar ages of red-giant stars for the purpose of Galactic archaeology is asteroseismic grid-based modelling (GBM). In asteroseismic GBM of red-giant stars with solar-like oscillations the large frequency separation ($\Delta\nu$) and the frequency of maximum oscillation power ($\nu_{\rm max}$) are commonly used asteroseismic observables, in addition to the usual spectroscopic parameters effective temperature ($T_{\rm eff}$) and metallicity ([Fe/H]). The precision with which $\Delta\nu$ and $\nu_{\rm max}$ can be determined largely depends on the length of the timeseries data (assuming the stars are bright enough that the oscillations can be detected). The question that we aim to answer here is: with what precision should [Fe/H] and $T_{\rm eff}$ be obtained to derive stellar ages of red-giant stars through asteroseismic GBM given the precision of $\Delta\nu$ and $\nu_{\rm max}$ that we can expect from the TESS data? In most cases presented here the value of $\sigma_{\rm age}$ shows only weak, or no, dependence on the uncertainties and biases in [Fe/H] and $T_{\rm eff}$ ingested in this study. Only in case a large ($>100$~K) $T_{\rm eff}$ biases is accompanied by a small (50 K) uncertainties in $T_{\rm eff}$ on the red-giant branch, $\Delta_{\rm age}$ is of the order of or larger than $\sigma_{\rm age}$ and in these cases we expect an impact of this bias on the resulting age estimate.

The goal of stellar evolution theory is to predict the structure of stars throughout their lifetimes. Usually, these predictions can be assessed only indirectly, for example by comparing predicted and observed effective temperatures and luminosities. Thanks now to asteroseismology, which can reveal the internal structure of stars, it becomes possible to compare the predictions from stellar evolution theory to actual stellar structures. In this work, we present an inverse analysis of the oscillation data from the solar-type star KIC 6225718, which was observed by the \emph{Kepler} space observatory during its nominal mission. As its mass is about 20% greater than solar, this star is predicted to transport energy by convection in its nuclear-burning core. We find significant differences between the predicted and actual structure of the star in the radiative interior near to the convective core. In particular, the predicted sound speed is higher than observed in the deep interior of the star, and too low at a fractional radius of 0.25 and beyond. The cause of these discrepancies is unknown, and is not remedied by known physics in the form of convective overshooting or elemental diffusion.

The strongest non-Gaussianity in single-scalar potential models of inflation is associated with features in the power spectrum. We stress the importance of accurately modelling the expected signal in order for the standard estimator to minimize contamination by random noise. We present explicit formulae which improve on the approximation introduced by Adshead, Hu, Dvorkin and Peiris. We also compare with a simple, analytic model of the first feature, and quantify our results using the correlators of Hung, Fergusson and Shellard.

The Transiting Exoplanet Survey Satellite (TESS) is performing a near all-sky survey for planets that transit bright stars. In addition, its excellent photometric precision enables asteroseismology of solar-type and red-giant stars, which exhibit convection-driven, solar-like oscillations. Simulations predict that TESS will detect solar-like oscillations in nearly 100 stars already known to host planets. In this paper, we present an asteroseismic analysis of the known red-giant host stars HD 212771 and HD 203949, both systems having a long-period planet detected through radial velocities. These are the first detections of oscillations in previously known exoplanet-host stars by TESS, further showcasing the mission's potential to conduct asteroseismology of red-giant stars. We estimate the fundamental properties of both stars through a grid-based modeling approach that uses global asteroseismic parameters as input. We discuss the evolutionary state of HD 203949 in depth and note the large discrepancy between its asteroseismic mass ($M_\ast = 1.23 \pm 0.15\,{\rm M}_\odot$ if on the red-giant branch or $M_\ast = 1.00 \pm 0.16\,{\rm M}_\odot$ if in the clump) and the mass quoted in the discovery paper ($M_\ast = 2.1 \pm 0.1\,{\rm M}_\odot$), implying a change $>30\,\%$ in the planet's mass. Assuming HD 203949 to be in the clump, we investigate the planet's past orbital evolution and discuss how it could have avoided engulfment at the tip of the red-giant branch. Finally, HD 212771 was observed by K2 during its Campaign 3, thus allowing for a preliminary comparison of the asteroseismic performances of TESS and K2. We estimate the ratio of the observed oscillation amplitudes for this star to be $A_{\rm max}^{\rm TESS}/A_{\rm max}^{\rm K2} = 0.75 \pm 0.14$, consistent with the expected ratio of $\sim0.85$ due to the redder bandpass of TESS.

In the asymptotic parameterisation of mode frequencies, the phase function $\epsilon(\nu)$ completely specifies the detailed structure of the frequency eigenvalues. In practice, however, this function of frequency is reduced to a single scalar $\epsilon$, defined, particularly by observers, as the intercept of a least-squares fit to the frequencies against radial order, or via the central value of this function. The procedure by which this is done is not unique. We derive a few simple expressions relating various observational estimators of $\epsilon$ for radial modes to each other, and to the underlying theoretical object. In particular we demonstrate that a ``reduced' functional parameterisation is both insensitive to mis-estimations of $\Delta\nu$, and easy to evaluate locally in terms of both observational and theoretical quantities. It has been shown previously that such a local definition of $\epsilon$ can distinguish between stars on the ascending part of the red giant branch and those in the red clump. We find that this sensitivity to evolutionary stage arises from differences in the local frequency derivative of the underlying phase function, a consequence of differences in internal structure. By constructing an HR-like diagram out of purely seismic observables, we provide a unified view of the \textit{Kepler} asteroseismic sample, as well as the initial results from \textit{TESS}. We investigate how various astrophysical quantities and modelling parameters affect the morphology of isochrones on this seismic diagram. We also show that $\epsilon$ can be used as an independent input when deriving stellar parameters from global asteroseismic quantities.

We explore the sensitivity of the frequencies of low-degree solar p-modes to near-surface magnetic flux on different spatial scales and strengths, specifically to active regions with strong magnetic fields and ephemeral regions with weak magnetic fields. We also use model reconstructions from the literature to calculate average frequency offsets back to the end of the Maunder minimum. We find that the p-mode frequencies are at least three times less sensitive (at 95% confidence) to the ephemeral-region field than they are to the active-region field. Frequency shifts between activity cycle minima and maxima are controlled predominantly by the change of active region flux. Frequency shifts at cycle minima (with respect to a magnetically quiet Sun) are determined largely by the ephemeral flux, and are estimated to have been $0.1\,\rm \mu Hz$ or less over the last few minima. We conclude that at epochs of cycle minimum, frequency shifts due to near-surface magnetic activity are negligible compared to the offsets between observed and model frequencies that arise from inaccurate modelling of the near-surface layers (the so-called surface term). The implication is that this will be the case for other Sun-like stars with similar activity, which has implications for asteroseismic modelling of stars.

We use helioseismic data from ground and space-based instruments to analyze how solar rotation has changed since the beginning of solar Cycle 23 with emphasis on studying the differences between Cycles 23 and 24. We find that the nature of solar rotation is indeed different for the two cycles. While the changes in the latitudinally independent component follows solar-cycle indices, some of the other components have a more complicated behavior. There is a substantial change in the behavior of the solar zonal flows and their spatial gradients too. While the zonal flows are in general weaker in Cycle 24 than those in Cycle 23, there are clear signs of the emergence of Cycle 25. We have also investigated the properties of the solar tachocline, in particular, its position, width, and the change (or jump) in the rotation rate across it. We find significant temporal variation in the change of the rotation rate across the tachocline. We also find that the changes in solar Cycle 24 were very different from those of Cycle 23. We do not find any statistically significant change in the position or the width of the tachocline.

We introduce a theory for the development of a transitional column density $\Sigma_{\rm TP}$ between the lognormal and the power-law forms of the probability distribution function (PDF) in a molecular cloud. Our turbulent magnetohydrodynamic simulations show that the value of $\Sigma_{\rm TP}$ increases as the strength of both the initial magnetic field and turbulence increases. We develop an analytic expression for $\Sigma_{\rm TP}$ based on the interplay of turbulence, a (strong) magnetic field, and gravity. The transition value $\Sigma_{\rm TP}$ scales with $\mathcal{M}^2_{\rm 0}$, the square of the initial sonic Mach number, and $\beta_{0}$, the initial ratio of gas pressure to magnetic pressure. We fit the variation of $\Sigma_{\rm TP}$ among different model clouds as a function of $\mathcal{M}^2_{\rm 0} \beta_{0}$, or equivalently the square of the initial Alfv\'enic Mach number $\mathcal{M}^2_{\rm A0}$. This implies that the transition value $\Sigma_{\rm TP}$ is an imprint of cloud initial conditions and is set by turbulent compression of a magnetic cloud. Physically, the value of $\Sigma_{\rm TP}$ denotes the boundary above which the mass-to-flux ratio becomes supercritical and gravity drives the evolution.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

One of the ideas to explain the existence of supermassive black holes (SMBH) that are in place by z~7 is that there was an earlier phase of very rapid accretion onto direct collapse black holes (DCBH) that started their lives with masses ~ 10^4-10^5 M_solar. Working in this scenario, we show that the mass function of SMBH after such a limited time period with growing formation rate paired with super-Eddington accretion can be described as a broken power-law with two characteristic features. There is a power-law at intermediate masses whose index is the dimensionless ratio {\alpha} = {\lambda}/{\gamma}, where {\lambda} is the growth rate of the number of DCBH during their formation era, and {\gamma} is the growth rate of DCBH masses by super-Eddington accretion during the DCBH growth era. A second feature is a break in the power law profile at high masses, above which the mass function declines rapidly. The location of the break is related to the dimensionless number \b{eta} = {\gamma} T, where T is the duration of the period of DCBH growth. If the SMBH continue to grow at later times at an Eddington-limited accretion rate, then the observed quasar luminosity function can be directly related to the tapered power-law function derived in this paper.

Asteroseismology of solar-like oscillators often relies on the comparisons between stellar models and stellar observations in order to determine the properties of stars. The values of the global seismic parameters, $\nu_\mathrm{max}$ (the frequency where the smoothed amplitude of the oscillations peak) and $\Delta \nu$ (the large frequency separation), are frequently used in grid-based modeling searches. However, the methods by which $\Delta \nu$ is calculated from observed data and how $\Delta \nu$ is calculated from stellar models are not the same. Typically for observed stars, especially for those with low signal-to-noise data, $\Delta \nu$ is calculated by taking the power spectrum of a power spectrum, or with autocorrelation techniques. However, for stellar models, the actual individual mode frequencies are calculated and the average spacing between them directly determined. In this work we try to determine the best way to combine model frequencies in order to obtain $\Delta \nu$ that can be compared with observations. For this we use stars with high signal-to-noise observations from Kepler as well as simulated TESS data of Ball et al. (2018). We find that when determining $\Delta \nu$ from individual mode frequencies the best method is to use the $\ell=0$ modes with either no weighting or with a Gaussian weighting around $\nu_\mathrm{max}$.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

Starting from a prestellar core with a size of $1.2\times10^4$ AU, we calculate the evolution of a gravitationally collapsing core until $\sim2000$ yr after protostar formation using a three-dimensional resistive magnetohydrodynamic simulation, in which the protostar is resolved with a spatial resolution of $5.6\times10^{-3}$ AU. Following protostar formation, a rotationally supported disk is formed. Although the disk size is as small as $\sim2-4$ AU, it remains present until the end of the simulation. Since the magnetic field dissipates and the angular momentum is then not effectively transferred by magnetic effects, the disk surface density gradually increases and spiral arms develop due to gravitational instability. The disk angular momentum is then transferred mainly by gravitational torques, which induce an episodic mass accretion onto the central protostar. The episodic accretion causes a highly time-variable mass ejection (the high-velocity jet) near the disk inner edge, where the magnetic field is well coupled with the neutral gas. As the mass of the central protostar increases, the jet velocity gradually increases and exceeds $\sim100$ km s$^{-1}$. The jet opening angle widens with time at its base, while the jet keeps a very good collimation on the large scale. In addition, a low-velocity outflow is driven from the disk outer edge. A cavity-like structure, a bow shock and several knots, all of which are usually observed in star-forming regions, are produced in the outflowing region.

We present the discovery of TOI-197.01, the first transiting planet identified by the Transiting Exoplanet Survey Satellite (TESS) for which asteroseismology of the host star is possible. TOI-197 (HIP116158) is a bright (V=8.2 mag), spectroscopically classified subgiant which oscillates with an average frequency of about 430 muHz and displays a clear signature of mixed modes. The oscillation amplitude confirms that the redder TESS bandpass compared to Kepler has a small effect on the oscillations, supporting the expected yield of thousands of solar-like oscillators with TESS 2-minute cadence observations. Asteroseismic modeling yields a robust determination of the host star radius (2.943+/-0.064 Rsun), mass (1.212 +/- 0.074 Msun) and age (4.9+/-1.1 Gyr), and demonstrates that it has just started ascending the red-giant branch. Combining asteroseismology with transit modeling and radial-velocity observations, we show that the planet is a "hot Saturn" (9.17+/-0.33 Rearth) with an orbital period of ~14.3 days, irradiance of 343+/-24 Fearth, moderate mass (60.5 +/- 5.7 Mearth) and density (0.431+/-0.062 gcc). The properties of TOI-197.01 show that the host-star metallicity - planet mass correlation found in sub-Saturns (4-8 Rearth) does not extend to larger radii, indicating that planets in the transition between sub-Saturns and Jupiters follow a relatively narrow range of densities. With a density measured to ~15%, TOI-197.01 is one of the best characterized Saturn-sized planets to date, augmenting the small number of known transiting planets around evolved stars and demonstrating the power of TESS to characterize exoplanets and their host stars using asteroseismology.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

The helium abundance of stars is a strong driver of evolutionary timescales, however it is difficult to measure in cool stars. We conduct an asteroseismic analysis of NGC 6791, an old, metal rich open cluster that previous studies have indicated also has a high helium abundance. The cluster was observed by Kepler and has unprecedented lightcurves for many of the red giant branch stars in the cluster. Previous asteroseismic studies with Kepler data have constrained the age through grid based modeling of the global asteroseismic parameters ($\Delta\nu$ and $\nu_{\mathrm{max}}$). However, with the precision of Kepler data, it is possible to do detailed asteroseismology of individual mode frequencies to better constrain the stellar parameters, something that has not been done for these cluster stars as yet. In this work, we use the observed mode frequencies in 27 hydrogen shell burning red giants to better constrain initial helium abundance ($Y_0$) and age of the cluster. The distributions of helium abundance and age for each individual red giant are combined to create a final probability distribution for age and helium abundance of the entire cluster. We find a helium abundance of $Y_0=0.297\pm0.003$ and a corresponding age of $8.2\pm0.3$ Gyr.

We present the target list of solar-type stars to be observed in short-cadence (2-min) for asteroseismology by the NASA Transiting Exoplanet Survey Satellite (TESS) during its 2-year nominal survey mission. The solar-like Asteroseismic Target List (ATL) is comprised of bright, cool main-sequence and subgiant stars and forms part of the larger target list of the TESS Asteroseismic Science Consortium (TASC). The ATL uses Gaia DR2 and the Extended Hipparcos Compilation (XHIP) to derive fundamental stellar properties, calculate detection probabilities and produce a rank-ordered target list. We provide a detailed description of how the ATL was produced and calculate expected yields for solar-like oscillators based on the nominal photometric performance by TESS. We also provide publicly available source code which can be used to reproduce the ATL, thereby enabling comparisons of asteroseismic results from TESS with predictions from synthetic stellar populations.

We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in the L1688 region of Ophiuchus. The CFS model predicts the profile of the relative field fluctuation, which is related to the observable dispersion in direction of the polarization vectors. Within the context of our model we find that all the cores have a transcritical mass-to-flux ratio.

The structural stratification of a solar-type main sequence star primarily depends on its mass and chemical composition. The surface heavy element abundances of the solar-type stars are reasonably well determined using conventional spectroscopy, however the second most abundant element helium is not. This is due to the fact that the envelope temperature of such stars is not high enough to excite helium. Since the helium abundance of a star affects its structure and subsequent evolution, the uncertainty in the helium abundance of a star makes estimates of its global properties (mass, radius, age etc.) uncertain as well. The detections of the signatures of the acoustic glitches from the precisely measured stellar oscillation frequencies provide an indirect way to estimate the envelope helium content. We use the glitch signature caused by the ionization of helium to determine the envelope helium abundance of 38 stars in the {\it Kepler} seismic LEGACY sample. Our results confirm that atomic diffusion does indeed take place in solar-type stars. We use the measured surface abundances in combination with the settling predicted by the stellar models to determine the initial abundances. The initial helium and metal mass fractions have subsequently been used to get the preliminary estimates of the primordial helium abundance, $Y_p = 0.244\pm0.019$, and the galactic enrichment ratio, $\Delta Y / \Delta Z = 1.226\pm0.849$. Although the current estimates have large errorbars due to the limited sample size, this method holds great promises to determine these parameters precisely in the era of upcoming space missions.

Asteroseismic large frequency separations possess great diagnostic value. However, their expressions as scaling relations are predicated on homology arguments which may not hold in general, resulting in mass- and temperature-dependent deviations. The first-order asymptotic expressions, which should in principle account for this structural evolution, also deviate more from fitted frequency-separation estimates than do the simple scaling relations, and exhibit qualitatively different behavior. We present a modified asymptotic estimator, and show that these discrepancies can be accounted for by the evolution of the acoustic turning points of the asteroseismic mode cavity, which is typically neglected in first-order asymptotic analysis. This permits us to use a single expression to accurately estimate the large frequency separations of main-sequence, ascending red giant branch, and red clump stellar models, except at transition points between two asymptotic regimes during the subgiant phase of evolution, where the WKB approach fails. The existence of such transition points provides theoretical justification for separately calibrated scaling relations for stars in different evolutionary stages.

Main sequence, solar-like stars (M < 1.5 Msun) have outer convective envelopes that are sufficiently thick to affect significantly their overall structure. The radii of these stars, in particular, are sensitive to the details of inefficient, super-adiabatic convection occurring in their outermost layers. The standard treatment of convection in stellar evolution models, based on the Mixing-Length Theory (MLT), provides only a very approximate description of convection in the super-adiabatic regime. Moreover, it contains a free parameter, alpha_MLT, whose standard calibration is based on the Sun, and is routinely applied to other stars ignoring the differences in their global parameters (e.g., effective temperature, gravity, chemical composition) and previous evolutionary history. In this paper, we present a calibration of alpha_MLT based on three-dimensional radiation-hydrodynamics (3D RHD) simulations of convection. The value of alpha_MLT is adjusted to match the specific entropy in the deep, adiabatic layers of the convective envelope to the corresponding value obtained from the 3D RHD simulations, as a function of the position of the star in the (log g, log T_eff) plane and its chemical composition. We have constructed a model of the present-day Sun using such entropy-based calibration. We find that its past luminosity evolution is not affected by the entropy calibration. The predicted solar radius, however, exceeds that of the standard model during the past several billion years, resulting in a lower surface temperature. This illustrative calculation also demonstrates the viability of the entropy approach for calibrating the radii of other late-type stars.

Oscillation frequencies of even the best stellar models differ from those of the stars they represent, and the difference is predominantly a function of frequency. This difference is caused by limitations of modeling the near-surface layers of a star. This frequency-dependent frequency error, usually referred to as the "surface term" can result in erroneous interpretation of the oscillation frequencies unless treated properly. Several techniques have been developed to minimize the effect of the surface term; it is either subtracted out, or frequency combinations insensitive to the surface terms are used, or the asteroseismic phase $\epsilon$ is used to determine a match between models and stars. In this paper we show that no matter what technique is used to account for the surface term, as long as the physics of the models is the same, the global parameters of a star --- mass, radius and, age --- obtained from frequency analyses are robust. This implies that even though we cannot model the internal structure of stars perfectly, we can have confidence in all results that use stellar global properties obtained through the analysis of stellar oscillation frequencies.

A model of magnetic field structure is presented to help test the prevalence of flux freezing in star-forming clouds of various shapes, orientations, and degrees of central concentration, and to estimate their magnetic field strength. The model is based on weak-field flux freezing in centrally condensed Plummer spheres and spheroids of oblate and prolate shape. For a spheroid of given density contrast, aspect ratio, and inclination, the model estimates the local field strength and direction, and the global field pattern of hourglass shape. Comparisons with a polarization simulation indicate typical angle agreement within 1 - 10 degrees. Scalable analytic expressions are given to match observed polarization patterns, and to provide inputs to radiative transfer codes for more accurate predictions. The model may apply to polarization observations of dense cores, elongated filamentary clouds, and magnetized circumstellar disks.

We present an asteroseismic analysis of the Kepler light curve of KIC 8263801, a red-giant star in the open cluster NGC 6866 that has previously been reported to be a helium-burning red-clump star. We extracted the frequencies of the radial and quadrupole modes from its frequency power spectrum and determined its properties using a grid of evolutionary models constructed with MESA. The oscillation frequencies were calculated using the GYRE code and the surface term was corrected using the Ball & Gizon(2014) prescription. We find that the star has a mass of $M/M_{\odot} = 1.793\pm 0.072$, age $t=1.48\pm 0.21$ Gyr and radius $R/R_{\odot} = 10.53\pm 0.28$. By analyzing the internal structure of the best-fitting model, we infer the evolutionary status of the star KIC 8263801 as being on the ascending part of the red giant branch, and not on the red clump. This result is verified using a purely asteroseismic diagnostic, the $\epsilon_{c}-\Delta\nu_{c}$ diagram which can distinguish red giant branch stars from red clump stars. Finally, by comparing its age with NGC 6866 ($t=0.65 \pm 0.1$ Gyr) we conclude that KIC 8263801 is not a member of this open cluster.

Low-degree solar p-mode observations from the long-lived Birmingham Solar Oscillations Network (BiSON) stretch back further than any other single helioseismic data set. Results from BiSON have suggested that the response of the mode frequency to solar activity levels may be different in different cycles. In order to check whether such changes can also be seen at higher degrees, we compare the response of medium-degree solar p-modes to activity levels across three solar cycles using data from Big Bear Solar Observatory (BBSO), Global Oscillation Network Group (GONG), Michelson Doppler Imager (MDI) and Helioseismic and Magnetic Imager (HMI), by examining the shifts in the mode frequencies and their sensitivity to solar activity levels. We compare these shifts and sensitivities with those from radial modes from BiSON. We find that the medium-degree data show small but significant systematic differences between the cycles, with solar cycle 24 showing a frequency shift about 10 per cent larger than cycle 23 for the same change in activity as determined by the 10.7 cm radio flux. This may support the idea that there have been changes in the magnetic properties of the shallow subsurface layers of the Sun that have the strongest influence on the frequency shifts.

We investigate how a single generation galactic mass function (SGMF) depends on the existence of variations in the initial stellar mass functions (IMF) of stellar clusters. We show that cluster-to-cluster variations of the IMF lead to a multicomponent SGMF where each component in a given mass range can be described by a distinct power-law function. We also show that a dispersion of $\approx 0.3$ M$_{\odot}$ in the characteristic mass of the IMF, as observed for young Galactic clusters, leads to a low mass slope of the SGMF that matches the observed Galactic stellar mass function even when the IMFs in the low mass end of individual clusters are much steeper.

We introduce a new dual power law (DPL) probability distribution function for the mass distribution of stellar and substellar objects at birth, otherwise known as the initial mass function (IMF). The model contains both deterministic and stochastic elements, and provides a unified framework within which to view the formation of brown dwarfs and stars resulting from an accretion process that starts from extremely low mass seeds. It does not depend upon a top down scenario of collapsing (Jeans) masses or an initial lognormal or otherwise IMF-like distribution of seed masses. Like the modified lognormal power law (MLP) distribution, the DPL distribution has a power law at the high mass end, as a result of exponential growth of mass coupled with equally likely stopping of accretion at any time interval. Unlike the MLP, a power law decay also appears at the low mass end of the IMF. This feature is closely connected to the accretion stopping probability rising from an initially low value up to a high value. This might be associated with physical effects of ejections sometimes (i.e., rarely) stopping accretion at early times followed by outflow driven accretion stopping at later times, with the transition happening at a critical time (therefore mass). Comparing the DPL to empirical data, the critical mass is close to the substellar mass limit, suggesting that the onset of nuclear fusion plays an important role in the subsequent accretion history of a young stellar object.

We present a catalog of stellar properties for a large sample of 6676 evolved stars with APOGEE spectroscopic parameters and \textit{Kepler} asteroseismic data analyzed using five independent techniques. Our data includes evolutionary state, surface gravity, mean density, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them. We employ a new empirical approach for combining asteroseismic measurements from different methods, calibrating the inferred stellar parameters, and estimating uncertainties. With high statistical significance, we find that asteroseismic parameters inferred from the different pipelines have systematic offsets that are not removed by accounting for differences in their solar reference values. We include theoretically motivated corrections to the large frequency spacing ($\Delta \nu$) scaling relation, and we calibrate the zero point of the frequency of maximum power ($\nu_{\rm max}$) relation to be consistent with masses and radii for members of star clusters. For most targets, the parameters returned by different pipelines are in much better agreement than would be expected from the pipeline-predicted random errors, but 22\% of them had at least one method not return a result and a much larger measurement dispersion. This supports the usage of multiple analysis techniques for asteroseismic stellar population studies. The measured dispersion in mass estimates for fundamental calibrators is consistent with our error model, which yields median random and systematic mass uncertainties for RGB stars of order 4\%. Median random and systematic mass uncertainties are at the 9\% and 8\% level respectively for RC stars.

Availability of reliable atomic and molecular opacity tables is essential in a wide variety of astronomical modeling: the solar and stellar interiors, stellar and planetary atmospheres, stellar evolution, pulsating stars, and protoplanetary disks, to name a few. With the advancement of powerful research techniques such as helio-seismology and asteroseismology, solar neutrino-flux measurements, exoplanet survey satellites, three-dimensional hydrodynamic atmospheric simulations (including non-LTE and granulation effects), high-performance computing of atomic and molecular data, and innovative plasma experiments the accuracy and completeness of opacity tables is being taken to an unprecedented level. The goal of the second Workshop on Astrophysical Opacities was to gather opacity data producers and consumers from both the atomic and molecular sectors to contribute to solving outstanding problems and to develop more effective and integrated interfaces. In this review we attempt to summa- rize the discussion at the workshop and propose future directions for opacity research.

HD 176465 is a binary system whose both components are solar-like pulsators and whose oscillation frequencies were observed by the \kepler\ mission. In this paper we have modeled the asteroseismic and spectroscopic data of the stars, and have determined their convection-zone helium abundances using the signatures left by the He{\sc ii} ionization zone on the mode frequencies. As expected we find that the components of the binary have the same age within uncertainties ($3.087 \pm 0.580$ Gyr and $3.569 \pm 0.912$ Gyr); they also have the same initial helium abundance (Y$_{\mathrm{init}}$=0.253 $\pm$ 0.006 and 0.254 $\pm$ 0.008). Their current metallicity ([Fe/H]=$-0.275 \pm 0.04$ and $-0.285 \pm 0.04$) is also the same within errors. Fits to the signature of the He{\sc ii} acoustic glitch yields current helium abundances of $Y_{\rm A} = 0.224 \pm 0.006$ and $Y_{\rm B} = 0.233 \pm 0.008$ for the two components. Analyzing the complete ensemble of models generated for this investigation we find that both the amplitude and acoustic depth of the glitch signature arising from the second helium ionization zone and the base of the convection zone (CZ) are functions of mass. We show that the acoustic depths of these glitches are positively correlated with each other. The analysis can help us to detect the internal structure and constrain the chemical compositions.

Stellar models typically use the mixing length approximation as a way to implement convection in a simplified manner. While conventionally the value of the mixing length parameter, $\alpha$, used is the solar calibrated value, many studies have shown that other values of $\alpha$ are needed to properly model stars. This uncertainty in the value of the mixing length parameter is a major source of error in stellar models and isochrones. Using asteroseismic data, we determine the value of the mixing length parameter required to properly model a set of about 450 stars ranging in $\log g$, $T_{\mathrm{eff}}$, and $\mathrm{[Fe/H]}$. The relationship between the value of $\alpha$ required and the properties of the star is then investigated. For Eddington atmosphere, non-diffusion models, we find that the value of $\alpha$ can be approximated by a linear model, in the form of $\alpha/\alpha_{\odot}=5.426 -0.101 \log (g) -1.071 \log (T_{\mathrm{eff}}) + 0.437 (\mathrm{[Fe/H]})$. This process is repeated using a variety of model physics as well as compared to previous studies and results from 3D convective simulations.

EarthFinder is a Probe Mission concept selected for study by NASA for input to the 2020 astronomy decadal survey. This study is currently active and a final white paper report is due to NASA at the end of calendar 2018. We are tasked with evaluating the scientific rationale for obtaining precise radial velocity (PRV) measurements in space, which is a two-part inquiry: What can be gained from going to space? What can't be done form the ground? These two questions flow down to these specific tasks for our study - Identify the velocity limit, if any, introduced from micro- and macro-telluric absorption in the Earth's atmosphere; Evaluate the unique advantages that a space-based platform provides to emable the identification and mitigation of stellar acitivity for multi-planet signal recovery.

We present a method for measuring internal stellar structure based on asteroseismology that we call "inversions for agreement." The method accounts for imprecise estimates of stellar mass and radius as well as the relatively limited oscillation mode sets that are available for distant stars. By construction, the results of the method are independent of stellar models. We apply this method to measure the isothermal sound speeds in the cores of the solar-type stars 16 Cyg A and B using asteroseismic data obtained from Kepler observations. We compare the asteroseismic structure that we deduce against best-fitting evolutionary models and find that the sound speeds in the cores of these stars exceed those of the models.

(Abridged) We present the first APOKASC catalog of spectroscopic and asteroseismic data for 415 dwarfs and subgiants. Asteroseismic data have been obtained by Kepler in short cadence. The spectroscopic parameters are based on spectra taken as part of APOGEE and correspond to DR13 of SDSS. We analyze our data using two Teff scales, the spectroscopic values from DR13 and those derived from SDSS griz photometry. We use the differences in our results arising from these choices as a test of systematic Teff, and find that they can lead to significant differences in the derived stellar properties. Determinations of surface gravity ($\log{g}$), mean density ($\rho$), radius ($R$), mass ($M$), and age ($\tau$) for the whole sample have been carried out with stellar grid-based modeling. We have assessed random and systematic error sources in the spectroscopic and seismic data, as well as in the grid-based modeling determination of the stellar quantities in the catalog. We provide stellar properties for both Teff scales. The median combined (random and systematic) uncertainties are 2% (0.01 dex; $\log{g}$), 3.4% ($\rho$), 2.6% ($R$), 5.1% ($M$), and 19% ($\tau$) for the photometric Teff scale and 2% ($\log{g}$), 3.5% ($\rho$), 2.7% ($R$), 6.3% ($M$), and 23% ($\tau$) for the spectroscopic scale. Comparisons with stellar quantities in the catalog by Chaplin et al.(2014) highlight the importance of metallicity measurements for determining stellar parameters accurately. We compare our results with those from other sources, including stellar radii determined from TGAS parallaxes and asteroseismic analyses based on individual frequencies. We find a very good agreement in all cases. Comparisons give strong support to the determination of stellar quantities based on global seismology, a relevant result for future missions such as TESS and PLATO. Table 5 corrected (wrongly listed SDSS Teff before).

We review star formation in molecular clouds and describe why magnetic fields may be important and how they can influence filamentary structure and the column density probability distribution function (PDF). We also comment on the origin of the stellar and substellar initial mass function (IMF), which may require explanations beyond a simple Jeans length argument in turbulent molecular clouds. A mathematical model of the modified lognormal power-law (MLP) distribution function provides a framework within which to connect accretion processes with the IMF.

Simulations generally show that non-self-gravitating clouds have a lognormal column density ($\Sigma$) probability distribution function (PDF), while self-gravitating clouds with active star formation develop a distinct power-law tail at high column density. Although the growth of the power law can be attributed to gravitational contraction leading to the formation of condensed cores, it is often debated if an observed lognormal shape is a direct consequence of supersonic turbulence alone, or even if it is really observed in molecular clouds. In this paper we run three-dimensional magnetohydrodynamic simulations including ambipolar diffusion with different initial conditions to see the effect of strong magnetic fields and nonlinear initial velocity perturbations on the evolution of the column density PDFs. Our simulations show that column density PDFs of clouds with supercritical mass-to-flux ratio, with either linear perturbations or nonlinear turbulence, quickly develop a power-law tail such that $dN/d \log \Sigma \propto \Sigma^{-\alpha}$ with index $\alpha \simeq 2$. Interestingly, clouds with subcritical mass-to-flux ratio also proceed directly to a power-law PDF, but with a much steeper index $\alpha \simeq 4$. This is a result of gravitationally-driven ambipolar diffusion. However, for nonlinear perturbations with a turbulent spectrum ($v_{k}^{2} \propto k^{-4}$), the column density PDFs of subcritical clouds do retain a lognormal shape for a major part of the cloud evolution, and only develop a distinct power-law tail with index $\alpha \simeq 2$ at greater column density when supercritical pockets are formed.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

There is a widespread belief among inflationary cosmologists that a local observer cannot sense super-horizon gravitons. The argument goes that a local observer would subsume super-horizon gravitons into a redefinition of his coordinate system. We show that adopting this view for pure gravity on de Sitter background leads to time variation in the Hubble parameter measured by a local observer. It also leads to a violation of the gravitational field equation $R = 4 \Lambda$ because that equation is obeyed by the full metric, rather than the one which has been cleansed of super-horizon modes.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

The Kepler space mission has made it possible to measure the rotational splittings of mixed modes in red giants, thereby providing an unprecedented opportunity to probe the internal rotation of these stars. Asymmetries have been detected in the rotational multiplets of several red giants. This is unexpected since all the red giants whose rotation have been measured thus far are found to rotate slowly, and low rotation, in principle, produces symmetrical multiplets. Our aim here is to explain these asymmetries and find a way of exploiting them to probe the internal rotation of red giants. We show that in the cases where asymmetrical multiplets were detected, near-degeneracy effects are expected to occur, because of the combined effects of rotation and mode mixing. Such effects have not been taken into account so far. By using both perturbative and non-perturbative approaches, we show that near-degeneracy effects produce multiplet asymmetries that are very similar to the observations. We then propose and validate a method based on the perturbative approach to probe the internal rotation of red giants using multiplet asymmetries. We successfully apply our method to the asymmetrical $l=2$ multiplets of the Kepler young red giant KIC7341231 and obtain precise estimates of its mean rotation in the core and the envelope. The observed asymmetries are reproduced with a good statistical agreement, which confirms that near-degeneracy effects are very likely the cause of the detected multiplet asymmetries. We expect near-degeneracy effects to be important for $l=2$ mixed modes all along the red giant branch (RGB). For $l=1$ modes, these effects can be neglected only at the base of the RGB. They must therefore be taken into account when interpreting rotational splittings and as shown here, they can bring valuable information about the internal rotation of red giants.

We examine the frequency shifts in low-degree helioseismic modes from the Birmingham Solar-Oscillations Network (BiSON) covering the period from 1985 - 2016, and compare them with a number of global activity proxies well as a latitudinally-resolved magnetic index. As well as looking at frequency shifts in different frequency bands, we look at a parametrization of the shift as a cubic function of frequency. While the shifts in the medium- and highfrequency bands are very well correlated with all of the activity indices (with the best correlation being with the 10.7 cm radio flux), we confirm earlier findings that there appears to have been a change in the frequency response to activity during solar cycle 23, and the low frequency shifts are less correlated with activity in the last two cycles than they were in Cycle 22. At the same time, the more recent cycles show a slight increase in their sensitivity to activity levels at medium and higher frequencies, perhaps because a greater proportion of activity is composed of weaker or more ephemeral regions. This lends weight to the speculation that a fundamental change in the nature of the solar dynamo may be in progress.

The scaling relations that relate the average asteroseismic parameters $\Delta \nu$ and $\nu_{\max}$ to the global properties of stars are used quite extensively to determine stellar properties. While the $\Delta \nu$ scaling relation has been examined carefully and the deviations from the relation have been well documented, the $\nu_{\max}$ scaling relation has not been examined as extensively. In this paper we examine the $\nu_{\max}$ scaling relation using a set of stellar models constructed to have a wide range of mass, metallicity, and age. We find that as with $\Delta \nu$, $\nu_{\max}$ does not follow the simple scaling relation. The most visible deviation is because of a mean molecular weight term and a $\Gamma_1$ term that are commonly ignored. The remaining deviation is more difficult to address. We find that the influence of the scaling relation errors on asteroseismically derived values of $\log g$ are well within uncertainties. The influence of the errors on mass and radius estimates is small for main sequence and subgiants, but can be quite large for red giants.

The masses and radii of solar-like oscillators can be estimated through the asteroseismic scaling relations. These relations provide a direct link between observables, i.e. effective temperature and characteristics of the oscillation spectra, and stellar properties, i.e. mean density and surface gravity (thus mass and radius). These scaling relations are commonly used to characterize large samples of stars. Usually, the Sun is used as a reference from which the structure is scaled. However, for stars that do not have a similar structure as the Sun, using the Sun as a reference introduces systematic errors as large as 10\% in mass and 5\% in radius. Several alternatives for the reference of the scaling relation involving the large frequency separation (typical frequency difference between modes of the same degree and consecutive radial order) have been suggested in the literature. In a previous paper, we presented a reference function with a dependence on both effective temperature and metallicity. The accuracy of predicted masses and radii improved considerably when using reference values calculated from our reference function. However, the residuals indicated that stars on the red-giant branch possess a mass dependence that was not accounted for. Here, we present a reference function for the scaling relation involving the large frequency separation that includes the mass dependence. This new reference function improves the derived masses and radii significantly by removing the systematic differences and mitigates the trend with $\nu_{\rm max}$ (frequency of maximum oscillation power) that exists when using the solar value as a reference.

We create isochrones of M67 using the Yale Rotating Stellar Evolution Code. In addition to metallicity, parameters that are traditionally held fixed, such as the mixing length parameter and initial helium abundance, also vary. The amount of convective overshoot is also changed in different sets of isochrones. Models are constructed both with and without diffusion. From the resulting isochrones that fit the cluster, the age range is between 3.6 and 4.8 Gyr and the distance is between 755 and 868 pc. We also confirm Michaud et al. (2004) claim that M67 can be fit without overshoot if diffusion is included.

With the advent of dedicated photometric space missions, the ability to rapidly process huge catalogues of stars has become paramount. Bellinger and Angelou et al. (2016) recently introduced a new method based on machine learning for inferring the stellar parameters of main-sequence stars exhibiting solar-like oscillations. The method makes precise predictions that are consistent with other methods, but with the advantages of being able to explore many more parameters while costing practically no time. Here we apply the method to 52 so-called "LEGACY" main-sequence stars observed by the Kepler space mission. For each star, we present estimates and uncertainties of mass, age, radius, luminosity, core hydrogen abundance, surface helium abundance, surface gravity, initial helium abundance, and initial metallicity as well as estimates of their evolutionary model parameters of mixing length, overshooting coefficient, and diffusion multiplication factor. We obtain median uncertainties in stellar age, mass, and radius of 14.8%, 3.6%, and 1.7%, respectively. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at: https://github.com/earlbellinger/asteroseismology

From its surface properties it can be difficult to determine whether a red-giant star is in its helium-core-burning phase or only burning hydrogen in a shell around an inert helium core. Stars in either of these stages can have similar effective temperatures, radii and hence luminosities, i.e. they can be located at the same position in the Hertzsprung-Russell diagram. Asteroseismology -- the study of the internal structure of stars through their global oscillations -- can provide the necessary additional constraints to determine the evolutionary states of red-giant stars. Here, we present a method that uses grid-based modelling based on global asteroseismic properties ($\nu_{\rm max}$, frequency of maximum oscillation power; and $\Delta\nu$, frequency spacing between modes of the same degree and consecutive radial orders) as well as effective temperature and metallicity to determine the evolutionary phases. This method is applicable even to timeseries data of limited length, although with a small fraction of miss-classifications.

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity dependent convective mixing length is used, with $ \Delta\alpha_{\rm ML, YREC} \sim 0.2$ per dex in metallicity, a trend inconsistent with the predictions of three dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of 2 even at modest deviations from solar metallicity ([Fe/H]=$-$0.5).

Astronomy is in an era where all-sky surveys are mapping the Galaxy. The plethora of photometric, spectroscopic, asteroseismic and astrometric data allows us to characterise the comprising stars in detail. Here we quantify to what extent precise stellar observations reveal information about the properties of a star, including properties that are unobserved, or even unobservable. We analyse the diagnostic potential of classical and asteroseismic observations for inferring stellar parameters such as age, mass and radius from evolutionary tracks of solar-like oscillators on the lower main sequence. We perform rank correlation tests in order to determine the capacity of each observable quantity to probe structural components of stars and infer their evolutionary histories. We also analyse the principal components of classic and asteroseismic observables to highlight the degree of redundancy present in the measured quantities and demonstrate the extent to which information of the model parameters can be extracted. We perform multiple regression using combinations of observable quantities in a grid of evolutionary simulations and appraise the predictive utility of each combination in determining the properties of stars. We identify the combinations that are useful and provide limits to where each type of observable quantity can reveal information about a star. We investigate the accuracy with which targets in the upcoming TESS and PLATO missions can be characterized. We demonstrate that the combination of observations from GAIA and PLATO will allow us to tightly constrain stellar masses, ages and radii with machine learning for the purposes of Galactic and planetary studies.

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, in rough agreement with the observations, and unlike the simple expectation based on a Jeans length argument.

I review some recent work on low-mass star formation, with an emphasis on theory, basic principles, and unresolved questions. Star formation is both a gravitational fragmentation problem as well as an accretion problem. Molecular cloud structure can be understood as a fragmentation process driven by the interplay of turbulence, magnetic fields, and gravity (acting on either a dynamical or ambipolar-diffusion time scale). This results in a natural way to understand filamentary structure as magnetic ribbons that have an apparent width that scales differently than the Jeans length. Recent work also shows that stellar mass accretion through a disk is episodic. We show through numerical simulations that bursts of FU Ori type may be clustered, since the clump that accretes to the center is tidally sheared apart in its last stage of infall. Finally, we utilize a simplified model of stellar mass accretion and accretion termination to derive an analytic form for the initial mass function that has a lognormal-like body and a power-law tail. This scenario is consistent with an expectation of a larger number of substellar objects than may have been previously detected.

We compute a new generation of standard solar models (SSMs) that includes recent updates on some important nuclear reaction rates and a more consistent treatment of the equation of state. Models also include a novel and flexible treatment of opacity uncertainties based on opacity kernels, required in the light of recent theoretical and experimental works on radiative opacity. Two large sets of SSMs, each based on a different canonical set of solar abundances with high and low metallicity (Z), are computed to determine model uncertainties and correlations among different observables. We present detailed comparisons of high- and low-Z models against different ensembles of solar observables including solar neutrinos, surface helium abundance, depth of convective envelope and sound speed profile. A global comparison, including all observables, yields a p-value of 2.7$\sigma$ for the high-Z model and 4.7$\sigma$ for the low-Z one. When the sound-speed differences in the narrow region of $0.65 < r/R_{sun} < 0.70$ are excluded from the analysis, results are 0.9$\sigma$ and 3.0$\sigma$ for high- and low-Z models respectively. These results show that: high-Z models agree well with solar data but have a systematic problem right below the bottom of the convective envelope linked to steepness of molecular weight and temperature gradients, and that low-Z models lead to a much more general disagreement with solar data. We also show that, while simple parametrizations of opacity uncertainties can strongly alleviate the solar abundance problem, they are insufficient to substantially improve the agreement of SSMs with helioseismic data beyond that obtained for high-Z models due to the intrinsic correlations of theoretical predictions. All the detailed information of the new solar models, can be found online: http://www.ice.cat/personal/aldos/Solar_Data.html

Acoustic glitches are regions inside a star where the sound speed or its derivatives change abruptly. These leave a small characteristic oscillatory signature in the stellar oscillation frequencies. With the precision achieved by {\it Kepler} seismic data, it is now possible to extract these small amplitude oscillatory signatures, and infer the locations of the glitches. We perform glitch analysis for all the 66 stars in the {\it Kepler} seismic LEGACY sample to derive the locations of the base of the envelope convection zone and the helium ionization zone. The signature from helium ionization zone is found to be robust for all stars in the sample, whereas the convection zone signature is found to be weak and problematic, particularly for relatively massive stars with large errorbars on the oscillation frequencies. We demonstrate that the helium glitch signature can be used to constrain the properties of the helium ionization layers and the helium abundance.

Previous studies show that the physical structures and kinematics of a region depend significantly on the ionisation fraction. In this paper, we extend our previous studies of the effect of ionisation fractions on star formation to clouds that include both non-ideal magnetohydrodynamics and turbulence. We aim to quantify the importance of a treatment of the ionisation fraction in turbulent magnetised media and investigate the effect of turbulence on shaping the clouds and filaments before star formation sets in. In particular, we investigate how the structure, mass and width of filamentary structures depend on the amount of turbulence in ionised media and the initial mass-to-flux ratio. We compare the resulting density and mass-to-flux ratio structures both qualitatively and quantitatively via filament and core masses and filament fitting techniques (Gaussian and Plummer profiles.) We find that even with almost no turbulence, filamentary structure still exists. Comparison of simulations show that for turbulent Mach numbers above 2, there is little structural difference between the Step-Like (SL) and Cosmic Ray only (CR-only) ionisation models, while below this threshold the ionisation structure significantly affects the formation of filaments. Analysis of the mass within cores and filaments show decrease in mass as the degree of turbulence increases. Finally, observed filaments within the Taurus L1495/B213 complex are best reproduced by models with supercritical mass-to-flux ratios and/or at least mildly supersonic turbulence, however, our models show that sterile fibres observed within Taurus may occur in highly ionised, subcritical environments. Based on this, we suggest that regions with fertile fibres likely indicate a trans- or supercritical mass-to-flux ratio within the region while sterile fibres are likely subcritical and transient.

Determining the ages of red-giant stars is a key problem in stellar astrophysics. One of the difficulties in this determination is to know the evolutionary state of the individual stars -- i.e. have they started to burn Helium in their cores? That is the topic of this paper. Asteroseismic data provide a route to achieving this information. What we present here is an highly autonomous way of determining the evolutionary state from an analysis of the power spectrum of the light curve. The method is fast and efficient and can provide results for a large number of stars. It uses the structure of the dipole-mode oscillations, which have a mixed character in red-giant stars, to determine some measures that are used in the categorisation. It does not require that all the individual components of any given mode be separately characterised. Some 6604 red giant stars have been classified. Of these 3566 are determined to be on the red-giant branch, 2077 are red-clump and 439 are secondary-clump stars. We do not specifically identify the low metallicity, horizontal-branch stars. The difference between red-clump and secondary-clump stars is dependent on the manner in which Helium burning is first initiated. We discuss that the way the boundary between these classifications is set may lead to mis-categorisation in a small number of stars. The remaining 522 stars were not classified either because they lacked sufficient power in the dipole modes (so-called depressed dipole modes) or because of conflicting values in the parameters.

We use asteroseismic data from the Kepler satellite to determine fundamental stellar properties of the 66 main-sequence targets observed for at least one full year by the mission. We distributed tens of individual oscillation frequencies extracted from the time series of each star among seven modelling teams who applied different methods to determine radii, masses, and ages for all stars in the sample. Comparisons among the different results reveal a good level of agreement in all stellar properties, which is remarkable considering the variety of codes, input physics and analysis methods employed by the different teams. Average uncertainties are of the order of $\sim$2\% in radius, $\sim$4\% in mass, and $\sim$10\% in age, making this the best-characterised sample of main-sequence stars available to date. Our predicted initial abundances and mixing-length parameters are checked against inferences from chemical enrichment laws $\Delta Y / \Delta Z$ and predictions from 3D atmospheric simulations. We test the accuracy of the determined stellar properties by comparing them to the Sun, angular diameter measurements, Gaia parallaxes, and binary evolution, finding excellent agreement in all cases and further confirming the robustness of asteroseismically-determined physical parameters of stars when individual frequencies of oscillation are available. Baptised as the Kepler dwarfs LEGACY sample, these stars are the solar-like oscillators with the best asteroseismic properties available for at least another decade. All data used in this analysis and the resulting stellar parameters are made publicly available for the community.

The advent of space-based missions like $Kepler$ has revolutionized the study of solar-type stars, particularly through the measurement and modeling of their resonant modes of oscillation. Here we analyze a sample of 66 $Kepler$ main-sequence stars showing solar-like oscillations as part of the $Kepler$ seismic LEGACY project. We use $Kepler$ short-cadence data, of which each star has at least 12 months, to create frequency power spectra optimized for asteroseismology. For each star we identify its modes of oscillation and extract parameters such as frequency, amplitude, and line width using a Bayesian Markov chain Monte Carlo `peak-bagging' approach. We report the extracted mode parameters for all 66 stars, as well as derived quantities such as frequency difference ratios, the large and small separations $\Delta\nu$ and $\delta\nu_{02}$; the behavior of line widths with frequency and line widths at $\nu_{\rm max}$ with $T_{\rm eff}$, for which we derive parametrizations; and behavior of mode visibilities. These average properties can be applied in future peak-bagging exercises to better constrain the parameters of the stellar oscillation spectra. The frequencies and frequency ratios can tightly constrain the fundamental parameters of these solar-type stars, and mode line widths and amplitudes can test models of mode damping and excitation.

The solar-cyle variation of acoustic mode frequencies has a frequency dependence related to the inverse mode inertia. The discrepancy between model predictions and measured oscillation frequencies for solar and solar-type stellar acoustic modes includes a significant frequency-dependent term known as the surface term that is also related to the inverse mode inertia. We parametrize both the surface term and the frequency variations for low-degree solar data from Birmingham Solar-Oscillations Network (BiSON) and medium-degree data from the Global Oscillations Network Group (GONG) using the mode inertia together with cubic and inverse frequency terms. We find that for the central frequency of rotationally split multiplets the cubic term dominates both the average surface term and the temporal variation, but for the medium-degree case the inverse term improves the fit to the temporal variation. We also examine the variation of the even-order splitting coefficients for the medium-degree data and find that, as for the central frequency, the latitude-dependent frequency variation, which reflects the changing latitudinal distribution of magnetic activity over the solar cycle, can be described by the combination of a cubic and an inverse function of frequency scaled by inverse mode inertia. The results suggest that this simple parametrization could be used to assess the activity-related frequency variation in solar-like asteroseismic targets.

It is widely believed that the leading secular loop corrections from quantum gravity can be subsumed into a coordinate redefinition. Hence the apparent infrared logarithm corrections to any quantity would be just the result of taking the expectation value of the tree order quantity at the transformed coordinates in the graviton vacuum. We term this the Transformation Ansatz and we compare its predictions against explicit one loop computations in Maxwell + Einstein and Dirac + Einstein on de Sitter background. In each case the ansatz fails.

We present an asteroseismic analysis of 33 solar-type stars observed in short cadence during Campaigns (C) 1-3 of the NASA K2 mission. We were able to extract both average seismic parameters and individual mode frequencies for stars with dominant frequencies up to ~3300{\mu}Hz, and we find that data for some targets are good enough to allow for a measurement of the rotational splitting. Modelling of the extracted parameters is performed by using grid-based methods using average parameters and individual frequencies together with spectroscopic parameters. For the target selection in C3, stars were chosen as in C1 and C2 to cover a wide range in parameter space to better understand the performance and noise characteristics. For C3 we still detected oscillations in 73% of the observed stars that we proposed. Future K2 campaigns hold great promise for the study of nearby clusters and the chemical evolution and age-metallicity relation of nearby field stars in the solar neighbourhood. We expect oscillations to be detected in ~388 short-cadence targets if the K2 mission continues until C18, which will greatly complement the ~500 detections of solar-like oscillations made for short-cadence targets during the nominal Kepler mission. For ~30-40 of these, including several members of the Hyades open cluster, we furthermore expect that inference from interferometry should be possible.

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We argue that this is a natural model for filaments and is based on the interplay between turbulence, strong magnetic fields, and gravitationally-driven ambipolar diffusion, rather than pure gravity and thermal pressure. An analytic model for the formation of magnetic ribbons that is based on numerical simulations is used to derive a lateral width of a magnetic ribbon. This differs from the thickness along the magnetic field direction, which is essentially the Jeans scale. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, similar to observations, and unlike the simple expectation based on a Jeans length argument.

We model the mass accretion rate $\dot{M}$ to stellar mass $M_*$ correlation that has been inferred from observations of intermediate to upper mass T Tauri stars---that is $\dot{M} \propto M_*^{1.3 \pm 0.3}$. We explain this correlation within the framework of quiescent disk evolution, in which accretion is driven largely by gravitational torques acting in the bulk of the mass and volume of the disk. Stresses within the disk arise from the action of gravitationally driven torques parameterized in our 1D model in terms of Toomre's $Q$ criterion. We do not model the hot inner sub-AU scale region of the disk that is likely stable according to this criterion, and appeal to other mechanisms to remove or redistribute angular momentum and allow accretion onto the star. Our model has the advantage of agreeing with large-scale angle-averaged values from more complex nonaxisymmetric calculations. The model disk transitions from an early phase (dominated by initial conditions inherited from the burst mode of accretion) into a later self-similar mode characterized by a steeper temporal decline in $\dot{M}$. The models effectively reproduce the spread in mass accretion rates that have been observed for protostellar objects of $0.2\,\mbox{M}_\odot \le M_* \le 3.0\,\mbox{M}_\odot$, such as those found in the $\rho$ Ophiuchus and Taurus star forming regions. We then compare realistically sampled populations of young stellar objects produced by our model to their observational counterparts. We find these populations to be statistically coincident, which we argue is evidence for the role of gravitational torques in the late time evolution of quiescent protostellar disks.

The Hyades open cluster was targeted during Campaign 4 (C4) of the NASA K2 mission, and short-cadence data were collected on a number of cool main-sequence stars. Here, we report results on two F-type stars that show detectable oscillations of a quality that allows asteroseismic analyses to be performed. These are the first ever detections of solar-like oscillations in main-sequence stars in an open cluster.

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

The Orion star-forming region is the nearest active high-mass star-forming region and has created a large superbubble, the Orion-Eridanus superbubble. Recent work by Ochsendorf et al. has extended the accepted boundary of the superbubble. We fit Kompaneets models of superbubbles expanding in exponential atmospheres to the new, larger shape of the Orion-Eridanus superbubble. We find that this larger morphology of the superbubble is consistent with the evolution of the superbubble being primarily controlled by expansion into the exponential Galactic disk ISM if the superbubble is oriented with the Eridanus side farther from the Sun than the Orion side. Unlike previous Kompaneets model fits that required abnormally small scale heights for the Galactic disk (<40 pc), we find morphologically consistent models with scale heights of 80 pc, similar to that expected for the Galactic disk.

We present the current status of the analytic theory of brown dwarf evolution and the lower mass limit of the hydrogen burning main sequence stars. In the spirit of a simplified analytic theory we also introduce some modifications to the existing models. We give an exact expression for the pressure of an ideal non-relativistic Fermi gas at a finite temperature, therefore allowing for non-zero values of the degeneracy parameter ($\psi = \frac{kT}{\mu_{F}}$, where $\mu_{F}$ is the Fermi energy). We review the derivation of surface luminosity using an entropy matching condition and the first-order phase transition between the molecular hydrogen in the outer envelope and the partially-ionized hydrogen in the inner region. We also discuss the results of modern simulations of the plasma phase transition, which illustrate the uncertainties in determining its critical temperature. Based on the existing models and with some simple modification we find the maximum mass for a brown dwarf to be in the range $0.064M_\odot-0.087M_\odot$. An analytic formula for the luminosity evolution allows us to estimate the time period of the non-steady state (i.e., non-main sequence) nuclear burning for substellar objects. Standard models also predict that stars that are just above the substellar mass limit can reach an extremely low luminosity main sequence after at least a few million years of evolution, and sometimes much longer. We estimate that $\simeq 11 \%$ of stars take longer than $10^7$ yr to reach the main-sequence, and $\simeq 5 \%$ of stars take longer than $10^8$ yr.

Owing to the remarkable photometric precision of space observatories like Kepler, stellar and planetary systems beyond our own are now being characterized en masse for the first time. These characterizations are pivotal for endeavors such as searching for Earth-like planets and solar twins, understanding the mechanisms that govern stellar evolution, and tracing the dynamics of our Galaxy. The volume of data that is becoming available, however, brings with it the need to process this information accurately and rapidly. While existing methods can constrain fundamental stellar parameters such as ages, masses, and radii from these observations, they require substantial computational efforts to do so. We develop a method based on machine learning for rapidly estimating fundamental parameters of main-sequence solar-like stars from classical and asteroseismic observations. We first demonstrate this method on a hare-and-hound exercise and then apply it to the Sun, 16 Cyg A & B, and 34 planet-hosting candidates that have been observed by the Kepler spacecraft. We find that our estimates and their associated uncertainties are comparable to the results of other methods, but with the additional benefit of being able to explore many more stellar parameters while using much less computation time. We furthermore use this method to present evidence for an empirical diffusion-mass relation. Our method is open source and freely available for the community to use. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at https://github.com/earlbellinger/asteroseismology

$\theta$ Cygni is an F3 spectral-type main-sequence star with visual magnitude V=4.48. This star was the brightest star observed by the original Kepler spacecraft mission. Short-cadence (58.8 s) photometric data using a custom aperture were obtained during Quarter 6 (June-September 2010) and subsequently in Quarters 8 and 12-17. We present analyses of the solar-like oscillations based on Q6 and Q8 data, identifying angular degree $l$ = 0, 1, and 2 oscillations in the range 1000-2700 microHz, with a large frequency separation of 83.9 plus/minus 0.4 microHz, and frequency with maximum amplitude 1829 plus/minus 54 microHz. We also present analyses of new ground-based spectroscopic observations, which, when combined with angular diameter measurements from interferometry and Hipparcos parallax, give T_eff = 6697 plus/minus 78 K, radius 1.49 plus/minus 0.03 solar radii, [Fe/H] = -0.02 plus/minus 0.06 dex, and log g = 4.23 plus/minus 0.03. We calculate stellar models matching the constraints using several methods, including using the Yale Rotating Evolution Code and the Asteroseismic Modeling Portal. The best-fit models have masses 1.35-1.39 solar masses and ages 1.0-1.6 Gyr. theta Cyg's T_eff and log g place it cooler than the red edge of the gamma Doradus instability region established from pre-Kepler ground-based observations, but just at the red edge derived from pulsation modeling. The pulsation models show gamma Dor gravity-mode pulsations driven by the convective-blocking mechanism, with frequencies of 1 to 3 cycles/day (11 to 33 microHz). However, gravity modes were not detected in the Kepler data, one signal at 1.776 cycles/day (20.56 microHz) may be attributable to a faint, possibly background, binary. Asteroseismic studies of theta Cyg and other A-F stars observed by Kepler and CoRoT, will help to improve stellar model physics and to test pulsation driving mechanisms.

The seismic study of the Sun and other stars offers a unique window into the interior of these stars. Thanks to helioseismology, we know the structure of the Sun to admirable precision. In fact, our knowledge is good enough to use the Sun as a laboratory. We have also been able to study the dynamics of the Sun in great detail. Helioseismic data also allow us to probe the changes that take place in the Sun as solar activity waxes and wanes. The seismic study of stars other than the Sun is a fairly new endeavour, but we are making great strides in this field. In this review I discuss some of the techniques used in helioseismic analyses and the results obtained using those techniques. In this review I focus on results obtained with global helioseismology, i.e., the study of the Sun using its normal modes of oscillation. I also briefly touch upon asteroseismology, the seismic study of stars other than the Sun, and discuss how seismic data of others stars are interpreted.

The scaling relations between global asteroseismic observables and stellar properties are widely used to estimate masses and radii of stars exhibiting solar-like oscillations. Since the mass and radius of the Sun are known independently, the Sun is commonly used as a reference to scale to. However, the validity of the scaling relations depends on the homology between the star under study and the reference star. Solar-like oscillators span a wide range of masses and metallicities, as well as evolutionary phases. Most of these stars are therefore not homologous to the Sun. This leads to errors of up to 10% (5%) in mass (radius) when using the asteroseismic scaling relations with the Sun as the reference. In this paper we derive a reference function to replace the solar-reference value used in the large-frequency-separation scaling relation. Our function is the first that depends on both effective temperature and metallicity, and is applicable from the end of the main sequence to just above the bump on the red giant branch. This reference function improves the estimates of masses and radii determined through scaling relations by a factor of 2, i.e. allows masses and radii to be recovered with an accuracy of 5% and 2%, respectively.

Context: Detailed oscillation spectra comprising individual frequencies for numerous solar-type stars and red giants are or will become available. These data can lead to a precise characterisation of stars. Aims: Our goal is to test and compare different methods for obtaining stellar properties from oscillation frequencies and spectroscopic constraints, in order to evaluate their accuracy and the reliability of the error bars. Methods: In the context of the SpaceInn network, we carried out a hare-and-hounds exercise in which one group produced "observed" oscillation spectra for 10 artificial solar-type stars, and various groups characterised these stars using either forward modelling or acoustic glitch signatures. Results: Results based on the forward modelling approach were accurate to 1.5 % (radius), 3.9 % (mass), 23 % (age), 1.5 % (surface gravity), and 1.8 % (mean density). For the two 1 Msun stellar targets, the accuracy on the age is better than 10 % thereby satisfying PLATO 2.0 requirements. The average accuracies for the acoustic radii of the base of the convection zone, the He II ionisation, and the Gamma_1 peak were 17 %, 2.4 %, and 1.9 %, respectively. Glitch fitting analysis seemed to be affected by aliasing problems for some of the targets. Conclusions: Forward modelling is the most accurate approach, but needs to be complemented by model-independent results from, e.g., glitch analysis. Furthermore, global optimisation algorithms provide more robust error bars.

One of the largest sources of uncertainty in stellar models is caused by the treatment of convection in stellar envelopes. One dimensional stellar models often make use of the mixing length or equivalent approximations to describe convection, all of which depend on various free parameters. There have been attempts to rectify this by using 3D radiative-hydrodynamic simulations of stellar convection, and in trying to extract an equivalent mixing length from the simulations. In this paper we show that the entropy of the deeper, adiabatic layers in these simulations can be expressed as a simple function of og g and log T_{eff} which holds potential for calibrating stellar models in a simple and more general manner.

Simulations predict that hot super-Earth sized exoplanets can have their envelopes stripped by photo-evaporation, which would present itself as a lack of these exoplanets. However, this absence in the exoplanet population has escaped a firm detection. Here we demonstrate, using asteroseismology on a sample of exoplanets and exoplanet candidates observed during the Kepler mission that, while there is an abundance of super-Earth sized exoplanets with low incident fluxes, none are found with high incident fluxes. We do not find any exoplanets with radii between 2.2 and 3.8 Earth radii with incident flux above 650 times the incident flux on Earth. This gap in the population of exoplanets is explained by evaporation of volatile elements and thus supports the predictions. The confirmation of a hot-super-Earth desert caused by evaporation will add an important constraint on simulations of planetary systems, since they must be able to reproduce the dearth of close-in super-Earths.

The angle $\psi$ between a planet's orbital axis and the spin axis of its parent star is an important diagnostic of planet formation, migration, and tidal evolution. We seek empirical constraints on $\psi$ by measuring the stellar inclination $i_{\rm s}$ via asteroseismology for an ensemble of 25 solar-type hosts observed with NASA's Kepler satellite. Our results for $i_{\rm s}$ are consistent with alignment at the 2-$\sigma$ level for all stars in the sample, meaning that the system surrounding the red-giant star Kepler-56 remains as the only unambiguous misaligned multiple-planet system detected to date. The availability of a measurement of the projected spin-orbit angle $\lambda$ for two of the systems allows us to estimate $\psi$. We find that the orbit of the hot-Jupiter HAT-P-7b is likely to be retrograde ($\psi=116.4^{+30.2}_{-14.7}\:{\rm deg}$), whereas that of Kepler-25c seems to be well aligned with the stellar spin axis ($\psi=12.6^{+6.7}_{-11.0}\:{\rm deg}$). While the latter result is in apparent contradiction with a statement made previously in the literature that the multi-transiting system Kepler-25 is misaligned, we show that the results are consistent, given the large associated uncertainties. Finally, we perform a hierarchical Bayesian analysis based on the asteroseismic sample in order to recover the underlying distribution of $\psi$. The ensemble analysis suggests that the directions of the stellar spin and planetary orbital axes are correlated, as conveyed by a tendency of the host stars to display large inclination values.

Kepler-454 (KOI-273) is a relatively bright (V = 11.69 mag), Sun-like starthat hosts a transiting planet candidate in a 10.6 d orbit. From spectroscopy, we estimate the stellar temperature to be 5687 +/- 50 K, its metallicity to be [m/H] = 0.32 +/- 0.08, and the projected rotational velocity to be v sin i <2.4 km s-1. We combine these values with a study of the asteroseismic frequencies from short cadence Kepler data to estimate the stellar mass to be 1.028+0:04-0:03 M_Sun, the radius to be 1.066 +/- 0.012 R_Sun and the age to be 5.25+1:41-1:39 Gyr. We estimate the radius of the 10.6 d planet as 2.37 +/- 0.13 R_Earth. Using 63 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 36 observations made with the HIRES spectrograph at Keck Observatory, we measure the mass of this planet to be 6.8 +/- 1.4M_Earth. We also detect two additional non-transiting companions, a planet with a minimum mass of 4.46 +/- 0.12 M_J in a nearly circular 524 d orbit and a massive companion with a period >10 years and mass >12.1M_J . The twelve exoplanets with radii <2.7 R_Earth and precise mass measurements appear to fall into two populations, with those <1.6 R_Earth following an Earth-like composition curve and larger planets requiring a significant fraction of volatiles. With a density of 2.76 +/- 0.73 g cm-3, Kepler-454b lies near the mass transition between these two populations and requires the presence of volatiles and/or H/He gas.

\Kepler has revolutionised our understanding of both exoplanets and their host stars. Asteroseismology is a valuable tool in the characterisation of stars and \Kepler is an excellent observing facility to perform asteroseismology. Here we select a sample of 35 \Kepler solar-type stars which host transiting exoplanets (or planet candidates) with detected solar-like oscillations. Using available \Kepler short cadence data up to Quarter 16 we create power spectra optimised for asteroseismology of solar-type stars. We identify modes of oscillation and estimate mode frequencies by ``peak bagging'' using a Bayesian MCMC framework. In addition, we expand the methodology of quality assurance using a Bayesian unsupervised machine learning approach. We report the measured frequencies of the modes of oscillation for all 35 stars and frequency ratios commonly used in detailed asteroseismic modelling. Due to the high correlations associated with frequency ratios we report the covariance matrix of all frequencies measured and frequency ratios calculated. These frequencies, frequency ratios, and covariance matrices can be used to obtain tight constraint on the fundamental parameters of these planet-hosting stars.

I review some steps in the conversion of molecular cloud gas into stars and planets, with an emphasis in this presentation on the early stage molecular cloud fragmentation that leads to elongated filaments/ribbons. Magnetic fields can play a crucial role in all stages and need to be invoked particularly for early stage fragmentation as well as in late core collapse where it may control disk formation. I also review some elements of hydrodynamic modeling of disk evolution.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

We present a new statistical analysis that combines helioseismology (sound speed, surface helium and convective radius) and solar neutrino observations (the $^8$B and $^7$Be fluxes) to place upper limits to the properties of non standard weakly interacting particles. Our analysis includes theoretical and observational errors, accounts for tensions between input parameters of solar models and can be easily extended to include other observational constraints. We present two applications to test the method: the well studied case of axions and axion-like particles and the more novel case of low mass hidden photons. For axions we obtain an upper limit at $3\sigma$ for the axion-photon coupling constant of $g_{a\gamma}\,<\,4.1 \cdot 10^{-10} \rm{GeV^{-1}}$. For hidden photons we obtain the most restrictive upper limit available accross a wide range of masses for the product of the kinetic mixing and mass of $\chi m < 1.8 \cdot 10^{-12} \rm{eV}$ at $3\sigma$. Both cases improve the previous solar constraints based on the Standard Solar Models showing the power of using a global statistical approach.

NASA's re-purposed Kepler mission -- dubbed K2 -- has brought new scientific opportunities that were not anticipated for the original Kepler mission. One science goal that makes optimal use of K2's capabilities, in particular its 360-degree ecliptic field of view, is galactic archaeology -- the study of the evolution of the Galaxy from the fossil stellar record. The thrust of this research is to exploit high-precision, time-resolved photometry from K2 in order to detect oscillations in red giant stars. This asteroseismic information can provide estimates of stellar radius (hence distance), mass and age of vast numbers of stars across the Galaxy. Here we present the initial analysis of a subset of red giants, observed towards the North Galactic Gap, during the mission's first full science campaign. We investigate the feasibility of using K2 data for detecting oscillations in red giants that span a range in apparent magnitude and evolutionary state (hence intrinsic luminosity). We demonstrate that oscillations are detectable for essentially all cool giants within the $\log g$ range $\sim 1.9-3.2$. Our detection is complete down to $\mathit{Kp}\sim 14.5$, which results in a seismic sample with little or no detection bias. This sample is ideally suited to stellar population studies that seek to investigate potential shortcomings of contemporary Galaxy models.

We present the first detections by the NASA K2 Mission of oscillations in solar-type stars, using short-cadence data collected during K2 Campaign\,1 (C1). We understand the asteroseismic detection thresholds for C1-like levels of photometric performance, and we can detect oscillations in subgiants having dominant oscillation frequencies around $1000\,\rm \mu Hz$. Changes to the operation of the fine-guidance sensors are expected to give significant improvements in the high-frequency performance from C3 onwards. A reduction in the excess high-frequency noise by a factor of two-and-a-half in amplitude would bring main-sequence stars with dominant oscillation frequencies as high as ${\simeq 2500}\,\rm \mu Hz$ into play as potential asteroseismic targets for K2.

We present a study of 33 {\it Kepler} planet-candidate host stars for which asteroseismic observations have sufficiently high signal-to-noise ratio to allow extraction of individual pulsation frequencies. We implement a new Bayesian scheme that is flexible in its input to process individual oscillation frequencies, combinations of them, and average asteroseismic parameters, and derive robust fundamental properties for these targets. Applying this scheme to grids of evolutionary models yields stellar properties with median statistical uncertainties of 1.2\% (radius), 1.7\% (density), 3.3\% (mass), 4.4\% (distance), and 14\% (age), making this the exoplanet host-star sample with the most precise and uniformly determined fundamental parameters to date. We assess the systematics from changes in the solar abundances and mixing-length parameter, showing that they are smaller than the statistical errors. We also determine the stellar properties with three other fitting algorithms and explore the systematics arising from using different evolution and pulsation codes, resulting in 1\% in density and radius, and 2\% and 7\% in mass and age, respectively. We confirm previous findings of the initial helium abundance being a source of systematics comparable to our statistical uncertainties, and discuss future prospects for constraining this parameter by combining asteroseismology and data from space missions. Finally we compare our derived properties with those obtained using the global average asteroseismic observables along with effective temperature and metallicity, finding an excellent level of agreement. Owing to selection effects, our results show that the majority of the high signal-to-noise ratio asteroseismic {\it Kepler} host stars are older than the Sun.

Asteroseismology is a powerful tool that can precisely characterize the mass, radius, and other properties of field stars. However, our inability to properly model the near-surface layers of stars creates a frequency-dependent frequency difference between the observed and the modeled frequencies, usually referred to as the "surface term". This surface term can add significant errors to the derived stellar properties unless removed properly. In this paper we simulate surface terms across a significant portion of the HR diagram, exploring four different masses ($M=0.8, 1.0, 1.2$, and $1.5$ M$_\odot$) at five metallicities ($[\rm{Fe/H}]=0.5, 0.0, -0.5 ,-1.0, and -1.5$) from main sequence to red giants for stars with $T_{\rm{eff}}<6500 K$ and explore how well the most common ways of fitting and removing the surface term actually perform. We find that the two-term model proposed by Ball & Gizon (2014) works much better than other models across a large portion of the HR diagram, including the red giants, leading us to recommend its use for future asteroseismic analyses.

We derive age constraints for 1639 red giants in the APOKASC sample for which seismic parameters from Kepler, as well as effective temperatures, metallicities and [alpha/Fe] values from APOGEE DR12 are available. We investigate the relation between age and chemical abundances for these stars, using a simple and robust approach to obtain ages. We first derive stellar masses using standard seismic scaling relations, then determine the maximum possible age for each star as function of its mass and metallicity, independently of its evolutionary stage. While the overall trend between maximum age and chemical abundances is a declining fraction of young stars with increasing [alpha/Fe], at least 14 out of 241 stars with [alpha/Fe]>0.13 are younger than 6 Gyr. Five stars with [alpha/Fe]>0.2 have ages below 4 Gyr. We examine the effect of modifications in the standard seismic scaling relations, as well as the effect of very low helium fractions, but these changes are not enough to make these stars as old as usually expected for alpha-rich stars (i.e., ages greater than 8-9 Gyr). Such unusual alpha-rich young stars have also been detected by other surveys, but defy simple explanations in a galaxy evolution context.

Ring-diagram analysis of acoustic waves observed at the photosphere can provide a relatively robust determination of the sub-surface flows at a particular time under a particular region. The depth of penetration of the waves is related to the size of the region, hence the depth extent of the measured flows is inversely proportional to the spatial resolution. Most ring-diagram analysis has focused on regions of extent ~15{\deg} (180 Mm) or more in order to provide reasonable mode sets for inversions. HMI data analysis also provides a set of ring fit parameters on a scale three times smaller. These provide flow estimates for the outer 1% (7 Mm) of the Sun only, with very limited depth resolution, but with spatial resolution adequate to map structures potentially associated with the belts and regions of magnetic activity. There are a number of systematic effects affecting the determination of flows from local helioseismic analysis of regions over different parts of the observable disk, not all well understood. In this study we characterize those systematic effects with higher spatial resolution, so that they may more effectively be accounted for in mapping temporal and spatial evolution of the flows. Leaving open the question of the mean structure of the global meridional circulation and the differential rotation, we describe the near-surface flow anomalies in time and latitude corresponding to the torsional oscillation pattern in differential rotation and analogous patterns in the meridional cell structure as observed by SDO/HMI.

Large-scale analyses of stellar samples comprised of cool, solar-like oscillators now commonly utilize the so-called asteroseismic scaling relations to estimate fundamental stellar properties. In this paper we present a test of the scaling relation for the global asteroseismic parameter $\nu_{\rm max}$, the frequency at which a solar-like oscillator presents its strongest observed pulsation amplitude. The classic relation assumes that this characteristic frequency scales with a particular combination of surface gravity and effective temperature that also describes the dependence of the cut-off frequency for acoustic waves in an isothermal atmosphere, i.e., $\nu_{\rm max} \propto gT_{\rm eff}^{-1/2}$. We test how well the oscillations of cool main-sequence and sub-giant stars adhere to this relation, using a sample of asteroseismic targets observed by the NASA \emph{Kepler} Mission. Our results, which come from a grid-based analysis, rule out departures from the classic $gT_{\rm eff}^{-1/2}$ scaling dependence at the level of $\simeq 1.5\,\rm per cent$ over the full $\simeq 1560\,\rm K$ range in $T_{\rm eff}$ that we tested. There is some uncertainty over the absolute calibration of the scaling. However, any variation with $T_{\rm eff}$ is evidently small, with limits similar to those above.

We report the discovery of Kepler-432b, a giant planet ($M_b = 5.41^{+0.32}_{-0.18} M_{\rm Jup}, R_b = 1.145^{+0.036}_{-0.039} R_{\rm Jup}$) transiting an evolved star $(M_\star = 1.32^{+0.10}_{-0.07} M_\odot, R_\star = 4.06^{+0.12}_{-0.08} R_\odot)$ with an orbital period of $P_b = 52.501129^{+0.000067}_{-0.000053}$ days. Radial velocities (RVs) reveal that Kepler-432b orbits its parent star with an eccentricity of $e = 0.5134^{+0.0098}_{-0.0089}$, which we also measure independently with asterodensity profiling (AP; $e=0.507^{+0.039}_{-0.114}$), thereby confirming the validity of AP on this particular evolved star. The well-determined planetary properties and unusually large mass also make this planet an important benchmark for theoretical models of super-Jupiter formation. Long-term RV monitoring detected the presence of a non-transiting outer planet (Kepler-432c; $M_c \sin{i_c} = 2.43^{+0.22}_{-0.24} M_{\rm Jup}, P_c = 406.2^{+3.9}_{-2.5}$ days), and adaptive optics imaging revealed a nearby (0\farcs87), faint companion (Kepler-432B) that is a physically bound M dwarf. The host star exhibits high signal-to-noise asteroseismic oscillations, which enable precise measurements of the stellar mass, radius and age. Analysis of the rotational splitting of the oscillation modes additionally reveals the stellar spin axis to be nearly edge-on, which suggests that the stellar spin is likely well-aligned with the orbit of the transiting planet. Despite its long period, the obliquity of the 52.5-day orbit may have been shaped by star-planet interaction in a manner similar to hot Jupiter systems, and we present observational and theoretical evidence to support this scenario. Finally, as a short-period outlier among giant planets orbiting giant stars, study of Kepler-432b may help explain the distribution of massive planets orbiting giant stars interior to 1 AU.

We demonstrate a new procedure to derive accurate and precise surface gravities from high resolution spectra without the use of external constraints. Our analysis utilizes Spectroscopy Made Easy (SME) with robust spectral line constraints and uses an iterative process to mitigate degeneracies in the fitting process. We adopt an updated radiative transfer code, a new treatment for neutral perturber broadening, a line list with multiple gravity constraints and separate fitting for global stellar properties and abundance determinations. To investigate the sources of temperature dependent trends in determining log g noted in previous studies, we obtained Keck HIRES spectra of 42 Kepler asteroseismic stars. In comparison to asteroseismically determined log g our spectroscopic analysis has a constant offset of 0.01 dex with a root mean square (RMS) scatter of 0.05 dex. We also analyzed 30 spectra which had published surface gravities determined using the $a/R_*$ technique from planetary transits and found a constant offset of 0.06 dex and RMS scatter of 0.07 dex. The two samples covered effective temperatures between 5000K and 6700K with log g between 3.7 and 4.6.

We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star formation. Using our numerical hydrodynamics model with improved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FU-Orionis-type eruptions. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods (~ 10^4 yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40\%--70\% of the star-forming cores can display bursts after forming a star-disk system.

We analyze the time evolution of the luminosity of a cluster of Population III protostars formed in the early universe. We argue from the Jeans criterion that primordial gas can collapse to form a cluster of first stars that evolve relatively independently of one another (i.e., with negligible gravitational interaction). We model the collapse of individual protostellar clumps using 2+1D nonaxisymmetric numerical hydrodynamics simulations. Each collapse produces a protostar surrounded by a massive disk (i.e., $M_{\rm disk} / M_{*} \gtrsim 0.1$), whose evolution we follow for a further 30--40 kyr. Gravitational instabilities result in the fragmentation and the formation of gravitationally bound clumps within the disk. The accretion of these fragments by the host protostar produces accretion and luminosity bursts on the order of $10^6\,\LSun$. Within the cluster, we show that a simultaneity of such events across several protostellar cluster members can elevate the cluster luminosity to 5--10${\times}$ greater than expected, and that the cluster spends $\sim15\%$ of it's star-forming history at these levels. This enhanced luminosity effect is particularly enabled in clusters of modest size with $\simeq$ 10--20 members. In one such instance, we identify a confluence of burst events that raise the luminosity to nearly $1000{\times}$ greater than the cluster mean luminosity, resulting in $L > 10^8\,\LSun$. This phenomenon arises solely through the gravitational-instability--driven episodic fragmentation and accretion that characterizes this early stage of protostellar evolution.

This work explores the mathematical properties of a distribution introduced by Basu & Jones (2004), and applies it to model the stellar initial mass function (IMF). The distribution arises simply from an initial lognormal distribution, requiring that each object in it subsequently undergoes exponential growth but with an exponential distribution of growth lifetimes. This leads to a modified lognormal with a power-law tail (MLP) distribution, which can in fact be applied to a wide range of fields where distributions are observed to have a lognormal-like body and a power-law tail. We derive important properties of the MLP distribution, like the cumulative distribution, the mean, variance, arbitrary raw moments, and a random number generator. These analytic properties of the distribution can be used to facilitate application to modeling the IMF. We demonstrate how the MLP function provides an excellent fit to the IMF compiled by Chabrier (2005) and how this fit can be used to quickly identify quantities like the mean, median, and mode, as well as number and mass fractions in different mass intervals.

The first discoveries of exoplanets around Sun-like stars have fueled efforts to find ever smaller worlds evocative of Earth and other terrestrial planets in the Solar System. While gas-giant planets appear to form preferentially around metal-rich stars, small planets (with radii less than four Earth radii) can form under a wide range of metallicities. This implies that small, including Earth-size, planets may have readily formed at earlier epochs in the Universe's history when metals were far less abundant. We report Kepler spacecraft observations of KOI-3158, a metal-poor Sun-like star from the old population of the Galactic thick disk, which hosts five planets with sizes between Mercury and Venus. We used asteroseismology to directly measure a precise age of 11.2+/-1.0 Gyr for the host star, indicating that KOI-3158 formed when the Universe was less than 20% of its current age and making it the oldest known system of terrestrial-size planets. We thus show that Earth-size planets have formed throughout most of the Universe's 13.8-billion-year history, providing scope for the existence of ancient life in the Galaxy.

The chemical composition of stars hosting small exoplanets (with radii less than four Earth radii) appears to be more diverse than that of gas-giant hosts, which tend to be metal-rich. This implies that small, including Earth-size, planets may have readily formed at earlier epochs in the Universe's history when metals were more scarce. We report Kepler spacecraft observations of Kepler-444, a metal-poor Sun-like star from the old population of the Galactic thick disk and the host to a compact system of five transiting planets with sizes between those of Mercury and Venus. We validate this system as a true five-planet system orbiting the target star and provide a detailed characterization of its planetary and orbital parameters based on an analysis of the transit photometry. Kepler-444 is the densest star with detected solar-like oscillations. We use asteroseismology to directly measure a precise age of 11.2+/-1.0 Gyr for the host star, indicating that Kepler-444 formed when the Universe was less than 20% of its current age and making it the oldest known system of terrestrial-size planets. We thus show that Earth-size planets have formed throughout most of the Universe's 13.8-billion-year history, leaving open the possibility for the existence of ancient life in the Galaxy. The age of Kepler-444 not only suggests that thick-disk stars were among the hosts to the first Galactic planets, but may also help to pinpoint the beginning of the era of planet formation.

We present direct radii measurements of the well-known transiting exoplanet host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and 0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD 189733 and HD 209458 are currently the only two transiting exoplanet systems where detection of the respective planetary companion's orbital motion from high resolution spectroscopy has revealed absolute masses for both star and planet. We use our new measurements together with the orbital information from radial velocity and photometric time series data, Hipparcos distances, and newly measured bolometric fluxes to determine the stellar effective temperatures (T_eff = 4875 +/- 43, 6093 +/- 103 K), stellar linear radii (R_* = 0.805 +/- 0.016, 1.203 +/- 0.061 R_sun), mean stellar densities (rho_* = 1.62 +/- 0.11, 0.58 +/- 0.14 rho_sun), planetary radii (R_p = 1.216 +/- 0.024, 1.451 +/- 0.074 R_Jup), and mean planetary densities (rho_p = 0.605 +/- 0.029, 0.196 +/- 0.033 rho_Jup) for HD 189733 b and HD 209458 b, respectively. The stellar parameters for HD 209458, a F9 dwarf, are consistent with indirect estimates derived from spectroscopic and evolutionary modeling. However, we find that models are unable to reproduce the observational results for the K2 dwarf, HD 189733. We show that, for stellar evolutionary models to match the observed stellar properties of HD 189733, adjustments lowering the solar-calibrated mixing length parameter from 1.83 to 1.34 need to be employed.

The solar analogs 16 Cyg A and 16 Cyg B are excellent asteroseismic targets in the \Kepler field of view and together with a red dwarf and a Jovian planet form an interesting system. For these more evolved Sun-like stars we cannot detect surface rotation with the current \Kepler data but instead use the technique of asteroseimology to determine rotational properties of both 16 Cyg A and B. We find the rotation periods to be $23.8^{+1.5}_{-1.8} \rm \, days$ and $23.2^{+11.5}_{-3.2} \rm \, days$, and the angles of inclination to be $56^{+6}_{-5} \, ^{\circ}$ and $36^{+17}_{-7} \, ^{\circ}$, for A and B respectively. Together with these results we use the published mass and age to suggest that, under the assumption of a solar-like rotation profile, 16 Cyg A could be used when calibrating gyrochronology relations. In addition, we discuss the known 16 Cyg B star-planet eccentricity and measured low obliquity which is consistent with Kozai cycling and tidal theory.

Bailey & Basu (2014) show analysis of density and mass-to-flux ratio maps for simulations with either an ionisation profile which takes into account photoionisation (step-like profile) or a cosmic ray only ionisation profile. We extend this study to analyse the effect of these ionisation profiles on velocity structures, kinematics, and synthetic spectra. Clump regions are found to occur at the convergence of two flows with a low velocity region and velocity direction transition occurring at the junction. Models with evident substructure show that core formation occurs on the periphery of these velocity valleys. Analysis of synthetic spectra reveals the presence of large non-thermal components within low-density gas, especially for models with the step-like ionisation profile. All cores show small, sub-thermal relative motions compared to background gas. Large deviations within this analysis are due to the line of sight intersecting low- and high-density regions across the velocity switch transition. Positive deviations correspond to a foreground core moving away from the observer while negative deviations correspond to a background core moving toward the observer. Comparison of velocities resulting from different ionisation profiles suggest that high ionisation fractions yield supersonic velocities, up to two times the sound speed, while regions with low ionisation fractions tend to be subsonic or mildly supersonic. This suggests that the transition to coherence within cores could be a transition between high and low ionisation fractions within the gas.

We present the first APOKASC catalog of spectroscopic and asteroseismic properties of 1916 red giants observed in the Kepler fields. The spectroscopic parameters provided from the Apache Point Observatory Galactic Evolution Experiment project are complemented with asteroseismic surface gravities, masses, radii, and mean densities determined by members of the Kepler Asteroseismology Science Consortium. We assess both random and systematic sources of error and include a discussion of sample selection for giants in the Kepler fields. Total uncertainties in the main catalog properties are of order 80 K in Teff , 0.06 dex in [M/H], 0.014 dex in log g, and 12% and 5% in mass and radius, respectively; these reflect a combination of systematic and random errors. Asteroseismic surface gravities are substantially more precise and accurate than spectroscopic ones, and we find good agreement between their mean values and the calibrated spectroscopic surface gravities. There are, however, systematic underlying trends with Teff and log g. Our effective temperature scale is between 0-200 K cooler than that expected from the Infrared Flux Method, depending on the adopted extinction map, which provides evidence for a lower value on average than that inferred for the Kepler Input Catalog (KIC). We find a reasonable correspondence between the photometric KIC and spectroscopic APOKASC metallicity scales, with increased dispersion in KIC metallicities as the absolute metal abundance decreases, and offsets in Teff and log g consistent with those derived in the literature. We present mean fitting relations between APOKASC and KIC observables and discuss future prospects, strengths, and limitations of the catalog data.

We present a first determination of distances and extinctions for individual stars in the first release of the APOKASC catalogue, built from the joint efforts of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and the Kepler Asteroseismic Science Consortium (KASC). Our method takes into account the spectroscopic constraints derived from the APOGEE Stellar Parameters and Chemical Abundances Pipeline, together with the asteroseismic parameters from KASC. These parameters are then employed to estimate intrinsic stellar properties, including absolute magnitudes, using the Bayesian tool PARAM. We then find the distance and extinction that best fit the observed photometry in SDSS, 2MASS, and WISE passbands. The first 1989 giants targeted by APOKASC are found at typical distances between 0.5 and 5 kpc, with individual uncertainties of just ~1.8 per cent. Our extinction estimates are systematically smaller than provided in the Kepler Input Catalogue and by the Schlegel, Finkbeiner and Davis maps. Distances to individual stars in the NGC 6791 and NGC 6819 star clusters agree to within their credible intervals. Comparison with the APOGEE red clump and SAGA catalogues provide another useful check, exhibiting agreement with our measurements to within a few percent. Overall, present methods seem to provide excellent distance and extinction determinations for the bulk of the APOKASC sample. Approximately one third of the stars present broad or multiple-peaked probability density functions and hence increased uncertainties. Uncertainties are expected to be reduced in future releases of the catalogue, when a larger fraction of the stars will have seismically-determined evolutionary status classifications.

We consider the prospects for detecting solar-like oscillations in the "super-Nyquist" regime of long-cadence (LC) Kepler photometry, i.e., above the associated Nyquist frequency of approximately 283 {\mu}Hz. Targets of interest are cool, evolved subgiants and stars lying at the base of the red-giant branch. These stars would ordinarily be studied using the short-cadence (SC) data, since the associated SC Nyquist frequency lies well above the frequencies of the detectable oscillations. However, the number of available SC target slots is quite limited. This imposes a severe restriction on the size of the ensemble available for SC asteroseismic study.We find that archival Kepler LC data from the nominal Mission may be utilized for asteroseismic studies of targets whose dominant oscillation frequencies lie as high as approximately 500 {\mu}Hz, i.e., about 1.75- times the LC Nyquist frequency. The frequency detection threshold for the shorter-duration science campaigns of the re-purposed Kepler Mission, K2, is lower. The maximum threshold will probably lie somewhere between approximately 400 and 450 {\mu}Hz. The potential to exploit the archival Kepler and K2 LC data in this manner opens the door to increasing significantly the number of subgiant and low-luminosity red-giant targets amenable to asteroseismic analysis, overcoming target limitations imposed by the small number of SC slots.We estimate that around 400 such targets are now available for study in the Kepler LC archive. That number could potentially be a lot higher for K2, since there will be a new target list for each of its campaigns.

We clarify the mechanism of accelerated core formation by large-scale nonlinear flows in subcritical magnetic clouds by finding a semi-analytical formula for the core formation time and describing the physical processes that lead to them. Recent numerical simulations show that nonlinear flows induce rapid ambipolar diffusion that leads to localized supercritical regions that can collapse. Here, we employ non-ideal magnetohydrodynamic simulations including ambipolar diffusion for gravitationally stratified sheets threaded by vertical magnetic fields. One of the horizontal dimensions is eliminated, resulting in a simpler two-dimensional simulation that can clarify the basic process of accelerated core formation. A parameter study of simulations shows that the core formation time is inversely proportional to the square of the flow speed when the flow speed is greater than the Alfv\'en speed. We find a semi-analytical formula that explains this numerical result. The formula also predicts that the core formation time is about three times shorter than that with no turbulence, when the turbulent speed is comparable to the Alfv\'en speed.

There are regions in stars, such as ionization zones and the interface between radiative and convective regions, that cause a localized sharp variation in the sound speed. These are known as "acoustic glitches". Acoustic glitches leave their signatures on the oscillation frequencies of stars, and hence these signature can be used as diagnostics of these regions. In particular, the signature of these glitches can be used as diagnostics of the position of the second helium ionization zone and that of the base of the envelope convection zone. With the help of stellar models we study the properties of these acoustic glitches in main-sequence stars. We find that the acoustic glitch due to the helium ionization zone does not correspond to the dip in the adiabatic index \Gamma_1 caused by the ionization of HeII, but to the peak in \Gamma_1 between the HeI and HeII ionization zones. We find that it is easiest to study the acoustic glitch due to the helium ionization zone in stars with masses in the range 0.9--1.2 M_\odot.

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar-evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5 to 10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc^3 over 5 kpc <~ R <~ 14 kpc, |Z| <~ 2 kpc, and -15 deg <~ Galactocentric azimuth <~ 30 deg. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from 2MASS, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ~45 degree-wide azimuthal region covered by the current sample to be <= 0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kpc from the Sun.

16 Cyg A and B are among the brightest stars observed by Kepler. What makes these stars more interesting is that they are solar analogs. 16 Cyg A and B exhibit solar-like oscillations. In this work we use oscillation frequencies obtained using 2.5 years of Kepler data to determine the current helium abundance of these stars. For this we use the fact that the helium ionization zone leaves a signature on the oscillation frequencies and that this signature can be calibrated to determine the helium abundance of that layer. By calibrating the signature of the helium ionization zone against models of known helium abundance, the helium abundance in the envelope of 16 Cyg A is found to lie in the range 0.231 to 0.251 and that of 16 Cyg B lies in the range 0.218 to 0.266.

We present the characterization of the Kepler-93 exoplanetary system, based on three years of photometry gathered by the Kepler spacecraft. The duration and cadence of the Kepler observations, in tandem with the brightness of the star, enable unusually precise constraints on both the planet and its host. We conduct an asteroseismic analysis of the Kepler photometry and conclude that the star has an average density of 1.652+/-0.006 g/cm^3. Its mass of 0.911+/-0.033 M_Sun renders it one of the lowest-mass subjects of asteroseismic study. An analysis of the transit signature produced by the planet Kepler-93b, which appears with a period of 4.72673978+/-9.7x10^-7 days, returns a consistent but less precise measurement of the stellar density, 1.72+0.02-0.28 g/cm^3. The agreement of these two values lends credence to the planetary interpretation of the transit signal. The achromatic transit depth, as compared between Kepler and the Spitzer Space Telescope, supports the same conclusion. We observed seven transits of Kepler-93b with Spitzer, three of which we conducted in a new observing mode. The pointing strategy we employed to gather this subset of observations halved our uncertainty on the transit radius ratio R_p/R_star. We find, after folding together the stellar radius measurement of 0.919+/-0.011 R_Sun with the transit depth, a best-fit value for the planetary radius of 1.481+/-0.019 R_Earth. The uncertainty of 120 km on our measurement of the planet's size currently renders it one of the most precisely measured planetary radii outside of the Solar System. Together with the radius, the planetary mass of 3.8+/-1.5 M_Earth corresponds to a rocky density of 6.3+/-2.6 g/cm^3. After applying a prior on the plausible maximum densities of similarly-sized worlds between 1--1.5 R_Earth, we find that Kepler-93b possesses an average density within this group.

Mixing length theory is the predominant treatment of convection in stellar models today. Usually described by a single free parameter, alpha, the common practice is to calibrate it using the properties of the Sun, and apply it to all other stellar models as well. Asteroseismic data from Kepler and CoRoT provide precise properties of other stars which can be used to determine alpha as well, and a recent study of stars in the Kepler field of view found alpha to vary with metallicity. Interpreting alpha obtained from calibrated stellar models, however, is complicated by the fact that the value for alpha depends on the surface boundary condition of the stellar model, or T-tau relation. Calibrated models that use typical T-tau relations, which are static and insensitive to chemical composition, do not include the complete effect of metallicity on alpha. We use 3D radiation-hydrodynamic simulations to extract metallicity-dependent T-tau relations and use them in calibrated stellar models. We find the previously reported alpha-metallicity trend to be robust, and not significantly affected by the surface boundary condition of the stellar models.

Fundamental stellar properties, such as mass, radius, and age, can be inferred using asteroseismology. Cool stars with convective envelopes have turbulent motions that can stochastically drive and damp pulsations. The properties of the oscillation frequency power spectrum can be tied to mass and radius through solar-scaled asteroseismic relations. Stellar properties derived using these scaling relations need verification over a range of metallicities. Because the age and mass of halo stars are well-constrained by astrophysical priors, they provide an independent, empirical check on asteroseismic mass estimates in the low-metallicity regime. We identify nine metal-poor red giants (including six stars that are kinematically associated with the halo) from a sample observed by both the Kepler space telescope and the Sloan Digital Sky Survey-III APOGEE spectroscopic survey. We compare masses inferred using asteroseismology to those expected for halo and thick-disk stars. Although our sample is small, standard scaling relations, combined with asteroseismic parameters from the APOKASC Catalog, produce masses that are systematically higher (<{\Delta}M>=0.17+/-0.05 Msun) than astrophysical expectations. The magnitude of the mass discrepancy is reduced by known theoretical corrections to the measured large frequency separation scaling relationship. Using alternative methods for measuring asteroseismic parameters induces systematic shifts at the 0.04 Msun level. We also compare published asteroseismic analyses with scaling relationship masses to examine the impact of using the frequency of maximum power as a constraint. Upcoming APOKASC observations will provide a larger sample of ~100 metal-poor stars, important for detailed asteroseismic characterization of Galactic stellar populations.

The NASA Kepler mission has been continuously monitoring the same field of the sky since the successful launch in March 2009, providing high-quality stellar lightcurves that are excellent data for asteroseismology, far superior to any other observations available at the present. In order to make a meaningful analysis and interpretation of the asteroseismic data, accurate fundamental parameters for the observed stars are needed. The currently available parameters are quite uncertain as illustrated by e.g. Thygesen et al. (A&A 543, A160, 2012), who found deviations as extreme as 2.0 dex in [Fe/H] and log g, compared to catalogue values. Thus, additional follow-up observations for these targets are needed in order to put firm limits on the parameter space investigated by the asteroseismic modellers. Here, we propose a metod for deriving accurate metallicities of main sequence and subgiant solar-like oscillators from medium resolution spectra with a moderate S/N. The method takes advantage of the additional constraints on the fundamental parameters, available from asteroseismology and multi-color photometry. The approach enables us to reduce the analysis overhead significantly when doing spectral synthesis, which in turn will increases the efficiency of follow-up observations.

We present revised properties for 196,468 stars observed by the NASA Kepler Mission and used in the analysis of Quarter 1-16 (Q1-Q16) data to detect and characterize transiting exoplanets. The catalog is based on a compilation of literature values for atmospheric properties (temperature, surface gravity, and metallicity) derived from different observational techniques (photometry, spectroscopy, asteroseismology, and exoplanet transits), which were then homogeneously fitted to a grid of Dartmouth stellar isochrones. We use broadband photometry and asteroseismology to characterize 11,532 Kepler targets which were previously unclassified in the Kepler Input Catalog (KIC). We report the detection of oscillations in 2,762 of these targets, classifying them as giant stars and increasing the number of known oscillating giant stars observed by Kepler by ~20% to a total of ~15,500 stars. Typical uncertainties in derived radii and masses are ~40% and ~20%, respectively, for stars with photometric constraints only, and 5-15% and ~10% for stars based on spectroscopy and/or asteroseismology, although these uncertainties vary strongly with spectral type and luminosity class. A comparison with the Q1-Q12 catalog shows a systematic decrease in radii for M dwarfs, while radii for K dwarfs decrease or increase depending on the Q1-Q12 provenance (KIC or Yonsei-Yale isochrones). Radii of F-G dwarfs are on average unchanged, with the exception of newly identified giants. The Q1-Q16 star properties catalog is a first step towards an improved characterization of all Kepler targets to support planet occurrence studies.

We present a novel method for estimating lower-limit surface gravities log g of Kepler targets whose data do not allow the detection of solar-like oscillations. The method is tested using an ensemble of solar-type stars observed in the context of the Kepler Asteroseismic Science Consortium. We then proceed to estimate lower-limit log g for a cohort of Kepler solar-type planet-candidate host stars with no detected oscillations. Limits on fundamental stellar properties, as provided by this work, are likely to be useful in the characterization of the corresponding candidate planetary systems. Furthermore, an important byproduct of the current work is the confirmation that amplitudes of solar-like oscillations are suppressed in stars with increased levels of surface magnetic activity.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

For the very best and brightest asteroseismic solar-type targets observed by Kepler, the frequency precision is sufficient to determine the acoustic depths of the surface convective layer and the helium ionization zone. Such sharp features inside the acoustic cavity of the star, which we call acoustic glitches, create small oscillatory deviations from the uniform spacing of frequencies in a sequence of oscillation modes with the same spherical harmonic degree. We use these oscillatory signals to determine the acoustic locations of such features in 19 solar-type stars observed by the Kepler mission. Four independent groups of researchers utilized the oscillation frequencies themselves, the second differences of the frequencies and the ratio of the small and large separation to locate the base of the convection zone and the second helium ionization zone. Despite the significantly different methods of analysis, good agreement was found between the results of these four groups, barring a few cases. These results also agree reasonably well with the locations of these layers in representative models of the stars. These results firmly establish the presence of the oscillatory signals in the asteroseismic data and the viability of several techniques to determine the location of acoustic glitches inside stars.

We model molecular cloud fragmentation with thin disk non-ideal magnetohydrodynamic simulations that include ambipolar diffusion and partial ionization that transitions from primarily ultraviolet dominated to cosmic ray dominated regimes. These simulations are used to determine the conditions required for star clusters to form through a two-stage fragmentation scenario. Recent linear analyses have shown that the fragmentation length and time scales can undergo a dramatic drop across the column density boundary that separates the ultraviolet and cosmic ray dominated ionization regimes. As found in earlier studies, the absence of an ionization drop and regular perturbations leads to a single-stage fragmentation on parsec scales in transcritical clouds, so that the nonlinear evolution yields the same fragment sizes as predicted by linear theory. However, we find that a combination of initial transcritical mass-to-flux ratio, evolution through a column density regime in which the ionization drop takes place, and regular small perturbations to the mass-to-flux ratio are sufficient to cause a second stage of fragmentation during the nonlinear evolution. Cores of size ~0.1 pc are formed within an initial fragment of ~ pc size. Regular perturbations to the mass-to-flux ratio also accelerate the onset of runaway collapse.

Stars hosting hot Jupiters are often observed to have high obliquities, whereas stars with multiple co-planar planets have been seen to have low obliquities. This has been interpreted as evidence that hot-Jupiter formation is linked to dynamical disruption, as opposed to planet migration through a protoplanetary disk. We used asteroseismology to measure a large obliquity for Kepler-56, a red giant star hosting two transiting co-planar planets. These observations show that spin-orbit misalignments are not confined to hot-Jupiter systems. Misalignments in a broader class of systems had been predicted as a consequence of torques from wide-orbiting companions, and indeed radial-velocity measurements revealed a third companion in a wide orbit in the Kepler-56 system.

We use asteroseismic data obtained by the NASA Kepler Mission to estimate the fundamental properties of more than 500 main-sequence and sub-giant stars. Data obtained during the first 10 months of Kepler science operations were used for this work, when these solar-type targets were observed for one month each in a survey mode. Stellar properties have been estimated using two global asteroseismic parameters and complementary photometric and spectroscopic data. Homogeneous sets of effective temperatures were available for the entire ensemble from complementary photometry; spectroscopic estimates of T_eff and [Fe/H] were available from a homogeneous analysis of ground-based data on a subset of 87 stars. [Abbreviated version... see paper for full abstract.]

We investigate the effect of helium abundance and $\alpha$-element enhancement on the properties of convection in envelopes of solar-like main-sequence stars stars using a grid of 3D radiation hydrodynamic simulations. Helium abundance increases the mean molecular weight of the gas, and alters opacity by displacing hydrogen. Since the scale of the effect of helium may depend on the metallicity, the grid consists of simulations with three helium abundances ($Y=0.1, 0.2, 0.3$), each with two metallicities ($Z=0.001, 0.020)$. We find that changing the helium mass fraction generally affects structure and convective dynamics in a way opposite to that of metallicity. Furthermore, the effect is considerably smaller than that of metallicity. The signature of helium differs from that of metallicity in the manner in which the photospheric velocity distribution is affected. \rev{We also find that helium abundance and surface gravity behave largely in similar ways, but differ in the way they affect the mean molecular weight}. A simple model for spectral line formation suggests that the bisectors and absolute Doppler shifts of spectral lines depends on the helium abundance. We look at the effect of $\alpha$-element enhancement and find that it has a considerably smaller effect on the convective dynamics in the SAL compared to that of helium abundance.

We show that the termination of the star formation process by winds from massive stars in protocluster forming clumps imposes dual constraints on the star formation efficiencies (SFEs) and stellar age spreads ($\Delta \tau_{*}$) in stellar clusters. We have considered two main classes of clump models. One class of models in one in which the core formation efficiency (CFE) per unit time and as a consequence the star formation rate (SFR) is constant in time and another class of models in which the CFE per unit time, and as a consequence the SFR, increases with time. Models with an increasing mode of star formation yield shorter age spreads (a few 0.1 Myrs) and typically higher SFEs than models in which star formation is uniform in time. We find that the former models reproduce remarkably well the SFE$-\Delta \tau_{*}$ values of starburst clusters such as NGC 3603 YC and Westerlund 1, while the latter describe better the star formation process in lower density environments such as in the Orion Nebula Cluster. We also show that the SFE and $\Delta \tau_{*}$ of massive clusters are expected to be higher in low metallicity environments. This could be tested with future large extragalactic surveys of stellar clusters. We advocate that placing a stellar cluster on the SFE-$\Delta \tau_{*}$ diagram is a powerful method to distinguish between different stellar clusters formation scenarios such as between generic gravitational instability of a gas cloud/clump or as the result of cloud-cloud collisions. It is also a very useful tool for testing star formation theories and numerical models versus the observations.

Supersonic turbulence in the interstellar medium (ISM) is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this general consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to take into account the roles of global structures of magnetic fields and the ISM. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in warm neutral medium. Early evolution of turbulent kinetic energy is consistent with previous results for decaying compressible MHD turbulence characterized by rapid energy decay with a power-law form of $E\propto t^{-1}$ and by short decay time compared to a flow crossing time. If initial magnetic fields are strong and perpendicular to the sheet, however, long-term evolution of the kinetic energy shows that a significant amount of turbulent energy ($\sim 0.2E_0$) still remains even after ten flow crossing times for models with periodic boundary conditions. The decay rate is also greatly reduced as the field strength increases for such initial and boundary conditions, but not if the boundary conditions are that for a completely isolated sheet. We analyze velocity power spectra of the remaining turbulence to show that in-plane, incompressible motions parallel to the sheet dominate at later times.

We comment on the potential for continuing asteroseismology of solar-type and red-giant stars in a 2-wheel Kepler Mission. Our main conclusion is that by targeting stars in the ecliptic it should be possible to perform high-quality asteroseismology, as long as favorable scenarios for 2-wheel pointing performance are met. Targeting the ecliptic would potentially facilitate unique science that was not possible in the nominal Mission, notably from the study of clusters that are significantly brighter than those in the Kepler field. Our conclusions are based on predictions of 2-wheel observations made by a space photometry simulator, with information provided by the Kepler Project used as input to describe the degraded pointing scenarios. We find that elevated levels of frequency-dependent noise, consistent with the above scenarios, would have a significant negative impact on our ability to continue asteroseismic studies of solar-like oscillators in the Kepler field. However, the situation may be much more optimistic for observations in the ecliptic, provided that pointing resets of the spacecraft during regular desaturations of the two functioning reaction wheels are accurate at the < 1 arcsec level. This would make it possible to apply a post-hoc analysis that would recover most of the lost photometric precision. Without this post-hoc correction---and the accurate re-pointing it requires---the performance would probably be as poor as in the Kepler-field case. Critical to our conclusions for both fields is the assumed level of pointing noise (in the short-term jitter and the longer-term drift). We suggest that further tests will be needed to clarify our results once more detail and data on the expected pointing performance becomes available, and we offer our assistance in this work.

The SDSS-III Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a three year survey that is collecting 100,000 high-resolution spectra in the near-IR across multiple Galactic populations. To derive stellar parameters and chemical compositions from this massive data set, the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP) has been developed. Here, we describe empirical calibrations of stellar parameters presented in the first SDSS-III APOGEE data release (DR10). These calibrations were enabled by observations of 559 stars in 20 globular and open clusters. The cluster observations were supplemented by observations of stars in NASA's Kepler field that have well determined surface gravities from asteroseismic analysis. We discuss the accuracy and precision of the derived stellar parameters, considering especially effective temperature, surface gravity, and metallicity; we also briefly discuss the derived results for the abundances of the alpha-elements, carbon, and nitrogen. Overall, we find that ASPCAP achieves reasonably accurate results for temperature and metallicity, but suffers from systematic errors in surface gravity. We derive calibration relations that bring the raw ASPCAP results into better agreement with independently determined stellar parameters. The internal scatter of ASPCAP parameters within clusters suggests that, metallicities are measured with a precision better than 0.1 dex, effective temperatures better than 150 K, and surface gravities better than 0.2 dex. The understanding provided by the clusters and Kepler giants on the current accuracy and precision will be invaluable for future improvements of the pipeline.

Asteroseismology, i.e. the study of the internal structures of stars via their global oscillations, is a valuable tool to obtain stellar parameters such as mass, radius, surface gravity and mean density. These parameters can be obtained using certain scaling relations which are based on an asymptotic approximation. Usually the observed oscillation parameters are assumed to follow these scaling relations. Recently, it has been questioned whether this is a valid approach, i.e., whether the order of the observed oscillation modes are high enough to be approximated with an asymptotic theory. In this work we use stellar models to investigate whether the differences between observable oscillation parameters and their asymptotic estimates are indeed significant. We compute the asymptotic values directly from the stellar models and derive the observable values from adiabatic pulsation calculations of the same models. We find that the extent to which the atmosphere is included in the models is a key parameter. Considering a larger extension of the atmosphere beyond the photosphere reduces the difference between the asymptotic and observable values of the large frequency separation. Therefore, we conclude that the currently suggested discrepancies in the scaling relations might have been overestimated. Hence, based on the results presented here we believe that the suggestions of Mosser et al. (2013) should not be followed without careful consideration.

Context: Asteroseismic surface gravity values can be of importance in determining spectroscopic stellar parameters. The independent log(g) value from asteroseismology can be used as a fixed value in the spectroscopic analysis to reduce uncertainties due to the fact that log(g) and effective temperature can not be determined independently from spectra. Since 2012, a combined analysis of seismically and spectroscopically derived stellar properties is ongoing for a large survey with SDSS/APOGEE and Kepler. Therefore, knowledge of any potential biases and uncertainties in asteroseismic log(g) values is now becoming important. Aims: The seismic parameter needed to derive log(g) is the frequency of maximum oscillation power (nu_max). Here, we investigate the influence of nu_max derived with different methods on the derived log(g) values. The large frequency separation between modes of the same degree and consecutive radial orders (Dnu) is often used as an additional constraint for the determination of log(g). Additionally, we checked the influence of small corrections applied to Dnu on the derived values of log(g). Methods We use methods extensively described in the literature to determine nu_max and Dnu together with seismic scaling relations and grid-based modeling to derive log(g). Results: We find that different approaches to derive oscillation parameters give results for log(g) with small, but different, biases for red-clump and red-giant-branch stars. These biases are well within the quoted uncertainties of ~0.01 dex (cgs). Corrections suggested in the literature to the Dnu scaling relation have no significant effect on log(g). However somewhat unexpectedly, method specific solar reference values induce biases of the order of the uncertainties, which is not the case when canonical solar reference values are used.

Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.

Asteroseismic data can be used to determine surface gravities with precisions of < 0.05 dex by using the global seismic quantities Deltanu and nu_max along with Teff and [Fe/H]. Surface gravity is also one of the four stellar properties to be derived by automatic analyses for 1 billion stars from Gaia data (workpackage GSP_Phot). We explore seismic data from MS F, G, K stars (solar-like stars) observed by Kepler as a potential calibration source for methods that Gaia will use for object characterisation (log g). We calculate log g for bright nearby stars for which radii and masses are known, and using their global seismic quantities in a grid-based method, we determine an asteroseismic log g to within 0.01 dex of the direct calculation, thus validating the accuracy of our method. We find that errors in Teff and mainly [Fe/H] can cause systematic errors of 0.02 dex. We then apply our method to a list of 40 stars to deliver precise values of surface gravity, i.e. sigma < 0.02 dex, and we find agreement with recent literature values. Finally, we explore the precision we expect in a sample of 400+ Kepler stars which have their global seismic quantities measured. We find a mean uncertainty (precision) on the order of <0.02 dex in log g over the full explored range 3.8 < log g < 4.6, with the mean value varying only with stellar magnitude (0.01 - 0.02 dex). We study sources of systematic errors in log g and find possible biases on the order of 0.04 dex, independent of log g and magnitude, which accounts for errors in the Teff and [Fe/H] measurements, as well as from using a different grid-based method. We conclude that Kepler stars provide a wealth of reliable information that can help to calibrate methods that Gaia will use, in particular, for source characterisation with GSP_Phot where excellent precision (small uncertainties) and accuracy in log g is obtained from seismic data.

Using asteroseismic data and stellar evolution models we make the first detection of a convective core in a Kepler field main-sequence star, putting a stringent constraint on the total size of the mixed zone and showing that extra mixing beyond the formal convective boundary exists. In a slightly less massive target the presence of a convective core cannot be conclusively discarded, and thus its remaining main-sequence life time is uncertain. Our results reveal that best-fit models found solely by matching individual frequencies of oscillations corrected for surface effects do not always properly reproduce frequency combinations. Moreover, slightly different criteria to define what the best-fit model is can lead to solutions with similar global properties but very different interior structures. We argue that the use of frequency ratios is a more reliable way to obtain accurate stellar parameters, and show that our analysis in field main-sequence stars can yield an overall precision of 1.5%, 4%, and 10% in radius, mass and age, respectively. We compare our results with those obtained from global oscillation properties, and discuss the possible sources of uncertainties in asteroseismic stellar modeling where further studies are still needed.

We review recent results of non-ideal magnetohydrodynamic models for the fragmentation of molecular clouds and the collapse of cloud cores to form protostar-disk systems. Thin-disk models can elucidate many aspects of the physical problem and allow the calculation of large dynamic range of time and length scales.

We have used asteroseismology to determine fundamental properties for 66 Kepler planet-candidate host stars, with typical uncertainties of 3% and 7% in radius and mass, respectively. The results include new asteroseismic solutions for four host stars with confirmed planets (Kepler-4, Kepler-14, Kepler-23 and Kepler-25) and increase the total number of Kepler host stars with asteroseismic solutions to 77. A comparison with stellar properties in the planet-candidate catalog by Batalha et al. shows that radii for subgiants and giants obtained from spectroscopic follow-up are systematically too low by up to a factor of 1.5, while the properties for unevolved stars are in good agreement. We furthermore apply asteroseismology to confirm that a large majority of cool main-sequence hosts are indeed dwarfs and not misclassified giants. Using the revised stellar properties, we recalculate the radii for 107 planet candidates in our sample, and comment on candidates for which the radii change from a previously giant-planet/brown-dwarf/stellar regime to a sub-Jupiter size, or vice versa. A comparison of stellar densities from asteroseismology with densities derived from transit models in Batalha et al. assuming circular orbits shows significant disagreement for more than half of the sample due to systematics in the modeled impact parameters, or due to planet candidates which may be in eccentric orbits. Finally, we investigate tentative correlations between host-star masses and planet candidate radii, orbital periods, and multiplicity, but caution that these results may be influenced by the small sample size and detection biases.

We study the formation and long-term evolution of primordial protostellar disks harbored by first stars using numerical hydrodynamics simulations in the thin-disk limit. The initial conditions are specified by pre-stellar cores with distinct mass, angular momentum, and temperature. This allows us to probe several tens of thousand years of the disk's initial evolution, during which we observe multiple episodes of fragmentation leading to the formation of gravitationally bound gaseous clumps within spiral arms. These fragments are torqued inward due to gravitational interaction with the spiral arms on timescales of 10^3 - 10^4 yr and accreted onto the growing protostar, giving rise to accretion and luminosity bursts. The burst phenomenon is fueled by continuing accretion of material falling onto the disk from the collapsing parent core, which replenishes the mass lost by the disk due to accretion, and triggers repetitive episodes of disk fragmentation. We show that the burst phenomenon is expected to occur for a wide spectrum of initial conditions in primordial pre-stellar cores and speculate on how the intense luminosities (~10^7 solar luminosities) produced by this mechanism may have important consequences for the disk evolution and subsequent growth of the protostar.

Starting with a mathematical boundary value problem for the magnetic vector potential in an axisymmetric cylindrical coordinate system, we derive a general solution for any arbitrary current distribution using the method of Green's functions. We use this to derive an analytic form for an hourglass magnetic field pattern created by electrical currents that are concentrated near (but not confined within) the equatorial plane of a cylindrical coordinate system. Our solution is not characterized by a cusp at the equatorial plane, as in previous solutions based on a current sheet. The pattern we derive provides a very good fit to hourglass magnetic field patterns emerging from three-dimensional numerical simulations of core formation, and can in principle be used for source-fitting of observed magnetic hourglass patterns.

We examine how metallicity affects convection and overshoot in the superadiabatic layer of main sequence stars. We present results from a grid of 3D radiation hydrodynamic simulations with four metallicities ($Z=0.040$, 0.020, 0.010, 0.001), and spanning a range in effective temperature ($4950 < \rm T_{\rm eff} < 6230$). We show that changing the metallicity alters properties of the convective gas dynamics, and the structure of the superadiabatic layer and atmosphere. Our grid of simulations show that the amount of superadiabaticity, which tracks the transition from efficient to inefficient convection, \rev{is sensitive to changes in metallicity. We find that increasing the metallicity forces the location of the transition region to lower densities and pressures, and results in larger mean and turbulent velocities throughout the superadiabatic region.} We also quantify the degree of convective overshoot in the atmosphere, and show that it increases with metallicity as well.

Results on the obliquity of exoplanet host stars -- the angle between the stellar spin axis and the planetary orbital axis -- provide important diagnostic information for theories describing planetary formation. Here we present the first application of asteroseismology to the problem of stellar obliquity determination in systems with transiting planets and Sun-like host stars. We consider two systems observed by the NASA Kepler Mission which have multiple transiting small (super-Earth sized) planets: the previously reported Kepler-50 and a new system, Kepler-65, whose planets we validate in this paper. Both stars show rich spectra of solar-like oscillations. From the asteroseismic analysis we find that each host has its rotation axis nearly perpendicular to the line of sight with the sines of the angles constrained at the 1-sigma level to lie above 0.97 and 0.91, respectively. We use statistical arguments to show that coplanar orbits are favoured in both systems, and that the orientations of the planetary orbits and the stellar rotation axis are correlated.

NASA's Kepler Mission has revealed two transiting planets orbiting Kepler-68. Follow-up Doppler measurements have established the mass of the innermost planet and revealed a third jovian-mass planet orbiting beyond the two transiting planets. Kepler-68b, in a 5.4 day orbit has mass 8.3 +/- 2.3 Earth, radius 2.31 +/- 0.07 Earth radii, and a density of 3.32 +/- 0.92 (cgs), giving Kepler-68b a density intermediate between that of the ice giants and Earth. Kepler-68c is Earth-sized with a radius of 0.953 Earth and transits on a 9.6 day orbit; validation of Kepler-68c posed unique challenges. Kepler-68d has an orbital period of 580 +/- 15 days and minimum mass of Msin(i) = 0.947 Jupiter. Power spectra of the Kepler photometry at 1-minute cadence exhibit a rich and strong set of asteroseismic pulsation modes enabling detailed analysis of the stellar interior. Spectroscopy of the star coupled with asteroseismic modeling of the multiple pulsation modes yield precise measurements of stellar properties, notably Teff = 5793 +/- 74 K, M = 1.079 +/- 0.051 Msun, R = 1.243 +/- 0.019 Rsun, and density 0.7903 +/- 0.0054 (cgs), all measured with fractional uncertainties of only a few percent. Models of Kepler-68b suggest it is likely composed of rock and water, or has a H and He envelope to yield its density of about 3 (cgs).

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds yields information about length scales involved in star formation. Combining these length scales with various distributions of other environmental variables, (i.e., column density and mass-to-flux ratio) and applying Monte Carlo methods allow us to produce synthetic core mass functions (CMFs) for different environmental conditions. Our analysis shows that the shape of the CMF is directly dependent on the physical conditions of the cloud. Specifically, magnetic fields act to broaden the mass function and develop a high-mass tail while ambipolar diffusion will truncate this high-mass tail. In addition, we analyze the effect of small number statistics on the shape and high-mass slope of the synthetic CMFs. We find that observed core mass functions are severely statistically limited, which has a profound effect on the derived slope for the high-mass tail.

A new view of disk evolution is emerging from self-consistent numerical simulation modeling of the formation of circumstellar disks from the direct collapse of prestellar cloud cores. This has implications for many aspects of star and planet formation, including the growth of dust and high-temperature processing of materials. A defining result is that the early evolution of a disk is crucially affected by the continuing mass loading from the core envelope, and is driven into recurrent phases of gravitational instability. Nonlinear spiral arms formed during these episodes fragment to form gaseous clumps in the disk. These clumps generally migrate inward due to gravitational torques arising from their interaction with a trailing spiral arm. Occasionally, a clump can open up a gap in the disk and settle into a stable orbit, revealing a direct pathway to the formation of companion stars, brown dwarfs, or giant planets. At other times, when multiple clumps are present, a low mass clump may even be ejected from the system, providing a pathway to the formation of free-floating brown dwarfs and giant planets in addition to low mass stars. Finally, it has been suggested that the inward migration of gaseous clumps can provide the proper conditions for the transport of high-temperature processed solids from the outer disk to the inner disk, and even possibly accelerate the formation of terrestrial planets in the inner disk. All of these features arising from clump formation and migration can be tied together conceptually in a Migrating Embryo model for disk evolution that can complement the well-known Core Accretion model for planet formation.

The active K2 dwarf epsilon Eri has been extensively characterized, both as a young solar analog and more recently as an exoplanet host star. As one of the nearest and brightest stars in the sky, it provides an unparalleled opportunity to constrain stellar dynamo theory beyond the Sun. We confirm and document the 3 year magnetic activity cycle in epsilon Eri originally reported by Hatzes and coworkers, and we examine the archival data from previous observations spanning 45 years. The data show coexisting 3 year and 13 year periods leading into a broad activity minimum that resembles a Maunder minimum-like state, followed by the resurgence of a coherent 3 year cycle. The nearly continuous activity record suggests the simultaneous operation of two stellar dynamos with cycle periods of 2.95+/-0.03 years and 12.7+/-0.3 years, which by analogy with the solar case suggests a revised identification of the dynamo mechanisms that are responsible for the so-called "active" and "inactive" sequences as proposed by Bohm-Vitense. Finally, based on the observed properties of epsilon Eri we argue that the rotational history of the Sun is what makes it an outlier in the context of magnetic cycles observed in other stars (as also suggested by its Li depletion), and that a Jovian-mass companion cannot be the universal explanation for the solar peculiarities.

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds shows that molecular clouds exhibit a preferred length scale for collapse that depends on the mass-to-flux ratio and neutral-ion collision time within the cloud. We extend this linear analysis to the context of clustered star formation. By combining the results of the linear analysis with a realistic ionization profile for the cloud, we find that a molecular cloud may evolve through two fragmentation events in the evolution toward the formation of stars. Our model suggests that the initial fragmentation into clumps occurs for a transcritical cloud on parsec scales while the second fragmentation can occur for transcritical and supercritical cores on subparsec scales. Comparison of our results with several star forming regions (Perseus, Taurus, Pipe Nebula) shows support for a two-stage fragmentation model.

Asteroseismology has been extremely successful in determining the properties of stars in different evolutionary stages with a remarkable level of precision. However, to fully exploit its potential, robust methods for estimating stellar parameters are required and independent verification of the results is needed. In this talk, I present a new technique developed to obtain stellar properties by coupling asteroseismic analysis with the InfraRed Flux Method. Using two global seismic observables and multi-band photometry, the technique determines masses, radii, effective temperatures, bolometric fluxes, and thus distances for field stars in a self-consistent manner. Applying our method to a sample of solar-like oscillators in the {\it Kepler} field that have accurate {\it Hipparcos} parallaxes, we find agreement in our distance determinations to better than 5%. Comparison with measurements of spectroscopic effective temperatures and interferometric radii also validate our results, and show that our technique can be applied to stars evolved beyond the main-sequence phase.

We employ ring-diagram analysis to study the sub-surface thermal structure of active regions. We present results using a large number of active regions over the course of Solar Cycle 23. We present both traditional inversions of ring-diagram frequency differences, with a total sample size of 264, and a statistical study using Principal Component Analysis. We confirm earlier results on smaller samples that sound speed and adiabatic index are changed below regions of strong magnetic field. We find that sound speed is decreased in the region between approximately r=0.99R_sun and r=0.995R_sun (depths of 3Mm to 7Mm), and increased in the region between r=0.97R_sun and r=0.985R_sun (depths of 11Mm to 21Mm). The adiabatic index is enhanced in the same deeper layers that sound-speed enhancement is seen. A weak decrease in adiabatic index is seen in the shallower layers in many active regions. We find that the magnitudes of these perturbations depend on the strength of the surface magnetic field, but we find a great deal of scatter in this relation, implying other factors may be relevant.

We examine the effect of different radiative transfer schemes on the properties of 3D simulations of near-surface stellar convection in the superadiabatic layer, where energy transport transitions from fully convective to fully radiative. We employ two radiative transfer schemes that fundamentally differ in the way they cover the 3D domain. The first solver approximates domain coverage with moments, while the second solver samples the 3D domain with ray integrations. By comparing simulations that differ only in their respective radiative transfer methods, we are able to isolate the effect that radiative efficiency has on the structure of the superadiabatic layer. We find the simulations to be in good general agreement, but they show distinct differences in the thermal structure in the superadiabatic layer and atmosphere.

Accurately determining the properties of stars is of prime importance for characterizing stellar populations in our Galaxy. The field of asteroseismology has been thought to be particularly successful in such an endeavor for stars in different evolutionary stages. However, to fully exploit its potential, robust methods for estimating stellar parameters are required and independent verification of the results is mandatory. With this purpose, we present a new technique to obtain stellar properties by coupling asteroseismic analysis with the InfraRed Flux Method. By using two global seismic observables and multi-band photometry, the technique allows us to obtain masses, radii, effective temperatures, bolometric fluxes, and hence distances for field stars in a self-consistent manner. We apply our method to 22 solar-like oscillators in the Kepler short-cadence sample, that have accurate Hipparcos parallaxes. Our distance determinations agree to better than 5%, while measurements of spectroscopic effective temperatures and interferometric radii also validate our results. We briefly discuss the potential of our technique for stellar population analysis and models of Galactic Chemical Evolution.

The solar magnetic activity cycle causes changes in the Sun on timescales that are relevant to human lifetimes. The minimum in solar activity that preceded the current solar cycle (cycle 24) was deeper and quieter than any other recent minimum. Using data from the Birmingham Solar-Oscillations Network (BiSON), we show that the structure of the solar sub-surface layers during the descending phase of the preceding cycle (cycle 23) was very different from that during cycle 22. This leads us to believe that a detailed examination of the data would have led to the prediction that the cycle-24 minimum would be out of the ordinary. The behavior of the oscillation frequencies allows us to infer that changes in the Sun that affected the oscillation frequencies in cycle 23 were localized mainly to layers above about 0.996Rsun, depths shallower than about 3000 km. In cycle 22, on the other hand, the changes must have also occurred in the deeper-lying layers.

Modeling of the self-consistent formation and evolution of disks as a result of prestellar core collapse reveals an intense early phase of recurrent gravitational instability and clump formation. These clumps generally migrate inward due to gravitational interaction with trailing spiral arms, and can be absorbed into the central object. However, in situations of multiple clump formation, gravitational scattering of clumps can result in the ejection of a low mass clump. These clumps can then give rise to free-floating low mass stars, brown dwarfs, or even giant planets. Detailed modeling of this process in the context of present-day star formation reveals that these clumps start out essentially as Larson first cores and grow subsequently by accretion. In the context of Pop III star formation, preliminary indications are that the disk clumps may also be of low mass. This mechanism of clump formation and possible ejection provides a channel for the formation of low mass objects in the first generation of stars.

We studied solar-like oscillations in 115 red giants in the three open clusters NGC 6791, NGC 6811, and NGC 6819, based on photometric data covering more than 19 months with NASA's Kepler space telescope. We present the asteroseismic diagrams of the asymptotic parameters \delta\nu_02, \delta\nu_01 and \epsilon, which show clear correlation with fundamental stellar parameters such as mass and radius. When the stellar populations from the clusters are compared, we see evidence for a difference in mass of the red giant branch stars, and possibly a difference in structure of the red clump stars, from our measurements of the small separations \delta\nu_02 and \delta\nu_01. Ensemble \'{e}chelle diagrams and upper limits to the linewidths of l = 0 modes as a function of \Delta\nu of the clusters NGC 6791 and NGC 6819 are also shown, together with the correlation between the l = 0 ridge width and the T_eff of the stars. Lastly, we distinguish between red giant branch and red clump stars through the measurement of the period spacing of mixed dipole modes in 53 stars among all the three clusters to verify the stellar classification from the color-magnitude diagram. These seismic results also allow us to identify a number of special cases, including evolved blue stragglers and binaries, as well as stars in late He-core burning phases, which can be potentially interesting targets for detailed theoretical modeling.

We present simulation results for the formation and long-term evolution of a primordial protostellar disk harbored by a first star. Using a 2+1D nonaxisymmetric thin disk numerical simulation, together with a barotropic relation for the gas, we are able to probe ~20 kyr of the disk's evolution. During this time period we observe fragmentation leading to loosely bound gaseous clumps within the disk. These are then torqued inward and accreted onto the growing protostar, giving rise to a burst phenomenon. The luminous feedback produced by this mechanism may have important consequences for the subsequent growth of the protostar.

Stellar models generally use simple parametrizations to treat convection. The most widely used parametrization is the so-called "Mixing Length Theory" where the convective eddy sizes are described using a single number, \alpha, the mixing-length parameter. This is a free parameter, and the general practice is to calibrate \alpha using the known properties of the Sun and apply that to all stars. Using data from NASA's Kepler mission we show that using the solar-calibrated \alpha is not always appropriate, and that in many cases it would lead to estimates of initial helium abundances that are lower than the primordial helium abundance. Kepler data allow us to calibrate \alpha for many other stars and we show that for the sample of stars we have studied, the mixing-length parameter is generally lower than the solar value. We studied the correlation between \alpha and stellar properties, and we find that \alpha increases with metallicity. We therefore conclude that results obtained by fitting stellar models or by using population-synthesis models constructed with solar values of \alpha are likely to have large systematic errors. Our results also confirm theoretical expectations that the mixing-length parameter should vary with stellar properties.

In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.

Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently were restricted to the Sun. In this paper, we report the detection of mixed modes - i.e. modes that behave both as g modes in the core and as p modes in the envelope - in the spectrum of the early red giant KIC7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last decades.

Solar-like oscillations have been observed by Kepler and CoRoT in several solar-type stars, thereby providing a way to probe the stars using asteroseismology. We provide the mode frequencies of the oscillations of various stars required to perform a comparison with those obtained from stellar modelling. We used a time series of nine months of data for each star. The 61 stars observed were categorised in three groups: simple, F-like and mixed-mode. The simple group includes stars for which the identification of the mode degree is obvious. The F-like group includes stars for which the identification of the degree is ambiguous. The mixed-mode group includes evolved stars for which the modes do not follow the asymptotic relation of low-degree frequencies. Following this categorisation, the power spectra of the 61 main sequence and subgiant stars were analysed using both maximum likelihood estimators and Bayesian estimators, providing individual mode characteristics such as frequencies, linewidths, and mode heights. We developed and describe a methodology for extracting a single set of mode frequencies from multiple sets derived by different methods and individual scientists. We report on how one can assess the quality of the fitted parameters using the likelihood ratio test and the posterior probabilities. We provide the mode frequencies of 61 stars (with their 1-sigma error bars), as well as their associated echelle diagrams.

Oscillations in solar-like oscillators tend to follow an approximately regular pattern in which oscillation modes of a certain degree and consecutive order appear at regular intervals in frequency, i.e. the so-called large frequency separation. This is true to first order approximation for acoustic modes. However, to a second order approximation it is evident that the large frequency separation changes as a function of frequency. This frequency dependence has been seen in the Sun and in other main-sequence stars. However, from observations of giant stars, this effect seemed to be less pronounced. We investigate the difference in frequency dependence of the large frequency separation between main-sequence and giant stars using YREC evolutionary models.

Context: The Kepler space mission is reaching continuous observing times long enough to start studying the fine structure of the observed p-mode spectra. Aims: In this paper, we aim to study the signature of stellar evolution on the radial and p-dominated l=2 modes in an ensemble of red giants that show solar-type oscillations. Results: We find that the phase shift of the central radial mode (eps_c) is significantly different for red giants at a given large frequency separation (Dnu_c) but which burn only H in a shell (RGB) than those that have already ignited core He burning. Even though not directly probing the stellar core the pair of local seismic observables (Dnu_c, eps_c) can be used as an evolutionary stage discriminator that turned out to be as reliable as the period spacing of the mixed dipole modes. We find a tight correlation between eps_c and Dnu_c for RGB stars and no indication that eps_c depends on other properties of these stars. It appears that the difference in eps_c between the two populations becomes if we use an average of several radial orders, instead of a local, i.e. only around the central radial mode, Dnu to determine the phase shift. This indicates that the information on the evolutionary stage is encoded locally, in the shape of the radial mode sequence. This shape turns out to be approximately symmetric around the central radial mode for RGB stars but asymmetric for core He burning stars. We computed radial modes for a sequence of RG models and find them to qualitatively confirm our findings. We also find that, at least in our models, the local Dnu is an at least as good and mostly better proxy for both the asymptotic spacing and the large separation scaled from the model density than the average Dnu. Finally, we investigate the signature of the evolutionary stage on the small frequency separation and quantify the mass dependency of this seismic parameter.

Context: Ideal MHD simulations have revealed catastrophic magnetic braking (MB) in the protostellar phase, which prevents the formation of a centrifugal disk around a nascent protostar. Aims: We determine if non-ideal MHD, including the effects of ambipolar diffusion and Ohmic dissipation determined from a detailed chemical network model, allows for disk formation at the earliest stages of star formation (SF). Methods: We employ the axisymmetric thin-disk approximation in order to resolve a dynamic range of 9 orders of magnitude in length and 16 in density, while also calculating partial ionization using up to 19 species in a detailed chemical equilibrium model. MB is applied using a steady-state approximation, and a barotropic relation is used to capture the thermal evolution. Results: We resolve the formation of the first and second cores, with expansion waves at the periphery of each, a magnetic diffusion shock, and prestellar infall profiles at larger radii. Power-law profiles in each region can be understood analytically. After the formation of the second core, centrifugal support rises rapidly and a low-mass disk of radius ~10 R_Sun is formed, when the second core has mass ~0.001 M_Sun. The mass-to-flux ratio is ~10,000 times the critical value in the central region. Conclusions: A small centrifugal disk can form in the earliest stage of SF, due to a shut-off of MB caused by magnetic field dissipation in the first core region. There is enough angular momentum loss to allow the second collapse to occur directly, and a low-mass stellar core to form with a surrounding disk. The disk mass and size will depend upon how the angular momentum transport mechanisms within the disk can keep up with mass infall onto the disk. We estimate that the disk will remain <~10 AU, undetectable even by ALMA, in the early Class 0 phase.

We present a detailed spectroscopic study of 93 solar-type stars that are targets of the NASA/Kepler mission and provide detailed chemical composition of each target. We find that the overall metallicity is well-represented by Fe lines. Relative abundances of light elements (CNO) and alpha-elements are generally higher for low-metallicity stars. Our spectroscopic analysis benefits from the accurately measured surface gravity from the asteroseismic analysis of the Kepler light curves. The log g parameter is known to better than 0.03 dex and is held fixed in the analysis. We compare our Teff determination with a recent colour calibration of V-K (TYCHO V magnitude minus 2MASS Ks magnitude) and find very good agreement and a scatter of only 80 K, showing that for other nearby Kepler targets this index can be used. The asteroseismic log g values agree very well with the classical determination using Fe1-Fe2 balance, although we find a small systematic offset of 0.08 dex (asteroseismic log g values are lower). The abundance patterns of metals, alpha elements, and the light elements (CNO) show that a simple scaling by [Fe/H] is adequate to represent the metallicity of the stars, except for the stars with metallicity below -0.3, where alpha-enhancement becomes important. However, this is only important for a very small fraction of the Kepler sample. We therefore recommend that a simple scaling with [Fe/H] be employed in the asteroseismic analyses of large ensembles of solar-type stars.

We present a calculation of protostellar disk formation and evolution in which gaseous clumps (essentially, the first Larson cores formed via disk fragmentation) are ejected from the disk during the early stage of evolution. This is a universal process related to the phenomenon of ejection in multiple systems of point masses. However, it occurs in our model entirely due to the interaction of compact, gravitationally-bound gaseous clumps and is free from the smoothing-length uncertainty that is characteristic of models using sink particles. Clumps that survive ejection span a mass range of 0.08--0.35 $M_\odot$, and have ejection velocities $0.8 \pm 0.35$ km s$^{-1}$, which are several times greater than the escape speed. We suggest that, upon contraction, these clumps can form substellar or low-mass stellar objects with notable disks, or even close-separation very-low-mass binaries. In this hybrid scenario, allowing for ejection of clumps rather than finished protostars/proto--brown-dwarfs, disk formation and the low velocity dispersion of low-mass objects are naturally explained, while it is also consistent with the observation of isolated low-mass clumps that are ejection products. We conclude that clump ejection and the formation of isolated low mass stellar and substellar objects is a common occurrence, with important implications for understanding the initial mass function, the brown dwarf desert, and the formation of stars in all environments and epochs.

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the stars are among the brightest in the Kepler field. We present an analysis of three months of nearly uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 46 and 41 oscillation frequencies for the two components respectively, including a clear detection of octupole (l=3) modes in both stars. We derive the properties of each star independently using the Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by fitting each component individually yield a common age (t=6.8+/-0.4 Gyr) and initial composition (Z_i=0.024+/-0.002, Y_i=0.25+/-0.01) within the uncertainties, as expected for the components of a binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data sets that will ultimately become available will allow future studies of differential rotation, convection zone depths, and long-term changes due to stellar activity cycles.

Asteroseismology with the Kepler space telescope is providing not only an improved characterization of exoplanets and their host stars, but also a new window on stellar structure and evolution for the large sample of solar-type stars in the field. We perform a uniform analysis of 22 of the brightest asteroseismic targets with the highest signal-to-noise ratio observed for 1 month each during the first year of the mission, and we quantify the precision and relative accuracy of asteroseismic determinations of the stellar radius, mass, and age that are possible using various methods. We present the properties of each star in the sample derived from an automated analysis of the individual oscillation frequencies and other observational constraints using the Asteroseismic Modeling Portal (AMP), and we compare them to the results of model-grid-based methods that fit the global oscillation properties. We find that fitting the individual frequencies typically yields asteroseismic radii and masses to \sim1% precision, and ages to \sim2.5% precision (respectively 2, 5, and 8 times better than fitting the global oscillation properties). The absolute level of agreement between the results from different approaches is also encouraging, with model-grid-based methods yielding slightly smaller estimates of the radius and mass and slightly older values for the stellar age relative to AMP, which computes a large number of dedicated models for each star. The sample of targets for which this type of analysis is possible will grow as longer data sets are obtained during the remainder of the mission.

A non-ideal MHD collapse calculation employing the axisymmetric thin-disk approximation is used to resolve cloud core collapse down to the scales of the second (stellar) core. Rotation and a magnetic braking torque are included in the model, and the partial ionization resulting in ambipolar diffusion and Ohmic dissipation is calculated from a detailed chemical network. We find that a centrifugal disk can indeed form in the earliest stage of star formation, due to a shut-off of magnetic braking caused by magnetic field diffusion in the first core region. Thus, there is no catastrophic magnetic braking in a model with realistic non-ideal MHD.

We present Kepler observations of the bright (V=8.3), oscillating star HD 179070. The observations show transit-like events which reveal that the star is orbited every 2.8 days by a small, 1.6 R_Earth object. Seismic studies of HD 179070 using short cadence Kepler observations show that HD 179070 has a frequencypower spectrum consistent with solar-like oscillations that are acoustic p-modes. Asteroseismic analysis provides robust values for the mass and radius of HD 179070, 1.34{\pm}0.06 M{\circ} and 1.86{\pm}0.04 R{\circ} respectively, as well as yielding an age of 2.84{\pm}0.34 Gyr for this F5 subgiant. Together with ground-based follow-up observations, analysis of the Kepler light curves and image data, and blend scenario models, we conservatively show at the >99.7% confidence level (3{\sigma}) that the transit event is caused by a 1.64{\pm}0.04 R_Earth exoplanet in a 2.785755{\pm}0.000032 day orbit. The exoplanet is only 0.04 AU away from the star and our spectroscopic observations provide an upper limit to its mass of ~10 M_Earth (2-{\sigma}). HD 179070 is the brightest exoplanet host star yet discovered by Kepler.

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.

Accurate and precise values of radii and masses of stars are needed to correctly estimate properties of extrasolar planets. We examine the effect of uncertainties in stellar model parameters on estimates of the masses, radii and average densities of solar-type stars. We find that in the absence of seismic data on solar-like oscillations, stellar masses can be determined to a greater accuracy than either stellar radii or densities; but to get reasonably accurate results the effective temperature, log g and metallicity must be measured to high precision. When seismic data are available, stellar density is the most well determined property, followed by radius, with mass the least well determined property. Uncertainties in stellar convection, quantified in terms of uncertainties in the value of the mixing length parameter, cause the most significant errors in the estimates of stellar properties.

We present an asteroseismic study of the solar-like stars KIC 11395018, KIC 10273246, KIC 10920273, KIC 10339342, and KIC 11234888 using short-cadence time series of more than eight months from the Kepler satellite. For four of these stars, we derive atmospheric parameters from spectra acquired with the Nordic Optical Telescope. The global seismic quantities (average large frequency separation and frequency of maximum power), combined with the atmospheric parameters, yield the mean density and surface gravity with precisions of 2% and ~0.03 dex, respectively. We also determine the radius, mass, and age with precisions of 2-5%, 7-11%, and ~35%, respectively, using grid-based analyses. We determine asteroseismic distances to these stars with a precision better than 10%, and constrain the stellar inclination for three of the stars. An Li abundance analysis yields an independent estimate of the age, but this is inconsistent with the asteroseismically determined age for one of the stars. We compare the results from five different grid-based analyses, and we find that they all provide radius and mass values to within 2.4sigma. The absence of a metallicity constraint when the average large frequency separation is measured with a precision of 1% biases the fitted radius and mass for the stars with non-solar metallicity (metal-rich KIC 11395018 and metal-poor KIC 10273246), while including a metallicity constraint reduces the uncertainties in both of these parameters by almost a factor of two. We found that including the average small frequency separation improves the determination of the age only for KIC 11395018 and KIC 11234888, and for the latter this improvement was due to the lack of strong atmospheric constraints. (Abridged).

In addition to its search for extra-solar planets, the NASA Kepler Mission provides exquisite data on stellar oscillations. We report the detections of oscillations in 500 solartype stars in the Kepler field of view, an ensemble that is large enough to allow statistical studies of intrinsic stellar properties (such as mass, radius and age) and to test theories of stellar evolution. We find that the distribution of observed masses of these stars shows intriguing differences to predictions from models of synthetic stellar populations in the Galaxy.

Mass loss of red giant branch (RGB) stars is still poorly determined, despite its crucial role in the chemical enrichment of galaxies. Thanks to the recent detection of solar-like oscillations in G-K giants in open clusters with Kepler, we can now directly determine stellar masses for a statistically significant sample of stars in the old open clusters NGC6791 and NGC6819. The aim of this work is to constrain the integrated RGB mass loss by comparing the average mass of stars in the red clump (RC) with that of stars in the low-luminosity portion of the RGB (i.e. stars with L <~ L(RC)). Stellar masses were determined by combining the available seismic parameters numax and Dnu with additional photometric constraints and with independent distance estimates. We measured the masses of 40 stars on the RGB and 19 in the RC of the old metal-rich cluster NGC6791. We find that the difference between the average mass of RGB and RC stars is small, but significant (Delta M=0.09 +- 0.03 (random) +- 0.04 (systematic) Msun). Interestingly, such a small DeltaM does not support scenarios of an extreme mass loss for this metal-rich cluster. If we describe the mass-loss rate with Reimers' prescription, a first comparison with isochrones suggests that the observed DeltaM is compatible with a mass-loss efficiency parameter in the range 0.1 <~ eta <~ 0.3. Less stringent constraints on the RGB mass-loss rate are set by the analysis of the ~ 2 Gyr-old NGC6819, largely due to the lower mass loss expected for this cluster, and to the lack of an independent and accurate distance determination. In the near future, additional constraints from frequencies of individual pulsation modes and spectroscopic effective temperatures, will allow further stringent tests of the Dnu and numax scaling relations, which provide a novel, and potentially very accurate, means of determining stellar radii and masses.

The large frequency separation between modes of the same degree and consecutive orders in a star is approximately proportional to the square root of its mean density. To determine the large frequency separation as accurately as possible a mean large frequency separation computed over several orders is often used. It is, however, known that the large frequency separation varies with frequency in a second order effect. From observations it has been shown that this frequency dependence is more important for main-sequence stars than it is for red-giant stars. Here we use YREC models to verify and explain this observational result. We find that for stars with R > 8 Rsun the effect of the Helium second ionisation zone is relatively small. For these stars the deep location of the He II zone induces a frequency modulation covering only a few large separations, while the amplitude of the modulation is low due to the relatively weak and extended He II layer, causing a shallow wide depression in the first adiabatic exponent. For less evolved stars the He II zone is located closer to the surface, and it is more confined, i.e. a deep narrow depression in the first adiabatic exponent. This causes frequency modulations with relatively high amplitudes covering up to about 20 large separations, inducing a relatively large frequency modulation. Additionally, we find that for less evolved stars the He II zone is stronger and more localised for more massive stars and for stars with low metallicities further increasing the amplitude of the frequency modulation.

We calculate precise stellar radii and surface gravities from the asteroseismic analysis of over 500 solar-type pulsating stars observed by the Kepler space telescope. These physical stellar properties are compared with those given in the Kepler Input Catalog (KIC), determined from ground-based multi-color photometry. For the stars in our sample, we find general agreement but we detect an average overestimation bias of 0.23 dex in the KIC determination of log (g) for stars with log (g)_KIC > 4.0 dex, and a resultant underestimation bias of up to 50% in the KIC radii estimates for stars with R_KIC < 2 R sun. Part of the difference may arise from selection bias in the asteroseismic sample; nevertheless, this result implies there may be fewer stars characterized in the KIC with R ~ 1 R sun than is suggested by the physical properties in the KIC. Furthermore, if the radius estimates are taken from the KIC for these affected stars and then used to calculate the size of transiting planets, a similar underestimation bias may be applied to the planetary radii.

The successful launches of the CoRoT and Kepler space missions have led to the detections of solar-like oscillations in large samples of red-giant stars. The large numbers of red giants with observed oscillations make it possible to investigate the properties of the sample as a whole: ensemble asteroseismology. In this article we summarise ensemble asteroseismology results obtained from data released by the Kepler Science Team (~150,000 field stars) as presented by Hekker et al. (2011b) and for the clusters NGC 6791, NGC 6811 and NGC 6819 (Hekker et al. 2011a) and we discuss the importance of such studies.

Asteroseismology of solar-type stars has entered a new era of large surveys with the success of the NASA \textit{Kepler} mission, which is providing exquisite data on oscillations of stars across the Hertzprung-Russell (HR) diagram. From the time-series photometry, the two seismic parameters that can be most readily extracted are the large frequency separation ($\Delta\nu$) and the frequency of maximum oscillation power ($\nu_\mathrm{max}$). After the survey phase, these quantities are available for hundreds of solar-type stars. By scaling from solar values, we use these two asteroseismic observables to identify for the first time an evolutionary sequence of 1-M$_\odot$ field stars, without the need for further information from stellar models. Comparison of our determinations with the few available spectroscopic results shows an excellent level of agreement. We discuss the potential of the method for differential analysis throughout the main-sequence evolution, and the possibility of detecting twins of very well-known stars.

Studying star clusters offers significant advances in stellar astrophysics due to the combined power of having many stars with essentially the same distance, age, and initial composition. This makes clusters excellent test benches for verification of stellar evolution theory. To fully exploit this potential, it is vital that the star sample is uncontaminated by stars that are not members of the cluster. Techniques for determining cluster membership therefore play a key role in the investigation of clusters. We present results on three clusters in the Kepler field of view based on a newly established technique that uses asteroseismology to identify fore- or background stars in the field, which demonstrates advantages over classical methods such as kinematic and photometry measurements. Four previously identified seismic non-members in NGC6819 are confirmed in this study, and three additional non-members are found -- two in NGC6819 and one in NGC6791. We further highlight which stars are, or might be, affected by blending, which needs to be taken into account when analysing these Kepler data.

The tachocline is believed to be the region where the solar dynamo operates. With over a solar cycle's worth of data available from the MDI and GONG instruments, we are in a position to investigate not merely the average structure of the solar tachocline, but also its time variations. We determine the properties of the tachocline as a function of time by fitting a two-dimensional model that takes latitudinal variations of the tachocline properties into account. We confirm that if we consider central position of the tachocline, it is prolate. Our results show that the tachocline is thicker at higher latitudes than the equator, making the overall shape of the tachocline more complex. Of the tachocline properties examined, the transition of the rotation rate across the tachocline, and to some extent the position of the tachocline, show some temporal variations.

Scaling relations that link asteroseismic quantities to global stellar properties are important for gaining understanding of the intricate physics that underpins stellar pulsation. The common notion that all stars in an open cluster have essentially the same distance, age, and initial composition, implies that the stellar parameters can be measured to much higher precision than what is usually achievable for single stars. This makes clusters ideal for exploring the relation between the mode amplitude of solar-like oscillations and the global stellar properties. We have analyzed data obtained with NASA's Kepler space telescope to study solar-like oscillations in 100 red giant stars located in either of the three open clusters, NGC 6791, NGC 6819, and NGC 6811. By fitting the measured amplitudes to predictions from simple scaling relations that depend on luminosity, mass, and effective temperature, we find that the data cannot be described by any power of the luminosity-to-mass ratio as previously assumed. As a result we provide a new improved empirical relation which treats luminosity and mass separately. This relation turns out to also work remarkably well for main-sequence and subgiant stars. In addition, the measured amplitudes reveal the potential presence of a number of previously unknown unresolved binaries in the red clump in NGC 6791 and NGC 6819, pointing to an interesting new application for asteroseismology as a probe into the formation history of open clusters.

Context: Four open clusters are present in the Kepler field of view and timeseries of nearly a year in length are now available. These timeseries allow us to derive asteroseismic global oscillation parameters of red-giant stars in the three open clusters NGC 6791, NGC 6819 and NGC 6811. From these parameters and effective temperatures, we derive mass, radii and luminosities for the clusters as well as field red giants. Aims: We study the influence of evolution and metallicity on the observed red-giant populations. Methods: The global oscillation parameters are derived using different published methods and the effective temperatures are derived from 2MASS colours. The observational results are compared with BaSTI evolution models. Results: We find that the mass has significant influence on the asteroseismic quantities delta_nu vs. nu_max relation, while the influence of metallicity is negligible, under the assumption that the metallicity does not affect the excitation / damping of the oscillations. The positions of the stars in the H-R diagram depend on both mass and metallicity. Furthermore, the stellar masses derived for the field stars are bracketed by those of the cluster stars. Conclusions: Both the mass and metallicity contribute to the observed difference in locations in the H-R diagram of the old metal-rich cluster NGC 6791 and the middle-aged solar-metallicity cluster NGC 6819. For the young cluster NGC 6811, the explanation of the position of the stars in the H-R diagram challenges the assumption of solar metallicity, and this open cluster might have significantly lower metallicity [Fe/H] in the range -0.3 to -0.7 dex. Also, nearly all the observed field stars seem to be older than NGC 6811 and younger than NGC 6791.

High degree solar mode frequencies as measured by ring diagrams are known to change in the presence of the strong magnetic fields found in active regions. We examine these changes in frequency for a large sample of active regions analyzed with data from the Michelson Doppler Imager (MDI) onboard the SoHO spacecraft, spanning most of solar cycle 23. We confirm that the frequencies increase with increasing magnetic field strength, and that this dependence is generally linear. We find that the dependence is slightly but significantly different for active regions with different sunspot types.

We use photometric observations of solar-type stars, made by the NASA Kepler Mission, to conduct a statistical study of the impact of stellar surface activity on the detectability of solar-like oscillations. We find that the number of stars with detected oscillations fall significantly with increasing levels of activity. The results present strong evidence for the impact of magnetic activity on the properties of near-surface convection in the stars, which appears to inhibit the amplitudes of the stochastically excited, intrinsically damped solar-like oscillations.

Asteroseismology of solar-type stars has an important part to play in the exoplanet program of the NASA Kepler Mission. Precise and accurate inferences on the stellar properties that are made possible by the seismic data allow very tight constraints to be placed on the exoplanetary systems. Here, we outline how to make an estimate of the detectability of solar-like oscillations in any given Kepler target, using rough estimates of the temperature and radius, and the Kepler apparent magnitude.

We present initial results on some of the properties of open clusters NGC 6791 and NGC 6819 derived from asteroseismic data obtained by NASA's Kepler mission. In addition to estimating the mass, radius and log g of stars on the red-giant branch of these clusters, we estimate the distance to the clusters and their ages. Our model-independent estimate of the distance modulus of NGC 6791 is (m-M)_0= 13.11\pm 0.06. We find (m-M)_0= 11.85\pm 0.05 for NGC 6819. The average mass of stars on the red-giant branch of NGC 6791 is 1.20 \pm 0.01 M_sun, while that of NGC 6819 is 1.68\pm 0.03M_sun. It should be noted that we do not have data that cover the entire red-giant branch and the actual mass will be somewhat lower. We have determined model-dependent estimates of ages of these clusters. We find ages between 6.8 and 8.6 Gyr for NGC 6791, however, most sets of models give ages around 7Gyr. We obtain ages between 2 and 2.4 Gyr for NGC 6819.

NASA's Kepler mission is providing basic asteroseismic data for hundreds of stars. One of the more common ways of determining stellar characteristics from these data is by so-called "grid based" modelling. We have made a detailed study of grid-based analysis techniques to study the errors (and error-correlations) involved. As had been reported earlier, we find that it is relatively easy to get very precise values of stellar radii using grid-based techniques. However, we find that there are small, but significant, biases that can result because of the grid of models used. The biases can be minimized if metallicity is known. Masses cannot be determined as precisely as the radii, and suffer from larger systematic effects. We also find that the errors in mass and radius are correlated. A positive consequence of this correlation is that log g can be determined both precisely and accurately with almost no systematic biases. Radii and log g can be determined with almost no model dependence to within 5% for realistic estimates of error in asteroseismic and conventional observations. Errors in mass can be somewhat higher unless accurate metallicity estimates are available. Age estimates of individual stars are the most model dependent. The errors are larger too. However, we find that for star-clusters, it is possible to get a relatively precise age if one assumes that all stars in a given cluster have the same age.

We employ the three-dimensional magnetohydrodynamic simulation including ambipolar diffusion to study the gravitationally-driven fragmentation of subcritical molecular clouds, in which the gravitational fragmentation is stabilized as long as magnetic flux-freezing applies. The simulations show that the cores in an initially subcritical cloud generally develop gradually over an ambipolar diffusion time, which is about a few $\times 10^7$ years in a typical molecular cloud. On the other hand, the formation of collapsing cores in subcritical clouds is accelerated by supersonic nonlinear flows. Our parameter study demonstrates that core formation occurs faster as the strength of the initial flow speed in the cloud increases. We found that the core formation time is roughly proportional to the inverse of the square root of the enhanced density created by the supersonic nonlinear flows. The density dependence is similar to that derived in quasistatically contracting magnetically supported clouds, although the core formation conditions are created by the nonlinear flows in our simulations. We have also found that the accelerated formation time is not strongly dependent on the initial strength of the magnetic field if the cloud is highly subcritical. Our simulation shows that the core formation time in our model subcritical clouds is several $\times 10^6$ years, due to the presence of large-scale supersonic flows ($\sim 3$ times sound speed). Once a collapsing core forms, the density, velocity, and magnetic field structure of the core does not strongly depend on the initial strength of the velocity fluctuation. The infall velocities of the cores are subsonic and the magnetic field lines show weak hourglass shapes.

In addition to the well-known 11-year solar cycle, the Sun's magnetic activity also shows significant variation on shorter time scales, e.g. between one and two years. We observe a quasi-biennial (2-year) signal in the solar p-mode oscillation frequencies, which are sensitive probes of the solar interior. The signal is visible in Sun-as-a-star data observed by different instruments and here we describe the results obtained using BiSON, GOLF, and VIRGO data. Our results imply that the 2-year signal is susceptible to the influence of the main 11-year solar cycle. However, the source of the signal appears to be separate from that of the 11-year cycle. We speculate as to whether it might be the signature of a second dynamo, located in the region of near-surface rotational shear.

Magnetic diffusion plays a vital role in star formation. We trace its influence from interstellar cloud scales down to star-disk scales. On both scales, we find that magnetic diffusion can be significantly enhanced by the buildup of strong gradients in magnetic field structure. Large scale nonlinear flows can create compressed cloud layers within which ambipolar diffusion occurs rapidly. However, in the flux-freezing limit that may be applicable to photoionized molecular cloud envelopes, supersonic motions can persist for long times if driven by an externally generated magnetic field that corresponds to a subcritical mass-to-flux ratio. In the case of protostellar accretion, rapid magnetic diffusion (through Ohmic dissipation with additional support from ambipolar diffusion) near the protostar causes dramatic magnetic flux loss. By doing so, it also allows the formation of a centrifugal disk, thereby avoiding the magnetic braking catastrophe.

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

Clear power excess in a frequency range typical for solar-type oscillations in red giants has been detected in more than 1000 stars, which have been observed during the first 138 days of the science operation of the NASA Kepler satellite. This sample includes stars in a wide mass and radius range with spectral types G and K, extending in luminosity from the bottom of the giant branch up to high-luminous red giants. The high-precision asteroseismic observations with Kepler provide a perfect source for testing stellar structure and evolutionary models, as well as investigating the stellar population in our Galaxy. We fit a global model to the observed frequency spectra, which allows us to accurately estimate the granulation background signal and the global oscillation parameters, such as the frequency of maximum oscillation power. We find regular patterns of radial and non-radial oscillation modes and use a new technique to automatically identify the mode degree and the characteristic frequency separations between consecutive modes of the same spherical degree. In most cases, we can also measure the small separation. The seismic parameters are used to estimate stellar masses and radii and to place the stars in an H-R diagram by using an extensive grid of stellar models that covers a wide parameter range. Using Bayesian techniques throughout our analysis allows us to determine reliable uncertainties for all parameters. We provide accurate seismic parameters and their uncertainties for a large sample of red giants and determine their asteroseismic fundamental parameters. We investigate the influence of the stars' metallicities on their positions in the H-R diagram. We study the red-giant populations in the red clump and bump and compare them to a synthetic population and find a mass and metallicity gradient in the red clump and clear evidence of a secondary-clump population.

The primary science goal of the Kepler Mission is to provide a census of exoplanets in the solar neighborhood, including the identification and characterization of habitable Earth-like planets. The asteroseismic capabilities of the mission are being used to determine precise radii and ages for the target stars from their solar-like oscillations. Chaplin et al. (2010) published observations of three bright G-type stars, which were monitored during the first 33.5 days of science operations. One of these stars, the subgiant KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that it has evolved significantly. We have derived asteroseismic estimates of the properties of KIC 11026764 from Kepler photometry combined with ground-based spectroscopic data. We present the results of detailed modeling for this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously. We determine both the radius and the age of KIC 11026764 with a precision near 1%, and an accuracy near 2% for the radius and 15% for the age. Continued observations of this star promise to reveal additional oscillation frequencies that will further improve the determination of its fundamental properties.

We perform axisymmetric resistive MHD calculations that demonstrate that centrifugal disks can indeed form around Class 0 objects despite magnetic braking. We follow the evolution of a prestellar core all the way to near-stellar densities and stellar radii. Under flux-freezing, the core is braked and disk formation is inhibited, while Ohmic dissipation renders magnetic braking ineffective within the first core. In agreement with observations that do not show evidence for large disks around Class 0 objects, the resultant disk forms in close proximity to the second core and has a radius of only $\approx 10~R_{\odot}$ early on. Disk formation does not require enhanced resistivity. We speculate that the disks can grow to the sizes observed around Class II stars over time under the influence of both Ohmic dissipation and ambipolar diffusion, as well as internal angular momentum redistribution.

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences in the rotational properties and the depths of the surface convection zones for stars with various masses and ages. In 2007 we initiated a long-term monitoring campaign of Ca II H and K emission for a sample of 57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when possible. We report the discovery of a 1.6-year magnetic activity cycle in the exoplanet host star iota Horologii, and we obtain an estimate of the rotation period that is consistent with Hyades membership. This is the shortest activity cycle so far measured for a solar-type star, and may be related to the short-timescale magnetic variations recently identified in the Sun and HD49933 from helio- and asteroseismic measurements. Future asteroseismic observations can be compared to those obtained near the magnetic minimum in 2006 to search for cycle-induced shifts in the oscillation frequencies. If such short activity cycles are common in F stars, then NASA's Kepler mission should observe their effects in many of its long-term asteroseismic targets.

The minimum of solar cycle 24 is significantly different from most other minima in terms of its duration as well as its abnormally low levels of activity. Using available helioseismic data that cover epochs from the minimum of cycle 23 to now, we study the differences in the nature of the solar rotation between the minima of cycles 23 and 24. We find that there are significant differences between the rotation rates during the two minima. There are differences in the zonal-flow pattern too. We find that the band of fast rotating region close to the equator bifurcated around 2005 and recombined by 2008. This behavior is different from that during the cycle 23 minimum. By auto-correlating the zonal-flow pattern with a time shift, we find that in terms of solar dynamics, solar cycle 23 lasted for a period of 11.7 years, consistent with the result of Howe et al. (2009). The autocorrelation coefficient also confirms that the zonal-flow pattern penetrates through the convection zone.

We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra, including the presence of a possible signature of faculae, and the presence of mixed modes in one of the three stars.

We revisit our original papers on the burst mode of accretion by incorporating a detailed energy balance equation into a thin-disk model for the formation and evolution of circumstellar disks around low-mass protostars.Our model includes the effect of radiative cooling, viscous and shock heating, and heating due to stellar and background irradiation. Following the collapse from the prestellar phase allows us to model the early embedded phase of disk formation and evolution. During this time, the disk is susceptible to fragmentation, depending upon the properties of the initial prestellar core. Globally, we find that higher initial core angular momentum and mass content favors more fragmentation, but higher levels of background radiation can moderate the tendency to fragment. A higher rate of mass infall onto the disk than that onto the star is a necessary but not sufficient condition for disk fragmentation. More locally, both the Toomre Q-parameter needs to be below a critical value _and_ the local cooling time needs to be shorter than a few times the local dynamical time. Fragments that form during the early embedded phase tend to be driven into the inner disk regions, and likely trigger mass accretion and luminosity bursts that are similar in magnitude to FU-Orionis-type or EX-Lupi-like events. Disk accretion is shown to be an intrinsically variable process, thanks to disk fragmentation, nonaxisymmetric structure, and the effect of gravitational torques. The additional effect of a generic \alpha-type viscosity acts to reduce burst frequency and accretion variability, and is likely to not be viable for values of \alpha significantly greater than 0.01.

We present a brief overview of the history of attempts to obtain a clear detection of solar-like oscillations in cluster stars, and discuss the results on the first clear detection, which was made by the Kepler Asteroseismic Science Consortium (KASC) Working Group 2.

Observations from the Kepler satellite were recently published for three bright G-type stars, which were monitored during the first 33.5d of science operations. One of these stars, KIC 11026764, exhibits a characteristic pattern of oscillation frequencies suggesting that the star has evolved significantly. We have derived initial estimates of the properties of KIC 11026764 from the oscillation frequencies observed by Kepler, combined with ground-based spectroscopic data. We present preliminary results from detailed modeling of this star, employing a variety of independent codes and analyses that attempt to match the asteroseismic and spectroscopic constraints simultaneously.

The Sun is a variable star whose magnetic activity varies most perceptibly on a timescale of approximately 11 years. However, significant variation is also observed on much shorter timescales. We observe a quasi-biennial (2 year) signal in the natural oscillation frequencies of the Sun. The oscillation frequencies are sensitive probes of the solar interior and so by studying them we can gain information about conditions beneath the solar surface. Our results point strongly to the 2 year signal being distinct and separate from, but nevertheless susceptible to the influence of, the main 11 year solar cycle.

We determine the dependence of the initial helium abundance and the present-day helium abundance in the convective envelope of solar models ($\yim$ and $\ysm$ respectively) on the parameters that are used to construct the models. We do so by using reference standard solar models to compute the power-law coefficients of the dependence of $\yim$ and $\ysm$ on the input parameters. We use these dependencies to determine the correlation between $\yim$ and $\ysm$ and use this correlation to eliminate uncertainties in $\yim$ from all solar model input parameters except the microscopic diffusion rate. We find an expression for $\yim$ that depends only on $\ysm$ and the diffusion rate. By adopting the helioseismic determination of solar surface helium abundance, $\ysms= 0.2485\pm0.0035$, and an uncertainty of 20% for the diffusion rate, we find that the initial solar helium abundance, $\yims$, is $0.278 \pm 0.006$ independently of the reference standard solar models (and particularly on adopted the solar abundances) used in the derivation of the correlation between $\yim$ and $\ysm$. When non-standard solar models with extra mixing are used, then we derive $\yims = 0.273 \pm 0.006$. In both cases, the derived $\yims$ value is higher than that directly derived from solar model calibrations when the low metalicity solar abundances (e.g. by Asplund et al.) are adopted in the models.

We have analyzed available full-disc data from the Michelson Doppler Imager (MDI) on board SoHO using the "ring diagram" technique to determine the behavior of solar meridional flows over solar cycle 23 in the outer 2% of the solar radius. We find that the dominant component of meridional flows during solar maximum was much lower than that during the minima at the beginning of cycles 23 and 24. There were differences in the flow velocities even between the two minima. The meridional flows show a migrating pattern with higher-velocity flows migrating towards the equator as activity increases. Additionally, we find that the migrating pattern of the meridional flow matches those of sunspot butterfly diagram and the zonal flows in the shallow layers. A high latitude band in meridional flow appears around 2004, well before the current activity minimum. A Legendre polynomial decomposition of the meridional flows shows that the latitudinal pattern of the flow was also different during the maximum as compared to that during the two minima. The different components of the flow have different time-dependences, and the dependence is different at different depths.

We calculate and analyze the longevity of magnetohydrodynamic (MHD) wave modes that occur in the plane of a magnetic thin sheet. Initial turbulent conditions applied to a magnetically subcritical cloud are shown to lead to relatively rapid energy decay if ambipolar diffusion is introduced at a level corresponding to partial ionization primarily by cosmic rays. However, in the flux-freezing limit, as may be applicable to photoionized molecular cloud envelopes, the turbulence persists at "nonlinear" levels in comparison with the isothermal sound speed $\cs$, with one-dimensional rms material motions in the range of $\approx 2\,\cs -5\,\cs$ for cloud sizes in the range of $\approx 2\,\pc - 16\,\pc$. These fluctuations persist indefinitely, maintaining a significant portion of the initial turbulent kinetic energy. We find the analytic explanation for these persistent fluctuations. They are magnetic-tension-driven modes associated with the interaction of the sheet with the external magnetic field. The phase speed of such modes is quite large, allowing residual motions to persist without dissipation in the flux-freezing limit, even as they are nonlinear with respect to the sound speed. We speculate that long-lived large-scale MHD modes such as these may provide the key to understanding observed supersonic motions in molecular clouds.

The Birmingham Solar-Oscillations Network (BiSON) has collected helioseismic data over three solar cycles. We use these data to determine how the internal properties of the Sun during this minimum differ from the previous two minima. The cycle 24 data show oscillatory differences with respect to the other two sets, indicating relatively localized changes in the solar interior. Analysis of MDI data from Cycle 23 and Cycle 24 also show significant signs of differences.

Motivated by the recent discovery of massive planets on wide orbits, we present a mechanism for the formation of such planets via disk fragmentation in the embedded phase of star formation. In this phase, the forming disk intensively accretes matter from the natal cloud core and undergoes several fragmentation episodes. However, most fragments are either destroyed or driven into the innermost regions (and probably onto the star) due to angular momentum exchange with spiral arms, leading to multiple FU-Ori-like bursts and disk expansion. Fragments that are sufficiently massive and form in the late embedded phase (when the disk conditions are less extreme) may open a gap and evolve into giant planets on typical orbits of several tens to several hundreds of AU. For this mechanism to work, the natal cloud core must have sufficient mass and angular momentum to trigger the burst mode and also form extended disks of the order of several hundreds of AU. When mass loading from the natal cloud core diminishes and the main fragmentation phase ends, such extended disks undergo a transient episode of contraction and density increase, during which they may give birth to a last and survivable set of giant planets on wide and relatively stable orbits.

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

Asteroseismology involves probing the interiors of stars and quantifying their global properties, such as radius and age, through observationsof normal modes of oscillation. The technical requirements for conducting asteroseismology include ultra-high precision measured in photometry in parts per million, as well as nearly continuous time series over weeks to years, and cadences rapid enough to sample oscillations with periods as shortas a few minutes. We report on results from the first 43 days of observations in which the unique capabilities of Kepler in providing a revolutionary advance in asteroseismology are already well in evidence. The Kepler asteroseismology program holds intrinsic importance in supporting the core planetary search program through greatly enhanced knowledge of host star properties, and extends well beyond this to rich applications in stellar astrophysics.

Asteroseismology of stars in clusters has been a long-sought goal because the assumption of a common age, distance and initial chemical composition allows strong tests of the theory of stellar evolution. We report results from the first 34 days of science data from the Kepler Mission for the open cluster NGC 6819 -- one of four clusters in the field of view. We obtain the first clear detections of solar-like oscillations in the cluster red giants and are able to measure the large frequency separation and the frequency of maximum oscillation power. We find that the asteroseismic parameters allow us to test cluster-membership of the stars, and even with the limited seismic data in hand, we can already identify four possible non-members despite their having a better than 80% membership probability from radial velocity measurements. We are also able to determine the oscillation amplitudes for stars that span about two orders of magnitude in luminosity and find good agreement with the prediction that oscillation amplitudes scale as the luminosity to the power of 0.7. These early results demonstrate the unique potential of asteroseismology of the stellar clusters observed by Kepler.

Understanding the reasons of the cyclic variation of the total solar irradiance is one of the most challenging targets of modern astrophysics. These studies prove to be essential also for a more climatologic issue, associated to the global warming. Any attempt to determine the solar components of this phenomenon must include the effects of the magnetic field, whose strength and shape in the solar interior are far from being completely known. Modelling the presence and the effects of a magnetic field requires a 2D approach, since the assumption of radial symmetry is too limiting for this topic. We present the structure of a 2D evolution code that was purposely designed for this scope; rotation, magnetic field and turbulence can be taken into account. Some preliminary results are presented and commented.

Normal modes of oscillation of the Sun are useful probes of the solar interior. In this work, we use the even-order splitting coefficients to study the evolution of magnetic fields in the convection zone over solar cycle 23, assuming that the frequency splitting is only due to rotation and a large scale magnetic field. We find that the data are best fit by a combination of a poloidal field and a double-peaked near-surface toroidal field. The toroidal fields are centered at r=0.999R_solar and r=0.996R_solar and are confined to the near-surface layers. The poloidal field is a dipole field. The peak strength of the poloidal field is 124 +/- 17G. The toroidal field peaks at 380 +/- 30G and 1.4 +/- 0.2kG for the shallower and deeper fields respectively. The field strengths are highly correlated with surface activity. The toroidal field strength shows a hysteresis-like effect when compared to the global 10.7 cm radio flux. The poloidal field strength shows evidence of saturation at high activity.

We construct updated solar models with different sets of solar abundances, including the most recent determinations by Asplund et al. (2009). The latter work predicts a larger ($\sim 10%$) solar metallicity compared to previous measurements by the same authors but significantly lower ($\sim 25%$) than the recommended value from a decade ago by Grevesse & Sauval (1998). We compare the results of our models with determinations of the solar structure inferred through helioseismology measurements. The model that uses the most recent solar abundance determinations predicts the base of the solar convective envelope to be located at $R_{\rm CZ}= 0.724{\rm R_\odot}$ and a surface helium mass fraction of $Y_{\rm surf}=0.231$. These results are in conflict with helioseismology data ($R_{\rm CZ}= 0.713\pm0.001{\rm R_\odot}$ and $Y_{\rm surf}=0.2485\pm0.0035$) at 5$-\sigma$ and 11$-\sigma$ levels respectively. Using the new solar abundances, we calculate the magnitude by which radiative opacities should be modified in order to restore agreement with helioseismology. We find that a maximum change of $\sim 15%$ at the base of the convective zone is required with a smooth decrease towards the core, where the change needed is $\sim 5%$. The required change at the base of the convective envelope is about half the value estimated previously. We also present the solar neutrino fluxes predicted by the new models. The most important changes brought about by the new solar abundances are the increase by $\sim 10%$ in the predicted $^{13}$N and $^{15}$O fluxes that arise mostly due to the increase in the C and N abundances in the newly determined solar composition.

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences in the rotational properties and the depths of the surface convection zones for stars with various masses and ages. Asteroseismic data will soon provide direct measurements of these quantities for individual stars, but many of the most promising targets are in the southern sky (e.g., alpha Cen A & B, beta Hyi, mu Ara, tau Cet, nu Ind), while long-term magnetic activity cycle surveys are largely confined to the north. In 2007 we began using the SMARTS 1.5-m telescope to conduct a long-term monitoring campaign of Ca II H & K emission for a sample of 57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when possible. This sample includes the most likely southern asteroseismic targets to be observed by the Stellar Oscillations Network Group (SONG), currently scheduled to begin operations in 2012. We present selected results from the first two years of the survey, and from the longer time baseline sampled by a single-epoch survey conducted in 1992.

Establishing relations between global stellar parameters and asteroseismic quantities can help improve our understanding of stellar astrophysics and facilitate the interpretation of observations. We present an observed relation between the large frequency separation, $\Delta\nu$, and the frequency of maximum power, $\nu_{max}$. We find that $\Delta\nu$ is proportional to $(\nu_{max})^0.77$, allowing prediction of $\Delta\nu$ to about 15 per cent given $\nu_{max}$. Our result is further supported by established scaling relations for $\Delta\nu$ and $\nu_{max}$ and by extended stellar model calculations, which confirm that $\Delta\nu$ can be estimated using this relation for basically any star showing solar-like oscillations in the investigated range (0.5<M/Msol<4.0).

The NASA Kepler mission is designed to find planets through transits. Accurate and precise radii of the detected planets depend on knowing the radius of the host star accurately, which is difficult unless the temperature and luminosity of the star are known precisely. Kepler, however, has an asteroseismology programme that will provide seismic variables that can characterise stellar radii easily, accurately, and extremely precisely. In this paper we describe the Yale-Birmingham (YB) method to determine stellar radii using a combination of seismic and conventional variables, and analyse the effect of these variables on the result. We find that for main-sequence stars, a knowledge of the parallax is not important to get accurate radii using the YB method: we can get results to an accuracy and precision of better than a few percent if we know the effective temperature and the seismic parameters for these stars. Metallicity does not make much difference either. However, good estimates of the effective temperature and metallicity, along with those of the seismic parameters, are essential to determine radii of sub giants properly. On the other hand, for red giants we find that determining radii properly is not possible without a good estimate of the parallax. We find that the so called "surface term" in the seismic data has minimal effect on the inferred radii. Uncertainties in the convective mixing length can matter under some circumstances and can cause a systematic shift in the inferred radii. Blind tests with data simulated to match those expected from the asteroseismic Survey Phase of Kepler show that it will be possible to infer stellar radii successfully using our method.

We present numerical solutions of the collapse of prestellar cores that lead to the formation and evolution of circumstellar disks. The disk evolution is then followed for up to three million years. A variety of models of different initial masses and rotation rates allows us to study disk accretion around brown dwarfs and low-mass T Tauri stars, with central object mass $M_st < 0.2 Msun$, as well as intermediate and upper-mass T Tauri stars (0.2 Msun < M_st < 3.0 Msun). Our models include self-gravity and allow for nonaxisymmetric motions. In addition to the self-consistently generated gravitational torques, we introduce an effective turbulent \alpha-viscosity with \alpha = 0.01, which allows us particularly to model accretion in the low-mass regime where disk self-gravity is diminishing. A range of models with observationally-motivated values of the initial ratio of rotational to gravitational energy yields a correlation between mass accretion rate \dot{M} and M_st that is relatively steep, as observed. Additionally, our modeling reveals evidence for a bimodality in the \dot{M}--M_st correlation, with a steeper slope at lower masses and a shallower slope at intermediate and upper masses, as also implied by observations. We show that the neglect of disk self-gravity leads to a much steeper \dot{M}--M_st relation for intermediate and upper-mass T Tauri stars. This demonstrates that an accurate treatment of global self-gravity is essential to understanding observations of circumstellar disks.

We present new results on the structure of the solar core, obtained with new sets of frequencies of solar low-degree p modes obtained from the BiSON network. We find that different methods used in extracting the different sets of frequencies cause shifts in frequencies, but the shifts are not large enough to affect solar structure results. We find that the BiSON frequencies show that the solar sound speed in the core is slightly larger than that inferred from data from MDI low-degree modes, and the uncertainties on the inversion results are smaller. Density results also change by a larger amount, and we find that solar models now tend to show smaller differences in density compared to the Sun. The result is seen at all radii, a result of the fact that conservation of mass implies that density differences in one region have to cancel out density differences in others, since our models are constructed to have the same mass as the Sun. The uncertainties on the density results are much smaller too. We attribute the change in results to having more, and lower frequency, low-degree mode frequencies available. These modes provide greater sensitivity to conditions in the core.

We analyze a suite of thin sheet magnetohydrodynamical simulations based on the formulation of Basu, Ciolek, Dapp & Wurster. These simulations allow us to examine the observational consequences to a star-forming region of varying the input level of turbulence (between thermal and a Mach number of 4) and the initial magnetic field strength corresponding to a range of mass to flux ratios between subcritical (mu_0=0.5) and supercritical (mu_0=10). The input turbulence is allowed to decay over the duration of the simulation. We compare the measured observable quantities with those found from surveying the Perseus molecular cloud. We find that only the most turbulent of simulations (high Mach number and weak magnetic field) have sufficient large-scale velocity dispersion (at ~1 pc) to match that observed across extinction regions in Perseus. Generally, the simulated core (~0.02 pc) and line of sight velocity dispersions provide a decent match to observations. The motion between the simulated core and its local environment, however, is far too large in simulations with high large-scale velocity dispersion.

We present a new analytical three-parameter formula to fit observed column density profiles of prestellar cores. It represents a line-of-sight integral through a spherically symmetric or disc-like isothermal cloud. The underlying model resembles the Bonnor-Ebert model in that it features a flat central region leading into a power-law decline \propto r^{-2} in density, and a well-defined outer radius. However, we do not assume that the cloud is in equilibrium, and can instead make qualitative statements about its dynamical state (expansion, equilibrium, collapse) using the size of the flat region as a proxy. Instead of having temperature as a fitting parameter, our model includes it as input, and thus avoids possible inconsistencies. It is significantly easier to fit to observational data than the Bonnor-Ebert sphere. We apply this model to L1689B and B68. We show that L1689B cannot be in equilibrium but instead appears to be collapsing, while our model verifies that B68 is not far from being a hydrostatic object.

We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field ("turbulence") is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway collapse (formation of the first star) in subcritical clouds is significantly accelerated by nonlinear flows in which a large-scale wave mode dominates the power spectrum. In addition, we find that a power spectrum with equal energy on all scales also accelerates collapse, but by a lesser amount. For a highly super-Alfvenic initial velocity field with most power on the largest scales, the collapse occurs promptly during the initial compression wave. However, for trans-Alfvenic perturbations, a subcritical magnetic field causes a rebound from the initial compression, and the system undergoes several oscillations before runaway collapse occurs. Models that undergo prompt collapse have highly supersonic infall motions at the core boundaries. Cores in magnetically subcritical models with trans-Alfvenic initial perturbations also pick up significant systematic speeds by inheriting motions associated with magnetically-driven oscillations. Core mass distributions are much broader than in models with small-amplitude initial perturbations, although the disturbed structure of cores that form due to nonlinear flows does not guarantee subsequent monolithic collapse. Our simulations also demonstrate that significant power can (if present initially) be maintained with negligible dissipation in large-scale compressive modes of a magnetic thin sheet, in the limit of perfect flux freezing.

We add the effect of turbulent viscosity via the \alpha-prescription to models of the self-consistent formation and evolution of protostellar discs. Our models are non-axisymmetric and carried out using the thin-disc approximation. Self-gravity plays an important role in the early evolution of a disc, and the later evolution is determined by the relative importance of gravitational and viscous torques. In the absence of viscous torques, a protostellar disc evolves into a self-regulated state with disk-averaged Toomre parameter Q \sim 1.5-2.0, non-axisymmetric structure diminishing with time, and maximum disc-to-star mass ratio \xi = 0.14. We estimate an effective viscosity parameter \alpha_eff associated with gravitational torques at the inner boundary of our simulation to be in the range 10^{-4}-10^{-3} during the late evolution. Addition of viscous torques with a low value \alpha = 10^{-4} has little effect on the evolution, structure, and accretion properties of the disc, and the self-regulated state is largely preserved. A sequence of increasing values of \alpha results in the discs becoming more axisymmetric in structure, being more gravitationally stable, having greater accretion rates, larger sizes, shorter lifetimes, and lower disc-to-star mass ratios. For \alpha=10^{-2}, the model is viscous-dominated and the self-regulated state largely disappears by late times. (Abridged)

In the second paper of this series we pursue two objectives. First, in order to make the code more sensitive to small effects, we remove many approximations made in Paper I. Second, we include turbulence and rotation in the two-dimensional framework. The stellar equilibrium is described by means of a set of five differential equations, with the introduction of a new dependent variable, namely the perturbation to the radial gravity, that is found when the non-radial effects are considered in the solution of the Poisson equation; following the scheme of the first paper, we write the equations in such a way that the two-dimensional effects can be easily disentangled. The key concept introduced in this series is the equipotential surface. We use the underlying cause-effect relation to develop a recurrence relation to calculate the equipotential surface functions for uniform rotation, differential rotation, rotation-like toroidal magnetic fields and turbulence. We also develop a more precise code to numerically solve the two-dimensional stellar structure and evolution equations based on the equipotential surface calculations. We have shown that with this formulation we can achieve the precision required by observations by appropriately selecting the convergence criterion. Several examples are presented to show that the method works well. Since we are interested in modeling the effects of a dynamo-type field on the detailed envelope structure and global properties of the Sun, the code has been optimized for short timescales phenomena (down to 1 yr). The time dependence of the code has so far been tested exclusively to address such problems.

We compute the p-mode oscillation frequencies and frequency splittings that arise in a two-dimensional model of the Sun that contains toroidal magnetic fields in its interior.

Helioseismic techniques such as ring-diagram analysis have often been used to determine the subsurface structural differences between solar active and quiet regions. Results obtained by inverting the frequency differences between the regions are usually interpreted as the sound-speed differences between them. These in turn are used as a measure of temperature and magnetic-field strength differences between the two regions. In this paper we first show that the "sound-speed" difference obtained from inversions is actually a combination of sound-speed difference and a magnetic component. Hence, the inversion result is not directly related to the thermal structure. Next, using solar models that include magnetic fields, we develop a formulation to use the inversion results to infer the differences in the magnetic and thermal structures between active and quiet regions. We then apply our technique to existing structure inversion results for different pairs of active and quiet regions. We find that the effect of magnetic fields is strongest in a shallow region above 0.985R_sun and that the strengths of magnetic-field effects at the surface and in the deeper (r < 0.98R_sun) layers are inversely related, i.e., the stronger the surface magnetic field the smaller the magnetic effects in the deeper layers, and vice versa. We also find that the magnetic effects in the deeper layers are the strongest in the quiet regions, consistent with the fact that these are basically regions with weakest magnetic fields at the surface. Because the quiet regions were selected to precede or follow their companion active regions, the results could have implications about the evolution of magnetic fields under active regions.

We present results from an extensive set of simulations of gravitational fragmentation in the presence of magnetic fields and ambipolar diffusion. The average fragmentation spacing in the nonlinear phase of evolution is in excellent agreement with the prediction of linear perturbation theory. The time scale for nonlinear growth and runaway of the first core is $\approx 10$ times the calculated growth time $\taugm$ of the eigenmode with minimum growth time, when starting from a uniform background state with small-amplitude white-noise perturbations. Subcritical and transcritical models typically evolve on a significantly longer time scale than the supercritical models. Infall motions in the nonlinear fully-developed contracting cores are subsonic on the core scale in subcritical and transcritical clouds, but are somewhat supersonic in supercritical clouds. Core mass distributions are sharply peaked with a steep decline to large masses, consistent with the existence of a preferred mass scale for each unique set of dimensionless free parameters. However, a sum total of results for various initial mass-to-flux ratios yields a broad distribution reminiscent of observed core mass distributions. Based on our results, we conclude that fragmentation spacings, magnitude of infall motions, core shapes, and, especially, the curvature of magnetic field morphology, may serve as indirect observational means of determining a cloud's ambient mass-to-flux ratio.

The frequencies of solar oscillations are known to change with solar activity. We use Principal Component Analysis to examine these changes with high precision. In addition to the well-documented changes in solar normal mode oscillations with activity as a function of frequency, which originate in the surface layers of the Sun, we find a small but statistically significant change in frequencies with an origin at and below the base of the convection zone. We find that at r=(0.712^{+0.0097}_{-0.0029})R_sun, the change in sound speed is \delta c^2 / c^2 = (7.23 +/- 2.08) x 10^{-5} between high and low activity. This change is very tightly correlated with solar activity. In addition, we use the splitting coefficients to examine the latitudinal structure of these changes. We find changes in sound speed correlated with surface activity for r >~ 0.9R_sun.

Available helioseismic data now span almost the entire solar activity cycle 23 making it possible to study solar-cycle related changes of the solar rotation rate in detail. In this paper we study how the solar rotation rate, in particular, the zonal flows change with time. In addition to the zonal flows that show a well known pattern in the solar convection zone, we also study changes in the radial and latitudinal gradients of the rotation rate, particularly in the shear layer that is present in the immediate sub-surface layers of the Sun. In the case of the zonal-flow pattern, we find that the band indicating fast rotating region close to the equator seems to have bifurcated around 2005. Our investigation of the rotation-rate gradients show that the relative variation in the rotation-rate gradients is about 20% or more of their average values, which is much larger than the relative variation in the rotation rate itself. These results can be used to test predictions of various solar dynamo models.

We demonstrate that the formation of collapsing cores in subcritical clouds is accelerated by nonlinear flows, by performing three-dimensional non-ideal MHD simulations. An initial random supersonic (and trans-Alfvenic) turbulent-like flow is input into a self-gravitating gas layer that is threaded by a uniform magnetic field (perpendicular to the layer) such that the initial mass-to-flux ratio is subcritical. Magnetic ambipolar diffusion occurs very rapidly initially due to the sharp gradients introduced by the turbulent flow. It subsequently occurs more slowly in the traditional near-quasistatic manner, but in regions of greater mean density than present in the initial state. The overall timescale for runaway growth of the first core(s) is several times, 10^6 yr, even though previous studies have found a timescale of several times, 10^7 yr when starting with linear perturbations and similar physical parameters. Large-scale supersonic flows exist in the cloud and provide an observationally testable distinguishing characteristic from core formation due to linear initial perturbations. However, the nonlinear flows have decayed sufficiently that the relative infall motions onto the first core are subsonic, as in the case of starting from linear initial perturbations. The ion infall motions are very similar to those of neutrals; however, they lag the neutral infall in directions perpendicular to the mean magnetic field direction and lead the neutral infall in the direction parallel to the mean magnetic field.

We have studied numerically the evolution of protostellar disks around intermediate and upper mass T Tauri stars (0.25 M_sun < M_st < 3.0 M_sun) that have formed self-consistently from the collapse of molecular cloud cores. In the T Tauri phase, disks settle into a self-regulated state, with low-amplitude nonaxisymmetric density perturbations persisting for at least several million years. Our main finding is that the global effect of gravitational torques due to these perturbations is to produce disk accretion rates that are of the correct magnitude to explain observed accretion onto T Tauri stars. Our models yield a correlation between accretion rate M_dot and stellar mass M_st that has a best fit M_dot \propto M_st^{1.7}, in good agreement with recent observations. We also predict a near-linear correlation between the disk accretion rate and the disk mass.

Improving our understanding of the initial conditions and earliest stages of star formation is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. We review the properties of low-mass dense cores as derived from recent millimeter/submillimeter observations of nearby molecular clouds and discuss them in the context of various contemporary scenarios for cloud core formation and evolution. None of the extreme scenarios can explain all observations. Pure laminar ambipolar diffusion has relatively long growth times for typical ionization levels and has difficulty satisfying core lifetime constraints. Purely hydrodynamic pictures have trouble accounting for the inefficiency of core formation and the detailed velocity structure of individual cores. A possible favorable scenario is a mixed model involving gravitational fragmentation of turbulent molecular clouds close to magnetic criticality. The evolution of the magnetic field and angular momentum in individual cloud cores after the onset of gravitational collapse is also discussed. In particular, we stress the importance of radiation-magnetohydrodynamical processes and resistive MHD effects during the protostellar phase. We also emphasize the role of the formation of the short-lived first (protostellar) core in providing a chance for sub-fragmentation into binary systems and triggering MHD outflows. Future submillimeter facilities such as Herschel and ALMA will soon provide major new observational constraints in this field. On the theoretical side, an important challenge for the future will be to link the formation of molecular clouds and prestellar cores in a coherent picture.

We review the impact of global helioseismology on key questions concerning the internal structure and dynamics of the Sun, and consider the exciting challenges the field faces as it enters a fourth decade of science exploitation. We do so with an eye on the past, looking at the perspectives global helioseismology offered in its earlier phases, in particular the mid-to-late 1970s and the 1980s. We look at how modern, higher-quality, longer datasets coupled with new developments in analysis, have altered, refined, and changed some of those perspectives, and opened others that were not previously available for study. We finish by discussing outstanding challenges and questions for the field.

Helioseismology has allowed us to study the structure of the Sun in unprecedented detail. One of the triumphs of the theory of stellar evolution was that helioseismic studies had shown that the structure of solar models is very similar to that of the Sun. However, this agreement has been spoiled by recent revisions of the solar heavy-element abundances. Heavy element abundances determine the opacity of the stellar material and hence, are an important input to stellar model calculations. The models with the new, low abundances do not satisfy helioseismic constraints. We review here how heavy-element abundances affect solar models, how these models are tested with helioseismology, and the impact of the new abundances on standard solar models. We also discuss the attempts made to improve the agreement of the low-abundance models with the Sun and discuss how helioseismology is being used to determine the solar heavy-element abundance. A review of current literature shows that attempts to improve agreement between solar models with low heavy-element abundances and seismic inference have been unsuccessful so far. The low-metallicity models that have the least disagreement with seismic data require changing all input physics to stellar models beyond their acceptable ranges. Seismic determinations of the solar heavy-element abundance yield results that are consistent with the older, higher values of the solar abundance, and hence, no major changes to the inputs to solar models are required to make higher-metallicity solar models consistent with helioseismic data.

Context: BL Lacertae objects have been the targets for numerous recent multiwavelength campaigns, continuum spectral variability studies, and theoretical spectral and variability modeling. A meaningful interpretation of the results of such studies requires a reliable knowledge of the objects' redshifts; however, the redshifts for many are still unknown or uncertain. Aims: Therefore, we hope to determine or constrain the redshifts of six BL Lac objects with unknown or poorly known redshifts. Methods: Observations were made of these objects with the MDM 2.4 m Hiltner telescope. Although no spectral features were detected, and thus no redshifts could be measured, lower redshift limits were assigned to the objects based on the expected equivalent widths of absorption features in their host galaxies. Redshifts were also estimated for some objects by assuming the host galaxies are standard candles and using host galaxy apparent magnitudes taken from the literature. Results: The commonly used redshift of $z=0.102$ for 1219+285 is almost certainly wrong, while the redshifts of the other objects studied remain undetermined.

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular nonaxisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual nonzero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q\approx 1.5-2.0. (Abstract abridged).

We use very precise frequencies of low-degree solar-oscillation modes measured from 4752 days of data collected by the Birmingham Solar-Oscillations Network (BiSON) to derive seismic information on the solar core. We compare these observations to results from a large Monte Carlo simulation of standard solar models, and use the results to constrain the mean molecular weight of the solar core, and the metallicity of the solar convection zone. We find that only a high value of solar metallicity is consistent with the seismic observations. We can determine the mean molecular weight of the solar core to a very high precision, and, dependent on the sequence of Monte Carlo models used, find that the average mean molecular weight in the inner 20% by radius of the Sun ranges from 0.7209 to 0.7231, with uncertainties of less than 0.5% on each value. Our lowest seismic estimate of solar metallicity is Z=0.0187 and our highest is Z=0.0239, with uncertainties in the range of 12--19%. Our results indicate that the discrepancies between solar models constructed with low metallicity and the helioseismic observations extend to the solar core and thus cannot be attributed to deficiencies in the modeling of the solar convection zone.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

Recent downward revision of solar heavy-element abundances using three-dimensional atmospheric model has introduced serious discrepancies between standard solar models and helioseismic inferences about solar structure. In this paper, we investigate the possibility of determining the heavy-element abundances using helioseismic inversion techniques with the hope of providing an independent estimate. We use the adiabatic index, Gamma_1 (logarithmic partial derivative of pressure with respect to density at constant entropy) as a probe to examine the effects of the total heavy-element abundance, as well as the effects due to the abundance of individual elements. Our inversion results show that the new, lower, abundance increases the discrepancy between the Sun and the solar models.

We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion-neutral friction in regulating gravitationally-driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infall speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (\ga 0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (rho) and the strength of the magnetic field (B) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both relations tend to B \propto \rho^{0.5}.

We have used 4752 days of data collected by the Birmingham Solar-Oscillations Network (BiSON) to determine very precise oscillation frequencies of acoustic low-degree modes that probe the solar core. We compare the fine (small frequency) spacings and frequency separation ratios formed from these data with those of different solar models. We find that models constructed with low metallicity are incompatible with the observations. The results provide strong support for lowering the theoretical uncertainties on the neutrino fluxes. These uncertainties had recently been raised due to the controversy over the solar abundances.

We present results on the structure of the near-surface layers of the Sun obtained by inverting frequencies of high-degree solar modes from "ring diagrams". We have results for eight epochs between June 1996 and October 2003. The frequencies for each epoch were obtained from ring diagrams constructed from MDI Dopplergrams spanning complete Carrington rotations. We find that there is a substantial latitudinal variation of both sound speed and the adiabatic index Gamma_1 in the outer 2% of the Sun. We find that both the sound-speed and Gamma_1 profiles change with changes in the level of solar activity. In addition, we also study differences between the northern and southern hemispheres of the Sun and find a small asymmetry that appears to reflect the difference in magnetic activity between the two hemispheres.

We present a detailed study of the small frequency separations as diagnostics of the mass of the convective core and evolutionary stage of solar-type stars. We demonstrate how the small separations can be combined to provide sensitive tests for the presence of convective overshoot at the edge of the core. These studies are focused on low degree oscillation modes, the only modes expected to be detected in distant stars. Using simulated data with realistic errors, we find that the mass of the convective core can be estimated to within 5% if the total stellar mass is known. Systematic errors arising due to uncertainty in the mass could be up to 20%. The evolutionary stage of the star, determined in terms of the central hydrogen abundance using our proposed technique, however, is much less sensitive to the mass estimate.

We formulate the problem of the formation and collapse of nonaxisymmetric protostellar cores in weakly ionized, self-gravitating, magnetic molecular clouds. In our formulation, molecular clouds are approximated as isothermal, thin (but with finite thickness) sheets. We present the governing dynamical equations for the multifluid system of neutral gas and ions, including ambipolar diffusion, and also a self-consistent treatment of thermal pressure, gravitational, and magnetic (pressure and tension) forces. The dimensionless free parameters characterizing model clouds are discussed. The response of cloud models to linear perturbations is also examined, with particular emphasis on length and time scales for the growth of gravitational instability in magnetically subcritical and supercritical clouds. We investigate their dependence on a cloud's initial mass-to-magnetic-flux ratio (normalized to the critical value for collapse), the dimensionless initial neutral-ion collision time, and also the relative external pressure exerted on a model cloud. Among our results, we find that nearly-critical model clouds have significantly larger characteristic instability lengthscales than do more distinctly sub- or supercritical models. Another result is that the effect of a greater external pressure is to reduce the critical lengthscale for instability. Numerical simulations showing the evolution of model clouds during the linear regime of evolution are also presented, and compared to the results of the dispersion analysis. They are found to be in agreement with the dispersion results, and confirm the dependence of the characteristic length and time scales on parameters such as the initial mass-to-flux ratio and relative external pressure.

We present new numerical simulations in the thin-disk approximation which characterize the burst mode of protostellar accretion. The burst mode begins upon the formation of a centrifugally balanced disk around a newly formed protostar. It is comprised of prolonged quiescent periods of low accretion rate (typically $\la 10^{-7} \Msun$ yr$^{-1}$) which are punctuated by intense bursts of accretion (typically $\ga 10^{-4} \Msun$ yr$^{-1}$, with duration $\la 100$ yr) during which most of the protostellar mass is accumulated. The accretion bursts are associated with the formation of dense protostellar/protoplanetary embryos, which are later driven onto the protostar by the gravitational torques that develop in the disk. Gravitational instability in the disk, driven by continuing infall from the envelope, is shown to be an effective means of transporting angular momentum outward, and mass inward to the protostar. We show that the disk mass always remains significantly less than the central protostar mass throughout this process. The burst phenomenon is robust enough to occur for a variety of initial values of rotation rate, frozen-in (supercritical) magnetic field, and density-temperature relations. Even in cases where the bursts are nearly entirely suppressed, a moderate increase in cloud size or rotation rate can lead to vigorous burst activity. We conclude that most (if not all) protostars undergo a burst mode of evolution during their early accretion history, as inferred empirically from observations of FU Orionis variables.

The recent downward revision of solar photospheric abundances of Oxygen and other heavy elements has resulted in serious discrepancies between solar models and solar structure as determined through helioseismology. In this work we investigate the possibility of determining the solar heavy-element abundance without reference to spectroscopy by using helioseismic data. Using the dimensionless sound-speed derivative in the solar convection zone, we find that the heavy element abundance, Z, of 0.0172 +/- 0.002, which is closer to the older, higher value of the abundances.

A high-precision two-dimensional stellar evolution code has been developed for studying solar variability due to structural changes produced by varying internal magnetic fields of arbitrary configurations. Specifically, we are interested in modeling the effects of a dynamo-type field on the detailed internal structure and on the global parameters of the Sun. The high precision is required both to model very small solar changes (of order of $10^{-4}$) and short time scales (or order of one year). It is accomplished by using the mass coordinate to replace the radial coordinate, by using fixed and adjustable time steps, a realistic stellar atmosphere, elements diffusion, and by adjusting the grid points. We have also built into the code the potential to subsequently include rotation and turbulence. The current code has been tested for several cases, including its ability to reproduce the 1-D results.

We use numerical simulations to study the effect of nonlinear MHD waves in a stratified, self-gravitating molecular cloud that is bounded by a hot and tenuous external medium. In a previous paper, we had shown the details of a standard model and studied the effect of varying the dimensionless amplitude. In this paper, we present the results of varying two other important free parameters: beta_0, the initial ratio of gas to magnetic pressure at the cloud midplane, and the dimensionless frequency of driving. Furthermore, we present the case of a temporally random driving force. Our results demonstrate that a very important consideration for the actual level of turbulent support against gravity is the ratio of driving wavelength lambda_0 to the the size of the initial non-turbulent cloud; maximum cloud expansion is achieved when this ratio is close to unity. The best consistency with the observational correlation of magnetic field strength, turbulent line width, and density is achieved by cloud models with beta_0 approx 1. We also calculate the spatial power spectra of the turbulent clouds, and show that significant power is developed on scales larger than the scale length H_0 of the initial cloud, even if the input wavelength of turbulence lambda_0 approx H_0. The cloud stratification and resulting increase of Alfven speed toward the cloud edge allows for a transfer of energy to wavelengths significantly larger than lambda_0. This explains why the relevant time scale for turbulent dissipation is the crossing time over the cloud scale rather than the crossing time over the driving scale.

We have evolved 10,000 solar models using 21 input parameters that are randomly drawn for each model from separate probability distributions for every parameter. We use the results of these models to determine the theoretical uncertainties in the predicted surface helium abundance, the profile of the sound speed versus radius, the profile of the density versus radius, the depth of the solar convective zone, the eight principal solar neutrino fluxes, and the fractions of nuclear reactions that occur in the CNO cycle or in the three branches of the p-p chains. We also determine the correlation coefficients of the neutrino fluxes for use in analysis of solar neutrino oscillations. Our calculations include the most accurate available input parameters, including radiative opacity, equation of state, and nuclear cross sections. We incorporate both the recently determined heavy element abundances recommended by Asplund, Grevesse & Sauval (2005) and the older (higher) heavy element abundances recommended by Grevesse & Sauval (1998). We present best-estimates of many characteristics of the standard solar model for both sets of recommended heavy element compositions.

We study numerically the evolution of rotating cloud cores, from the collapse of a magnetically supercritical core to the formation of a protostar and the development of a protostellar disk during the main accretion phase. We find that the disk quickly becomes unstable to the development of a spiral structure similar to that observed recently in AB Aurigae. A continuous infall of matter from the protostellar envelope makes the protostellar disk unstable, leading to spiral arms and the formation of dense protostellar/protoplanetary clumps within them. The growing strength of spiral arms and ensuing redistribution of mass and angular momentum creates a strong centrifugal disbalance in the disk and triggers bursts of mass accretion during which the dense protostellar/protoplanetary clumps fall onto the central protostar. These episodes of clump infall may manifest themselves as episodes of vigorous accretion rate (\ge 10^{-4} M_sun/yr) as is observed in FU Orionis variables. Between these accretion bursts, the protostar is characterized by a low accretion rate (< 10^{-6} M_sun/yr). During the phase of episodic accretion, the mass of the protostellar disk remains less than or comparable to the mass of the protostar.

We show in this paper that the changes of the solar diameter in response to variations of large scale magnetic fields and turbulence are not homologous. For the best current model, the variation at the photospheric level is over 1000 times larger than the variation at a depth of 5 Mm, which is about the level at which f-mode solar oscillations determine diameter variations. This model is supported by observations that indicate larger diameter changes for high degree f-modes than for low degree f-modes, since energy of the former are concentrated at shallower layers than the latter.

We investigate the role of non-isothermality in gravitational collapse and protostellar accretion by explicitly including the effects of molecular radiative cooling, gas-dust energy transfer, and cosmic ray heating in models of spherical hydrodynamic collapse. Isothermal models have previously shown an initial decline in the mass accretion rate \dot{M}, due to a gradient of infall speed that develops in the prestellar phase. Our results show that: (1) in the idealized limit of optically thin cooling, a positive temperature gradient is present in the prestellar phase which effectively cancels out the effect of the velocity gradient, producing a near constant \dot{M} in the early accretion phase; (2) in the more realistic case including cooling saturation at higher densities, \dot{M} may initially be either weakly increasing or weakly decreasing with time, for low (T_d ~ 6 K) and high dust temperature (T_d ~ 10 K) cases, respectively. Hence, our results show that the initial decline in \dot{M} seen in isothermal models is definitely not enhanced by non-isothermal effects, and is often suppressed by them. In all our models, \dot{M} does eventually decline rapidly due to the finite mass condition on our cores and a resulting inward propagating rarefaction wave. Thus, any explanation for a rapid decline of $\dot{M}$ in the accretion phase likely needs to appeal to the global molecular cloud structure and possible envelope support, which results in a finite mass reservoir for cores.

In this article we analyze the radiation loss from a high energy cosmic ray proton propagating in a spacetime with non-systematic Lorentz violation. From an effective field theory perspective we illuminate flaws in previous attempts that use threshold approaches to analyze this problem. We argue that in general such approaches are of rather limited use when dealing with non-systematic Lorentz violating scenarios. The main issues we raise are a) the limited applicability of threshold energy conservation rules when translation invariance is broken and b) the large amounts of proton particle production due to the time dependence of the fluctuations. Ignoring particle production, we derive a constraint on the magnitude of velocity fluctuation $|v_f|<10^{-6.5}$, much weaker than has been previously argued. However, we show that in fact particle production makes any such constraint completely unreliable.

We have evolved a series of thirteen complete solar models that utilize different assumed heavy element compositions. Models that are based upon the heavy element abundances recently determined by Asplund, Grevesse, and Sauval (2005) are inconsistent with helioseismological measurements. However, models in which the neon abundance is increased by 0.4-0.5 dex to log N(Ne) = 8.29 +- 0.05 (on the scale in which log N(H) = 12) are consistent with the helioseismological measurements even though the other heavy element abundances are in agreement with the determinations of Asplund et al. (2005). These results sharpen and strengthen an earlier study by Antia and Basu (2005). The predicted solar neutrino fluxes are affected by the uncertainties in the composition by less than their 1sigma theoretical uncertainties.

Motivated by recent observations which detect an outer boundary for starless cores, and evidence for time-dependent mass accretion in the Class 0 and Class I protostellar phases, we reexamine the case of spherical isothermal collapse in the case of a finite mass reservoir. The presence of a core boundary results in the generation of an inward propagating rarefaction wave. This steepens the gas density profile from r^{-2} to r^{-3} or steeper. After a protostar forms, the mass accretion rate \dot{M} evolves through three distinct phases: (1) an early phase of decline in \dot{M}, which is a non-self-similar effect due to spatially nonuniform infall in the prestellar phase; (2) for large cores, an intermediate phase of near-constant \dot{M} from the infall of the outer part of the self-similar density profile; (3) a late phase of rapid decline in \dot{M} when accretion occurs from the region affected by the inward propagating rarefaction wave. Our model clouds of small to intermediate size make a direct transition from phase (1) to phase (3) above. Both the first and second phase are characterized by a temporally increasing bolometric luminosity L_bol, while L_bol is decreasing in the third (final) phase. We identify the period of temporally increasing L_bol with the Class 0 phase, and the later period of terminal accretion and decreasing L_bol with the Class I phase. The peak in L_bol corresponds to the evolutionary time when 50% \pm 10% of the cloud mass has been accreted by the protostar. This is in agreement with the classification scheme proposed by Andre et al. (1993). We show how our results can be used to explain tracks of envelope mass M_env versus L_bol for protostars in Taurus and Ophiuchus. We also develop an analytic formalism which reproduces the protostellar accretion rate.

We construct solar models with the newly calculated radiative opacities from the Opacity Project (OP) and recently determined (lower) heavy element abundances. We compare results from the new models with predictions of a series of models that use OPAL radiative opacities, older determinations of the surface heavy element abundances, and refinements of nuclear reaction rates. For all the variations we consider, solar models that are constructed with the newer and lower heavy element abundances advocated by Asplund et al. (2005) disagree by much more than the estimated measuring errors with helioseismological determinations of the depth of the solar convective zone, the surface helium composition, the internal sound speeds, and the density profile. Using the new OP radiative opacities, the ratio of the 8B neutrino flux calculated with the older and larger heavy element abundances (or with the newer and lower heavy element abundances) to the total neutrino flux measured by the Sudbury Neutrino Observatory is 1.09 (0.87) with a 9% experimental uncertainty and a 16% theoretical uncertainty, 1 sigma errors.

There have been recent downward revisions of the solar photospheric abundances of Oxygen and other heavy elements. These revised abundances along with OPAL opacities are not consistent with seismic constraints. In this work we show that the recently released OP opacity tables cannot resolve this discrepancy either. While the revision in opacities does not seem to resolve this conflict, an upward revision of Neon abundance in solar photosphere offers a possible solution to this problem.

We perform numerical hydrodynamic modeling of various physical processes that can form an HI ring as is observed in Holmberg I. Three energetic mechanisms are considered: multiple supernova explosions (SNe), a hypernova explosion associated with a gamma ray burst (GRB), and the vertical impact of a high velocity cloud (HVC). The total released energy has an upper limit of 10^54 ergs. We find that multiple SNe are in general more effective in producing shells that break out of the disk than a hypernova explosion of the same total energy. As a consequence, multiple SNe form rings with a high ring-to-center contrast K<100 in the HI column density, whereas single hypernova explosions form rings with K<10. Only multiple SNe can reproduce both the size (diameter \~1.7 kpc) and the ring-to-center contrast (K ~ 15-20) of the HI ring in Hoolmberg I. High velocity clouds create HI rings that are much smaller in size (< 0.8 kpc) and contrast (K < 4.5) than seen in Holmberg I. We construct model position-velocity (pV) diagrams and find that they can be used to distinguish among different HI ring formation mechanisms. The observed pV-diagrams of Holmberg I are best reproduced by multiple SNe. We conclude that the giant HI ring in Holmberg I is most probably formed by multiple SNe. We also find that the appearance of the SNe-driven shell in the integrated HI image depends on the inclination angle of the galaxy. In nearly face-on galaxies, the integrated HI image shows a ring of roughly constant HI column density surrounding a deep central depression, whereas in considerably inclined galaxies (i > 45 deg) the HI image is characterized by two kidney-shaped density enhancements and a mild central depression.

We discuss several categories of models which may explain the IMF, including the possible role of turbulence and magnetic fields.

We present a minimum hypothesis model for an IMF that resembles a lognormal distribution at low masses but has a distinct power-law tail. Even if the central limit theorem ensures a lognormal distribution of condensation masses at birth, a power-law tail in the distribution arises due to accretion from the ambient cloud, coupled with a non-uniform (exponential) distribution of accretion times.

We present evidence that solar structure changes with changes in solar activity. We find that the adiabatic index, Gamma_1, changes near the second helium ionization, i.e., at a depth of about 0.98 R_sun. We believe that this change is a result of the change in the effective equation of state caused by magnetic fields. Inversions should be able to detect the changes in Gamma_1 if mode sets with reliable and precise high-degree modes are available.

Observations of magnetic field strengths imply that molecular cloud fragments are individually close to being in a magnetically critical state, even though both magnetic field and column density measurements range over two orders of magnitude. The turbulent pressure also approximately balances the self-gravitational pressure. These results together mean that the one-dimensional velocity dispersion $\sigv$ is proportional to the mean \Alf speed of a cloud $\va$. Global models of MHD turbulence in a molecular cloud show that this correlation is naturally satisfied for a range of different driving strengths of the turbulence. For example, an increase of turbulent driving causes a cloud expansion which also increases $\va$. Clouds are in a time averaged balance but exhibit large oscillatory motions, particularly in their outer rarefied regions. We also discuss models of gravitational fragmentation in a sheet-like region in which turbulence has already dissipated, including the effects of magnetic fields and ion-neutral friction. Clouds with near-critical mass-to-flux ratios lead to subsonic infall within cores, consistent with some recent observations of motions in starless cores. Conversely, significantly supercritical clouds are expected to produce extended supersonic infall.

We show that standard solar models are in good agreement with the helioseismologically determined sound speed and density as a function of solar radius, the depth of the convective zone, and the surface helium abundance, as long as those models do not incorporate the most recent heavy element abundance determinations. However, sophisticated new analyses of the solar atmosphere infer lower abundances of the lighter metals (like C, N, O, Ne, and Ar) than the previously widely used surface abundances. We show that solar models that include the lower heavy element abundances disagree with the solar profiles of sound speed and density as well as the depth of the convective zone and the helium abundance. The disagreements for models with the new abundances range from factors of several to many times the quoted uncertainties in the helioseismological measurements. The disagreements are at temperatures below what is required for solar interior fusion reactions and therefore do not significantly affect solar neutrino emission. If errors in thecalculated OPAL opacities are solely responsible for the disagreements, then the corrections in the opacity must extend from 2 times 10^6 K (R = 0.7R_Sun)to 5 times 10^6 K (R = 0.4 R_Sun), with opacity increases of order 10%.

We propose that the mushroom-shaped structure of the Galactic worm GW 123.4--1.5 is created by a cloud collision with the Galactic gas disk. A hydrodynamic simulation shows that a mushroom-shaped structure is created after the cloud crosses the Galactic midplane. The lifetime of the mushroom-shaped structure is of order the dynamical time scale of the disk, \sim 10^7 years. We find that the velocities across the cap of the mushroom-shaped structure in the simulation are consistent with the observed values. The simulation also predicts a structure on the opposite side of the Galactic plane which is created by the Kelvin-Helmholtz instability after the cloud passes through the disk.

We extend our earlier work on ambipolar diffusion induced formation of protostellar cores in isothermal sheet-like magnetic interstellar clouds, by studying nonaxisymmetric collapse for the physically interesting regime of magnetically critical and supercritical model clouds ($\mui \geq 1$, where $\mui$ is the initial mass-to-magnetic flux ratio in units of the critical value for gravitational collapse). Cores that form in model simulations are effectively triaxial, with shapes that are typically closer to being oblate, rather than prolate. Infall velocities in the critical model ($\mui = 1$) are subsonic; in contrast, a supercritical model ($\mui = 2$) has extended supersonic infall that may be excluded by observations. For the magnetically critical model, ambipolar diffusion forms cores that are supercritical ($\muc > 1$) and embedded within subcritical envelopes ($\muenv < 1$). Cores in our models have density profiles that eventually merge into a near-uniform background, which is suggestive of observed properties of cloud cores.

Recent analyses of solar photospheric abundances suggest that the oxygen abundance in the solar atmosphere needs to be revised downwards. In this study we investigate the consequence of this revision on helioseismic analyses of the depth of the solar convection zone and the helium abundance in the solar envelope and find no significant effect. We also find that the revised abundances along with the current OPAL opacity tables are not consistent with seismic data. A significant upward revision of the opacity tables is required to make solar models with lower oxygen abundance consistent with seismic observations.

Intermediate degree modes of the solar oscillations have previously been used to determine the solar helium abundance to a high degree of precision. However, we cannot expect to observe such modes in other stars. In this work we investigate whether low degree modes that should be available from space-based asteroseismology missions can be used to determine the helium abundance, Y, in stellar envelopes with sufficient precision. We find that the oscillatory signal in the frequencies caused by the depression in \Gamma_1 in the second helium ionisation zone can be used to determine the envelope helium abundance of low mass main sequence stars. For frequency errors of 1 part in 10^4, we expect errors \sigma_Y in the estimated helium abundance to range from 0.03 for 0.8M_sun stars to 0.01 for 1.2M_sun stars. The task is more complicated in evolved stars, such as subgiants, but is still feasible if the relative errors in the frequencies are less than 10^{-4}.

We explore the idea that the power-law tail in the mass function of protostellar condensations and stars arises from the accretion of ambient cloud material on to a condensation, coupled with a nonuniform (exponential) distribution of accretion lifetimes. This model allows for the generation of power-law distributions in all star-forming regions, even if condensations start with a lognormal mass distribution, as may be expected from the central limit theorem, and supported by some recent numerical simulations of turbulent molecular clouds. For a condensation mass $m$ with growth rate $dm/dt \propto m$, an analytic three-parameter probability density function is derived; it resembles a lognormal at low mass and has a pure power-law high-mass tail. An approximate power-law tail is also expected for other growth laws, and we calculate the distribution for the plausible case $dm/dt \propto m^{2/3}$. Furthermore, any single time snapshot of the masses of condensations that are still accreting (and are of varying ages) also yields a distribution with a power-law tail similar to that of the IMF.

We perform numerical simulations of nonlinear MHD waves in a gravitationally stratified molecular cloud that is bounded by a hot and tenuous external medium. We study the relation between the strength of the turbulence and various global properties of a molecular cloud, within a 1.5-dimensional approximation. Under the influence of a driving source of Alfvenic disturbances, the cloud is lifted up by the pressure of MHD waves and reaches a steady-state characterized by oscillations about a new time-averaged equilibrium state. The nonlinear effect results in the generation of longitudinal motions and many shock waves; however, the wave kinetic energy remains predominantly in transverse, rather than longitudinal, motions. There is an approximate equipartition of energy between the transverse velocity and fluctuating magnetic field (aspredicted by small-amplitude theory) in the region of the stratified cloud which contains most of the mass; however, this relation breaks down in the outer regions, particularly near the cloud surface, where the motions have a standing-wave character. This means that the Chandrasekhar-Fermi formula applied to molecular clouds must be significantly modified in such regions. Models of an ensemble of clouds show that, for various strengths of the input energy, the velocity dispersion in the cloud $\sigma \propto Z^{0.5}$, where $Z$ is a characteristic size of the cloud.Furthermore, $\sigma$ is always comparable to the mean Alfven velocity of the cloud, consistent with observational results.

Data obtained by the GONG and MDI instruments over the last seven years are used to study how solar dynamics -- both rotation and other large scale flows -- have changed with time. In addition to the well known phenomenon of bands of faster and slower rotation moving towards the equator and pole, we find that the zonal flow pattern rises upwards with time. Like the zonal flows, the meridional flows also show distinct solar activity related changes. In particular, the anti-symmetric component of the meridional flow shows a decrease in speed with activity. We do not see any significant temporal variations in the dynamics of the tachocline region where the solar dynamo is believed to be operating.

The Sun is believed to have been the recipient of a substantial amount of metal-rich material over the course of its evolution, particularly in the early stages of the Solar System. With a long diffusion timescale, the majority of this accreted matter should still exist in the solar convection zone, enhancing its observed surface abundance, and implying a lower-abundance core. While helioseismology rules out solar models with near-zero metallicity cores, some solar models with enhanced metallicity in the convection zone might be viable, as small perturbations to the standard model. Because of the reduced interior opacity and core temperature, the neutrino flux predicted for such models is lower than that predicted by the standard solar model. This paper examines how compatible inhomogeneous solar models of this kind are with the observed low and intermediate degree p-mode oscillation data, and with the solar neutrino data from the SNO Collaboration. We set an upper limit on how much metal-rich accretion took place during the early evolution of the Sun at about 2 Earth masses of iron (or about 40 Earth masses of meteoric material).