We present an analysis of a slightly eccentric ($e=0.05$), partially eclipsing long-period ($P = 69.73$ d) main sequence binary system (WOCS 12009, Sanders 1247) in the benchmark old open cluster M67. Using Kepler K2 and ground-based photometry along with a large set of new and reanalyzed spectra, we derived highly precise masses ($1.111\pm0.015$ and $0.748\pm0.005 M_\odot$) and radii ($1.071\pm0.008\pm0.003$ and $0.713\pm0.019\pm0.026 R_\odot$, with statistical and systematic error estimates) for the stars. The radius of the secondary star is in agreement with theory. The primary, however, is approximately $15\%$ smaller than reasonable isochrones for the cluster predict. Our best explanation is that the primary star was produced from the merger of two stars, as this can also account for the non-detection of photospheric lithium and its higher temperature relative to other cluster main sequence stars at the same $V$ magnitude. To understand the dynamical characteristics (low measured rotational line broadening of the primary star and the low eccentricity of the current binary orbit), we believe that the most probable (but not the only) explanation is the tidal evolution of a close binary within a primordial triple system (possibly after a period of Kozai-Lidov oscillations), leading to merger approximately 1Gyr ago. This star appears to be a future blue straggler that is being revealed as the cluster ages and the most massive main sequence stars die out.

We use the Athena hydrodynamics code to model propagating two-dimensional relativistic jets as approximations to the growth of radio-loud active galactic nuclei for various input jet velocities and jet-to-ambient matter density ratios. Using results from these simulations we estimate the changing synchrotron emission by summing the fluxes from a vertical strip of zones behind the reconfinement shock, which is nearly stationary, and from which a substantial portion of the flux variability should arise. We explore a wide range of time scales by considering two light curves from each simulation; one uses a relativistic turbulence code with bulk velocities taken from our simulations as input, while the other uses the bulk velocity data to compute fluctuations caused by variations in the Doppler boosting due to changes in the direction and the speed of the flow through all zones in the strip. We then calculate power spectral densities (PSDs) from the light curves for both turbulent and bulk velocity origins for variability. The range of the power-law slopes of the PSDs for the turbulence induced variations is -1.8 to -2.3, while for the bulk velocity produced variations this range is -2.1 to -2.9; these are in agreement with most observations. When superimposed, these power spectra span a very large range in frequency (about five decades), with the turbulent fluctuations yielding most of the shorter timescale variations and the bulk flow changes dominating the longer periods.