Fully-Coupled Simulation of Cosmic Reionization. I: Numerical Methods and Tests

We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. We do, however, employ a simple subgrid model for star formation which we calibrate to observations. Radiation transport is done in the grey flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. We illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with $3200^3$ Eulerian grid cells and dark matter particles.Comment: 32 pages, 23 figures. ApJ Supp accepted. New title and substantial revisions re. v

Gamma-ray transfer and energy deposition in supernovae

Solutions to the energy-independent (gray) radiative transfer equations are compared to results of Monte Carlo simulations of the \Ni\ and \Co\ radioactive decay \GR\ energy deposition in supernovae. The comparison shows that an effective, purely absorptive, gray opacity, \KG\ \sim (0.06 \pm 0.01)Y_e cm^2 g^{-1}, where Y_e is the total number of electrons per baryon, accurately describes the interaction of \GRs\ with the cool supernova gas and the local \GR\ energy deposition within the gas. The nature of the \GR\ interaction process (dominated by Compton scattering in the relativistic regime) creates a weak dependence of \KG\ on the optical thickness of the (spherically symmetric) supernova atmosphere: The maximum value of \KG\ applies during optically thick conditions when individual \GRs\ undergo multiple scattering encounters and the lower bound is reached at the phase characterized by a total Thomson optical depth to the center of the atmosphere \te\ \LA\ 1. Our results quantitatively confirm that the quick and efficient solution to the gray transfer problem provides an accurate representation of \GR\ energy deposition for a broad range of supernova conditions

Development of a robust mating system for use in the autonomous assembly of planetary drill strings

Volume-constrained robotic missions seeking to obtain samples from beneath a planetary subsurface may wish to use a rigid drill string consisting of multiple, individual drill bit sections connected together, as opposed to a single, lengthy drill bit. To ensure that drill strings can be assembled and disassembled reliably, it is essential that a robust connection system be used. The authors propose a geometry that seeks to address the requirements of such a mating interface. The proposed solution is based on the bayonet interface, using L- and T-shaped so-called female grooves and male studs connected and disconnected together through a series of clockwise and counterclockwise rotations and single-point clamping events. This routine allows the transfer of both percussion through the drill string and torque in both directions of rotation, while permitting the accurate disconnection of individual drills bits at the required location. Sustained laboratory and field drilling operations suggest that bayonet-style connections offer a reliable solution to the problem of autonomous assembly and disassembly of drill strings in a planetary exploration setting. This paper discusses the development of such a connection system, based on the bayonet connection, which has been implemented in the overall architecture of the Ultrasonic Planetary Core Drill (UPCD). The design trade-off study, which sought to evaluate the use of the bayonet system in comparison with the more conventional screw thread interface, will be discussed, alongside experimental results from percussion transmission testing and drill string assembly testing

Turbulent Mixing in the Interstellar Medium -- an application for Lagrangian Tracer Particles

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.Comment: 10 pages, 5 figures, conference proceedings "Turbulent Mixing and Beyond 2007

Characterizing and Predicting Canadian Adolescents’ Internalizing Symptoms In The First Year Of The COVID-19 Pandemic

To date, most longitudinal studies of adolescents’ internalizing symptoms during the COVID-19 pandemic include few time points, limiting knowledge about the long-term course of adolescents’ mental health during the pandemic. Moreover, examining intraindividual variability in symptoms, which may have important implications for adolescents’ adjustment beyond mean or “typical” symptoms, requires multiple time points. We examined the course of internalizing symptoms in 271 Ontario adolescents (mean n = 193 across time points) during the first year of the pandemic (March 2020–April 2021) via mixed-effect location scale models, drawing upon established internalizing symptom risk factors as predictors of mean trends and intraindividual variability. Adolescents’ internalizing symptoms were relatively stable and generally low over the first year of the pandemic, with severity peaking in February and April 2021. Girls showed more symptoms on average and greater intraindividual variability in symptoms. Parents’ depressive symptoms predicted intraindividual variability in adolescents’ anxious and depressive symptoms. Adolescents’ symptoms were stable and generally below clinical cutoffs. However, female adolescents and those whose parents experienced more depressive symptoms were most vulnerable to the stress of the pandemic. Implications for intervention and prevention efforts are discussed

Studying stellar binary systems with the Laser Interferometer Space Antenna using Delayed Rejection Markov chain Monte Carlo methods

Bayesian analysis of LISA data sets based on Markov chain Monte Carlo methods has been shown to be a challenging problem, in part due to the complicated structure of the likelihood function consisting of several isolated local maxima that dramatically reduces the efficiency of the sampling techniques. Here we introduce a new fully Markovian algorithm, a Delayed Rejection Metropolis-Hastings Markov chain Monte Carlo method, to efficiently explore these kind of structures and we demonstrate its performance on selected LISA data sets containing a known number of stellar-mass binary signals embedded in Gaussian stationary noise.Comment: 12 pages, 4 figures, accepted in CQG (GWDAW-13 proceedings

Manipulation and removal of defects in spontaneous optical patterns

Defects play an important role in a number of fields dealing with ordered structures. They are often described in terms of their topology, mutual interaction and their statistical characteristics. We demonstrate theoretically and experimentally the possibility of an active manipulation and removal of defects. We focus on the spontaneous formation of two-dimensional spatial structures in a nonlinear optical system, a liquid crystal light valve under single optical feedback. With increasing distance from threshold, the spontaneously formed hexagonal pattern becomes disordered and contains several defects. A scheme based on Fourier filtering allows us to remove defects and to restore spatial order. Starting without control, the controlled area is progressively expanded, such that defects are swept out of the active area.Comment: 4 pages, 4 figure

Evidence for Supernova Signatures in the Spectrum of the Late-time Bump of the Optical Afterglow of GRB 021211

We present photometric and spectroscopic observations of the gamma-ray burst GRB 021211 obtained during the late stages of its afterglow. The light curve shows a rebrightening occurring ~25 days after the GRB. The analysis of a VLT spectrum obtained during the bump (27 days after the GRB) reveals a suggestive resemblance with the spectrum of the prototypical type-Ic SN 1994I, obtained about ~10 days past maximum light. Particularly we have measured a strong, broad absorption feature at 3770 A, which we have identified with Ca II blueshifted by ~14400 km/s, thus indicating that a supernova (SN) component is indeed powering the `bump' in the afterglow decay. Assuming SN 1994I as a template, the spectroscopic and photometric data together indicate that the SN and GRB explosions were at most separated by a few days. Our results suggest that GRBs might be associated also to standard type-Ic supernovae.Comment: 6 pages, 4 color figures. Accepted for publication in A&A Letters. Fig. 4 does not appair in the A&A version due to space restrictions. Includes aa.cls and txfonts.st