Parameter Optimisation for Woodwind Single-Reed Models

Time-domain modelling of single-reed woodwind instruments usually involves a lumped model of the excitation mechanism. The parameters of this lumped model have to be estimated for use in numerical simulations. Several attempts have been made to estimate these parameters, including observations of the mechanics of isolated reeds, measurements under artificial or real playing conditions and estimations based on numerical simulations. In this study an optimisation routine is presented, that can estimate reed-model parameters, given the pressure and flow signals in the mouthpiece. The method is validated, tested on a series of numerically synthesised data. In order to incorporate the actions of the player in the parameter estimation process, the optimisation routine has to be applied to signals obtained under real playing conditions. The estimated parameters can then be used to resynthesise the pressure and flow signals in the mouthpiece. In the case of measured data, as opposed to numerically synthesised data, special care needs to be taken while modelling the bore of the instrument. In fact, a careful study of various experimental datasets revealed that for resynthesis to work, the bore termination impedance should be known very precisely from theory. An example is given, where the above requirement is satisfied, and the resynthesised signals closely match the original signals generated by the player.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

Atmospheric parameters and rotational velocities for a sample of Galactic B-type supergiants

High-resolution optical spectra of 57 Galactic B-type supergiant stars have been analysed to determine their rotational and macroturbulent velocities. In addition, their atmospheric parameters (effective temperature, surface gravity and microturbulent velocity) and surface nitrogen abundances have been estimated using a non-local thermodynamic equilibrium grid of model atmospheres. Comparisons of the projected rotational velocities have been made with the predictions of stellar evolutionary models and in general good agreement was found. However, for a small number of targets, their observed rotational velocities were significantly larger than predicted, although their nitrogen abundances were consistent with the rest of the sample. We conclude that binarity may have played a role in generating their large rotational velocities. No correlation was found between nitrogen abundances and the current projected rotational velocities. However, a correlation was found with the inferred projected rotational velocities of the main-sequence precursors of our supergiant sample. This correlation is again in agreement with the predictions of single star evolutionary models that incorporate rotational mixing. The origin of the macroturbulence and microturbulent velocity fields is discussed and our results support previous theoretical studies that link the former to subphotospheric convection and the latter to non-radial gravity mode oscillations. In addition, we have attempted to identify differential rotation in our most rapidly rotating targets.

Optical and near-infrared coverage of SN 2004et: physical parameters and comparison with other Type IIP supernovae

We present new optical and near-infrared (NIR) photometry and spectroscopy of the Type IIP supernova (SN), SN 2004et. In combination with already published data, this provides one of the most complete studies of optical and NIR data for any Type IIP SN from just after explosion to +500 d. The contribution of the NIR flux to the bolometric light curve is estimated to increase from 15 per cent at explosion to around 50 per cent at the end of the plateau and then declines to 40 per cent at 300 d. SN 2004et is one of the most luminous IIP SNe which has been well studied and characterized, and with a luminosity of log L = 42.3 erg s-1 and a 56Ni mass of 0.06 +/- 0.04 M-circle dot, it is two times brighter than SN 1999em. We provide parametrized bolometric corrections as a function of time since explosion for SN 2004et and three other IIP SNe that have extensive optical and NIR data. These can be used as templates for future events in optical and NIR surveys without full wavelength coverage. We compare the physical parameters of SN 2004et with those of other well-studied IIP SNe and find that the kinetic energies span a range of 1050-1051 erg. We compare the ejected masses calculated from hydrodynamic models with the progenitor masses and limits derived from pre-discovery images. Some of the ejected mass estimates are significantly higher than the progenitor mass estimates, with SN 2004et showing perhaps the most serious mass discrepancy. With the current models, it appears difficult to reconcile 100 d plateau lengths and high expansion velocities with the low ejected masses of 5-6 M-circle dot implied from 7-8 M-circle dot progenitors. The nebular phase is studied using very late-time Hubble Space Telescope photometry, along with optical and NIR spectroscopy. The light curve shows a clear flattening at 600 d in the optical and the NIR, which is likely due to the ejecta impacting on circumstellar material. We further show that the [O i] 6300, 6364 A line strengths in the nebular spectra of four Type IIP SNe imply ejected oxygen masses of 0.5-1.5 M-circle dot.

DiveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors

