From analytical studies of tidal heating, eclipses and planetary illumination, it is clear that the exomoon habitable zone (EHZ) - the set of moon and host planet orbits that permit liquid water on an Earthlike moon's surface - is a manifold of higher dimension than the planetary HZ. This paper outlines the first attempt to produce climate models of exomoons which possess all the above sources and sinks of energy. We expand on our previous 1D latitudinal energy balance models (LEBMs), which follow the evolution of the temperature on an Earthlike moon orbiting a Jupiterlike planet, by adding planetary illumination. We investigate the EHZ in four dimensions, running two separate suites of simulations. The first investigates the EHZ by varying the planet's orbit, keeping the moon's orbit fixed, to compare the EHZ with planetary habitable zones. In general, planetary illumination pushes EHZs slightly further away from the star. Secondly, we fix the planet's orbit and vary the moon's orbit, to investigate the circumplanetary inner habitable edge. We demonstrate that an outer edge can exist due to eclipses (rather than merely orbital stability), but this edge may be pushed outwards when the effect of the carbonate-silicate cycle is taken into account.

We report our investigation of the first transiting planet candidate from the YETI project in the young (~4 Myr old) open cluster Trumpler 37. The transit-like signal detected in the lightcurve of the F8V star 2M21385603+5711345 repeats every 1.364894+/-0.000015 days, and has a depth of 54.5+/-0.8 mmag in R. Membership to the cluster is supported by its mean radial velocity and location in the color-magnitude diagram, while the Li diagnostic and proper motion are inconclusive in this regard. Follow-up photometric monitoring and adaptive optics imaging allow us to rule out many possible blend scenarios, but our radial-velocity measurements show it to be an eclipsing single-lined spectroscopic binary with a late-type (mid-M) stellar companion, rather than one of planetary nature. The estimated mass of the companion is 0.15-0.44 solar masses. The search for planets around very young stars such as those targeted by the YETI survey remains of critical importance to understand the early stages of planet formation and evolution.

We present a new solution of the direct problem of planet transits based on transformation of double integrals to single ones. On the basis of our direct problem solution we created the code TAC-maker for rapid and interactive calculation of synthetic planet transits by numerical computations of the integrals. The validation of our approach was made by comparison with the results of the wide-spread Mandel & Agol (2002) method for the cases of linear, quadratic and squared root limb-darkening laws and various combinations of model parameters. For the first time our approach allows the use of arbitrary limb-darkening law of the host star. This advantage together with the practically arbitrary precision of the calculations make the code a valuable tool that faces the challenges of the continuously increasing photometric precision of the ground-based and space observations.