SN 1954J in NGC 2403 and SN 1961V in NGC 1058 were two luminous transients whose definitive classification as either non-terminal eruptions or supernovae remains elusive. A critical question is whether a surviving star can be significantly obscured by dust formed from material ejected during the transient. We use three lines of argument to show that the candidate surviving stars are not significantly optically extincted ($\tau \lesssim 1$) by dust formed in the transients. First, we use SED fits to new HST optical and near-IR photometry. Second, neither source is becoming brighter as required by absorption from an expanding shell of ejected material. Third, the ejecta masses implied by the H$\alpha$ luminosities are too low to produce significant dust absorption. The latter two arguments hold independent of the dust properties. The H$\alpha$ fluxes should also be declining with time as $t^{-3}$, and this seems not to be observed. As a result, it seems unlikely that recently formed dust can be responsible for the present faintness of the sources compared to their progenitors, although this can be verified with \textit{JWST}. This leaves three possibilities: 1) the survivors were misidentified; 2) they are intrinsically less luminous; 3) SN 1954J and SN 1961V were true supernovae.

A crucial ingredient in population synthesis studies involving massive stars is the determination of whether they explode or implode in the end. While the final fate of a massive star is sensitive to its core structure at the onset of collapse, the existing binary population synthesis studies do not reach core-collapse. Instead, they employ simple prescriptions to infer their final fates without knowing the presupernova core structure. We explore a potential solution to this problem by treating the carbon-oxygen (CO) core independently from the rest of the star. Using the implicit hydrodynamics code KEPLER, we have computed an extensive grid of 3496 CO-core models from a diverse range of initial conditions, each evolved from carbon ignition until core-collapse. The final core structure, and thus the explodability, varies non-monotonically and depends sensitively on both the mass and initial composition of the CO-core. Although bare CO-cores are not perfect substitutes for cores embedded in massive stars, our models compare well both with MESA and full hydrogenic and helium star calculations. Our results can be used to infer the presupernova core structures from population synthesis estimates of CO-core properties, thus to determine the final outcomes based on the results of modern neutrino-driven explosion simulations. A sample application is presented for a population of Type-IIb supernova progenitors. All of our models are available at https://doi.org/10.5281/zenodo.3785377.

A model for terrestrial planets, inclusive of viscous fluid behavior and featuring finite normal stress differences, is developed. This work offers new insights for the interpretation of planetary survey data. Evolution equations for poloidal and toroidal motions include gradients of density, viscosity, and two normal stress moduli. The poloidal field exhibits gradients in the cubic dilation, which couple non-isotropic pressures to the combined deformation field. In contrast, the toroidal field exhibits vorticity gradients with magnitudes proportional to the natural time. This holds even in the absence of material gradients. Consequently, viscosity gradients are not required to drive toroidal motions. The toroidal field is governed by an inhomogeneous diharmonic equation, exhibiting dynamic shear localization. The strain-energy density for this model, as a function of temperature, is found via thermodynamics. Assuming heat transfer with characteristic diffusivity, a radial model parameterized by thermomechanical competence is found. Shear dislocations, admissible for thermomechanical competence greater than 1/2, are found to coincide with supershear rupture speeds for in-plane (Mode II) cracks. This range of thermomechanical competence coincides with depths in the crust, upper mantle and transition zone where earthquake foci are observed. Consequently, all seismic sources must exhibit some supershearing component. Observed variations in Earth's gravity-topography admittance and correlation spectra, and earthquake moment-depth release are interpreted in light of this hypothetical structure.

Analysis of the global centroid-moment tensor catalog reveals significant regional variations of seismic energy release to 290 km depth. These variations reflect radial and lateral contrasts in thermomechanical competence, consistent with a shear-dominated non-adiabatic boundary layer some 700-km thick, capped by denser oceanic lithosphere as much as 100 km thick, or lighter continental tectosphere 170 to 260 km thick. Thus, isobaric shearing at fractally-distributed depths likely facilitates toroidal plate rotations while minimizing global energy dissipation. Shear localization in the shallow crust occurs as dislocations at finite angles with respect to the shortening direction, with a 30 degree angle being the most likely. Consequently, relatively low-angle reverse faults, steep normal faults, and triple junctions with orthogonal or hexagonal symmetry are likely to form in regions of crustal shortening, extension, and transverse motion, respectively. Thermomechanical theory also predicts adiabatic conditions in the mantle below about 1000-km depth, consistent with observed variations in bulk sound speed.