An update is given on the exoplanet research collaboration between Nielsen (a marketing research company), Brigham Young University, and NZ universities with the National University of Singapore, which has been expanded to include a community college in the US. Key achievements from the past year are outlined, including density estimates for HD 209458 and Kepler 1 from radial velocity and transit fits. A comparison between the WinFitter optimizer and other techniques is outlined, showing that WinFitter estimated statistical errors are essentially in line (bar a scaling proportion) with those estimated via Markov Chain Monte Carlo techniques.

We have carried out an intensive study of photometric (Kepler Mission) and spectroscopic data on the system Kepler-2 (HAT-P-7A) using the dedicated software WinFitter. The mean individual data-point error of the normalized flux values for this system is 0.00015, leading to the model's specification for the mean reference flux to an accuracy of $\sim$0.5 ppm. This testifies to the remarkably high accuracy of the binned data-set, derived from over 1.8 million individual observations. Spectroscopic data are reported with the similarly high-accuracy radial velocity amplitude measure of $\sim$2 m s$^{-1}$. The analysis includes discussion of the fitting quality and model adequacy. Our derived absolute parameters for Kepler-2 are as follows: $M_p$ (Jupiter) 1.80 $\pm$ 0.13; $R_{\star}$ 1.46 $\pm 0.08 \times 10^6$ km; $R_p$ 1.15 $\pm 0.07 \times 10^5$ km. These values imply somewhat larger and less condensed bodies than previously catalogued, but within reasonable error estimates of such literature parameters. We find also tidal, reflection and Doppler effect parameters, showing that the optimal model specification differs slightly from a `cleaned' model that reduces the standard deviation of the $\sim$3600 binned light curve points to less than 0.9 ppm. We consider these slight differences, making comparisons with the hot-jupiter systems Kepler-1 (TrES-2) and 13.

We present an updated O-C diagram of the light-time variations of the eclipsing binary (component B) in the system QZ Carinae as it moves in the long-period orbit around the non-eclipsing pair (component A). This includes new Variable Stars South members' measures from 2017 to 2019, BRITE satellite observations in 2017 and 2018, and 100 previously unpublished measures made at Auckland Observatory from 1974 to 1978. We conclude that QZ Carinae has not yet completed one orbit of the two pairs since discovery in 1971. The duration of totality of primary eclipses was measured to be 0.295 +/- 0.02 day (7.08 +/-0.48 hours), rather longer than earlier values from light curve models. Other observational findings include the shape of primary and secondary eclipses and small-scale short-term brightness changes.

A brief review is given of a university outreach programme by a commercial organisation, which uses the Kepler exoplanet data. Key insights derived from this research are presented, along with discussion of the benefi ts and challenges of such a collaboration between industry and academia. It is hoped that this account will be an inspiring example for others to emulate.

This paper employs a simple model, considering just geometry and linear or quadratic limb darkening, to fit Kepler transit data via a Markov Chain Monte Carlo (MCMC) methodology for Kepler-1b, 5b, 8b, 12b, 77b, 428b, 491b, 699b, 706b, and 730b. Additional fits were made of the systems using the more sophisticated modeller Winfitter, which gives results in general agreement with the simpler model. Analysis of data with longer integration times showed biasing of the derived parameters, as expected from the literature,} leading to larger estimates for radii and reducing estimates of the system inclination.

A simple transit model is described along with tests of this model against published results for 4 exoplanet systems (Kepler-1, 2, 8, and 77). Data from the Kepler mission are used. The Markov Chain Monte Carlo (MCMC) method is applied to obtain realistic error estimates. Optimization of limb darkening coefficients is subject to data quality. It is more likely for MCMC to derive an empirical limb darkening coefficient for light curves with S/N (signal to noise) above 15. Finally, the model is applied to Kepler data for 4 Kepler candidate systems (KOI 760.01, 767.01, 802.01, and 824.01) with previously unpublished results. Error estimates for these systems are obtained via the MCMC method.

New photometric follow-up observations of exoplanet transits of 'hot Jupiters' TrES-3b and Qatar-1b are presented. Weighted mean values of the solutions of light curves in R-fiter for both planetary systems are reported and compared with the previous results. The transit light curves were analysed using the WINFITTER code. The physical properties of the planets were estimated. The planet radii are found to be $R_p = 1.381 \pm 0.033 R_J$ for TrES-3b and $R_p = 1.142 \pm 0.025 R_J$ for Qatar-1b. Transit times and their uncertainties were also determined and a new linear ephemeris was computed for both systems. Analysis of transit times showed that a significant signal could not be determined for TrES-3b, while weak evidence was found for Qatar-1b, which might be tested using more precise future transit times.

