The planar Green`s function in an infinite multiplying medium

Throughout the history of neutron transport theory, the study of simplified problems that have analytical or semi-analytical solutions has been a foundation for more complicated analyses. Analytical transport results are often used as benchmarks or in pedagogical settings. Benchmark problems in infinite homogeneous media have been studied continually, beginning with the monograph by Case, DeHoffmann, and Placzek. A fundamental problem considered in this work is the Green`s function in an infinite medium. The Green`s function problem considers an infinite planar source emitting neutral particles in the single directions`. Recently, this Green`s function has been used to obtain solutions for finite media. These solutions, which hinge on accurate and fast evaluation of the infinite medium Green`s function, use Fourier and Laplace transform inversion techniques for the evaluation. Until now, only absorbing media have been considered, and applications were therefore limited to non-multiplying media. In an effort to relax this limitation, the infinite medium Green`s function is numerically evaluated for an infinite multiplying medium using the double-sided Laplace transform inversion. Of course, no steady-state mathematical solution exists for an infinite multiplying medium with a source present; however, the non-physical solution in an infinite medium can be used in finite media problems where the solution is physically realizable

Analytical three-dimensional neutron transport benchmarks for verification of nuclear engineering codes. Final report

Because of the requirement of accountability and quality control in the scientific world, a demand for high-quality analytical benchmark calculations has arisen in the neutron transport community. The intent of these benchmarks is to provide a numerical standard to which production neutron transport codes may be compared in order to verify proper operation. The overall investigation as modified in the second year renewal application includes the following three primary tasks. Task 1 on two dimensional neutron transport is divided into (a) single medium searchlight problem (SLP) and (b) two-adjacent half-space SLP. Task 2 on three-dimensional neutron transport covers (a) point source in arbitrary geometry, (b) single medium SLP, and (c) two-adjacent half-space SLP. Task 3 on code verification, includes deterministic and probabilistic codes. The primary aim of the proposed investigation was to provide a suite of comprehensive two- and three-dimensional analytical benchmarks for neutron transport theory applications. This objective has been achieved. The suite of benchmarks in infinite media and the three-dimensional SLP are a relatively comprehensive set of one-group benchmarks for isotropically scattering media. Because of time and resource limitations, the extensions of the benchmarks to include multi-group and anisotropic scattering are not included here. Presently, however, enormous advances in the solution for the planar Green`s function in an anisotropically scattering medium have been made and will eventually be implemented in the two- and three-dimensional solutions considered under this grant. Of particular note in this work are the numerical results for the three-dimensional SLP, which have never before been presented. The results presented were made possible only because of the tremendous advances in computing power that have occurred during the past decade

Ring Formation from an Oscillating Black Hole

Massive black hole (BH) mergers can result in the merger remnant receiving a "kick", of order 200 km s$^{-1}$ or more, which will cause the remnant to oscillate about the galaxy centre. Here we analyze the case where the BH oscillates through the galaxy centre perpendicular or parallel to the plane of the galaxy for a model galaxy consisting of an exponential disk, a Plummer model bulge, and an isothermal dark matter halo. For the perpendicular motion we find that there is a strong resonant forcing of the disk radial motion near but somewhat less than the "resonant radii" $r_R$ where the BH oscillation frequency is equal one-half, one-fourth, (1/6, etc.) of the radial epicyclic frequency in the plane of the disk. Near the resonant radii there can be a strong enhancement of the radial flow and disk density which can lead to shock formation. In turn the shock may trigger the formation of a ring of stars near $r_R$. As an example, for a BH mass of $10^8 M_\odot$ and a kick velocity of 150 km s$^{-1}$, we find that the resonant radii lie between 0.2 and 1 kpc. For BH motion parallel to the plane of the galaxy we find that the BH leaves behind it a supersonic wake where star formation may be triggered. The shape of the wake is calculated as well as the slow-down time of the BH. The differential rotation of the disk stretches the wake into ring-like segments.Comment: 7 pages, 7 figure

A Photometric Method for Quantifying Asymmetries in Disk Galaxies

A photometric method for quantifying deviations from axisymmetry in optical images of disk galaxies is applied to a sample of 32 face-on and nearly face-on spirals. The method involves comparing the relative fluxes contained within trapezoidal sectors arranged symmetrically about the galaxy center of light, excluding the bulge and/or barred regions. Such a method has several advantages over others, especially when quantifying asymmetry in flocculent galaxies. Specifically, the averaging of large regions improves the signal-to-noise in the measurements; the method is not strongly affected by the presence of spiral arms; and it identifies the kinds of asymmetry that are likely to be dynamically important. Application of this "method of sectors" to R-band images of 32 disk galaxies indicates that about 30% of spirals show deviations from axisymmetry at the 5-sigma level.Comment: 17 pages, 2 tables and 6 figures, uses psfig and AAS LaTex; to appear in A

