-
UCLCHEMCMC: A MCMC Inference tool for Physical Parameters of Molecular Clouds
Authors:
Marcus Keil,
Serena Viti,
Jonathan Holdship
Abstract:
We present the publicly available, open source code UCLCHEMCMC, designed to estimate physical parameters of an observed cloud of gas by combining Monte Carlo Markov Chain (MCMC) sampling with chemical and radiative transfer modeling. When given the observed values of different emission lines, UCLCHEMCMC runs a Bayesian parameter inference, using a MCMC algorithm to sample the likelihood and produc…
▽ More
We present the publicly available, open source code UCLCHEMCMC, designed to estimate physical parameters of an observed cloud of gas by combining Monte Carlo Markov Chain (MCMC) sampling with chemical and radiative transfer modeling. When given the observed values of different emission lines, UCLCHEMCMC runs a Bayesian parameter inference, using a MCMC algorithm to sample the likelihood and produce an estimate of the posterior probability distribution of the parameters. UCLCHEMCMC takes a full forward modeling approach, generating model observables from the physical parameters via chemical and radiative transfer modeling. While running UCLCHEMCMC, the created chemical models and radiative transfer code results are stored in an SQL database, preventing redundant model calculations in future inferences. This means that the more UCLCHEMCMC is used, the more efficient it becomes. Using UCLCHEM and RADEX, the increase of efficiency is nearly two orders of magnitude, going from 5185.33 \pm 1041.96 s for ten walkers to take one thousand steps when the database is empty, to 68.89 \pm 45.39 s when nearly all models requested are in the database. In order to demonstrate its usefulness we provide an example inference of UCLCHEMCMC to estimate the physical parameters of mock data, and perform two inferences on the well studied prestellar core, L1544, one of which show that it is important to consider the substructures of an object when determining which emission lines to use.
△ Less
Submitted 4 February, 2022;
originally announced February 2022.
-
The Complex Circumstellar and Circumbinary Environment of V356 Sgr
Authors:
Jamie R. Lomax,
Andrew G. Fullard,
Michael A. Malatesta,
Brian Babler,
Daniel Bednarski,
Jodi R. Berdis,
Karen S. Bjorkman,
Jon E. Bjorkman,
Alex C. Carciofi,
James W. Davidson Jr.,
Marcus Keil,
Marilyn R. Meade,
Kenneth Nordsieck,
Matt Scheffler,
Jennifer L. Hoffman,
John P. Wisniewski
Abstract:
We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable…
▽ More
We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable eclipses in the polarization phase curves indicates the presence of a substantial pool of scatterers not occulted by either star. We suggest that these scatterers form either a circumbinary disk coplanar with the gainer's accretion disk or an elongated structure perpendicular to the orbital plane of V356 Sgr, possibly formed by bipolar outflows. We also observe small-scale, cycle-to-cycle variations in the magnitude of intrinsic polarization at individual phases, which we interpret as evidence of variability in the amount of scattering material present within and around the system. This may indicate a mass transfer or mass loss rate that varies on the time-scale of the system's orbital period. Finally, we compare the basic polarimetric properties of V356 Sgr with those of the well studied beta Lyr system; the significant differences observed between the two systems suggests diversity in the basic circumstellar geometry of Roche-lobe overflow systems
△ Less
Submitted 23 September, 2016;
originally announced September 2016.
-
Kepler Flares II: The Temporal Morphology of White-Light Flares on GJ 1243
Authors:
James R. A. Davenport,
Suzanne L. Hawley,
Leslie Hebb,
John P. Wisniewski,
Adam F. Kowalski,
Emily C. Johnson,
Michael Malatesta,
Jesus Peraza,
Marcus Keil,
Steven M. Silverberg,
Tiffany C. Jansen,
Matthew S. Scheffler,
Jodi R. Berdis,
Daniel M. Larsen,
Eric J. Hilton
Abstract:
We present the largest sample of flares ever compiled for a single M dwarf, the active M4 star GJ 1243. Over 6100 individual flare events, with energies ranging from $10^{29}$ to $10^{33}$ erg, are found in 11 months of 1-minute cadence data from Kepler. This sample is unique for its completeness and dynamic range. We have developed automated tools for finding flares in short-cadence Kepler light…
▽ More
We present the largest sample of flares ever compiled for a single M dwarf, the active M4 star GJ 1243. Over 6100 individual flare events, with energies ranging from $10^{29}$ to $10^{33}$ erg, are found in 11 months of 1-minute cadence data from Kepler. This sample is unique for its completeness and dynamic range. We have developed automated tools for finding flares in short-cadence Kepler light curves, and performed extensive validation and classification of the sample by eye. From this pristine sample of flares we generate a median flare template. This template shows that two exponential cooling phases are present during the white-light flare decay, providing fundamental constraints for models of flare physics. The template is also used as a basis function to decompose complex multi-peaked flares, allowing us to study the energy distribution of these events. Only a small number of flare events are not well fit by our template. We find that complex, multi-peaked flares occur in over 80% of flares with a duration of 50 minutes or greater. The underlying distribution of flare durations for events 10 minutes and longer appears to follow a broken power law. Our results support the idea that sympathetic flaring may be responsible for some complex flare events.
△ Less
Submitted 13 November, 2014;
originally announced November 2014.
-
Supernova Neutrino Spectra and Applications to Flavor Oscillations
Authors:
Mathias Th. Keil
Abstract:
We study the flavor-dependent neutrino spectra formation in the core of a supernova (SN) by means of Monte Carlo simulations. A high-statistics neutrino signal from a galactic SN may contain information that severely constrains the parameter space for neutrino oscillations. Therefore, reliable predictions for flavor-dependent fluxes and spectra are urgently needed. In all traditional hydrodynami…
▽ More
We study the flavor-dependent neutrino spectra formation in the core of a supernova (SN) by means of Monte Carlo simulations. A high-statistics neutrino signal from a galactic SN may contain information that severely constrains the parameter space for neutrino oscillations. Therefore, reliable predictions for flavor-dependent fluxes and spectra are urgently needed. In all traditional hydrodynamic simulations the nu_mu,tau and nu_mu,tau-bar interactions commonly included are rather schematic. With our Monte Carlo simulations we find that the most relevant sources for nu_mu,tau and nu_mu,tau-bar are traditionally not included. In comparing our numerical results for all flavors we find the standard hierarchy of mean energies nu_e < nu_e-bar < nu_mu,tau, with, however, very similar values for nu_mu,tau and nu_e-bar. The luminosities of nu_mu,tau and nu_mu,tau-bar can differ by up to a factor of 2 from L_nue-bar and L_nue, the latter two are very similar. The Garching Group obtains similar results from their self-consistent simulation with the full set of interactions. These results are almost orthogonal to the previous standard picture of exactly equal luminosities of all flavors and differences in mean energies of up to a factor of 2. Existing concepts for identifying oscillation effects in a SN neutrino signal need to be revised. We present two methods for detecting the earth-matter effect that are rather independent of predictions from SN simulations.
△ Less
Submitted 13 August, 2003;
originally announced August 2003.
-
Identifying Earth matter effects on supernova neutrinos at a single detector
Authors:
Amol S. Dighe,
Mathias Th. Keil,
Georg G. Raffelt
Abstract:
The neutrino oscillations in Earth matter introduce modulations in the supernova neutrino spectra. These modulations can be exploited to identify the presence of Earth effects on the spectra, which would enable us to put a limit on the value of the neutrino mixing angle $θ_{13}$ and to identify whether the mass hierarchy is normal or inverted. We demonstrate how the Earth effects can be identifi…
▽ More
The neutrino oscillations in Earth matter introduce modulations in the supernova neutrino spectra. These modulations can be exploited to identify the presence of Earth effects on the spectra, which would enable us to put a limit on the value of the neutrino mixing angle $θ_{13}$ and to identify whether the mass hierarchy is normal or inverted. We demonstrate how the Earth effects can be identified at a single detector without prior assumptions about the flavor-dependent source spectra, using the Fourier transform of the ``inverse-energy'' spectrum of the signal. We explore the factors affecting the efficiency of this method, and find that the energy resolution of the detector is the most crucial one. In particular, whereas water Cherenkov detectors may need a few ten thousand events to identify the Earth effects, a few thousand may be enough at scintillation detectors, which generically have a much better energy resolution. A successful identification of the Earth effects through this method can also provide $Δm^2_\odot$ to a good accuracy. The relative strength of the detected Earth effects as a function of time provides a test for supernova models.
△ Less
Submitted 17 June, 2003; v1 submitted 15 April, 2003;
originally announced April 2003.
-
Detecting the Neutrino Mass Hierarchy with a Supernova at IceCube
Authors:
A. S. Dighe,
M. T. Keil,
G. G. Raffelt
Abstract:
IceCube, a future km^3 antarctic ice Cherenkov neutrino telescope, is highly sensitive to a galactic supernova (SN) neutrino burst. The Cherenkov light corresponding to the total energy deposited by the SN neutrinos in the ice can be measured relative to background fluctuations with a statistical precision much better than 1%. If the SN is viewed through the Earth, the matter effect on neutrino…
▽ More
IceCube, a future km^3 antarctic ice Cherenkov neutrino telescope, is highly sensitive to a galactic supernova (SN) neutrino burst. The Cherenkov light corresponding to the total energy deposited by the SN neutrinos in the ice can be measured relative to background fluctuations with a statistical precision much better than 1%. If the SN is viewed through the Earth, the matter effect on neutrino oscillations can change the signal by more than 5%, depending on the flavor-dependent source spectra and the neutrino mixing parameters. Therefore, IceCube together with another high-statistics experiment like Hyper-Kamiokande can detect the Earth effect, an observation that would identify specific neutrino mixing scenarios that are difficult to pin down with long-baseline experiments. In particular, the normal mass hierarchy can be clearly detected if the third mixing angle is not too small, sin^2 theta_13 < 10^-3. The small flavor-dependent differences of the SN neutrino fluxes and spectra that are found in state-of-the-art simulations suffice for this purpose. Although the absolute calibration uncertainty at IceCube may exceed 5%, the Earth effect would typically vary by a large amount over the duration of the SN signal, obviating the need for a precise calibration. Therefore, IceCube with its unique geographic location and expected longevity can play a decisive role as a "co-detector" to measure SN neutrino oscillations. It is also a powerful stand-alone SN detector that can verify the delayed-explosion scenario.
△ Less
Submitted 28 May, 2003; v1 submitted 25 March, 2003;
originally announced March 2003.
-
Supernova neutrinos: Flavor-dependent fluxes and spectra
Authors:
G. G. Raffelt,
M. Th. Keil,
R. Buras,
H. -T. Janka,
M. Rampp
Abstract:
Transporting nu_mu and nu_tau in a supernova (SN) core involves several processes that have been neglected in traditional simulations. Based on a Monte Carlo study we find that the flavor-dependent spectral differences are much smaller than is often stated in the literature. A full-scale SN simulation using a Boltzmann solver and including all relevant neutrino reactions confirms these results.…
▽ More
Transporting nu_mu and nu_tau in a supernova (SN) core involves several processes that have been neglected in traditional simulations. Based on a Monte Carlo study we find that the flavor-dependent spectral differences are much smaller than is often stated in the literature. A full-scale SN simulation using a Boltzmann solver and including all relevant neutrino reactions confirms these results. The flavor-dependent flux differences are largest during the initial accretion phase.
△ Less
Submitted 11 March, 2003;
originally announced March 2003.
-
Monte Carlo Study of Supernova Neutrino Spectra Formation
Authors:
Mathias Th. Keil,
Georg G. Raffelt,
Hans-Thomas Janka
Abstract:
The neutrino flux and spectra formation in a supernova core is studied by using a Monte Carlo code. The dominant opacity contribution for nu_mu and nu_tau is elastic scattering on nucleons. In addition we switch on or off a variety of processes which allow for the exchange of energy or the creation and destruction of neutrino pairs, notably nucleon bremsstrahlung, the e^+ e^- pair annihilation p…
▽ More
The neutrino flux and spectra formation in a supernova core is studied by using a Monte Carlo code. The dominant opacity contribution for nu_mu and nu_tau is elastic scattering on nucleons. In addition we switch on or off a variety of processes which allow for the exchange of energy or the creation and destruction of neutrino pairs, notably nucleon bremsstrahlung, the e^+ e^- pair annihilation process and nu_e-bar nu_e -> nu_{mu,tau} nu_{mu,tau}-bar, recoil and weak magnetism in elastic nucleon scattering, elastic scattering on electrons and positrons and elastic scattering on electron neutrinos and anti-neutrinos. The least important processes are neutrino-neutrino scattering and e^+ e^- annihilation. The formation of the spectra and fluxes of nu_mu is dominated by the nucleonic processes, i.e. bremsstrahlung and elastic scattering with recoil, but also nu_e nu_e-bar annihilation and nu_mu e^\pm scattering contribute significantly. When all processes are included, the spectral shape of the emitted neutrino flux is always ``pinched,'' i.e. the width of the spectrum is smaller than that of a thermal spectrum with the same average energy. In all of our cases we find that the average nu_mu-bar energy exceeds the average nu_e-bar energy by only a small amount, 10% being a typical number. Weak magnetism effects cause the opacity of nu_mu to differ slightly from that of nu_mu-bar, translating into differences of the luminosities and average energies of a few percent. Depending on the density, temperature, and composition profile, the flavor-dependent luminosities L_{nu_e}$, L_{nu_e-bar}, and L_{nu_mu} can mutually differ from each other by up to a factor of two in either direction.
△ Less
Submitted 21 February, 2003; v1 submitted 1 August, 2002;
originally announced August 2002.
-
Electron Neutrino Pair Annihilation: A New Source for Muon and Tau Neutrinos in Supernovae
Authors:
Robert Buras,
Hans-Thomas Janka,
Mathias Th. Keil,
Georg G. Raffelt,
Markus Rampp
Abstract:
We show that in a supernova core the annihilation process nu_e nu_e-bar -> nu_{mu,tau} nu_{mu,tau}-bar is always more important than the traditional reaction e^+ e^- -> nu_{mu,tau} nu_{mu,tau}-bar as a source for muon and tau neutrino pairs. We study the impact of the new process by means of a Monte Carlo transport code with a static stellar background model and by means of a self-consistent hyd…
▽ More
We show that in a supernova core the annihilation process nu_e nu_e-bar -> nu_{mu,tau} nu_{mu,tau}-bar is always more important than the traditional reaction e^+ e^- -> nu_{mu,tau} nu_{mu,tau}-bar as a source for muon and tau neutrino pairs. We study the impact of the new process by means of a Monte Carlo transport code with a static stellar background model and by means of a self-consistent hydrodynamical simulation with Boltzmann neutrino transport. Nucleon bremsstrahlung NN -> NN nu_{mu,tau} nu_{mu,tau}-bar is also included as another important source term. Taking into account nu_e nu_e-bar -> nu_{mu,tau} nu_{mu,tau}-bar increases the nu_mu and nu_tau luminosities by as much as 20% while the spectra remain almost unaffected. In our hydrodynamical simulation the shock was somewhat weakened. Elastic nu_{mu,tau} nu_e and nu_{mu,tau} nu_e scattering is not negligible but less important than nu_{mu,tau} e^+ or e^- scattering. Its influence on the nu_{mu,tau} fluxes and spectra is small after all other processes have been included.
△ Less
Submitted 1 May, 2002;
originally announced May 2002.
-
Gauged Inflation
Authors:
Ralf Hofmann,
Mathias Th. Keil
Abstract:
We propose a model for cosmic inflation which is based on an effective description of strongly interacting, nonsupersymmetric matter within the framework of dynamical Abelian projection and centerization. The underlying gauge symmetry is assumed to be $SU(N+1)$ with $N \gg 1$. Appealing to a thermodynamical treatment, the ground-state structure of the model is classically determined by a potenti…
▽ More
We propose a model for cosmic inflation which is based on an effective description of strongly interacting, nonsupersymmetric matter within the framework of dynamical Abelian projection and centerization. The underlying gauge symmetry is assumed to be $SU(N+1)$ with $N \gg 1$. Appealing to a thermodynamical treatment, the ground-state structure of the model is classically determined by a potential for the inflaton field (dynamical monopole condensate) which allows for nontrivially BPS saturated and thereby stable solutions. For $T<M_P$ this leads to decoupling of gravity from the inflaton dynamics. The ground state dynamics implies a heat capacity for the vacuum leading to inflation for temperatures comparable to the mass scale $M$ of the potential. The dynamics has an attractor property. In contrast to the usual slow-roll paradigm we have $m\gg H$ during inflation. As a consequence, density perturbations generated from the inflaton are irrelevant for the formation of large-scale structure, and the model has to be supplemented with an inflaton independent mechanism for the generation of spatial curvature perturbations. Within a small fraction of the Hubble time inflation is terminated by a transition of the theory to its center symmetric phase. The spontaneously broken $Z_{N+1}$ symmetry stabilizes relic vector bosons in the epochs following inflation. These heavy relics contribute to the cold dark matter of the universe and potentially originate the UHECRs beyond the GZK bound.
△ Less
Submitted 28 May, 2002; v1 submitted 7 November, 2001;
originally announced November 2001.