Summary: We present a new R package, diveRsity, for the calculation of various diversity statistics, including common diversity partitioning statistics (?, G) and population differentiation statistics (D, GST ', ? test for population heterogeneity), among others. The package calculates these estimators along with their respective bootstrapped confidence intervals for loci, sample population pairwise and global levels. Various plotting tools are also provided for a visual evaluation of estimated values, allowing users to critically assess the validity and significance of statistical tests from a biological perspective. diveRsity has a set of unique features, which facilitate the use of an informed framework for assessing the validity of the use of traditional F-statistics for the inference of demography, with reference to specific marker types, particularly focusing on highly polymorphic microsatellite loci. However, the package can be readily used for other co-dominant marker types (e.g. allozymes, SNPs). Detailed examples of usage and descriptions of package capabilities are provided. The examples demonstrate useful strategies for the exploration of data and interpretation of results generated by diveRsity. Additional online resources for the package are also described, including a GUI web app version intended for those with more limited experience using R for statistical analysis. © 2013 British Ecological Society.

Test-retest variability in structural and functional parameters of glaucoma damage in the glaucoma imaging longitudinal study.

PURPOSE:

To determine the test-retest variability in perimetric, optic disc, and macular thickness parameters in a cohort of treated patients with established glaucoma.

PATIENTS AND METHODS:

In this cohort study, the authors analyzed the imaging studies and visual field tests at the baseline and 6-month visits of 162 eyes of 162 participant in the Glaucoma Imaging Longitudinal Study (GILS). They assessed the difference, expressed as the standard error of measurement, of Humphrey field analyzer II (HFA) Swedish Interactive Threshold Algorithm fast, Heidelberg retinal tomograph (HRT) II, and retinal thickness analyzer (RTA) parameters between the two visits and assumed that this difference was due to measurement variability, not pathologic change. A statistically significant change was defined as twice the standard error of measurement.

RESULTS:

In this cohort of treated glaucoma patients, it was found that statistically significant changes were 3.2 dB for mean deviation (MD), 2.2 for pattern standard deviation (PSD), 0.12 for cup shape measure, 0.26 mm for rim area, and 32.8 microm and 31.8 microm for superior and inferior macular thickness, respectively. On the basis of these values, it was estimated that the number of potential progression events detectable in this cohort by the parameters of MD, PSD, cup shape measure, rim area, superior macular thickness, and inferior macular thickness was 7.5, 6.0, 2.3, 5.7, 3.1, and 3.4, respectively.

CONCLUSIONS:

The variability of the measurements of MD, PSD, and rim area, relative to the range of possible values, is less than the variability of cup shape measure or macular thickness measurements. Therefore, the former measurements may be more useful global measurements for assessing progressive glaucoma damage.

Influence of corneal structure, corneal responsiveness, and other ocular parameters on tonometric measurement of intraocular pressure.

PURPOSE: To estimate the relationships between ocular parameters and tonometrically measured intraocular pressure (IOP), to determine the influence of ocular parameters on different instrument measurements of IOP, and to evaluate the association of ocular parameters with a parameter called hysteresis. METHODS: Patients presenting at a glaucoma clinic were recruited for this study. Subjects underwent IOP measurement with the Goldmann applanation tonometer (GAT), the TonoPen, and the Reichert Ocular Response Analyzer (ORA), and also measurements of central corneal thickness (CCT), axial length, corneal curvature, corneal astigmatism, central visual acuity, and refractive error. Chart information was reviewed to determine glaucoma treatment history. The ORA instrument provided a measurement called corneal hysteresis. The association between measured IOP and the other ocular characteristics was estimated using generalized estimating equations. RESULTS: Among 230 patients, IOP measurements from the TonoPen read lowest, and ORA read highest, and GAT measurements were closest to the mean IOP of the 3 instruments. In a multiple regression model adjusting for age, sex, race, and other ocular characteristics, a 10 microm increase in CCT was associated with an increase of 0.79 mm Hg measured IOP in untreated eyes (P<0.0001). Of the 3 tonometers, GAT was the least affected by CCT (0.66 mm Hg/10 mum, P<0.0001). Hysteresis was significantly correlated with CCT with a modest correlation coefficient (r=0.20, P<0.0007). CONCLUSIONS: Among parameters related to measured IOP, features in addition to CCT, such as hysteresis and corneal curvature, may also be important. Tonometric instruments seem to be affected differently by various physiologic characteristics.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Context: Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. 

Aims: Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLT-FLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. 

Methods: tlusty non-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The effects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s^-1 with 65% in the range 30 km s^-1 ≤ ν_e sin i ≤ 60 km s^-1. About ten per cent have larger ve sin i (≥100 km s^-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km s^-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. 

Conclusions: The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to T_eff = 20 000 K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.