Mining for diagnostic information in body surface potential maps: A comparison of feature selection techniques

BACKGROUND: In body surface potential mapping, increased spatial sampling is used to allow more accurate detection of a cardiac abnormality. Although diagnostically superior to more conventional electrocardiographic techniques, the perceived complexity of the Body Surface Potential Map (BSPM) acquisition process has prohibited its acceptance in clinical practice. For this reason there is an interest in striking a compromise between the minimum number of electrocardiographic recording sites required to sample the maximum electrocardiographic information. METHODS: In the current study, several techniques widely used in the domains of data mining and knowledge discovery have been employed to mine for diagnostic information in 192 lead BSPMs. In particular, the Single Variable Classifier (SVC) based filter and Sequential Forward Selection (SFS) based wrapper approaches to feature selection have been implemented and evaluated. Using a set of recordings from 116 subjects, the diagnostic ability of subsets of 3, 6, 9, 12, 24 and 32 electrocardiographic recording sites have been evaluated based on their ability to correctly asses the presence or absence of Myocardial Infarction (MI). RESULTS: It was observed that the wrapper approach, using sequential forward selection and a 5 nearest neighbour classifier, was capable of choosing a set of 24 recording sites that could correctly classify 82.8% of BSPMs. Although the filter method performed slightly less favourably, the performance was comparable with a classification accuracy of 79.3%. In addition, experiments were conducted to show how (a) features chosen using the wrapper approach were specific to the classifier used in the selection model, and (b) lead subsets chosen were not necessarily unique. CONCLUSION: It was concluded that both the filter and wrapper approaches adopted were suitable for guiding the choice of recording sites useful for determining the presence of MI. It should be noted however that in this study recording sites have been suggested on their ability to detect disease and such sites may not be optimal for estimating body surface potential distributions

The AMIGA sample of isolated galaxies. IV. A catalogue of neighbours around isolated galaxies

Studies of the effects of environment on galaxy properties and evolution require well defined control samples. Such isolated galaxy samples have up to now been small or poorly defined. The AMIGA project (Analysis of the interstellar Medium of Isolated GAlaxies) represents an attempt to define a statistically useful sample of the most isolated galaxies in the local (z < 0.05) Universe. A suitable large sample for the AMIGA project already exists, the Catalogue of Isolated Galaxies (CIG, Karachentseva 1973; 1050 galaxies), and we use this sample as a starting point to refine and perform a better quantification of its isolation properties. Digitised POSS-I E images were analysed out to a minimum projected radius R > 0.5 Mpc around 950 CIG galaxies (those within Vr = 1500 km s-1 were excluded). We identified all galaxy candidates in each field brighter than B = 17.5 with a high degree of confidence using the LMORPHO software. We generated a catalogue of approximately 54 000 potential neighbours (redshifts exist for 30% of this sample). Six hundred sixty-six galaxies pass and two hundred eighty-four fail the original CIG isolation criterion. The available redshift data confirm that our catalogue involves a largely background population rather than physically associated neighbours. We find that the exclusion of neighbours within a factor of four in size around each CIG galaxy, employed in the original isolation criterion, corresponds to Delta Vr ~ 18000 km s-1 indicating that it was a conservative limit. Galaxies in the CIG have been found to show different degrees of isolation. We conclude that a quantitative measure of this is mandatory. It will be the subject of future work based on the catalogue of neighbours obtained here.Comment: Accepted by A&A, 10 pages, 8 figures, 4 table

Simulation of human factors for material safety

The D-1 (Statistical Sciences) and D-2 (Stockpile Complex Modeling and Analysis) groups frequently collaborate to analyze production capabilities at Los Alamos National Laboratory. The facilities in question run the gamut from traditional machining to the fabrication of Plutonium components. This paper documents our efforts to extend our modeling capabilities from traditional discrete event simulation modeling to include agent based models

A Runaway Black Hole in COSMOS: Gravitational Wave or Slingshot Recoil?

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.Comment: 13 figures; submitted to ApJ. Sent back to the referee after the first interaction and awaiting the final comment

Local Star formation triggered by SN shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a SNR with a neutral magnetized cloud we derived analytically a set of conditions which are favorable for driving gravitational instability in the cloud and thus star formation. We have built diagrams of the SNR radius, versus the cloud density, that constrain a domain in the parameter space where star formation is allowed. The diagrams are also tested with fully 3-D MHD simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity ~1 $mu$G results only a small shrinking of the star formation triggering zone in the diagrams, a larger magnetic field (~10 $mu$G) causes a significant shrinking, as expected. Applications of the diagrams to a few regions of our own galaxy have revealed that star formation in those sites could have been triggered by shock waves from SNRs. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is smaller than the observed values in our own Galaxy (sfe ~0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star forming galaxies, not even when the magnetic field in the neutral clouds is neglected. (abridged)Comment: 19 pages, 13 figures, accepted for pubblication in MNRA