Can Vertical Profiles of Tropospheric Methane on Titan Be Derived from Radio-Occultation Soundings?

The intensity of the received signal at Earth in the radio occultations of Titan is attenuated both by refractive defocusing and pressure-induced absorption from N2-N2 and CH4-N2 pairs. Because the absorption strength is different for the two sets of pairs, matching the retrieved absorptivity profile can in principle yield the vertical variation in gaseous methane in the troposphere. There are two factors that make this difficult. The first is the propagation of noise in the phase and amplitude of the received signal in the absorption retrieval. The phase data is first inverted to retrieve vertical profiles of refractivity, from which the refractive defocusing is calculated. This is then subtracted from the observed. intensity attenuation of the received signal to generate a profile of atmospheric absorption. The second problem is the uncertainty in the pressure-induced absorption coefficients. Laboratory data at radio wavelengths is only available near room temperature (see, e.g., [1] for N2-N2), and the extrapolation to the low temperatures in Titan's troposphere is not well established. Ab initio calculations by Borysow et al. [2, 3] provide absorption coefficients at low temperatures and long wavelengths, but their accuracy has come into question. We present examples from Cassini radio occultations of Titan to illustrate the difficulties. For methane mole fractions in the lower troposphere comparable to that inferred from the Huygens probe (approximately 0.05), it will be difficult to separate the contributions of N2-N2 collisions from those of N2-CH4, collisions to the retrieved absorption. However, higher concentrations of CH4 and/or a higher signal-to-noise ratio from a future uplink experiment could result in a successful separation of the two components. However, key to this are highly accurate estimates of the absorption from a combination of laboratory measurements at love temperatures and long wavelengths, and possibly improved theoretical calculations

The Structure of Titan's Atmosphere from Cassini Radio Occultations

We present results from the two radio occultations of the Cassini spacecraft by Titan in 2006, which probed mid-southern latitudes. Three of the ingress and egress soundings occurred within a narrow latitude range, 31.34 deg S near the surface, and the fourth at 52.8 deg S. Temperature - altitude profiles for all four occultation soundings are presented, and compared with the results of the Voyager 1 radio occultation (Lindal et al., 1983), the HASI instrument on the Huygens descent probe (Fulchignoni et al., 2005), and Cassini CIRS results (Flasar et al., 2005; Achterberg et al., 2008b). Sources of error in the retrieved temperature - altitude profiles are also discussed, and a major contribution is from spacecraft velocity errors in the reconstructed ephemeris. These can be reduced by using CIRS data at 300 km to make along-track adjustments of the spacecraft timing. The occultation soundings indicate that the temperatures just above the surface at 31-34 deg S are about 93 K, while that at 53 deg S is about 1 K colder. At the tropopause, the temperatures at the lower latitudes are all about 70 K, while the 53 deg S profile is again 1 K colder. The temperature lapse rate in the lowest 2 km for the two ingress (dawn) profiles at 31 and 33 deg S lie along a dry adiabat except within approximately 200m of the surface, where a small stable inversion occurs. This could be explained by turbulent mixing with low viscosity near the surface. The egress profile near 34 deg S shows a more complex structure in the lowest 2 km, while the egress profile at 53 deg S is more stable

Pulsar kicks from neutrino oscillations

Neutrino oscillations can explain the observed motion of pulsars. We show that two different models of neutrino emission from a cooling neutron star are in good quantitative agreement and predict the same order of magnitude for the pulsar kick velocity, consistent with the data.Comment: revtex; 4 page

scX: A user-friendly tool for scRNA-seq exploration

Single-cell RNA sequencing (scRNA-seq) has transformed our ability to explore biological systems. Nevertheless, proficient expertise is essential for handling and interpreting the data. In this paper, we present scX, an R package built on the Shiny framework that streamlines the analysis, exploration, and visualization of single-cell experiments. With an interactive graphic interface, implemented as a web application, scX provides easy access to key scRNAseq analyses, including marker identification, gene expression profiling, and differential gene expression analysis. Additionally, scX seamlessly integrates with commonly used single-cell Seurat and SingleCellExperiment R objects, resulting in efficient processing and visualization of varied datasets. Overall, scX serves as a valuable and user-friendly tool for effortless exploration and sharing of single-cell data, simplifying some of the complexities inherent in scRNAseq analysis.Comment: 10 pages, 2 figures. Source code can be downloaded from https://github.com/chernolabs/scX. User manual available at https://chernolabs.github.io/scX/. Docker image available from dockerhub as chernolabs/sc

Neutrino - nucleon reaction rates in the supernova core in the relativistic random phase approximation

