The Longitudinal Study of Astronomy Graduate Students (LSAGS) arose from the 2003 Women in Astronomy Conference, where it was noted that a majority of young members of the American Astronomical Society were women. The astronomy community wishes to make every effort to retain young women in astronomy, so they commissioned a longitudinal study to be conducted that would pinpoint the factors that contribute to retention in general, with a focus on differences between women and men. The LSAGS follows a cohort of people who were graduate students in astronomy or astrophysics during 2006-07. The first survey was conducted during 2007-08, the second during 2012-13, and the third during 2015. The analysis presented in this paper, which is an update to our previous paper on this topic, used a subset of the respondents, all of whom had PhDs in astronomy, astrophysics, or a related field at the time of the third survey. We tested the effects of four major concepts on attrition from physics and astronomy. These concepts included: the imposter syndrome, mentoring and advising during graduate school, the so-called "two-body problem" that occurs when a couple needs to find two jobs in the same geographic area, and gender of the respondent. Having a mentor in grad school did not contribute to working outside of physics or astronomy. Showing characteristics of the imposter syndrome and gender of the respondent had indirect effects on working outside the field. Encouragement of the graduate advisor, the two-body problem, and completing a postdoc, had significant direct effects on working in physics or astronomy. This research identifies specific areas of concern that can be addressed by the scientific community to increase the retention of all people, but especially women, in astronomy and astrophysics.

Building on the framework of Zhang \& Shu \cite{zhangShu_2010a,zhangShu_2010b}, we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function $f$. The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky \& Bludman \cite{cernohorskyBludman_1994} and Banach \& Larecki \cite{banachLarecki_2017a}. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., $f\in[0,1]$). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. Numerical results demonstrate the realizability-preserving properties of the scheme. We also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.