We have applied the graphical user interfaced close binary system analysis program WinFitter to an intensive study of Kepler 91 using all the available photometry from the NASA Exoplanet Archive (NEA) at the Caltech website: http://exoplanetarchive.ipac.caltech.edu . Our fitting function for the tidal distortion derives from the relevant Radau equation and includes terms up to the fifth power of the fractional radius. This results in a systematic improvement in the mass ratio estimation over that of Lillo Box et al (2014a) and our derived value for the mass ratio is in close agreement with that inferred from recent high-resolution spectroscopic data. It is clear that the data analysis in terms of simply an eclipsing binary system is compromised by the presence of significant other causes of light variation, in particular non-radial pulsations. We apply a low-frequency filtering procedure to separate out some of this additional light variation. Whilst the derived eccentricity appears then reduced, an eccentric effect remains in the light curve. We consider how this may be maintained in spite of likely frictional effects operating over a long time. There are also indications of the possible effects of Trojan or other period-resonant mass concentrations. Suggestions of a possible secular period variation are briefly discussed.

New high-resolution spectroscopy and $BVR$ photometry, together with literature data, on the Gould's Belt close binary systems GG Lup and $\mu^1$ Sco are presented and analysed. In the case of GG Lup, light and radial velocity curve fittings confirm a near-Main-Sequence picture of a pair of close stars. Absolute parameters are found, to within a few percent, thus: $M_1$ = 4.16$\pm$0.12, $M_2$ = 2.64$\pm$0.12, $R_{1}$ = 2.42$\pm$0.05, $R_2$ = 1.79$\pm$0.04, ($\odot$); $T_{1}$ $\sim$13000, $T_2$ $\sim$10600 (K); photometric distance $\sim$ 160 (pc). The high eccentricity and relatively short period (105 y) of apsidal revolution may be related to an apparent `slow B-type pulsator' (SPB) oscillation. Disturbances of the outer envelope of at least one of the components then compromise comparisons to standard evolutionary models, at least regarding the age of the system. A rate of apsidal advance is derived, which allows a check on the mean internal structure constant $\bar{k_2} = 0.0058 \pm 0.0004$. This is in agreement with values recently derived for young stars of solar composition and mass $\sim$3${\odot}$. For $\mu^1$ Sco, we agree with previous authors that the secondary component is considerably oversized for its mass, implying binary (interactive) stellar evolution, probably of the `Case A' type. The primary appears relatively little affected by this evolution, however. Its parameters show consistency with a star of its derived mass at age about 13 My, consistent with the star's membership of the Sco-Cen OB2 Association. The absolute parameters are as follows: $M_1$ = 8.3$\pm$1.0, $M_2$ = 4.6$\pm$1.0, $R_{1}$ = 3.9$\pm$0.3, $R_2$ = 4.6$\pm$0.4, ($\odot$); $T_{1}$ $\sim$24000, $T_2$ $\sim$17000 (K); photometric distance $\sim$ 135 (pc).

We review the results of the 1988 multi-wavelength campaign on the late-type eclipsing binary YY Geminorum. Observations include: broad-band optical and near infra-red photometry, simultaneous optical and ultraviolet (IUE) spectroscopy, X-ray (Ginga) and radio (VLA) data. From models fitted to the optical light curves, fundamental physical parameters have been determined together with evidence for transient maculations (spots) located near quadrature longitudes and intermediate latitudes. Eclipses were observed at optical, ultraviolet and radio wavelengths. Significant drops in 6cm radio emission near the phases of both primary and secondary eclipse indicate relatively compact radio emitting volumes that may lie between the binary components. IUE observations during secondary eclipse are indicative of a uniform chromosphere saturated with MgII plage-type emission and an extended volume of Ly$\alpha$ emission. Profile fitting of high-dispersion H alpha spectra confirms the chromospheric saturation and indicates significant H$\alpha$ opacity to heights of a few percent of the photospheric radius. There is evidence for an enhanced H alpha emission region visible near phase 0.25-0.35 which may be associated with a large spot on the primary and with two small optical flares which were also observed at other wavelengths: one in microwave radiation and the other in X-rays. For both flares, L_X/L_opt is consistent with energy release in closed magnetic structures.

We have observed the bright, magnetically active multiple star AB Doradus in a multiwavelength campaign centring around two large facility allocations in November 2006 and January, 2007. Our observations have covered at least three large flares. These flares were observed to produce significant hardening of the X-ray spectra during their very initial stages. We monitored flare-related effects using the Suzaku X-ray satellite and the Australia Telescope Compact Array at 3.6 and 6 cm. Observations at 11 and 21 cm were also included, but they were compromised by interference. From our multiwavelength coverage we find that the observed effects can be mainly associated with a large active region near longitude zero. The second major X-ray and microwave flare of Jan 8, 2007 was observed with a favourable geometry that allowed its initial high-energy impulsive phase to be observed in the higher frequency range of Suzaku's XIS detectors. The fractional circular polarisation was measured for the complete runs, for 25 min integrations and, at 4.80 GHz, for 5 min integrations. Most of the full data sets showed circular polarisation fractions from AB Dor B that were significant at greater than the 3 sigma level. In several of the 5 min integrations at 4.80 and 8.64 GHz this fraction reached a significance level between 3 and 9 sigma. Lack of angular resolution prevented identification of these high V/I values with one or other of the two low-mass red-dwarf components of AB Dor B.

Aims: We aim to derive the absolute parameters of the components of AD And, AL Cam, and V338 Her, interpret their orbital period changes and discuss their evolutionary status. Methods: New and complete multi-filter light curves of the eclipsing binaries AD And, AL Cam, and V338 Her were obtained and analysed with modern methods. Using all reliably observed times of minimum light, we examined orbital period irregularities using the least squares method. In addition, we acquired new spectroscopic observations during the secondary eclipses for AL Cam and V338 Her. Results: For AL Cam and V338 Her, we derive reliable spectral types for their primary stars. Statistical checks of orbital period analysis for all systems are very reassuring in the cases of V338 Her and AD And, although less so for AL Cam. The LIght-Time Effect (LITE) results are checked by inclusion of a third light option in the photometric analyses. Light curve solutions provide the means to calculate the absolute parameters of the components of the systems and reliably estimate their present evolutionary status. Conclusions: AL Cam and V338 Her are confirmed as classical Algols of relatively low mass in similar configurations. Unlike AL Cam, however, V338 Her is still transferring matter between its components, raising interest in the determinability of the evolutionary histories of Algols. AD And is found to be a detached system in which both close stars are of age ~109 yr and is probably a "non-classical" young triple, at an interesting stage of its dynamical evolution.

We have modified the graphical user interfaced close binary system analysis program CurveFit to the form WinKepler and applied it to 16 representative planetary candidate light curves found in the NASA Exoplanet Archive (NEA) at the Caltech website http://exoplanetarchive.ipac.caltech.edu, with an aim to compare different analytical approaches. WinKepler has parameter options for a realistic physical model, including gravity-brightening and structural parameters derived from the relevant Radau equation. We tested our best-fitting parameter-sets for formal determinacy and adequacy. A primary aim is to compare our parameters with those listed in the NEA. Although there are trends of agreement, small differences in the main parameter values are found in some cases, and there may be some relative bias towards a 90 degrees value for the NEA inclinations. These are assessed against realistic error estimates. Photometric variability from causes other than planetary transits affects at least 6 of the data-sets studied; with small pulsational behaviour found in 3 of those. For the false positive KOI 4.01, we found that the eclipses could be modelled by a faint background classical Algol as effectively as by a transiting exoplanet. Our empirical checks of limb-darkening, in the cases of KOI 1.01 and 12.01, revealed that the assigned stellar temperatures are probably incorrect. For KOI 13.01, our empirical mass-ratio differs by about 7% from that of Mislis and Hodgkin (2012), who neglected structural effects and higher order terms in the tidal distortion. Such detailed parameter evaluation, additional to the usual main geometric ones, provides an additional objective for this work.

CCD photometric observations of the Algol-type eclipsing binary AT Peg have been obtained. The light curves are analyzed with modern techniques and new geometric and photometric elements are derived. A new orbital period analysis of the system, based on the most reliable timings of minima found in the literature, is presented and apparent period modulations are discussed with respect to the Light-Time effect (LITE) and secular changes in the system. The results of these analyses are compared and interpreted in order to obtain a coherent view of the system's behaviour.

A search for ultra-high energy gamma-rays emitted by the young, nearby supernova remnant that was discovered recently by the COMPTEL and ROSAT satellites was made using the JANZOS database for the period 1987-1993. A 95% confidence upper limit on the flux above 100 TeV of 3x10^(-13)cm^(-2)sec^(-1) was obtained. This is an order of magnitude below the expected flux based on the standard model of cosmic ray acceleration in supernova shocks. An optical survey of the region that has been commenced is also reported. This uses UK and ESO Schmidt plates, and CCD images by a NZ/Japan microlensing group.