In view of the application to supernova simulations, we calculate neutrino reaction rates with nucleons via the neutral and charged currents in the supernova core in the relativistic random phase approximation (RPA) and study their effects on the opacity of the supernova core. The formulation is based on the Lagrangian employed in the calculation of nuclear equation of state (EOS) in the relativistic mean field theory (RMF). The nonlinear meson terms are treated appropriately so that the consistency of the density correlation derived in RPA with the thermodynamic derivative obtained from EOS by RMF is satisfied in the static and long wave length limit. We employ pion and rho meson exchange interactions together with the phenomenological Landau-Migdal parameters for the isospin-dependent nuclear interactions. We find that both the charged and neutral current reaction rates are suppressed from the standard Bruenn's approximate formula considerably in the high density regime. In the low density regime, on the other hand, the vector current contribution to the neutrino-nucleon scattering rate is enhanced in the vicinity of the boundary of the liquid-gas phase transition, while the other contributions are moderately suppressed there also. In the high temperature regime or in the regime where electrons have a large chemical potential, the latter of which is important only for the electron capture process and its inverse process, the recoil of nucleons cannot be neglected and further reduces the reaction rates with respect to the standard approximate formula which discards any energy transfer in the processes. These issues could have a great impact on the neutrino heating mechanism of collapse-driven supernovae.Comment: 16pages, 19figures, submitted to PR

General Relativistic Effects in the Core Collapse Supernova Mechanism

We apply our recently developed code for spherically symmetric, fully general relativistic (GR) Lagrangian hydrodynamics and multigroup flux-limited diffusion neutrino transport to examine the effects of GR on the hydrodynamics and transport during collapse, bounce, and the critical shock reheating phase of core collapse supernovae. Comparisons of models computed with GR versus Newtonian hydrodynamics show that collapse to bounce takes slightly less time in the GR limit, and that the shock propagates slightly farther out in radius before receding. After a secondary quasistatic rise in the shock radius, the shock radius declines considerably more rapidly in the GR simulations than in the corresponding Newtonian simulations. During the shock reheating phase, core collapse computed with GR hydrodynamics results in a substantially more compact structure from the center out to the stagnated shock. The inflow speed of material behind the shock is also increased. Comparisons also show that the luminosity and rms energy of any neutrino flavor during the shock reheating phase increases when switching from Newtonian to GR hydrodynamics, and decreases when switching from Newtonian to GR transport. This latter decrease in neutrino luminosities and rms energies is less in magnitude than the increase that arise when switching from Newtonian to GR hydrodynamics, with the result that a fully GR simulation gives higher neutrino luminosities and harder neutrino spectra than a fully Newtonian simulation of the same precollapse model.Comment: 35 pages, 23 figure

Neutrino Interactions in Hot and Dense Matter

We study the charged and neutral current weak interaction rates relevant for the determination of neutrino opacities in dense matter found in supernovae and neutron stars. We establish an efficient formalism for calculating differential cross sections and mean free paths for interacting, asymmetric nuclear matter at arbitrary degeneracy. The formalism is valid for both charged and neutral current reactions. Strong interaction corrections are incorporated through the in-medium single particle energies at the relevant density and temperature. The effects of strong interactions on the weak interaction rates are investigated using both potential and effective field-theoretical models of matter. We investigate the relative importance of charged and neutral currents for different astrophysical situations, and also examine the influence of strangeness-bearing hyperons. Our findings show that the mean free paths are significantly altered by the effects of strong interactions and the multi-component nature of dense matter. The opacities are then discussed in the context of the evolution of the core of a protoneutron star.Comment: 41 pages, 25 figure

Three Dimensional Numerical General Relativistic Hydrodynamics I: Formulations, Methods, and Code Tests

This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and convergence of our general relativistic hydrodynamic treatment and its coupling to the spacetime evolutions described by the full set of Einstein equations with a perfect fluid source. The numerical treatment of the general relativistic hydrodynamic equations is based on high resolution shock capturing schemes. These schemes rely on the characteristic information of the system. A spectral decomposition for general relativistic hydrodynamics suitable for a general spacetime metric is presented. Evolutions based on three different approximate Riemann solvers coupled to four different discretizations of the Einstein equations are studied and compared. The coupling between the hydrodynamics and the spacetime (the right and left hand side of the Einstein equations) is carried out in a treatment which is second order accurate in {\it both} space and time. Convergence tests for all twelve combinations with a variety of test beds are studied, showing consistency with the differential equations and correct convergence properties. The test-beds examined include shocktubes, Friedmann-Robertson-Walker cosmology tests, evolutions of self-gravitating compact (TOV) stars, and evolutions of relativistically boosted TOV stars. Special attention is paid to the numerical evolution of strongly gravitating objects, e.g., neutron stars, in the full theory of general relativity, including a simple, yet effective treatment for the surface region of the star (where the rest mass density is abruptly dropping to zero).Comment: 45 pages RevTeX, 34 figure

Light curve and neutrino spectrum emitted during the collapse of a nonrotating, supermassive star

The formation of a neutrino pulse emitted during the relativistic collapse of a spherical supermassive star is considered. The free collapse of a body with uniform density in the absence of rotation and with the free escape of the emitted neutrinos can be solved analytically by quadrature. The light curve of the collapsing star and the spectrum of the emitted neutrinos at various times are calculated.Comment: 17 pages, 2 figures, published in Astronomy Report