HW3-Problem 3 (Jocelyn Quispe)

The inserted parameters:

  • Event Midpoint JD UT (days) =  2458186.69173
  • Period (days) = 1.327347000
Formula used to find next known transit:
T_transit (JD UT) = JD + P *n 

I used a online converter to change the JD back to Calendar date and time. The boxes highlighted in green at the night where it would be visible for us to see. So out of the 50 transits, only 19 would occur at night on campus. 

Planet
Event Midpoint JD UT (days)
Period (days)
n
T_transit (JD UT days)
T_transit (Calendar UT)
HAT-P 36
2458186.69173
1.327347000
1
2458188.019077
2018-Mar-10 12:27:28.25



2
2458189.346424
2018-Mar-11 20:18:51.03



3
2458190.673771
2018-Mar-13 04:10:13.81



4
2458192.001118
2018-Mar-14 12:01:36.6



5
2458193.328465
2018-Mar-15 19:52:59.38



6
2458194.655812
2018-Mar-17 03:44:22.16



7
2458195.983159
2018-Mar-18 11:35:44.94



8
2458197.310506
2018-Mar-19 19:27:07.72



9
2458198.637853
2018-Mar-21 03:18:30.5



10
2458199.9652
2018-Mar-22 11:09:53.28



11
2458201.292547
2018-Mar-23 19:01:16.06



12
2458202.619894
2018-Mar-25 02:52:38.84



13
2458203.947241
2018-Mar-26 10:44:01.62



14
2458205.274588
2018-Mar-27 18:35:24.4



15
2458206.601935
2018-Mar-29 02:26:47.18



16
2458207.929282
2018-Mar-30 10:18:09.96



17
2458209.256629
2018-Mar-31 18:09:32.75



18
2458210.583976
2018-Apr-02 02:00:55.53



19
2458211.911323
2018-Apr-03 09:52:18.31



20
2458213.23867
2018-Apr-04 17:43:41.09



21
2458214.566017
2018-Apr-06 01:35:03.87



22
2458215.893364
2018-Apr-07 09:26:26.65



23
2458217.220711
2018-Apr-08 17:17:49.43



24
2458218.548058
2018-Apr-10 01:09:12.21



25
2458219.875405
2018-Apr-11 09:00:34.99



26
2458221.202752
2018-Apr-12 16:51:57.77



27
2458222.530099
2018-Apr-14 00:43:20.55



28
2458223.857446
2018-Apr-15 08:34:43.33



29
2458225.184793
2018-Apr-16 16:26:06.12



30
2458226.51214
2018-Apr-18 00:17:28.9



31
2458227.839487
2018-Apr-19 08:08:51.68



32
2458229.166834
2018-Apr-20 16:00:14.46



33
2458230.494181
2018-Apr-21 23:51:37.24



34
2458231.821528
2018-Apr-23 07:43:00.02



35
2458233.148875
2018-Apr-24 15:34:22.8



36
2458234.476222
2018-Apr-25 23:25:45.58



37
2458235.803569
2018-Apr-27 07:17:08.36



38
2458237.130916
2018-Apr-28 15:08:31.14



39
2458238.458263
2018-Apr-29 22:59:53.92



40
2458239.78561
2018-May-01 06:51:16.7



41
2458241.112957
2018-May-02 14:42:39.48



42
2458242.440304
2018-May-03 22:34:02.27



43
2458243.767651
2018-May-05 06:25:25.05



44
2458245.094998
2018-May-06 14:16:47.83



45
2458246.422345
2018-May-07 22:08:10.61



46
2458247.749692
2018-May-09 05:59:33.39



47
2458249.077039
2018-May-10 13:50:56.17



48
2458250.404386
2018-May-11 21:42:18.95



49
2458251.731733
2018-May-13 05:33:41.73



50
2458253.05908
2018-May-14 13:25:04.51

Comments

Post a Comment

Popular posts from this blog

HW5 Q1 (William Matzko)

Image
1.  For radial velocities, we have the relation k = 28.4 (m/s) p^-(1/3)(M_p*sin(i)) (M_s)^(-2/3), where M_p is the planet mass in Jupiter masses, and M_s is the mass of the star in solar masses. For astrometry, we have  α  = (M_p/M_s)(a)(d)^(-1), where a is the semimajor axis in AU, and d is the distance to the system in parsecs. Just to make things easier to look at, we can make some simplifying assumptions. We are already told to assume sin(i) = 1 in all instances and M_s is the mass of the Sun. We might as well further assume that our target is one Jupiter mass. Further, we know P is proportional to a ^(3/2). The two above equations then reduce to α = 0.001 (a) (d)^-1 k = 28.4 (m/s) (a)^(-1/2) We are given that astrometry is sensitive to  α = 10 micro arcseconds and above, and radial velocities are sensitive to k = 1 m/s and above. Note that only the astrometry method depends on d. We must now find when the astrometry method becomes more sensitive than the rad...
Read more

HW5 Q2, Q3, Q5 (William Matzko)

2. We know the phase amplitude is  ε = A*g*(R/r)^2, where A is the geometric albedo, R is the radius of the planet, r = a(1-e^2)/(1+ecos(f)), and g = 10^(-0.4*m). Here, m is just  m = 0.09( α/100) + 2.39( α/100)^2 - 0.65( α/100), where  α is in degrees (the 100's in the denominator are likewise in degrees). We are working with a Venus analog, so we can assume, in this part, that e = 0, so that r = a. The phase amplitude will be maximum when  α = 0. So, our maximum phase amplitude is then  ε = (0.67) (1)(6050km/108E6km)^2, where I have plugged in the geometric albedo of Venus, the radius of Venus, and Venus' semimajor axis. This yeilds  ε = 2E-9.  Now, we are to assume i = 40 degrees and w = 270 degrees. The effect of w is to orient the orbit such that we have a maximum phase angle at f = 0 . We know the relation between  α and these orbital elements: cos( α) = -sin(i)sin(w+f). Plugging in the numbers and solving for  α yields  α = 50 d...
Read more

HW 5 Q6 (William Matzko)

6. (a)Usually, the temperature on a planet decreases with altitude. However, in inversion layers, this is not the case. In an inversion layer, temperature increases with altitude. Such a thing is important in the atmosphere of exoplanets because it can give insight into what that planet's atmosphere could be composed of. These inversion layers are usually composed of molecules that absorb UV radiation, so finding these inversion layers on exoplanets is desirable if we want to try and find life on them. (b) We know that P = P_0 exp(-z/H), where P is the pressure, P_0 is the pressure at the surface level, z is the altitude above surface level and H is the scale height. We are given that, on Earth z = 8.8 km and H = 8 km. This corresponds to a pressure of ~0.333*P_0. On a modified planet, to find the z at which we have the same pressure, we first must compute the scale height. We can express the scale height equation as H = k*T/(u*m*g), where k is the Boltzmann constant, T is the te...
Read more